Nghiên cứu phương pháp xác định vị trí sự cố trên đường dây truyền tải điện thuộc lưới điện phức tạp700

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI NGUYỄN XUÂN VINH NGHIÊN CỨU PHƯƠNG PHÁP XÁC ĐỊNH VỊ TRÍ SỰ CỐ TRÊN ĐƯỜNG DÂY TRUYỀN TẢI ĐIỆN THUỘC LƯỚI ĐIỆN PHỨC TẠP LUẬN ÁN TIẾN SĨ KỸ THUẬT ĐIỆN Hà Nội – 2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI NGUYỄN XUÂN VINH NGHIÊN CỨU PHƯƠNG PHÁP XÁC ĐỊNH VỊ TRÍ SỰ CỐ TRÊN ĐƯỜNG DÂY TRUYỀN TẢI ĐIỆN THUỘC LƯỚI ĐIỆN PHỨC TẠP Chuyên ngành: Kỹ thuật điện Mã số: 62520202 LUẬN ÁN TIẾN SĨ KỸ THUẬT ĐIỆN NGƯỜI HƯỚNG DẪN KHOA HỌC: TS NGUYỄN XUÂN TÙNG TS NGUYỄN ĐỨC HUY Hà Nội – 2017 LỜI CAM ĐOAN Tôi xin cam đoan kết khoa học trình bày luận án thành nghiên cứu thân suốt thời gian làm nghiên cứu sinh chưa xuất công bố tác giả khác Các kết đạt xác trung thực Hà Nội, ngày 30 tháng 10 năm 2017 XÁC NHẬN CỦA TẬP THỂ HƯỚNG DẪN GV HƯỚNG DẪN GV HƯỚNG DẪN TÁC GIẢ LUẬN ÁN TS Nguyễn Xuân Tùng TS Nguyễn Đức Huy Nguyễn Xuân Vinh i LỜI CẢM ƠN Đầu tiên, xin bày tỏ lời cảm ơn sâu sắc kính trọng đến hai thầy hướng dẫn khoa học trực tiếp, TS Nguyễn Xuân Tùng TS Nguyễn Đức Huy trực tiếp hướng dẫn, định hướng khoa học trình nghiên cứu Hai thầy dành nhiều thời gian tâm huyết, hỗ trợ mặt để tác giả hoàn thành luận án Tác giả xin trân trọng cảm ơn Lãnh đạo trường Đại học Bách Khoa Hà Nội, Viện Đào tạo Sau Đại học, Viện Điện Bộ môn Hệ thống Điện tạo điều kiện thuận lợi cho nghiên cứu sinh suốt trình học tập nghiên cứu Chân thành cảm ơn Giảng viên cán Bộ mơn Hệ thống điện, hỗ trợ tận tình giúp đỡ q trình thực luận án Tơi xin trân trọng cảm ơn Ban Giám hiệu trường Đại học Sư phạm Kỹ thuật Vĩnh Long, Ban Chủ nhiệm khoa Điện – Điện tử tạo điều kiện thuận lợi cho tác giả tập trung nghiên cứu Hà Nội suốt thời gian qua Xin chân thành cảm ơn quan tâm, giúp đỡ động viên đồng nghiệp, nhóm nghiên cứu sinh Viện Điện Cuối cùng, tác giả thực cảm động từ đáy lịng xin bày tỏ lịng biết ơn đến Ông bà nội ngoại hai bên người vợ yêu quý trai thân yêu bên tác giả lúc khó khăn nhất, lúc mệt mỏi nhất, để động viên, để hỗ trợ tài tinh thần, giúp tác giả đứng vững trình nghiên cứu hồn thiện luận án Tác giả luận án Nguyễn Xuân Vinh ii MỤC LỤC Trang LỜI CAM ĐOAN LỜI CẢM ƠN DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT vi DANH MỤC CÁC BẢNG vi DANH MỤC CÁC HÌNH ẢNH, ĐỒ THỊ viii MỞ ĐẦU 1 Tính cấp thiết đề tài Mục tiêu, đối tượng, phương pháp phạm vi nghiên cứu Các đóng góp luận án Cấu trúc nội dung luận án TỔNG QUAN .7 1.1 Tình hình nghiên cứu nước 1.2 Tình hình nghiên cứu ngồi nước .7 Thuật toán định vị cố áp dụng đường dây đơn đồng với u cầu biết xác thơng số đường dây Thuật tốn định vị cố khơng u cầu thơng số đường dây áp dụng với đường dây đơn đồng .10 Thuật toán định vị cố áp dụng cho đường dây đơn không đồng với yêu cầu biết xác thơng số đường dây 12 Thuật toán định vị cố áp dụng cho đường dây truyền tải điện rẽ nhánh với u cầu biết xác thơng số đường dây 14 Đánh giá ảnh hưởng sai số đo lường đến kết định vị cố .16 1.3 Những vấn đề tồn hướng nghiên cứu .17 THUẬT TOÁN ĐỊNH VỊ SỰ CỐ CẢI TIẾN ÁP DỤNG VỚI ĐƯỜNG DÂY ĐƠN KHÔNG ĐỒNG NHẤT 18 2.1 Đặt vấn đề .1 2.2 Thuật toán xác định thông số phân đoạn thuộc đường dây khơng đồng sử dụng tín hiệu đo lường từ hai phía đầu đường dây 18 Mạng hai cửa tương đương đường dây khơng đồng .18 Thuật tốn xác định thơng số phân đoạn thuộc đường dây đơn không đồng .2 2.3 Thuật toán định vị cố sử dụng tín hiệu đo lường từ hai phía đầu đường dây .2 Bước 1: Giả thiết cố xảy phân đoạn SC 26 Bước 2: Giả thiết cố xảy phân đoạn RC 27 iii Bước 3: Xác định vị trí cố 28 2.4 Mơ kiểm chứng thuật tốn định vị cố xảy đường dây truyền tải điện không đồng 29 Kết định vị cố 30 So sánh kết định vị cố thuật toán đề xuất với kết định vị cố thuật toán [48] .30 2.5 Kết luận 3 PHƯƠNG PHÁP ĐỊNH VỊ SỰ CỐ KHÔNG BIẾT TRƯỚC THÔNG SỐ ĐƯỜNG DÂY ÁP DỤNG VỚI ĐƯỜNG DÂY DÂY TRUYỀN TẢI ĐIỆN RẼ NHÁNH KHÔNG ĐỒNG NHẤT CHỦNG LOẠI DÂY 33 3.1 Đặt vấn đề .3 3.2 Thuật tốn xác định thơng số đường dây đường dây rẽ nhánh khơng đồng sử dụng tín hiệu đo lường từ ba đầu đường dây 33 Mạng hai cửa tương đương đường dây rẽ nhánh 33 Trường hợp phân đoạn đường dây rẽ nhánh có thông số đường dây khác 3.3 Thuật tốn xác định vị trí cố cho đường dây rẽ nhánh khơng đồng sử dụng tín hiệu đo lường từ ba đầu đường dây 49 Bước 1: Giả thiết cố xảy phân đoạn SJ 49 Bước 2: Giả thiết cố xảy phân đoạn RJ .49 Bước 3: Giả thiết cố xảy phân đoạn TJ 49 Bước 4: Xác định vị trí cố 50 3.4 Thuật tốn xác định thơng số đường dây đồng thời định vị cố áp dụng với đường dây rẽ nhánh, phân đoạn có chủng loại dây, sử dụng ghi cố .5 Bước 1: Giả thiết cố xảy phân đoạn SJ 51 Bước 2: Giả thiết cố xảy phân đoạn RJ .51 Bước 3: Giả thiết cố xảy phân đoạn TJ 51 Bước 4: Xác định thơng số đường dây vị trí cố 52 3.5 Mơ kiểm chứng thuật tốn định vị cố xảy đường dây truyền tải điện rẽ nhánh Trường hợp đường dây truyền tải điện rẽ nhánh có phân đoạn sử dụng chủng loại dây khác 53 Trường hợp đường dây rẽ nhánh có hai phân đoạn có thông số đường dây khác thông số phân đoạn thứ ba 55 Trường hợp phân đoạn đường dây rẽ nhánh có thơng số đường dây giống 3.6 Kết luận ĐÁNH GIÁ ẢNH HƯỞNG CỦA SAI SỐ ĐO LƯỜNG ĐẾN KẾT QUẢ ĐỊNH VỊ SỰ CỐ .6 iv 4.1 Cơ sở đề xuất phương pháp 62 4.2 Mô Monte Carlo 63 4.3 Đánh giá ảnh hưởng sai số đo lường đến kết định vị cố thuật tốn xác định vị trí cố xảy đường dây truyền tải điện không đồng 64 Kịch mô .65 Đặc tính thống kê 65 4.4 Đánh giá ảnh hưởng sai số đo lường đến kết định vị cố thuật toán định vị cố xảy đường dây truyền tải điện rẽ nhánh 69 Trường hợp phân đoạn thuộc đường dây rẽ nhánh có thơng số đường dây khác Trường hợp hai phân đoạn có thơng số đường dây khác thông số phân đoạn thứ ba 73 Trường hợp phân đoạn thuộc đường dây rẽ nhánh sử dụng chủng loại dây .7 4.5 Kết luận KẾT LUẬN VÀ KIẾN NGHỊ 83 Đóng góp khoa học luận án 83 Kiến nghị nghiên cứu .84 TÀI LIỆU THAM KHẢO 85 DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA LUẬN ÁN 91 PHỤ LỤC A.1 Kết mô định vị cố xảy đường dây truyền tải điện không đồng A.2 Kết mô định vị cố xảy đường dây truyền tải điện rẽ nhánh: Trường hợp phân đoạn SJ, RJ, TJ có thơng số có thơng số đường dây khác 96 A.3 Kết định vị sư cố xảy đường dây truyền tải điện rẽ nhánh: Trường hợp phân đoạn SJ, TJ có thơng số đường dây khác thông số đường dây phân đoạn RJ A.4 Kết định vị cố xảy đường dây truyền tải điện rẽ nhánh: Trường hợp đường dây rẽ nhánh đồng nhất, phân đoạn SJ, RJ, TJ có thơng số đường dây 10 A.5 Phương trình (2.48) 102 A.6 Phương trình (3.103) 106 A.7 Phương trình (3.105) 118 A.8 Phương trình (3.111) 126 A.9 Phương trình (3.112) 128 A.10 Phương trình (3.113) 129 v DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT ABC Ngắn mạch pha AG Ngắn mạch pha A với đất BC Ngắn mạch hai pha BC BCG Ngắn mạch hai pha BC với đất AC Dòng điện xoay chiều (Alternating Current) ANFIS Mạng nơ ron thích nghi mờ (Adaptive Neuro-Fuzzy training of Sugeno-type) COMTRADE Tiêu chuẩn định dạng tập tin để trao đổi tín hiệu cố (COMmon format for TRAnsient Data Exchange ) CT Máy biến dòng điện (Current transformer) d Vị trí cố DC Dịng điện chiều (Direct Current) DFT Biến đổi Fourier rời rạc (Discrete Fourier Transform) GPS Hệ thống định vị toàn cầu (Global Positioning System) IEC Ủy ban Kỹ thuật Điện Quốc tế (International Electrotechnical Commission) HT Hệ thống điện IEEE Viện kỹ thuật điện điện tử (Institute of Electrical and Electronics Engineers) imag Phần ảo số phức l Chiều dài đường dây MLP Mạng nơ ron truyền thẳng nhiều lớp (Multilayer Perceptron) real Phần thực số phức VT Máy biến áp (Voltage Transformer) Zc Tổng trở sóng đường dây γ Hệ số truyền sóng đường dây vi DANH MỤC CÁC BẢNG Bảng 1.1 Định nghĩa Vsf , I sf , ΔI sf cho loại cố Bảng 2.1 Thông số phân đoạn SC, l1=100km 29 Bảng 2.2 Thông số phân đoạn CR, l2=20km 29 Bảng 2.3 Thông số nguồn 29 Bảng 2.4 Kết định vị cố xảy đường dây không đồng 30 Bảng 3.1 Thông số nguồn 53 Bảng 3.2 Thông số phân đoạn SJ 53 Bảng 3.3 Thông số phân đoạn RJ 53 Bảng 3.4 Thông số phân đoạn TJ 53 Bảng 3.5 Kết định vị cố xảy đường dây rẽ nhánh không đồng 54 Bảng 3.6 Kết định vị cố đường dây có thông số SJ, TJ khác RJ 56 Bảng 3.7 Kết định vị cố xảy đường dây rẽ nhánh đồng 58 Bảng 4.1 Sai số giới hạn máy biến dòng điện [34] 64 Bảng 4.2 Sai số giới hạn máy biến điện áp [35] 64 Bảng 4.3 Chỉ số thống kê loại cố AG – đường dây không đồng 66 Bảng 4.4 Chỉ số thống kê loại cố BC 67 Bảng 4.5 Chỉ số thống kê loại cố BCG 68 Bảng 4.6 Chỉ số thống kê loại cố ABC 68 Bảng 4.7 Chỉ số thống kê loại cố AG – đường dây rẽ nhánh không đồng 71 Bảng 4.8 Chỉ số thống kê loại cố BC 71 Bảng 4.9 Chỉ số thống kê loại cố BCG 72 Bảng 4.10 Chỉ số thống kê loại cố ABC 73 Bảng 4.11 Chỉ số thống kê loại cố AG – SJ, TJ khác thông số RJ 75 Bảng 4.12 Chỉ số thống kê loại cố BC 76 Bảng 4.13 Chỉ số thống kê loại cố BCG 76 Bảng 4.14 Chỉ số thống kê loại cố ABC 77 Bảng 4.15 Chỉ số thống kê loại cố AG – đường dây rẽ nhánh đồng 79 Bảng 4.16 Chỉ số thống kê loại cố BC 80 Bảng 4.17 Chỉ số thống kê loại cố BCG 80 Bảng 4.18 Chỉ số thống kê loại cố ABC 81 vii DANH MỤC CÁC HÌNH ẢNH, ĐỒ THỊ Hình 1.1 Đường dây truyền tải đồng Hình 1.2 Đường dây không đồng 12 Hình 1.3 Đường dây rẽ nhánh 14 Hình 2.1 Đường dây truyền tải khơng đồng 18 Hình 2.2 Nối tiếp mạng hai cửa 19 Hình 2.3 Bản ghi cố xảy đường dây 220kV Hưng Đông – Hà Tĩnh 20 Hình 2.4 Sự cố xảy phân đoạn SC vị trí F1 21 Hình 2.5 Mạng hai cổng đường dây không đồng cố F1 xảy SC 21 Hình 2.6 Sự cố xảy phân đoạn SC vị trí F2 23 Hình 2.7 Mạng hai cổng đường dây không đồng cố F2 xảy SC 23 Hình 2.8 Sự cố xảy phân đoạn RC 24 Hình 2.9 Mạng hai cổng đường dây khơng đồng cố F2 xảy RC 24 Hình 2.10 Thuật tốn xác định thơng số đường dây không đồng 26 Hình 2.11 Sự cố xảy phân đoạn SC 26 Hình 2.12 Sự cố xảy phân đoạn RC 27 Hình 2.13 Thuật tốn định vị cố xảy đường dây không đồng 28 Hình 2.14 Mơ hình mơ đường dây khơng đồng 29 Hình 2.15 So sánh kết định vị cố thuật toán đề xuất với thuật toán [48] 31 Hình 3.1 Đường dây truyền tải rẽ nhánh 33 Hình 3.2 Mơ hình mạng hai cửa đường dây rẽ nhánh 34 Hình 3.3 Sự cố phân đoạn SJ 39 Hình 3.4 Mơ hình mạng hai cửa cố phân đoạn SJ 39 Hình 3.5 Sự cố phân đoạn RJ 41 Hình 3.6 Mơ hình mạng hai cửa cố phân đoạn RJ 41 Hình 3.7 Sự cố phân đoạn TJ 44 Hình 3.8 Mơ hình mạng hai cửa cố phân đoạn TJ 44 Hình 3.9 Thuật tốn xác định thơng số đường dây rẽ nhánh không đồng 47 Hình 3.10 Thuật tốn xác định thơng số đường dây rẽ nhánh – SJ, TJ thông số khác RJ 48 Hình 3.11 Thuật tốn xác định vị trí cố xảy đường dây truyền tải điện rẽ nhánh50 Hình 3.12 Thuật tốn định vị cố áp dụng với đường dây truyền tải rẽ nhánh đồng nhất52 Hình 3.13 Mơ hình mơ đường dây rẽ nhánh 53 Hình 3.14 So sánh kết định vị cố thuật toán đề xuất với thuật toán [47] 55 Hình 3.15 So sánh kết định vị cố thuật toán đề xuất với thuật toán [47] 57 viii lb*x10)*x4)*cosh(la*x7)*sin(la*x8))*x1/(x1^2+x2^2)+((cosh(db*x9-lb*x9)*sin(db*x10lb*x10)*x3+sinh(db*x9-lb*x9)*cos(db*x10lb*x10)*x4)*sinh(la*x7)*cos(la*x8)+(sinh(db*x9-lb*x9)*cos(db*x10-lb*x10)*x3cosh(db*x9-lb*x9)*sin(db*x10lb*x10)*x4)*cosh(la*x7)*sin(la*x8))*x2/(x1^2+x2^2))*v14+((-cosh(db*x9lb*x9)*sin(db*x10-lb*x10)*x3-sinh(db*x9-lb*x9)*cos(db*x10lb*x10)*x4)*cosh(la*x7)*cos(la*x8)+(-sinh(db*x9-lb*x9)*cos(db*x10lb*x10)*x3+cosh(db*x9-lb*x9)*sin(db*x10-lb*x10)*x4)*sinh(la*x7)*sin(la*x8))*i13+((sinh(db*x9-lb*x9)*cos(db*x10-lb*x10)*x3+cosh(db*x9-lb*x9)*sin(db*x10lb*x10)*x4)*cosh(la*x7)*cos(la*x8)-(-cosh(db*x9-lb*x9)*sin(db*x10-lb*x10)*x3sinh(db*x9-lb*x9)*cos(db*x10-lb*x10)*x4)*sinh(la*x7)*sin(la*x8))*i14=0 f111(x1,x 2,x3 ,x4,x5,x 6,x7 ,x8,x 9,x10 ,x11,x12 ) = real(PT6) = f11 = cosh(dc*x11)*cos(dc*x12)*v23-sinh(dc*x11)*sin(dc*x12)*v24+(-sinh(dc*x11) *cos(dc*x12)*x5+cosh(dc*x11)*sin(dc*x12)*x6)*i23-(-cosh(dc*x11)*sin(dc*x12)*x5sinh(dc*x11)*cos(dc*x12)*x6)*i24+(-cosh(la*x7)*cos(la*x8)*cosh(dc*x11lc*x11)*cos(dc*x12-lc*x12)+sinh(la*x7)*sin(la*x8)*sinh(dc*x11-lc*x11)*sin(dc*x12lc*x12)-((-sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5+cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x6)*sinh(la*x7)*cos(la*x8)-(-cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x1/(x1^2+x2^2)-((-cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(la*x7)*cos(la*x8)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x2/(x1^2+x2^2))*v19-(sinh(la*x7)*sin(la*x8)*cosh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)cosh(la*x7)*cos(la*x8)*sinh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)-((-cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(la*x7)*cos(la*x8)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x1/(x1^2+x2^2)+((-sinh(dc*x11lc*x11)*cos(dc*x12-lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(la*x7)*cos(la*x8)-(-cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x2/(x1^2+x2^2))*v20+((cosh(dc*x11lc*x11)*cos(dc*x12-lc*x12)*x1-sinh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x2)*sinh(la*x7)*cos(la*x8)-(sinh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x1+cosh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x2)*cosh(la*x7)*sin(la*x8)+(sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(la*x7)*cos(la*x8)-(-cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x6)*sinh(la*x7)*sin(la*x8))*i19-((sinh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x1+cosh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x2)*sinh(la*x7)*cos(la*x8)+(cosh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x1sinh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x2)*cosh(la*x7)*sin(la*x8)+(-cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(la*x7)*cos(la*x8)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(la*x7)*sin(la*x8))*i20+(((sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5116 cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)-(cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x3/(x3^2+x4^2)+((cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)+(sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x4/(x3^2+x4^2))*v21-(((cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)+(sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x3/(x3^2+x4^2)((sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5-cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)-(cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x4/(x3^2+x4^2))*v22+((-sinh(dc*x11lc*x11)*cos(dc*x12-lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(lb*x9)*cos(lb*x10)-(-cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x6)*sinh(lb*x9)*sin(lb*x10))*i21-((cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(lb*x9)*cos(lb*x10)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(lb*x9)*sin(lb*x10))*i22=0 f12(x1,x 2,x3 ,x4 ,x5,x6 ,x7,x8 ,x 9,x 10,x 11,x 12) = imag(PT6) = f12 = sinh(dc*x11)*sin(dc*x12)*v23+cosh(dc*x11)*cos(dc*x12)*v24+(-cosh(dc*x11) *sin(dc*x12)*x5-sinh(dc*x11)*cos(dc*x12)*x6)*i23+(sinh(dc*x11)*cos(dc*x12)*x5+cosh(dc*x11)*sin(dc*x12)*x6)*i24+(sinh(la*x7)*sin(la*x8)*cosh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)cosh(la*x7)*cos(la*x8)*sinh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)-((-cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(la*x7)*cos(la*x8)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x1/(x1^2+x2^2)+((-sinh(dc*x11lc*x11)*cos(dc*x12-lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(la*x7)*cos(la*x8)-(-cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x2/(x1^2+x2^2))*v19+(cosh(la*x7)*cos(la*x8)*cosh(dc*x11-lc*x11)*cos(dc*x12lc*x12)+sinh(la*x7)*sin(la*x8)*sinh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)-((-sinh(dc*x11lc*x11)*cos(dc*x12-lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(la*x7)*cos(la*x8)-(-cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x1/(x1^2+x2^2)((-cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(la*x7)*cos(la*x8)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(la*x7)*sin(la*x8))*x2/(x1^2+x2^2))*v20+((sinh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x1+cosh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x2)*sinh(la*x7)*cos(la*x8)+(cosh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x1sinh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x2)*cosh(la*x7)*sin(la*x8)+(-cosh(dc*x11117 lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(la*x7)*cos(la*x8)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(la*x7)*sin(la*x8))*i19+((cosh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x1sinh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x2)*sinh(la*x7)*cos(la*x8)-(sinh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x1+cosh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x2)*cosh(la*x7)*sin(la*x8)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x6)*cosh(la*x7)*cos(la*x8)-(cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(la*x7)*sin(la*x8))*i20+(((cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)+(sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x3/(x3^2+x4^2)((sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5-cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)-(cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x4/(x3^2+x4^2))*v21+(((sinh(dc*x11lc*x11)*cos(dc*x12-lc*x12)*x5-cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)-(cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x3/(x3^2+x4^2)+((cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5+sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(lb*x9)*cos(lb*x10)+(sinh(dc*x11-lc*x11)*cos(dc*x12-lc*x12)*x5cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*cosh(lb*x9)*sin(lb*x10))*x4/(x3^2+x4^2))*v22+((-cosh(dc*x11lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*cosh(lb*x9)*cos(lb*x10)+(-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12lc*x12)*x6)*sinh(lb*x9)*sin(lb*x10))*i21+((-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x5+cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x6)*cosh(lb*x9)*cos(lb*x10)-(cosh(dc*x11-lc*x11)*sin(dc*x12-lc*x12)*x5-sinh(dc*x11-lc*x11)*cos(dc*x12lc*x12)*x6)*sinh(lb*x9)*sin(lb*x10))*i22=0 A.7 Phương trình (3.105) f1(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = real(PT1) = f1 = v1+(-cosh(la*x5)*cos(la*x6)*cosh(lb*x7)*cos(lb*x8)+sinh(la*x5)*sin(la*x6) *sinh(lb*x7)*sin(lb*x8)-((x1*sinh(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x3/(x3^ 2+x4^2)((x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)+(x1*sin h(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x4/(x3^2+x4^2))*v3-(sinh(la*x5)*sin(la*x6)*cosh(lb*x7)*cos(lb*x8)cosh(la*x5)*cos(la*x6)*sinh(lb*x7)*sin(lb*x8)((x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)+(x1*sin 118 h(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x3/(x3^2+x4^2)+((x1*sinh(la*x5)*c os(la*x6)-x2*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x4/(x3^ 2+x4^2))*v4+((cosh(la*x5)*cos(la*x6)*x3sinh(la*x5)*sin(la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)(sinh(la*x5)*sin(la*x6)*x3+cosh(la*x5)*cos(la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(x1*sin h(la*x5)*cos(la*x6)-x2*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*cos(lb*x8)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*sin(lb*x8))*i3((sinh(la*x5)*sin(la*x6)*x3+cosh(la*x5)*cos(la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)+(cosh(l a*x5)*cos(la*x6)*x3sinh(la*x5)*sin(la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(x2*sinh(la*x5)*cos(la*x6)+x1*cosh (la*x5)*sin(la*x6))*cosh(lb*x7)*cos(lb*x8)+(x1*sinh(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*sin(lb*x8))*i4+(sinh(la*x5)*cos(la*x6)*sinh(lc*x5)*cos(lc*x6)+cosh(la*x5)*sin(la*x6)*cosh(lc*x5)*sin(l c*x6))*v5-(-cosh(la*x5)*sin(la*x6)*sinh(lc*x5)*cos(lc*x6)sinh(la*x5)*cos(la*x6)*cosh(lc*x5)*sin(lc*x6))*v6+((x1*sinh(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*cosh(lc*x5)*cos(lc*x6)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lc*x5)*sin(lc*x6))*i5((x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*cosh(lc*x5)*cos(lc*x6)+(x1*sin h(la*x5)*cos(la*x6)-x2*cosh(la*x5)*sin(la*x6))*sinh(lc*x5)*sin(lc*x6))*i6=0 f2(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = imag (PT1) = f2 = v2+(-sinh(la*x5)*sin(la*x6)*cosh(lb*x7)*cos(lb*x8)-cosh(la*x5)*cos(la*x6) *sinh(lb*x7)*sin(lb*x8)((x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)+(x1*sin h(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x3/(x3^2+x4^2)+((x1*sinh(la*x5)*c os(la*x6)-x2*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x4/(x3^ 2+x4^2))*v3+(cosh(la*x5)*cos(la*x6)*cosh(lb*x7)*cos(lb*x8)+sinh(la*x5)*sin(la*x6)*sinh(lb*x7)*sin(l b*x8)-((x1*sinh(la*x5)*cos(la*x6)-x2*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x3/(x3^ 2+x4^2)((x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*cos(lb*x8)+(x1*sin h(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*sin(lb*x8))*x4/(x3^2+x4^2))*v4+((sinh(la*x5)*s in(la*x6)*x3+cosh(la*x5)*cos(la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)+(cosh(la*x5)*cos(la*x 6)*x3sinh(la*x5)*sin(la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(x2*sinh(la*x5)*cos(la*x6)+x1*cosh (la*x5)*sin(la*x6))*cosh(lb*x7)*cos(lb*x8)+(x1*sinh(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*sin(lb*x8))*i3+((cosh(la*x5)*cos(la*x6)*x3sinh(la*x5)*sin(la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)(sinh(la*x5)*sin(la*x6)*x3+cosh(la*x5)*cos(la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(x1*sin h(la*x5)*cos(la*x6)-x2*cosh(la*x5)*sin(la*x6))*cosh(lb*x7)*cos(lb*x8)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lb*x7)*sin(lb*x8))*i4+(cosh(la*x5)*sin(la*x6)*sinh(lc*x5)*cos(lc*x6)119 sinh(la*x5)*cos(la*x6)*cosh(lc*x5)*sin(lc*x6))*v5+(sinh(la*x5)*cos(la*x6)*sinh(lc*x5)*cos(lc*x6)+cosh(la*x5)*sin(la*x6)*cosh(lc*x5)*sin(l c*x6))*v6+((x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*cosh(lc*x5)*cos(lc* x6)+(x1*sinh(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*sinh(lc*x5)*sin(lc*x6))*i5+((x1*sinh(la*x5)*cos(la*x6)x2*cosh(la*x5)*sin(la*x6))*cosh(lc*x5)*cos(lc*x6)(x2*sinh(la*x5)*cos(la*x6)+x1*cosh(la*x5)*sin(la*x6))*sinh(lc*x5)*sin(lc*x6))*i6=0 f3(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = real(PT2) = f3 = v3+(-cosh(lb*x7)*cos(lb*x8)*cosh(lc*x5)*cos(lc*x6)+sinh(lb*x7)*sin(lb*x8) *sinh(lc*x5)*sin(lc*x6)-((x3*sinh(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x1/(x1^ 2+x2^2)((x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)+(x3*sin h(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x2/(x1^2+x2^2))*v5-(sinh(lb*x7)*sin(lb*x8)*cosh(lc*x5)*cos(lc*x6)cosh(lb*x7)*cos(lb*x8)*sinh(lc*x5)*sin(lc*x6)((x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)+(x3*sin h(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x1/(x1^2+x2^2)+((x3*sinh(lb*x7)*c os(lb*x8)-x4*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x2/(x1^ 2+x2^2))*v6+((cosh(lb*x7)*cos(lb*x8)*x1sinh(lb*x7)*sin(lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)(sinh(lb*x7)*sin(lb*x8)*x1+cosh(lb*x7)*cos(lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(x3*sinh (lb*x7)*cos(lb*x8)-x4*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*cos(lc*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*sin(lc*x6))*i5((sinh(lb*x7)*sin(lb*x8)*x1+cosh(lb*x7)*cos(lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)+(cosh(l b*x7)*cos(lb*x8)*x1sinh(lb*x7)*sin(lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh( lb*x7)*sin(lb*x8))*cosh(lc*x5)*cos(lc*x6)+(x3*sinh(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*sin(lc*x6))*i6+(((x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)-(x4*sinh(lb*x7)*cos(lb*x8)x3*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x1/(x1^2+x2^2)+((x4*sinh(lb*x7)*cos(lb*x8)-x3*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)+(x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x2/(x1^2 +x2^2))*v1-(((-x4*sinh(lb*x7)*cos(lb*x8)x3*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)+(x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x1/(x1^2 +x2^2)-((x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)-(x4*sinh(lb*x7)*cos(lb*x8)x3*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x2/(x1^2+x2^2))*v2+((x3*sinh(lb*x 7)*cos(lb*x8)-x4*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*cos(la*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*sin(la*x6))*i1120 ((x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*cos(la*x6)+(x3*si nh(lb*x7)*cos(lb*x8)-x4*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*sin(la*x6))*i2=0 f4(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = imag (PT2) = f4 = v4+(-sinh(lb*x7)*sin(lb*x8)*cosh(lc*x5)*cos(lc*x6)-cosh(lb*x7)*cos(lb*x8) *sinh(lc*x5)*sin(lc*x6)((x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)+(x3*sin h(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x1/(x1^2+x2^2)+((x3*sinh(lb*x7)*c os(lb*x8)-x4*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x2/(x1^ 2+x2^2))*v5+(cosh(lb*x7)*cos(lb*x8)*cosh(lc*x5)*cos(lc*x6)+sinh(lb*x7)*sin(lb*x8)*sinh(lc*x5)*sin(l c*x6)-((x3*sinh(lb*x7)*cos(lb*x8)-x4*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x1/(x1^ 2+x2^2)((x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*cos(lc*x6)+(x3*sin h(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*sin(lc*x6))*x2/(x1^2+x2^2))*v6+((sinh(lb*x7)*s in(lb*x8)*x1+cosh(lb*x7)*cos(lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)+(cosh(lb*x7)*cos(lb*x 8)*x1sinh(lb*x7)*sin(lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh( lb*x7)*sin(lb*x8))*cosh(lc*x5)*cos(lc*x6)+(x3*sinh(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*sin(lc*x6))*i5+((cosh(lb*x7)*cos(lb*x8)*x1sinh(lb*x7)*sin(lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)(sinh(lb*x7)*sin(lb*x8)*x1+cosh(lb*x7)*cos(lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(x3*sinh (lb*x7)*cos(lb*x8)-x4*cosh(lb*x7)*sin(lb*x8))*cosh(lc*x5)*cos(lc*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(lc*x5)*sin(lc*x6))*i6+(((x4*sinh(lb*x7)*cos(lb*x8)-x3*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)+(x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x1/(x1^2 +x2^2)-((x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)-(x4*sinh(lb*x7)*cos(lb*x8)x3*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x2/(x1^2+x2^2))*v1+(((x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)-(x4*sinh(lb*x7)*cos(lb*x8)x3*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x1/(x1^2+x2^2)+((x4*sinh(lb*x7)*cos(lb*x8)-x3*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*cos(la*x6)+(x3*sinh(lb*x7)*cos(lb*x8)+x4*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*sin(la*x6))*x2/(x1^2 +x2^2))*v2+((x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*cos(la *x6)+(x3*sinh(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*sin(la*x6))*i1+((x3*sinh(lb*x7)*cos(lb*x8)x4*cosh(lb*x7)*sin(lb*x8))*cosh(la*x5)*cos(la*x6)(x4*sinh(lb*x7)*cos(lb*x8)+x3*cosh(lb*x7)*sin(lb*x8))*sinh(la*x5)*sin(la*x6))*i2=0 f5(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = real(PT4) = f5 = cosh(da*x5)*cos(da*x6)*v7-sinh(da*x5)*sin(da*x6)*v8+(-sinh(da*x5)*cos(da*x6) *x1+cosh(da*x5)*sin(da*x6)*x2)*i7-(-cosh(da*x5)*sin(da*x6)*x1121 sinh(da*x5)*cos(da*x6)*x2)*i8+(-cosh(da*x5-la*x5)*cos(da*x6la*x6)*cosh(lb*x7)*cos(lb*x8)+sinh(da*x5-la*x5)*sin(da*x6la*x6)*sinh(lb*x7)*sin(lb*x8)-((-sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x1+cosh(da*x5la*x5)*sin(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)-(-cosh(da*x5-la*x5)*sin(da*x6la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x3/(x3^2+x4^2)-((-cosh(da*x5-la*x5)*sin(da*x6la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)+(sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x4/(x3^2+x4^2))*v9-(-sinh(da*x5-la*x5)*sin(da*x6la*x6)*cosh(lb*x7)*cos(lb*x8)-cosh(da*x5-la*x5)*cos(da*x6la*x6)*sinh(lb*x7)*sin(lb*x8)-((-cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5la*x5)*cos(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)+(-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x3/(x3^2+x4^2)+((-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)-(cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x4/(x3^2+x4^2))*v10+((cosh(da*x5la*x5)*cos(da*x6-la*x6)*x3-sinh(da*x5-la*x5)*sin(da*x6la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)-(sinh(da*x5-la*x5)*sin(da*x6-la*x6)*x3+cosh(da*x5la*x5)*cos(da*x6-la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x2)*cosh(lb*x7)*cos(lb*x8)-(cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*sinh(lb*x7)*sin(lb*x8))*i9-((sinh(da*x5-la*x5)*sin(da*x6la*x6)*x3+cosh(da*x5-la*x5)*cos(da*x6la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)+(cosh(da*x5-la*x5)*cos(da*x6-la*x6)*x3sinh(da*x5-la*x5)*sin(da*x6-la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(-cosh(da*x5la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*cosh(lb*x7)*cos(lb*x8)+(-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*sinh(lb*x7)*sin(lb*x8))*i10+(sinh(da*x5-la*x5)*cos(da*x6la*x6)*sinh(lc*x5)*cos(lc*x6)-cosh(da*x5-la*x5)*sin(da*x6la*x6)*cosh(lc*x5)*sin(lc*x6))*v11-(cosh(da*x5-la*x5)*sin(da*x6la*x6)*sinh(lc*x5)*cos(lc*x6)+sinh(da*x5-la*x5)*cos(da*x6la*x6)*cosh(lc*x5)*sin(lc*x6))*v12+((-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x2)*cosh(lc*x5)*cos(lc*x6)-(cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*sinh(lc*x5)*sin(lc*x6))*i11-((-cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x2)*cosh(lc*x5)*cos(lc*x6)+(-sinh(da*x5la*x5)*cos(da*x6-la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*sinh(lc*x5)*sin(lc*x6))*i12=0 f6(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = imag (PT4) = f6 = sinh(da*x5)*sin(da*x6)*v7+cosh(da*x5)*cos(da*x6)*v8+(-cosh(da*x5)*sin(da*x6) *x1-sinh(da*x5)*cos(da*x6)*x2)*i7+(sinh(da*x5)*cos(da*x6)*x1+cosh(da*x5)*sin(da*x6)*x2)*i8+(-sinh(da*x5la*x5)*sin(da*x6-la*x6)*cosh(lb*x7)*cos(lb*x8)-cosh(da*x5-la*x5)*cos(da*x6la*x6)*sinh(lb*x7)*sin(lb*x8)-((-cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5la*x5)*cos(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)+(-sinh(da*x5-la*x5)*cos(da*x6122 la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x3/(x3^2+x4^2)+((-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)-(cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x4/(x3^2+x4^2))*v9+(-cosh(da*x5la*x5)*cos(da*x6-la*x6)*cosh(lb*x7)*cos(lb*x8)+sinh(da*x5-la*x5)*sin(da*x6la*x6)*sinh(lb*x7)*sin(lb*x8)-((-sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x1+cosh(da*x5la*x5)*sin(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)-(-cosh(da*x5-la*x5)*sin(da*x6la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x3/(x3^2+x4^2)-((-cosh(da*x5-la*x5)*sin(da*x6la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x2)*sinh(lb*x7)*cos(lb*x8)+(sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*cosh(lb*x7)*sin(lb*x8))*x4/(x3^2+x4^2))*v10+((sinh(da*x5la*x5)*sin(da*x6-la*x6)*x3+cosh(da*x5-la*x5)*cos(da*x6la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)+(cosh(da*x5-la*x5)*cos(da*x6-la*x6)*x3sinh(da*x5-la*x5)*sin(da*x6-la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(-cosh(da*x5la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*cosh(lb*x7)*cos(lb*x8)+(-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*sinh(lb*x7)*sin(lb*x8))*i9+((cosh(da*x5-la*x5)*cos(da*x6-la*x6)*x3sinh(da*x5-la*x5)*sin(da*x6-la*x6)*x4)*sinh(lb*x7)*cos(lb*x8)-(sinh(da*x5la*x5)*sin(da*x6-la*x6)*x3+cosh(da*x5-la*x5)*cos(da*x6la*x6)*x4)*cosh(lb*x7)*sin(lb*x8)+(-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x2)*cosh(lb*x7)*cos(lb*x8)-(cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*sinh(lb*x7)*sin(lb*x8))*i10+(cosh(da*x5-la*x5)*sin(da*x6la*x6)*sinh(lc*x5)*cos(lc*x6)+sinh(da*x5-la*x5)*cos(da*x6la*x6)*cosh(lc*x5)*sin(lc*x6))*v11+(sinh(da*x5-la*x5)*cos(da*x6la*x6)*sinh(lc*x5)*cos(lc*x6)-cosh(da*x5-la*x5)*sin(da*x6la*x6)*cosh(lc*x5)*sin(lc*x6))*v12+((-cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1sinh(da*x5-la*x5)*cos(da*x6-la*x6)*x2)*cosh(lc*x5)*cos(lc*x6)+(-sinh(da*x5la*x5)*cos(da*x6-la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6la*x6)*x2)*sinh(lc*x5)*sin(lc*x6))*i11+((-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x1+cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x2)*cosh(lc*x5)*cos(lc*x6)-(cosh(da*x5-la*x5)*sin(da*x6-la*x6)*x1-sinh(da*x5-la*x5)*cos(da*x6la*x6)*x2)*sinh(lc*x5)*sin(lc*x6))*i12=0 f7(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = real(PT5) = f7 = cosh(db*x7)*cos(db*x8)*v15-sinh(db*x7)*sin(db*x8)*v16+(-sinh(db*x7) *cos(db*x8)*x3+cosh(db*x7)*sin(db*x8)*x4)*i15-(-cosh(db*x7)*sin(db*x8)*x3sinh(db*x7)*cos(db*x8)*x4)*i16+(-cosh(db*x7-lb*x7)*cos(db*x8lb*x8)*cosh(lc*x5)*cos(lc*x6)+sinh(db*x7-lb*x7)*sin(db*x8lb*x8)*sinh(lc*x5)*sin(lc*x6)-((-sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7lb*x7)*sin(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)-(-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x1/(x1^2+x2^2)-((-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8123 lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x2/(x1^2+x2^2))*v17-(-sinh(db*x7-lb*x7)*sin(db*x8lb*x8)*cosh(lc*x5)*cos(lc*x6)-cosh(db*x7-lb*x7)*cos(db*x8lb*x8)*sinh(lc*x5)*sin(lc*x6)-((-cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7lb*x7)*cos(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)+(-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x1/(x1^2+x2^2)+((-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)-(cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x2/(x1^2+x2^2))*v18+((cosh(db*x7lb*x7)*cos(db*x8-lb*x8)*x1-sinh(db*x7-lb*x7)*sin(db*x8lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)-(sinh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x1+cosh(db*x7lb*x7)*cos(db*x8-lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*cosh(lc*x5)*cos(lc*x6)-(cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*sinh(lc*x5)*sin(lc*x6))*i17-((sinh(db*x7-lb*x7)*sin(db*x8lb*x8)*x1+cosh(db*x7-lb*x7)*cos(db*x8lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)+(cosh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x1-sinh(db*x7lb*x7)*sin(db*x8-lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x4)*cosh(lc*x5)*cos(lc*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*sinh(lc*x5)*sin(lc*x6))*i18+(((sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*sinh(la*x5)*cos(la*x6)-(cosh(db*x7lb*x7)*sin(db*x8-lb*x8)*x3+sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x1/(x1^2+x2^2)+((cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3+sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*sinh(la*x5)*cos(la*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x2/(x1^2+x2^2))*v13(((cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3+sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*sinh(la*x5)*cos(la*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x1/(x1^2+x2^2)((sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*sinh(la*x5)*cos(la*x6)-(cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3+sinh(db*x7lb*x7)*cos(db*x8-lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x2/(x1^2+x2^2))*v14+((sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*cosh(la*x5)*cos(la*x6)-(-cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x4)*sinh(la*x5)*sin(la*x6))*i13-((-cosh(db*x7lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*cosh(la*x5)*cos(la*x6)+(-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*sinh(la*x5)*sin(la*x6))*i14=0 f8(x1,x 2,x3,x4,x 5,x6,x 7,x 8) = imag(PT5) = f8 = sinh(db*x7)*sin(db*x8)*v15+cosh(db*x7)*cos(db*x8)*v16+(-cosh(db*x7) *sin(db*x8)*x3-sinh(db*x7)*cos(db*x8)*x4)*i15+(sinh(db*x7)*cos(db*x8)*x3+cosh(db*x7)*sin(db*x8)*x4)*i16+(-sinh(db*x7lb*x7)*sin(db*x8-lb*x8)*cosh(lc*x5)*cos(lc*x6)-cosh(db*x7-lb*x7)*cos(db*x8lb*x8)*sinh(lc*x5)*sin(lc*x6)-((-cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7lb*x7)*cos(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)+(-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8124 lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x1/(x1^2+x2^2)+((-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)-(cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x2/(x1^2+x2^2))*v17+(-cosh(db*x7lb*x7)*cos(db*x8-lb*x8)*cosh(lc*x5)*cos(lc*x6)+sinh(db*x7-lb*x7)*sin(db*x8lb*x8)*sinh(lc*x5)*sin(lc*x6)-((-sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7lb*x7)*sin(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)-(-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x1/(x1^2+x2^2)-((-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x4)*sinh(lc*x5)*cos(lc*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*cosh(lc*x5)*sin(lc*x6))*x2/(x1^2+x2^2))*v18+((sinh(db*x7lb*x7)*sin(db*x8-lb*x8)*x1+cosh(db*x7-lb*x7)*cos(db*x8lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)+(cosh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x1-sinh(db*x7lb*x7)*sin(db*x8-lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x4)*cosh(lc*x5)*cos(lc*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*sinh(lc*x5)*sin(lc*x6))*i17+((cosh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x1sinh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x2)*sinh(lc*x5)*cos(lc*x6)-(sinh(db*x7lb*x7)*sin(db*x8-lb*x8)*x1+cosh(db*x7-lb*x7)*cos(db*x8lb*x8)*x2)*cosh(lc*x5)*sin(lc*x6)+(-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*cosh(lc*x5)*cos(lc*x6)-(cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*sinh(lc*x5)*sin(lc*x6))*i18+(((cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x3+sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*sinh(la*x5)*cos(la*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x1/(x1^2+x2^2)((sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*sinh(la*x5)*cos(la*x6)-(cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3+sinh(db*x7lb*x7)*cos(db*x8lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x2/(x1^2+x2^2))*v13+(((sinh(db*x7lb*x7)*cos(db*x8-lb*x8)*x3-cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*sinh(la*x5)*cos(la*x6)-(cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3+sinh(db*x7lb*x7)*cos(db*x8-lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x1/(x1^2+x2^2)+((cosh(db*x7lb*x7)*sin(db*x8-lb*x8)*x3+sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*sinh(la*x5)*cos(la*x6)+(sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*cosh(la*x5)*sin(la*x6))*x2/(x1^2+x2^2))*v14+((-cosh(db*x7lb*x7)*sin(db*x8-lb*x8)*x3-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x4)*cosh(la*x5)*cos(la*x6)+(-sinh(db*x7-lb*x7)*cos(db*x8lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x4)*sinh(la*x5)*sin(la*x6))*i13+((sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x3+cosh(db*x7-lb*x7)*sin(db*x8lb*x8)*x4)*cosh(la*x5)*cos(la*x6)-(-cosh(db*x7-lb*x7)*sin(db*x8-lb*x8)*x3sinh(db*x7-lb*x7)*cos(db*x8-lb*x8)*x4)*sinh(la*x5)*sin(la*x6))*i14=0 125 A.8 Phương trình (3.111) f1(x1,x 2,x3,x4,x 5) = real(PT1) = = v1-(cosh(la*x3)*cos(la*x4)*cosh(lb*x3)*cos(lb*x4)f1 sinh(la*x3)*sin(la*x4)*sinh(lb*x3)*sin(lb*x4)+sinh(la*x3)*cos(la*x4)*sinh(lb*x3)*cos(l b*x4)cosh(la*x3)*sin(la*x4)*cosh(lb*x3)*sin(lb*x4))*v3+(sinh(la*x3)*sin(la*x4)*cosh(lb*x3) *cos(lb*x4)+cosh(la*x3)*cos(la*x4)*sinh(lb*x3)*sin(lb*x4)+cosh(la*x3)*sin(la*x4)*sinh (lb*x3)*cos(lb*x4)+sinh(la*x3)*cos(la*x4)*cosh(lb*x3)*sin(lb*x4))*v4((cosh(la*x3)*cos(la*x4)*x1-sinh(la*x3)*sin(la*x4)*x2)*sinh(lb*x3)*cos(lb*x4)(sinh(la*x3)*sin(la*x4)*x1+cosh(la*x3)*cos(la*x4)*x2)*cosh(lb*x3)*sin(lb*x4)+(x1*sin h(la*x3)*cos(la*x4)-x2*cosh(la*x3)*sin(la*x4))*cosh(lb*x3)*cos(lb*x4)(x2*sinh(la*x3)*cos(la*x4)+x1*cosh(la*x3)*sin(la*x4))*sinh(lb*x3)*sin(lb*x4))*i3+((sin h(la*x3)*sin(la*x4)*x1+cosh(la*x3)*cos(la*x4)*x2)*sinh(lb*x3)*cos(lb*x4)+(cosh(la*x3 )*cos(la*x4)*x1sinh(la*x3)*sin(la*x4)*x2)*cosh(lb*x3)*sin(lb*x4)+(x2*sinh(la*x3)*cos(la*x4)+x1*cosh (la*x3)*sin(la*x4))*cosh(lb*x3)*cos(lb*x4)+(x1*sinh(la*x3)*cos(la*x4)x2*cosh(la*x3)*sin(la*x4))*sinh(lb*x3)*sin(lb*x4))*i4(sinh(la*x3)*cos(la*x4)*sinh(lc*x3)*cos(lc*x4)cosh(la*x3)*sin(la*x4)*cosh(lc*x3)*sin(lc*x4))*v5+(cosh(la*x3)*sin(la*x4)*sinh(lc*x3)* cos(lc*x4)+sinh(la*x3)*cos(la*x4)*cosh(lc*x3)*sin(lc*x4))*v6((x1*sinh(la*x3)*cos(la*x4)-x2*cosh(la*x3)*sin(la*x4))*cosh(lc*x3)*cos(lc*x4)(x2*sinh(la*x3)*cos(la*x4)+x1*cosh(la*x3)*sin(la*x4))*sinh(lc*x3)*sin(lc*x4))*i5+((x2 *sinh(la*x3)*cos(la*x4)+x1*cosh(la*x3)*sin(la*x4))*cosh(lc*x3)*cos(lc*x4)+(x1*sinh(la *x3)*cos(la*x4)-x2*cosh(la*x3)*sin(la*x4))*sinh(lc*x3)*sin(lc*x4))*i6 =0 f2(x1,x 2,x3,x4,x 5) = imag(PT1) = = v2f2 (sinh(la*x3)*sin(la*x4)*cosh(lb*x3)*cos(lb*x4)+cosh(la*x3)*cos(la*x4)*sinh(lb*x3)*sin( lb*x4)+cosh(la*x3)*sin(la*x4)*sinh(lb*x3)*cos(lb*x4)+sinh(la*x3)*cos(la*x4)*cosh(lb* x3)*sin(lb*x4))*v3-(cosh(la*x3)*cos(la*x4)*cosh(lb*x3)*cos(lb*x4)sinh(la*x3)*sin(la*x4)*sinh(lb*x3)*sin(lb*x4)+sinh(la*x3)*cos(la*x4)*sinh(lb*x3)*cos(l b*x4)-cosh(la*x3)*sin(la*x4)*cosh(lb*x3)*sin(lb*x4))*v4((sinh(la*x3)*sin(la*x4)*x1+cosh(la*x3)*cos(la*x4)*x2)*sinh(lb*x3)*cos(lb*x4)+(cosh(l a*x3)*cos(la*x4)*x1sinh(la*x3)*sin(la*x4)*x2)*cosh(lb*x3)*sin(lb*x4)+(x2*sinh(la*x3)*cos(la*x4)+x1*cosh (la*x3)*sin(la*x4))*cosh(lb*x3)*cos(lb*x4)+(x1*sinh(la*x3)*cos(la*x4)x2*cosh(la*x3)*sin(la*x4))*sinh(lb*x3)*sin(lb*x4))*i3-((cosh(la*x3)*cos(la*x4)*x1sinh(la*x3)*sin(la*x4)*x2)*sinh(lb*x3)*cos(lb*x4)(sinh(la*x3)*sin(la*x4)*x1+cosh(la*x3)*cos(la*x4)*x2)*cosh(lb*x3)*sin(lb*x4)+(x1*sin h(la*x3)*cos(la*x4)-x2*cosh(la*x3)*sin(la*x4))*cosh(lb*x3)*cos(lb*x4)(x2*sinh(la*x3)*cos(la*x4)+x1*cosh(la*x3)*sin(la*x4))*sinh(lb*x3)*sin(lb*x4))*i4(cosh(la*x3)*sin(la*x4)*sinh(lc*x3)*cos(lc*x4)+sinh(la*x3)*cos(la*x4)*cosh(lc*x3)*sin( lc*x4))*v5-(sinh(la*x3)*cos(la*x4)*sinh(lc*x3)*cos(lc*x4)cosh(la*x3)*sin(la*x4)*cosh(lc*x3)*sin(lc*x4))*v6((x2*sinh(la*x3)*cos(la*x4)+x1*cosh(la*x3)*sin(la*x4))*cosh(lc*x3)*cos(lc*x4)+(x1*sin h(la*x3)*cos(la*x4)-x2*cosh(la*x3)*sin(la*x4))*sinh(lc*x3)*sin(lc*x4))*i5126 ((x1*sinh(la*x3)*cos(la*x4)-x2*cosh(la*x3)*sin(la*x4))*cosh(lc*x3)*cos(lc*x4)(x2*sinh(la*x3)*cos(la*x4)+x1*cosh(la*x3)*sin(la*x4))*sinh(lc*x3)*sin(lc*x4))*i6 =0 f3(x1,x 2,x3,x4,x 5) = real(PT2) = f3 = v3+(cosh(lb*x3)*cos(lb*x4)*cosh(lc*x3)*cos(lc*x4)+sinh(lb*x3)*sin(lb*x4)*sinh(lc*x3)*sin(l c*x4)sinh(lb*x3)*cos(lb*x4)*sinh(lc*x3)*cos(lc*x4)+cosh(lb*x3)*sin(lb*x4)*cosh(lc*x3)*sin(l c*x4))*v5-(-sinh(lb*x3)*sin(lb*x4)*cosh(lc*x3)*cos(lc*x4)cosh(lb*x3)*cos(lb*x4)*sinh(lc*x3)*sin(lc*x4)cosh(lb*x3)*sin(lb*x4)*sinh(lc*x3)*cos(lc*x4)sinh(lb*x3)*cos(lb*x4)*cosh(lc*x3)*sin(lc*x4))*v6+(-(cosh(lb*x3)*cos(lb*x4)*x1sinh(lb*x3)*sin(lb*x4)*x2)*sinh(lc*x3)*cos(lc*x4)+(sinh(lb*x3)*sin(lb*x4)*x1+cosh(lb* x3)*cos(lb*x4)*x2)*cosh(lc*x3)*sin(lc*x4)-(x1*sinh(lb*x3)*cos(lb*x4)x2*cosh(lb*x3)*sin(lb*x4))*cosh(lc*x3)*cos(lc*x4)+(x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh (lb*x3)*sin(lb*x4))*sinh(lc*x3)*sin(lc*x4))*i5-((sinh(lb*x3)*sin(lb*x4)*x1+cosh(lb*x3)*cos(lb*x4)*x2)*sinh(lc*x3)*cos(lc*x4)(cosh(lb*x3)*cos(lb*x4)*x1-sinh(lb*x3)*sin(lb*x4)*x2)*cosh(lc*x3)*sin(lc*x4)(x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh(lb*x3)*sin(lb*x4))*cosh(lc*x3)*cos(lc*x4)(x1*sinh(lb*x3)*cos(lb*x4)-x2*cosh(lb*x3)*sin(lb*x4))*sinh(lc*x3)*sin(lc*x4))*i6+(sinh(la*x3)*cos(la*x4)*sinh(lb*x3)*cos(lb*x4)+cosh(la*x3)*sin(la*x4)*cosh(lb*x3)*sin(l b*x4))*v1-(-sinh(la*x3)*cos(la*x4)*cosh(lb*x3)*sin(lb*x4)cosh(la*x3)*sin(la*x4)*sinh(lb*x3)*cos(lb*x4))*v2+((x1*sinh(lb*x3)*cos(lb*x4)x2*cosh(lb*x3)*sin(lb*x4))*cosh(la*x3)*cos(la*x4)(x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh(lb*x3)*sin(lb*x4))*sinh(la*x3)*sin(la*x4))*i1((x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh(lb*x3)*sin(lb*x4))*cosh(la*x3)*cos(la*x4)+(x1*si nh(lb*x3)*cos(lb*x4)-x2*cosh(lb*x3)*sin(lb*x4))*sinh(la*x3)*sin(la*x4))*i2 =0 f4(x1,x 2,x3,x4,x 5) = imag(PT2) = = v4+(-sinh(lb*x3)*sin(lb*x4)*cosh(lc*x3)*cos(lc*x4)f4 cosh(lb*x3)*cos(lb*x4)*sinh(lc*x3)*sin(lc*x4)cosh(lb*x3)*sin(lb*x4)*sinh(lc*x3)*cos(lc*x4)sinh(lb*x3)*cos(lb*x4)*cosh(lc*x3)*sin(lc*x4))*v5+(cosh(lb*x3)*cos(lb*x4)*cosh(lc*x3)*cos(lc*x4)+sinh(lb*x3)*sin(lb*x4)*sinh(lc*x3)*sin(l c*x4)sinh(lb*x3)*cos(lb*x4)*sinh(lc*x3)*cos(lc*x4)+cosh(lb*x3)*sin(lb*x4)*cosh(lc*x3)*sin(l c*x4))*v6+((sinh(lb*x3)*sin(lb*x4)*x1+cosh(lb*x3)*cos(lb*x4)*x2)*sinh(lc*x3)*cos(lc*x4)(cosh(lb*x3)*cos(lb*x4)*x1-sinh(lb*x3)*sin(lb*x4)*x2)*cosh(lc*x3)*sin(lc*x4)(x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh(lb*x3)*sin(lb*x4))*cosh(lc*x3)*cos(lc*x4)(x1*sinh(lb*x3)*cos(lb*x4)-x2*cosh(lb*x3)*sin(lb*x4))*sinh(lc*x3)*sin(lc*x4))*i5+((cosh(lb*x3)*cos(lb*x4)*x1sinh(lb*x3)*sin(lb*x4)*x2)*sinh(lc*x3)*cos(lc*x4)+(sinh(lb*x3)*sin(lb*x4)*x1+cosh(lb* x3)*cos(lb*x4)*x2)*cosh(lc*x3)*sin(lc*x4)-(x1*sinh(lb*x3)*cos(lb*x4)x2*cosh(lb*x3)*sin(lb*x4))*cosh(lc*x3)*cos(lc*x4)+(x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh (lb*x3)*sin(lb*x4))*sinh(lc*x3)*sin(lc*x4))*i6+(sinh(la*x3)*cos(la*x4)*cosh(lb*x3)*sin(lb*x4)cosh(la*x3)*sin(la*x4)*sinh(lb*x3)*cos(lb*x4))*v1+(127 sinh(la*x3)*cos(la*x4)*sinh(lb*x3)*cos(lb*x4)+cosh(la*x3)*sin(la*x4)*cosh(lb*x3)*sin(l b*x4))*v2+((x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh(lb*x3)*sin(lb*x4))*cosh(la*x3)*cos(la* x4)+(x1*sinh(lb*x3)*cos(lb*x4)x2*cosh(lb*x3)*sin(lb*x4))*sinh(la*x3)*sin(la*x4))*i1+((x1*sinh(lb*x3)*cos(lb*x4)x2*cosh(lb*x3)*sin(lb*x4))*cosh(la*x3)*cos(la*x4)(x2*sinh(lb*x3)*cos(lb*x4)+x1*cosh(lb*x3)*sin(lb*x4))*sinh(la*x3)*sin(la*x4))*i2=0 f5(x1,x 2,x3,x4,x 5) = real(PT4) = f5 = cosh(x5*x3)*cos(x5*x4)*v7-sinh(x5*x3)*sin(x5*x4)*v8+(sinh(x5*x3)*cos(x5*x4)*x1+cosh(x5*x3)*sin(x5*x4)*x2)*i7-(cosh(x5*x3)*sin(x5*x4)*x1-sinh(x5*x3)*cos(x5*x4)*x2)*i8+(-cosh(la*x3x3*x5)*cos(la*x4-x4*x5)*cosh(lb*x3)*cos(lb*x4)+sinh(la*x3-x3*x5)*sin(la*x4x4*x5)*sinh(lb*x3)*sin(lb*x4)-sinh(la*x3-x3*x5)*cos(la*x4x4*x5)*sinh(lb*x3)*cos(lb*x4)+cosh(la*x3-x3*x5)*sin(la*x4x4*x5)*cosh(lb*x3)*sin(lb*x4))*v9-(-sinh(la*x3-x3*x5)*sin(la*x4x4*x5)*cosh(lb*x3)*cos(lb*x4)-cosh(la*x3-x3*x5)*cos(la*x4x4*x5)*sinh(lb*x3)*sin(lb*x4)-cosh(la*x3-x3*x5)*sin(la*x4x4*x5)*sinh(lb*x3)*cos(lb*x4)-sinh(la*x3-x3*x5)*cos(la*x4x4*x5)*cosh(lb*x3)*sin(lb*x4))*v10+((cosh(la*x3-x3*x5)*cos(la*x4-x4*x5)*x1sinh(la*x3-x3*x5)*sin(la*x4-x4*x5)*x2)*sinh(lb*x3)*cos(lb*x4)-(sinh(la*x3x3*x5)*sin(la*x4-x4*x5)*x1+cosh(la*x3-x3*x5)*cos(la*x4x4*x5)*x2)*cosh(lb*x3)*sin(lb*x4)+(sinh(la*x3-x3*x5)*cos(la*x4-x4*x5)*x1-cosh(la*x3x3*x5)*sin(la*x4-x4*x5)*x2)*cosh(lb*x3)*cos(lb*x4)-(cosh(la*x3-x3*x5)*sin(la*x4x4*x5)*x1+sinh(la*x3-x3*x5)*cos(la*x4-x4*x5)*x2)*sinh(lb*x3)*sin(lb*x4))*i9((sinh(la*x3-x3*x5)*sin(la*x4-x4*x5)*x1+cosh(la*x3-x3*x5)*cos(la*x4x4*x5)*x2)*sinh(lb*x3)*cos(lb*x4)+(cosh(la*x3-x3*x5)*cos(la*x4-x4*x5)*x1-sinh(la*x3x3*x5)*sin(la*x4-x4*x5)*x2)*cosh(lb*x3)*sin(lb*x4)+(cosh(la*x3-x3*x5)*sin(la*x4x4*x5)*x1+sinh(la*x3-x3*x5)*cos(la*x4x4*x5)*x2)*cosh(lb*x3)*cos(lb*x4)+(sinh(la*x3-x3*x5)*cos(la*x4-x4*x5)*x1cosh(la*x3-x3*x5)*sin(la*x4-x4*x5)*x2)*sinh(lb*x3)*sin(lb*x4))*i10+(-sinh(la*x3x3*x5)*cos(la*x4-x4*x5)*sinh(lc*x3)*cos(lc*x4)+cosh(la*x3-x3*x5)*sin(la*x4x4*x5)*cosh(lc*x3)*sin(lc*x4))*v11-(-cosh(la*x3-x3*x5)*sin(la*x4x4*x5)*sinh(lc*x3)*cos(lc*x4)-sinh(la*x3-x3*x5)*cos(la*x4x4*x5)*cosh(lc*x3)*sin(lc*x4))*v12+((sinh(la*x3-x3*x5)*cos(la*x4-x4*x5)*x1cosh(la*x3-x3*x5)*sin(la*x4-x4*x5)*x2)*cosh(lc*x3)*cos(lc*x4)-(cosh(la*x3x3*x5)*sin(la*x4-x4*x5)*x1+sinh(la*x3-x3*x5)*cos(la*x4x4*x5)*x2)*sinh(lc*x3)*sin(lc*x4))*i11-((cosh(la*x3-x3*x5)*sin(la*x4x4*x5)*x1+sinh(la*x3-x3*x5)*cos(la*x4x4*x5)*x2)*cosh(lc*x3)*cos(lc*x4)+(sinh(la*x3-x3*x5)*cos(la*x4-x4*x5)*x1-cosh(la*x3x3*x5)*sin(la*x4-x4*x5)*x2)*sinh(lc*x3)*sin(lc*x4))*i12 =0 A.9 Phương trình (3.112) f1, f2, f 3, f4 tương tự phụ lục A.9 f5(x1,x 2,x3,x4,x 5) = real(PT4) = 128 f5 = cosh(x5*x3)*cos(x5*x4)*v9-sinh(x5*x3)*sin(x5*x4)*v10(sinh(x5*x3)*cos(x5*x4)*x1cosh(x5*x3)*sin(x5*x4)*x2)*i9+(cosh(x5*x3)*sin(x5*x4)*x1+sinh(x5*x3)*cos(x5*x4)*x2 )*i10-(cosh(lb*x3-x3*x5)*cos(lb*x4-x4*x5)*cosh(lc*x3)*cos(lc*x4)-sinh(lb*x3x3*x5)*sin(lb*x4-x4*x5)*sinh(lc*x3)*sin(lc*x4)+sinh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*sinh(lc*x3)*cos(lc*x4)-cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*cosh(lc*x3)*sin(lc*x4))*v11+(sinh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*cosh(lc*x3)*cos(lc*x4)+cosh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*sinh(lc*x3)*sin(lc*x4)+cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*sinh(lc*x3)*cos(lc*x4)+sinh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*cosh(lc*x3)*sin(lc*x4))*v12-(-(cosh(lb*x3-x3*x5)*cos(lb*x4-x4*x5)*x1sinh(lb*x3-x3*x5)*sin(lb*x4-x4*x5)*x2)*sinh(lc*x3)*cos(lc*x4)+(sinh(lb*x3x3*x5)*sin(lb*x4-x4*x5)*x1+cosh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*x2)*cosh(lc*x3)*sin(lc*x4)-(sinh(lb*x3-x3*x5)*cos(lb*x4-x4*x5)*x1-cosh(lb*x3x3*x5)*sin(lb*x4-x4*x5)*x2)*cosh(lc*x3)*cos(lc*x4)+(cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*x1+sinh(lb*x3-x3*x5)*cos(lb*x4-x4*x5)*x2)*sinh(lc*x3)*sin(lc*x4))*i11+((sinh(lb*x3-x3*x5)*sin(lb*x4-x4*x5)*x1+cosh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*x2)*sinh(lc*x3)*cos(lc*x4)-(cosh(lb*x3-x3*x5)*cos(lb*x4-x4*x5)*x1-sinh(lb*x3x3*x5)*sin(lb*x4-x4*x5)*x2)*cosh(lc*x3)*sin(lc*x4)-(cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*x1+sinh(lb*x3-x3*x5)*cos(lb*x4-x4*x5)*x2)*cosh(lc*x3)*cos(lc*x4)-(sinh(lb*x3x3*x5)*cos(lb*x4-x4*x5)*x1-cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*x2)*sinh(lc*x3)*sin(lc*x4))*i12-(sinh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*sinh(la*x3)*cos(la*x4)-cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*cosh(la*x3)*sin(la*x4))*v7+(cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*sinh(la*x3)*cos(la*x4)+sinh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*cosh(la*x3)*sin(la*x4))*v8-((-sinh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*x1+cosh(lb*x3-x3*x5)*sin(lb*x4-x4*x5)*x2)*cosh(la*x3)*cos(la*x4)-(cosh(lb*x3-x3*x5)*sin(lb*x4-x4*x5)*x1-sinh(lb*x3-x3*x5)*cos(lb*x4x4*x5)*x2)*sinh(la*x3)*sin(la*x4))*i7+((-cosh(lb*x3-x3*x5)*sin(lb*x4-x4*x5)*x1sinh(lb*x3-x3*x5)*cos(lb*x4-x4*x5)*x2)*cosh(la*x3)*cos(la*x4)+(-sinh(lb*x3x3*x5)*cos(lb*x4-x4*x5)*x1+cosh(lb*x3-x3*x5)*sin(lb*x4x4*x5)*x2)*sinh(la*x3)*sin(la*x4))*i8=0 A.10 Phương trình (3.113) f1, f2, f 3, f4 tương tự phụ lục A.9 f5(x1,x 2,x3,x4,x 5) = real(PT4) = = cosh(x5*x3)*cos(x5*x4)*v11-sinh(x5*x3)*sin(x5*x4)*v12+(f5 sinh(x5*x3)*cos(x5*x4)*x1+cosh(x5*x3)*sin(x5*x4)*x2)*i11-(cosh(x5*x3)*sin(x5*x4)*x1-sinh(x5*x3)*cos(x5*x4)*x2)*i12+(cosh(la*x3)*cos(la*x4)*cosh(lc*x3-x3*x5)*cos(lc*x4x4*x5)+sinh(la*x3)*sin(la*x4)*sinh(lc*x3-x3*x5)*sin(lc*x4-x4*x5)-sinh(lc*x3x3*x5)*cos(lc*x4-x4*x5)*sinh(la*x3)*cos(la*x4)+cosh(lc*x3-x3*x5)*sin(lc*x4x4*x5)*cosh(la*x3)*sin(la*x4))*v7-(-sinh(la*x3)*sin(la*x4)*cosh(lc*x3x3*x5)*cos(lc*x4-x4*x5)-cosh(la*x3)*cos(la*x4)*sinh(lc*x3-x3*x5)*sin(lc*x4-x4*x5)129 cosh(lc*x3-x3*x5)*sin(lc*x4-x4*x5)*sinh(la*x3)*cos(la*x4)-sinh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*cosh(la*x3)*sin(la*x4))*v8+((cosh(lc*x3-x3*x5)*cos(lc*x4-x4*x5)*x1sinh(lc*x3-x3*x5)*sin(lc*x4-x4*x5)*x2)*sinh(la*x3)*cos(la*x4)-(sinh(lc*x3x3*x5)*sin(lc*x4-x4*x5)*x1+cosh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*x2)*cosh(la*x3)*sin(la*x4)+(sinh(lc*x3-x3*x5)*cos(lc*x4-x4*x5)*x1-cosh(lc*x3x3*x5)*sin(lc*x4-x4*x5)*x2)*cosh(la*x3)*cos(la*x4)-(cosh(lc*x3-x3*x5)*sin(lc*x4x4*x5)*x1+sinh(lc*x3-x3*x5)*cos(lc*x4-x4*x5)*x2)*sinh(la*x3)*sin(la*x4))*i7((sinh(lc*x3-x3*x5)*sin(lc*x4-x4*x5)*x1+cosh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*x2)*sinh(la*x3)*cos(la*x4)+(cosh(lc*x3-x3*x5)*cos(lc*x4-x4*x5)*x1-sinh(lc*x3x3*x5)*sin(lc*x4-x4*x5)*x2)*cosh(la*x3)*sin(la*x4)+(cosh(lc*x3-x3*x5)*sin(lc*x4x4*x5)*x1+sinh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*x2)*cosh(la*x3)*cos(la*x4)+(sinh(lc*x3-x3*x5)*cos(lc*x4-x4*x5)*x1-cosh(lc*x3x3*x5)*sin(lc*x4-x4*x5)*x2)*sinh(la*x3)*sin(la*x4))*i8+(-sinh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*sinh(lb*x3)*cos(lb*x4)+cosh(lc*x3-x3*x5)*sin(lc*x4x4*x5)*cosh(lb*x3)*sin(lb*x4))*v9-(-cosh(lc*x3-x3*x5)*sin(lc*x4x4*x5)*sinh(lb*x3)*cos(lb*x4)-sinh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*cosh(lb*x3)*sin(lb*x4))*v10+((sinh(lc*x3-x3*x5)*cos(lc*x4-x4*x5)*x1cosh(lc*x3-x3*x5)*sin(lc*x4-x4*x5)*x2)*cosh(lb*x3)*cos(lb*x4)-(cosh(lc*x3x3*x5)*sin(lc*x4-x4*x5)*x1+sinh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*x2)*sinh(lb*x3)*sin(lb*x4))*i9-((cosh(lc*x3-x3*x5)*sin(lc*x4x4*x5)*x1+sinh(lc*x3-x3*x5)*cos(lc*x4x4*x5)*x2)*cosh(lb*x3)*cos(lb*x4)+(sinh(lc*x3-x3*x5)*cos(lc*x4-x4*x5)*x1-cosh(lc*x3x3*x5)*sin(lc*x4-x4*x5)*x2)*sinh(lb*x3)*sin(lb*x4))*i10=0 130 ... đầu đường dây xác định thơng số đường dây; xác định vị trí cố xảy đường dây truyền tải điện đồng Nghiên cứu [15] giới thiệu thuật toán định vị cố mà không yêu cầu sử dụng thông số đường dây Thuật... chọn luận án “ Nghiên cứu phương pháp xác định vị trí cố đường dây truyền tải điện thuộc lưới điện phức tạp” cần thiết có ý nghĩa quan trọng cơng ty vận hành, quản lý lưới điện truyền tải Mục tiêu,... XUÂN VINH NGHIÊN CỨU PHƯƠNG PHÁP XÁC ĐỊNH VỊ TRÍ SỰ CỐ TRÊN ĐƯỜNG DÂY TRUYỀN TẢI ĐIỆN THUỘC LƯỚI ĐIỆN PHỨC TẠP Chuyên ngành: Kỹ thuật điện Mã số: 62520202 LUẬN ÁN TIẾN SĨ KỸ THUẬT ĐIỆN NGƯỜI