luận án tiến sĩ kỹ thuật 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

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 q trình nghiên cứu hồn thiện luận án Tác giả luận án Nguyễn Xuân Vinh ii MỤC LỤC Trang LỜI CAM ĐOAN 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 ngồi nước .7 Thuật toán định vị cố áp dụng đường dây đơn đồng với u cầu biết xác thơng số đường dây Thuật tốn định vị cố khơng 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 u cầu biết xác thơng số đường dây 12 Thuật toán định vị cố áp dụng cho đường dây truyền tải điện rẽ nhánh với u cầu biết xác thơng số đường dây 14 Đánh giá ảnh hưởng sai số đo lường đến kết định vị cố .16 1.3 Những vấn đề tồn hướng nghiên cứu .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 tốn xác định thơng số phân đoạn thuộc đường dây khơng đồng sử dụng tín hiệu đo lường từ hai phía đầu đường dây 18 Mạng hai cửa tương đương đường dây khơng đồng .18 Thuật tốn xác định thông số phân đoạn thuộc đường dây đơn không đồng .20 2.3 Thuật tố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 tốn xác định thơng số đường dây đường dây rẽ nhánh không đồng sử dụng tín hiệu đo lường từ ba đầu đường dây 33 Mạng hai cửa tương đương đường dây rẽ nhánh 33 Trường hợp phân đoạn đường dây rẽ nhánh có thơng số đường dây khác 34 3.3 Thuật tốn xác định vị trí cố cho đường dây rẽ nhánh không đồng sử dụng tín hiệu đo lường từ ba đầu đường dây 49 Bước 1: Giả thiết cố xảy phân đoạn SJ 49 Bước 2: Giả thiết cố xảy phân đoạn RJ .49 Bước 3: Giả thiết cố xảy phân đoạn TJ .49 Bước 4: Xác định vị trí cố 50 3.4 Thuật tốn xác định thơng số đường dây đồng thời định vị cố áp dụng với đường dây rẽ nhánh, phân đoạn có chủng loại dây, sử dụng ghi cố .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 tố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 tốn xác định vị trí cố xảy đường dây truyền tải điện không đồng 64 Kịch mô .65 Đặc tính thống kê 65 4.4 Đánh giá ảnh hưởng sai số đo lường đến kết định vị cố thuật toán định vị cố xảy đường dây truyền tải điện rẽ nhánh 69 Trường hợp phân đoạn thuộc đường dây rẽ nhánh có thơng số đường dây khác 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 tốn xác định thơng số đường dây không đồng 26 Hình 2.11 Sự cố xảy phân đoạn SC .26 Hình 2.12 Sự cố xảy phân đoạn RC 27 Hình 2.13 Thuật tốn định vị cố xảy đường dây không đồng 28 Hình 2.14 Mơ hình mơ đường dây khơng đồng 29 Hình 2.15 So sánh kết định vị cố thuật tốn đề xuất với thuật tốn [48] 31 Hình 3.1 Đường dây truyền tải rẽ nhánh 33 Hình 3.2 Mơ hình mạng hai cửa đường dây rẽ nhánh 34 Hình 3.3 Sự cố phân đoạn SJ 39 Hình 3.4 Mơ hình mạng hai cửa cố phân đoạn SJ .39 Hình 3.5 Sự cố phân đoạn RJ 41 Hình 3.6 Mơ hình mạng hai cửa cố phân đoạn RJ 41 Hình 3.7 Sự cố phân đoạn TJ .44 Hình 3.8 Mơ hình mạng hai cửa cố phân đoạn TJ 44 Hình 3.9 Thuật tốn xác định thông số đường dây rẽ nhánh không đồng 47 Hình 3.10 Thuật tốn xác định thơng số đường dây rẽ nhánh – SJ, TJ thông số khác RJ 48 Hình 3.11 Thuật tốn xác định vị trí cố xảy đường dây truyền tải điện rẽ nhánh50 Hình 3.12 Thuật tốn định vị cố áp dụng với đường dây truyền tải rẽ nhánh đồng nhất52 Hình 3.13 Mơ hình mơ đường dây rẽ nhánh 53 Hình 3.14 So sánh kết định vị cố thuật toán đề xuất với thuật toán [47] 55 Hình 3.15 So sánh kết định vị cố thuật toán đề xuất với thuật toán [47] 57 viii Hình 3.16 So sánh kết định vị cố thuật toán đề xuất với thuật tốn [47] 59 Hình 4.1 Thuật tốn đánh giá ảnh hưởng sai số đo lường đến kết định vị cố 63 Hình 4.2 Hàm mật độ phân bố sai số định vị cố - đường dây khơng đồng 65 Hình 4.3 Phân tích độ nhạy kết định vị cố với biến ngõ vào 66 Hình 4.4 Hàm mật độ phân bố sai số định vị cố ˗ SJ, RJ, TJ khác thông số 69 Hình 4.5 Phân tích độ nhạy kết định vị cố với biến ngõ vào .70 Hình 4.6 Hàm mật độ phân bố sai số định vị cố - SJ, TJ khác thông số RJ 74 Hình 4.7 Phân tích độ nhạy kết định vị cố với biến ngõ vào .74 Hình 4.8 Hàm mật độ phân bố sai số định vị cố - đường dây rẽ nhánh đồng 78 Hình 4.9 Phân tích độ nhạy kết định vị cố với biến ngõ vào 78 ix 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ố thố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í 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 vị trí cố thực chụp ảnh ghi nhận báo cáo b Các khó khăn cị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 ngun 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 tố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 có tính chất 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