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

139 264 0
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

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

1. Tính cấp thiết của đề tài Việc xác định chính xác điểm sự cố trên đƣờng dây truyền tải điện mang một ý nghĩa thiết thực đối với hệ thống truyền tải điện; định vị chính xác điểm sự cố sẽ làm giảm thời gian ngừng cấp điện với các sự cố duy trì, giảm nhân công huy động tìm kiếm vị trí sự cố và góp phần nâng cao độ tin cậy của lƣới điện truyền tải. a. Qui trình tìm kiếm vị trí sự cố xảy ra trên đường dây truyền tải đang được áp dụng tại các công ty truyền tải điện Việt Nam  Khi có sự cố trên đƣờng dây, hệ thống rơle bảo vệ tại trạm biến áp sẽ cắt điện đƣờng dây đó và thực hiện tự đóng lại; đồng thời chức năng báo khoảng cách của các rơle bảo vệ khoảng cách hoặc bảo vệ so lệch sẽ báo vị trí sự cố tính theo km từ vị trí đặt điểm đo của rơle.  Sự cố có thể là thoáng qua hoặc duy trì, tuy nhiên theo qui định hiện hành các đội quản lý đƣờng dây đều phải tìm kiếm và xác định vị trí cũng nhƣ loại sự cố để có các giải pháp tránh sự cố tƣơng tự có thể lặp lại hoặc có thể nhanh chóng thay thế các trang thiết bị hỏng hóc.  Thông thƣờng các đơn vị quản lý đƣờng dây cần phải cử nhân công đến vị trí cột tƣơng ứng với khoảng cách sự cố rơle đã báo và chia thành hai nhóm để tìm kiếm về hai phía. Tại các vị trí cột đều cần phải kiểm tra bằng mắt thƣờng xem có hƣ hỏng hay bất thƣờng. Quá trình tìm kiếm tiếp diễn cho đến khi nào tìm đƣợc vị trí sự cố thực thì sẽ chụp ảnh ghi nhận và báo cáo. b. Các khó khăn còn tồn tại  Chức năng báo khoảng cách sự cố của các rơle bảo vệ khoảng cách thƣờng báo vị trí sự cố với sai số tƣơng đối lớn. Sai số của vị trí sự cố trong nhiều trƣờng hợp có thể lên tới nhiều kilômét. Số liệu thống kê ở một số nƣớc cho thấy cho thấy sai số gặp phải có thể từ 0,5÷2% [57], với đƣờng dây dài 300km thì sai số ±1% tƣơng đƣơng với việc phải đi tìm kiếm trong phạm vi 6 kilômét (khoảng 20 khoảng cột).  Các đƣờng dây truyền tải điện thƣờng đi qua các địa hình đồi núi hoặc xa dân cƣ, xa đƣờng giao thông (lƣới điện truyền tải khu vực phía Tây Bắc). Do vậy rất nhiều trƣờng hợp phải mất tới cả ngày để tìm chính xác một vị trí sự cố với lƣợng nhân công từ 4÷6 ngƣời. c. Các nguyên nhân kỹ thuật có thể liên quan tới việc định vị sự cố không chính xác Phân tích các yếu tố đầu vào của các thuật toán định xác định điểm sự cố cho thấy sai số về đo lƣờng và mô hình đƣờng dây sử dụng trong các thuật toán là các yếu tố chính ảnh hƣởng tới sai số của kết quả định vị:  Các thuật toán định vị sự cố hiện nay đều yêu cầu dữ liệu đầu vào là tín hiệu điện áp và d ng điện đo lƣờng từ các đầu đƣờng dây truyền tải điện, vì vậy khi có sai số đo lƣờng thì kết quả định vị sự cố sẽ gặp phải sai số. Trong đó, sai số đo lƣờng bao gồm:  Sai số của thiết bị biến đổi: theo tiêu chuẩn IEC60044_1_1996 về sai số của máy biến điện áp (VT) và IEC60044_2_1997 cho máy biến d ng điện (CT), thiết bị biến đổi có thể có sai số lớn nhất ±10% đối với CT và ±6% đối với VT. Do đó sai số của các thiết bị này có ảnh hƣởng đáng kể đến độ chính xác của phƣơng pháp định vị, độ không đồng nhất của sai số cũng là yếu tố có thể có ảnh hƣởng lớn đến kết quả định vị sự cố. Tuy nhiên, sai số của thiết bị biến đổi

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 MỤC LỤC Tr ang 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 vii DANH MỤC CÁC BẢNG viii DANH MỤC CÁC HÌNH ẢNH, ĐỒ THỊ ix 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 1.1 Tình hình nghiên cứu nƣớc 1.2 Tình hình nghiên cứu nƣớc 1.2.1 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 1.2.2 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 1.2.3 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 1.2.4 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 1.2.5 Đánh giá ảnh hƣởng sai số đo lƣờng đến kết định vị cố 16 1.3 Những vấn đề c 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 đề 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 2.2.1 Mạng hai cửa tƣơng đƣơng đƣờng dây không đồng 18 2.2.2 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 27 2.3.1 Bƣớc 1: Giả thiết cố xảy phân đoạn SC 27 iii 2.3.2 Bƣớc 2: Giả thiết cố xảy phân đoạn RC 28 2.3.3 Bƣớc 3: Xác định vị trí cố 29 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 30 2.4.1 Kết định vị cố 31 2.4.2 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] 31 2.5 Kết luận 33 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 34 3.1 Đặt vấn đề 34 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 34 3.2.1 Mạng hai cửa tƣơng đƣơng đƣờng dây rẽ nhánh 34 3.2.2 Trƣờng hợp phân đoạn đƣờng dây rẽ nhánh có thông số đƣờng dây khác 35 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 50 3.3.1 Bƣớc 1: Giả thiết cố xảy phân đoạn SJ 50 3.3.2 Bƣớc 2: Giả thiết cố xảy phân đoạn RJ 50 3.3.3 Bƣớc 3: Giả thiết cố xảy phân đoạn TJ 51 3.3.4 Bƣớc 4: Xác định vị trí cố 51 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ố 52 3.4.1 Bƣớc 1: Giả thiết cố xảy phân đoạn SJ 52 3.4.2 Bƣớc 2: Giả thiết cố xảy phân đoạn RJ 52 3.4.3 Bƣớc 3: Giả thiết cố xảy phân đoạn TJ 53 3.4.4 Bƣớc 4: Xác định thông số đƣờng dây vị trí cố 53 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 54 3.5.1 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 55 3.5.2 Trƣờng hợp đƣờng dây rẽ nhánh có hai phân đoạn có thông số đƣờng dây nhƣng khác thông số phân đoạn thứ ba 57 3.5.3 Trƣờng hợp phân đoạn đƣờng dây rẽ nhánh có thông số đƣờng dây giống 59 3.6 Kết luận 61 ĐÁNH GIÁ ẢNH HƢỞNG CỦA SAI SỐ ĐO LƢỜNG ĐẾN KẾT QUẢ iv ĐỊNH VỊ SỰ CỐ 63 4.1 Cơ sở đề xuất phƣơng pháp 63 4.2 Mô Monte Carlo 64 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 65 4.3.1 Kịch mô 66 4.3.2 Đặc tính thống kê 66 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 70 4.4.1 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 70 4.4.2 Trƣờng hợp hai phân đoạn có thông số đƣờng dây nhƣng khác thông số phân đoạn thứ ba 74 4.4.3 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 78 4.5 Kết luận 82 KẾT LUẬN VÀ KIẾN NGHỊ 84 Đóng góp khoa học luận án 84 Kiến nghị nghiên cứu 85 TÀI LIỆU THAM KHẢO 86 DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ CỦA LUẬN ÁN 91 PHỤ LỤC 92 A.1 Kết mô định vị cố xảy đƣờng dây truyền tải điện không đồng 92 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 nhƣng khác thông số đƣờng dây phân đoạn RJ 98 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 100 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 vi 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 vii 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 30 Bảng 2.2 Thông số phân đoạn CR, l2=20km 30 Bảng 2.3 Thông số nguồn 30 Bảng 2.4 Kết định vị cố xảy đƣờng dây không đồng 31 Bảng 3.1 Thông số nguồn 54 Bảng 3.2 Thông số phân đoạn SJ 54 Bảng 3.3 Thông số phân đoạn RJ 54 Bảng 3.4 Thông số phân đoạn TJ 55 Bảng 3.5 Kết định vị cố xảy đƣờng dây rẽ nhánh không đồng 55 Bảng 3.6 Kết định vị cố đƣờng dây có thông số SJ, TJ khác RJ 57 Bảng 3.7 Kết định vị cố xảy đƣờng dây rẽ nhánh đồng 59 Bảng 4.1 Sai số giới hạn máy biến d ng điện [34] 65 Bảng 4.2 Sai số giới hạn máy biến điện áp [35] 65 Bảng 4.3 Chỉ số thống kê loại cố AG – đƣờng dây không đồng 67 Bảng 4.4 Chỉ số thống kê loại cố BC 68 Bảng 4.5 Chỉ số thống kê loại cố BCG 69 Bảng 4.6 Chỉ số thống kê loại cố ABC 69 Bảng 4.7 Chỉ số thống kê loại cố AG – đƣờng dây rẽ nhánh không đồng 72 Bảng 4.8 Chỉ số thống kê loại cố BC 72 Bảng 4.9 Chỉ số thống kê loại cố BCG 73 Bảng 4.10 Chỉ số thống kê loại cố ABC 74 Bảng 4.11 Chỉ số thống kê loại cố AG – SJ, TJ khác thông số RJ 76 Bảng 4.12 Chỉ số thống kê loại cố BC 77 Bảng 4.13 Chỉ số thống kê loại cố BCG 77 Bảng 4.14 Chỉ số thống kê loại cố ABC 78 Bảng 4.15 Chỉ số thống kê loại cố AG – đƣờng dây rẽ nhánh đồng 80 Bảng 4.16 Chỉ số thống kê loại cố BC 81 Bảng 4.17 Chỉ số thống kê loại cố BCG 81 Bảng 4.18 Chỉ số thống kê loại cố ABC 82 viii 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 19 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 22 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 27 Hình 2.11 Sự cố xảy phân đoạn SC 27 Hình 2.12 Sự cố xảy phân đoạn RC 28 Hình 2.13 Thuật toán định vị cố xảy đƣờng dây không đồng 29 Hình 2.14 Mô hình mô đƣờng dây không đồng 30 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] 32 Hình 3.1 Đƣờng dây truyền tải rẽ nhánh 34 Hình 3.2 Mô hình mạng hai cửa đƣờng dây rẽ nhánh 35 Hình 3.3 Sự cố phân đoạn SJ 40 Hình 3.4 Mô hình mạng hai cửa cố phân đoạn SJ 40 Hình 3.5 Sự cố phân đoạn RJ 42 Hình 3.6 Mô hình mạng hai cửa cố phân đoạn RJ 42 Hình 3.7 Sự cố phân đoạn TJ 45 Hình 3.8 Mô hình mạng hai cửa cố phân đoạn TJ 45 Hình 3.9 Thuật toán xác định thông số đƣờng dây rẽ nhánh không đồng 48 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 49 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ánh 51 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 53 Hình 3.13 Mô hình mô đƣờng dây rẽ nhánh 54 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] 56 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] 58 ix Hình 3.16 So sánh kết định vị cố thuật toán đề xuất với thuật toán [47] 60 Hình 4.1 Thuật toán đánh giá ảnh hƣởng sai số đo lƣờng đến kết định vị cố 64 Hình 4.2 Hàm mật độ phân bố sai số định vị cố - đƣờng dây không đồng 66 Hình 4.3 Phân tích độ nhạy kết định vị cố với biến ngõ vào 67 Hình 4.4 Hàm mật độ phân bố sai số định vị cố ˗ SJ, RJ, TJ khác thông số 70 Hình 4.5 Phân tích độ nhạy kết định vị cố với biến ngõ vào 71 Hình 4.6 Hàm mật độ phân bố sai số định vị cố - SJ, TJ khác thông số RJ 75 Hình 4.7 Phân tích độ nhạy kết định vị cố với biến ngõ vào 75 Hình 4.8 Hàm mật độ phân bố sai số định vị cố - đƣờng dây rẽ nhánh đồng 79 Hình 4.9 Phân tích độ nhạy kết định vị cố với biến ngõ vào 79 x MỞ ĐẦU Tính cấp thiết đề tài Việc xác định xác điểm cố đƣờng dây truyền tải điện mang ý nghĩa thiết thực hệ thống truyền tải điện; định vị xác điểm cố làm giảm thời gian ngừng cấp điện với cố trì, giảm nhân công huy động tìm kiếm vị trí cố góp phần nâng cao độ tin cậy lƣới điện truyền tải a Qui trình tìm kiếm vị trí cố xảy đường dây truyền tải áp dụng công ty truyền tải điện Việt Nam  Khi có cố đƣờng dây, hệ thống rơle bảo vệ trạm biến áp cắt điện đƣờng dây thực tự đóng lại; đồng thời chức báo khoảng cách rơle bảo vệ khoảng cách bảo vệ so lệch báo vị trí cố tính theo km từ vị trí đặt điểm đo rơle  Sự cố thoáng qua trì, nhiên theo qui định hành đội quản lý đƣờng dây phải tìm kiếm xác định vị trí nhƣ loại cố để có giải pháp tránh cố tƣơng tự lặp lại nhanh chóng thay trang thiết bị hỏng hóc  Thông thƣờng đơn vị quản lý đƣờng dây cần phải cử nhân công đến vị trí cột tƣơng ứng với khoảng cách cố rơle báo chia thành hai nhóm để tìm kiếm hai phía Tại vị trí cột cần phải kiểm tra mắt thƣờng xem có hƣ hỏng hay bất thƣờng Quá trình tìm kiếm tiếp diễn tìm đƣợc vị trí cố thực chụp ảnh ghi nhận báo cáo b Các khó khăn tồn  Chức báo khoảng cách cố rơle bảo vệ khoảng cách thƣờng báo vị trí cố với sai số tƣơng đối lớn Sai số vị trí cố nhiều trƣờng hợp lên tới nhiều kilômét Số liệu thống kê số nƣớc cho thấy cho thấy sai số gặp phải từ 0,5÷2% [57], với đƣờng dây dài 300km sai số ±1% tƣơng đƣơng với việc phải tìm kiếm phạm vi kilômét (khoảng 20 khoảng cột)  Các đƣờng dây truyền tải điện thƣờng qua địa hình đồi núi xa dân cƣ, xa đƣờng giao thông (lƣới điện truyền tải khu vực phía Tây Bắc) Do nhiều trƣờng hợp phải tới ngày để tìm xác vị trí cố với lƣợng nhân công từ 4÷6 ngƣời c Các nguyên nhân kỹ thuật liên quan tới việc định vị cố không xác Phân tích yếu tố đầu vào thuật toán định xác định điểm cố cho thấy sai số đo lƣờng mô hình đƣờng dây sử dụng thuật toán yếu tố ảnh hƣởng tới sai số kết định vị:  Các thuật toán định vị cố yêu cầu liệu đầu vào tín hiệu điện áp d ng điện đo lƣờng từ đầu đƣờng dây truyền tải điện, có sai số đo lƣờng kết định vị cố gặp phải sai số Trong đó, sai số đo lƣờng bao gồm:  Sai số thiết bị biến đổi: theo tiêu chuẩn IEC60044_1_1996 sai số máy biến điện áp (VT) IEC60044_2_1997 cho máy biến d ng điện (CT), thiết bị biến đổi có sai số lớn ±10% CT ±6% VT Do sai số thiết bị có ảnh hƣởng đáng kể đến độ xác phƣơng pháp định vị, độ không đồng sai số yếu tố có ảnh hƣởng lớn đến kết định vị cố Tuy nhiên, sai số thiết bị biến đổi 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)*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,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*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,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 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,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)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,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))*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,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)*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,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*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,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*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,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*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,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))*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,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+(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,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ự nhƣ phụ lục A.9 f5(x1,x2,x3,x4,x5) = 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, f3, f4 tƣơng tự nhƣ 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)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

Ngày đăng: 05/09/2017, 08:58

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan