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 29 tháng 08 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 đỡ 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 hoà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 i LỜI CẢM ƠN ii DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT vi DANH MỤC CÁC BẢNG vii 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 nước .7 Thuật toán định vị cố áp dụng đường dây đơn đồng với yêu cầu biết xác thông số đường dây Thuật toán định vị cố không yê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 yê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 .16 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 đề .18 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 toán xác định thông số phân đoạn thuộc đường dây đơn không đồng .20 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 .26 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 toá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 32 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 đề .33 3.2 Thuật toá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 34 3.3 Thuật toá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 toá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ố .51 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 toán định vị cố xảy đường dây truyền tải điện rẽ nhánh 53 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 57 3.6 Kết luận 60 ĐÁNH GIÁ ẢNH HƯỞNG CỦA SAI SỐ ĐO LƯỜNG ĐẾN KẾT QUẢ ĐỊNH VỊ SỰ CỐ .62 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 toá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 69 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 .77 4.5 Kết luận 81 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 90 PHỤ LỤC 91 A.1 Kết mô định vị cố xảy đường dây truyền tải điện không đồng 91 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 95 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 97 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 99 A.5 Phương trình (2.48) 101 A.6 Phương trình (3.103) 105 A.7 Phương trình (3.105) 117 A.8 Phương trình (3.111) 125 A.9 Phương trình (3.112) 127 A.10 Phương trình (3.113) 128 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 , Isf , ΔIsf 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 .8 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 toá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 toá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 toán xác định thông số đường dây rẽ nhánh không đồng 47 Hình 3.10 Thuật toá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 toá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 toá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,x2,x3,x4,x5,x6,x7,x8,x9,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)*x5115 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,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12) = 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*x11116 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,x2,x3,x4,x5,x6,x7,x8) = 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 117 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,x2,x3,x4,x5,x6,x7,x8) = 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)118 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,x2,x3,x4,x5,x6,x7,x8) = 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))*i1119 ((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,x2,x3,x4,x5,x6,x7,x8) = 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,x2,x3,x4,x5,x6,x7,x8) = 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)*x1120 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,x2,x3,x4,x5,x6,x7,x8) = 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*x6121 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,x2,x3,x4,x5,x6,x7,x8) = 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*x8122 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,x2,x3,x4,x5,x6,x7,x8) = 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*x8123 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 124 A.8 Phương trình (3.111) f1(x1,x2,x3,x4,x5) = real(PT1) = f1 = v1-(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))*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,x2,x3,x4,x5) = imag(PT1) = f2 = v2(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))*i5125 ((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,x2,x3,x4,x5) = 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,x2,x3,x4,x5) = imag(PT2) = f4 = v4+(-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))*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+(126 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,x2,x3,x4,x5) = 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, f3, f4 tương tự phụ lục A.9 f5(x1,x2,x3,x4,x5) = real(PT4) = 127 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, f3, f4 tương tự phụ lục A.9 f5(x1,x2,x3,x4,x5) = real(PT4) = f5 = cosh(x5*x3)*cos(x5*x4)*v11-sinh(x5*x3)*sin(x5*x4)*v12+(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)128 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 129 ... 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... 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,... 01: cố xảy TJ, Dt vị trí cố xác Nghiên cứu [10], [32] đề xuất thuật toán cho toán định vị cố xảy đường dây truyền