Đăng ký
  1. Trang chủ
  2. » Luận Văn - Báo Cáo

Một phương pháp đánh giá mức độ anờtoàn của kết cấu không chịu tải trọng động theo lý thuyết tập mờ

175 3 0

Đ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

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC XÂY DỰNG -⁄ fo NÀ_ ` + ®, THNG DẠI HỌC XÂY DUNG N KHOA ĐÀO TẠO „ SAU ĐẠI HỌC LE CONG DUY NG MỘT PHƯƠNG PHÁP ĐÁNH GIÁ MỨC ĐỘ AN TOÀN CỦA KẾT CẤU KHUNG CHỊU TẢI TRỌNG ĐỘNG THEO LÝ THUYẾT TẬP MỜ Chuyên ngành : Kỹ thuật xây dựng cơng trình dân dụng cơng nghiệp Mã số 62.58.02.08 LUẬN ÁN TIẾN SỸ KỸ THUẬT NGƯỜI HƯỚNG DẪN KHOA HỌC GS.TS LÊ XUÂN HUỚM Ha Noi — 2014 // LOI CAM DOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết nêu luận án trung thực chưa công bố cơng trình khác Tác giả luận án “) Lê Công uy li LOI CAM ON Tôi xin bày tỏ lòng biết ơn chân thành GS.TS.Lê Xuân Huỳnh, người thầy tận tình hướng dẫn dạy bảo nhiều, thường xuyên động viên, tạo điều kiện thuận lợi, giúp đỡ tơi hồn thành luận án nâng cao kiến thức khoa học Tơi xin chân thành cảm ơn PGS.TS.Trần Văn Liên, GS.TS.Nguyễn Văn Phó có ý kiến đóng góp quý báu giúp tơi hồn thành luận án Tơi xin trân trọng cảm ơn TS.Nguyễn Xuân Thành, TS.Nguyén Tiến Dũng thầy, cô giáo Bộ môn Cơ Học Kết cấu cán Khoa Đào tạo sau Đại học Trường Đại Học Xây Dựng Hà Nội ln tạo điều kiện thuận lợi để tơi hồn thành luận án Tơi xin chân thành cảm ơn thầy, cô giáo, cán Khoa Xây dựng Trường Đại Học Duy Tân Đà Nắng, bạn đồng nghiệp tạo điều kiện thuận lợi, giúp đỡ, động viên cho tơi q trình học tập nghiên cứu Cuối xin tỏ lòng biết ơn sâu sắc đến bố, mẹ hai bên gia đình đặc biệt người bạn đời yêu quý hỗ trợ, động viên tinh thần lẫn vật chất giúp tơi hồn thành tốt luận án Nghiên Cứu Sinh Lê Công Duy 11 MỤC LỤC Lời cam đOaI - - 0v; 1:0 00000 00 11000000 TT 000 0000 00 I li MUC LUC 0.0 eee ill Danh mục ký hiệu luận án 55c S1 se Vil Danh mục hình vẽ đỒ thị - - -GG QC SH ve ree 1X Danh mục bảng biểu -(G26 SG 3E 1131 3K kg ng ri xH 6062700355 - ,ÔỎ 1 Ý nghĩa khoa học thực tiễn đề tầi - se csEsEsEx+e sex cee J Mục tiêu nội dung nghiên cứu để tài -555555 NLPSolve(N*(56-N/(2*Rb*b))}+-52*Fs*Rs, Rs = 25.2 30.8, 106.043); > NLPSolve(N*(S56-N/(2*Rb*b))+52*FS*Rs, Rs = 25.9 30.1, 106.043 106.043); >NLPSolve(N*(56-N/(2*Rb*b))+52*Fs*Rs, Rs = 26.6 29.4, 106.043 106.043); >NLPSolve(N*(56-N/(2*Rb*b))+52*Fs*Rs, Rs = 27.3 28.7, Fs = 29.454 29.454, Rb = 1.305 1.595, b= 25 25, N = 106.043 Fs = 29.454 29.454, Rb = 1.3413 1.5588, b = 25 25, N= Fs = 29.454 29.454, Rb = 1.3775 1.5225, b= 25 25, N= Fs = 29.454 29.454, Rb = 1.4138 1.4863, b=25 25, N= 106.043 106.043); >NLPSolve(N*(56-N/(2*Rb*b))+52*Fs*Rs, Rs = 25.2 30.8, Fs = 29.454 29.454, Rb = 1.305 1.595, b= 25 25, N= 106.043 106.043, maximize); >NLPSolve(N*(56-N/(2*Rb*b))+52*Fs*Rs, Rs = 25.9 30.1, Fs = 29.454 29.454, Rb = 1.3413 1.5588, b= 25 25, N= 106.043 106.043, maximize); >NLPSolve(N*(56-N/(2*Rb*b))+52*Fs*Rs, Rs = 26.6 29.4, Fs = 29.454 29.454, Rb = 1.3775 1.5225, b= 25 25,N= 106.043 106.043, maximize); >NLPSolve(N*¥(56-N/(2*Rb*b))+52*Fs*Rs, Rs = 27.3 28.7, Fs = 29.454 29.454, Rb = 1.4138 1.4863, b= 25 25, N= 106.043 106.043, maximize) ; >NLPSolve(N*(56-N/(2*Rb*b))+52*Fs*Rs, Rs = 28 28, Fs = 29.454 29.454, Rb = 1.45 1.45, b = 25 25, N = 106.043 106.043); Khao sat d6 tin cay theo biên độ mờ tải > with(Optimization); tinh*chuyen*vi-x7max- bdo0%>; bdo5%; bdol0%; bdol 5%; bdo20%; bdo25%; print( output redirected -); # input placeholder tinh chuyen vi - x7max - bdo0%- bdo5%- bdol10%- bdolS%- bdo20%- bdo25% >%bdo0%; > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi=0.le-) .1,P=15 15,t=0 3); > NLPSolve(x7, m = 5454 .634, E = 245] 2848.8, xi = 0.2e-1 0.875e-1, P= 15 15, t =1.372 1.372); > NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 15 15, t =1.372 1.372); > NLPSolve(x7, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-] 0.625e-1, P=15 > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi =0.le-1 1, P=15 15, t=1.372 1.372); 15, t = 1.372 1.372, maximize); > NLPSolve(x7, m = 5454 .634, E = 2451] 2848.8, xi = 0.2e-1 0.875e-1, P= 15 1S, t= 1.372 1.372, maximize); > NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, xi=0.2e-1 0.75e-1, P=15 18, t = 1.372 1.372, maximize); > NLPSolve(x7, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-] 0.62S5e-1, P= 15 15, t = 1.372 1.372,maximize); > NLPSolve(x7, m = 5862 .5862, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 15 15, t =1.372 1.372, maximize); >%bdo5%; > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi =0.le-1 .1, P= 14.25 15.75, t =1.372 1.372); > NLPSolve(x7, m = 5454 .634, E= 2451 2848.8, xi = 0.2e-1 0.875e-1, P = 14.438 15.563, t = 1.372 1.372); > NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 14.625 15.375, t = 1.372 1.372); > NLPSolve(x7, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-], P= 14.813 15.188, t = 1.372 1.372); > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi=0.le-1 .1, P= 14.25 15.75, t= 3, maximize); > NLPSolve(x7, m = 5454 .634, E = 245] 2848.8, xi = 0.2e-1 0.875e-1, P= 14.438 15.563, t= 1.372 1.372, maximize); > NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 14.625 15.375, t= 1.372 1.372, maximize); > NLPSolve(x7, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-] 0.625e-1, P= 14.813 15.188, t = 1.372 1.372,maximize); > NLPSolve(x7, m = 5862 .5862, E = 2650 2650, x1 = 0.5e-1 0.5e-1, P=15 15, t= 1.372 1.372, maximize); >bdol0%; > NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi =0.le-] 1,P= 13.5 16.5,t=1.372 1.372); > NLPSolve(x7, m = 5454 .634, E= 245] 2848.8, xi = 0.2e-1 0.875e-1, P = 13.875 16.125, t= 1.372 1.372); > NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 14.25 15.75, t = 1.372 1.372); > NLPSolve(x7, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.625 15.375, t = 1.372 1.372); > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi = O0.le-1 1, P= 13.5 16.5, t= 1.372 1.372, maximize); > NLPSolve(x7, m = 5454 .634, E = 245] 2848.8, xi = 0.2e-1 0.875e-1, P = 13.875 16.125, t= 1.372 1.372, maximize); 12 > NLPSolve(x7, m = 5601 619, E = 2517.5 2782.5, x1 = 0.2e-1 0.75e-1, P= 14.25 15.75, t = 1.372 1.372, maximize); > NLPSolve(x7, m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.625 15.375, t= 1.372 1.372,maximize); > NLPSolve(x7, m = 5862 5862, E = 2650 2650, xi = 0.5e-1 O.5e-1, P= 15 15, t= 1.372 1.372, maximize); >bdol15%; > NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi = 0.le-1 .1, P= 12.75 17.25, t = 1.372 1.372); > NLPSolve(x7, m = 5454 .634, E = 2451 2848.8, xi = 0.2e-1 0.875e-1, P = 13.313 16.688, t = 1.372 1.372); > NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 13.875 16.125, t= 1.372 1.372); > NLPSolve(x7, m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.438 15.563, t = 1.372 1.372); > NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi =0.1le-1] .1, P= 12.75 17.25, t =1.372 1.372); > NLPSolve(x7, m = 5454 .634, E = 245] 2848.8, xi = 0.2e-1 0.875e-1, P= 13.313 16.688, t = 1.372 1.372, maximize); > NLPSolve(x7, m = 5601 619, E= 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P = 13.875 16.125, t= 1.372 1.372, maximize); > NLPSolve(x7, m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-] 0.625e-1, P= 14.438 15.563, t= 1.372 1.372,maximize); > NLPSolve(x7, m = 5862 5862, E = 2650 2650, xi = 0.Se-1 0.Se-1, P= 15 15, t= 1.372 1.372, maximize); >bdo20%; > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi=O.Je-] .1, P= 12 18, t=1.372 1.372); > NLPSolve(x7, m = 5454 634, E= 2451 2848.8, xi = 0.2e-1 0.875e-1, P= 12.75 17.25, t = 1.372 1.372); > NLPSolve(x7, m = 5601 619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 13.5 16.5, t= 1.372 1.372); > NLPSolve(x7, m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P = 14.25 15.75, t= 1.372 1.372); > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi =0.le-1 1, P= 12 18, t =0 3,smaximize); > NLPSolve(x7, m = 5454 634, E = 245] 2848.8, xi = 0.2e-1 0.875e-1, P= 12.75 17.25, t = 1.372 1.372, maximize); > NLPSolve(x7, m = 5601 619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 13.5 16.5, t= 1.372 1.372, maximize); > NLPSolve(x7, m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-] 0.625e-1, P= 14.25 15.75, t= 1.372 1.372,maximize); > NLPSolve(x7, m = 5862 5862, E = 2650 2650, xi = 0.Se-1 0.5e-1, P= 15 15, t = 1.372 1.372, maximize); > wIith(Optimization); tình-trangthaiQmax-`bdoO0%`; bdo5%; bdo10%; bdo1l 5%; bdo20%; bdo25%; pnnt( output redirected `); # input placeholder tinh-trangthaiQmax - bdo0%- bdo5%- bdol10%- bdo15%- bdo20%- bdo25% >%bdo0%; > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 53066 64858, E= 2385 2915, xi = 0.1e-] 1,P 15, t =1.301 1.301); > NLPSolve((13.946* E*x1/N+0.26)*N, N=106.043 106.043, m = 5454 15 634, E = 2451 2848.8, xi = 0.2e-1 0.875e-1, P= 15 15, t=1.301 1.301); > NLPSolve(((13.946*E*x 1/N+0.26)*N, N=106.043 106.043, m = 5601 619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 15 15, t =1.30) 1.301); > NLPSolve((13.946* E*x1/N+0.26)*N, N=106.043 106.043, m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 15 15, t =1.301 1.301); > NLPSolve((1 3.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 53066 64858, E = 2385 2915, xi =0.le-1 1,P=15 15, t = 1.301 1.301, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = S454 634, E = 2451 2848.8, xi = 0.2e-1 0.875e-1, P= 15 15, t= 1.301 1.301, maximize); > NLPSolve((J 3.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5601 619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 15 15, t=].301 1.301,maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5749 15 15, t=1.301 1.301, maximize); 604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= > NLPSolve((13.946* E*x 1/N+0.26)*N, N=106.043 106.043, m = 5862 5862, E = 2650 2650, xi = 0.5e-1 O.5e-1, P= 15 15, t=1.301 1.301, maximize); >%bdo05%; > NLPSolve((13.946*E*x 1/N+0.26)*N, N=106.043 106.043, m = 53066 64858, E = 2385 2915, x1 = O.le-1 .1, P= 14.25 15.75, t=1.301 1.301); > NLPSolve((13.946* E*x 1/N+0.26)*N, N=106.043 106.043, m = 5454 634, E = 2451 2848.8, xi = 0.2e-1 0.875e-1, P= 14.438 15.563, t=].301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5601 619, E = 2517.5 2782.5, xi = 0.2e-] 0.75e-1, P= 14.625 15.375, t=1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043,m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.813 15.188, t=1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 53066 64858, E = 2385 2915, x1 =O.le-1 .1, P= 14.25 15.75, t =1.30] 1.301 maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5454 .634, E = 245] 2848.8, xi = 0.2e-1 0.875e-1, P= 14.438 15.563, t=1.301 1.301, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5601 619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 14.625 15.375, t=1.30] 1.301, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5749 604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.813 15.188, t =1.301 1.301,maximize); > NLPSolve((13.946* E*x1/N+0.26)*N, N=106.043 106.043, m = 5862 5862, E = 2650 2650, xi = 0.5e-1 0.Se-1, P=15 15, t =1.301 1.301, maximize); >bdo10%; 13 > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 53066 .64858, E = 2385 2915, xi=0.le-1 1, P= 13.5 16.5, t=1.372 1.372); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5454 .634, E = 2451 2848.8, xi = 0.2e-] 0.875e-1, P= 13.875 16.125, t =1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 14.25 15.75, t=1.301 1.301); > NLPSolve((13.946*E*x 1/N+0.26)*N, N=106.043 106.043, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.625 15.375, t=1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=]106.043 106.043, m = 53066 .64858, E = 2385 2915, xi = O.le-1] .1, P=13.5 16.5, t= 1.372 1.372, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5454 .634, E= 245] 2848.8, xi = 0.2e-1] 0.875e-1, P= 13.875 16.125, t =1.301 1.301, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 14.25 15.75, t=1.301 1.301, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.625 15.375, t =1.301 1.301,maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5862 .5862, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 15 15, t=1.301 1.301, maximize); >bdol 5%; > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 53066 .64858, E = 2385 2915, xi = O.le-1} .1, P= 12.75 17.25, t=1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5454 .634, E = 2451 2848.8, xi = 0.2e-1 0.875e-1, P= 13.313 16.688, t =1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5601 .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 13.875 16.125, t= 1.30] 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.438 1§.563, t=1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 53066 .64858, E = 2385 2915, xi=0.Je-1] .1, P=12.75 17.25, t =1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5454 .634, E= 245] 2848.8, xi = 0.2e-1 0.875e-1, P= 13.313 16.688, t =1.301 1.301, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5601 .619, E = 2517.5 2782.5, x1 = 0.2e-1 0.75e-1, P= 13.875 16.125, t =1.301 1.301, maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.438 15.563, t=1.301 1.301 ,maximize); > NLPSolve((13.946* E*x1/N+0.26)*N, N=106.043 106.043, m = 5862 .5862, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 15 15, t=1.301 1.301, maximize); >bdo20%; > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5366 64858, E = 238§5 2915, xi =0.1e-1 1,P= 12 18, t =].301 1.301); > NLPSolve((13.946*E*x 1/N+0.26)*N, N=106.043 106.043, m = 5454 .634, E = 245] 2848.8, xi = 0.2e-1 0.875e-1, P= 12.75 17.25, t =1.301 1.301); > NLPSolve((13.946*E*x1, N=106.043 106.043, m = 5601 .619, E = 2517.5 2782.5, x1 = 0.2e-] 0.7Se-1, P= 13.5 16.5, t =].301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5749 .604, E = 2583.8 2716.25, x1 = 0.4e-1 0.625e-1, P= 14.25 15.75, t = 1.301 1.301); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 53066 64858, E = 2385 2915, xi = O.}e-] .1, P= 12 18, t =1.301 1.301,maximize); > NLPSolve(13.946*E*x 1/N+0.26)*N, N=106.043 106.043, m = 5454 .634, E= 2451 2848.8, xi = 0.2e-1 0.87Se-1, P= 12.75 17.25, t =1.301 1.301, maximize); > NLPSolve((13.946*E*x 1/N+0.26)*N, N=106.043 106.043, m = 560] .619, E = 2517.5 2782.5, xi = 0.2e-1 0.75e-1, P= 13.5 16.5, t= 1.301 1.301, maximize); > NLPSolve((13.946*E*x /N+0.26)*N, N=106.043 106.043, m = 5749 .604, E = 2583.8 2716.25, xi = 0.4e-1 0.625e-1, P= 14.25 15.75, t= 1.301 1.301 ,maximize); > NLPSolve((13.946*E*x1/N+0.26)*N, N=106.043 106.043, m = 5862 .5862, E = 2650 2650, xi = 0.5e-1 0.5e-1, P=15 15, t= 1.301 1.301, maximize); Tải trọng ngắn hạn dạng hình chữ nhật > with(linalg); ddong-cnhat; print( output redirected °); # input placeholder ddong - cnhat > K :=matrix([.404*E, -.202*E, 0, 0, 0, 0, 0], [-.202*E, 404*E, -.202*E, 0, 0, 0, 0], [0, -.202*E, 404*E, -.202*E, 0, 0, 0], {0, 0, 202*E, 298*E, -0.96e-1*E, 0, 0], [0, 0, 0, -0.96e-1*E, 192*E, -0.96e-1*E, 0], [0, 0, 0, 0, -0.96e-1*E, 192*E, -0.96e-1*E], [0, 0, 0, 0, 0, -0.96e-1*E, 0.96e-1*E]]}); M := diag(1.033*m, m, m, m, m, m, 737*m); B := evalm(K-M*f); pttanso := det(B); tanso := solve(pttanso, f); > f1 := 0.8192537790e-2*E/m; f2 := 0.5246674796e-1*E/m; f3 := 1509264763*E/m; f4 := 2496215899*E/m; f5 := 3339001967*E/m; f6 := 4986868863*E/m; f7 := 7175572683*E/m; b] := evalm(K-M*f]); b2 := evalm(K-M*f2); b3 := evalm(K- 14 M*f3); b4 := evalm(K-M*f4); b5 := evalm(K-M*f5); b6 := evalm(K-M*f6); b7 := evalm(K-M*f7); BJ := delcols(delrows(bl, 1), 1); B11 := delcols(delrows(b1, 1), 7); phis] := -evalm( &* (inverse(B1), B11)); phil := stackmatrix(vector([{1]), phis1); B2 := delcols(delrows(b2, 1), 1); B21 := delcols(delrows(b2, 1), 7); phis2 := -evalm( &* (inverse(B2), B21)); phi2 := stackmatrix(vector([1]), phis2); B3 := delcols(delrows(b3, 1), 1); B31 := delcols(delrows(b3, 1), 7); phis3 := evalm( &**(inverse(B3), B31)); phi3 := stackmatrix(vector((1]), phis3); B4 := delcols(delrows(b4, 1), 1); B41 := delcols(delrows(b4, 1), 7); phis4 := -evalm( &**(inverse(B4), B41)); phi4 := stackmatnx(vector({ ]), phis4); B5 := delcols(delrows(bŠ, 1), 1); B51 := delcols(delrows(b5, 1), 7); phis5S := -evalm( &**(inverse(B5), B51)); phi5 := stackmatrix(vector((]]), phis5); B6 := delcols(delrows(b6, 1), 1); B61 := delcols(delrows(b6, 1), 7); phis6 := evalm( &**(inverse(B6), B61)); phi6 := stackmatrix(vector([1]), phis6); B7 := delcols(delrows(b7, 1), 1); B71] := delcols(delrows(b7, 1), 7); phis7 := -evalm( &**(inverse(B7), B71)); phi? := stackmatrix(vector([1]), phis7); phi := concat(phil, phi2, phi3, phi4, phi5, phi6, phi7); > F := matrix(7, 1, [P, P, P, P, P, P, P]); Fl := evalm( &**(transpose(phil), F)); F2 := evalm(C &*`(transpose(phi2), F)); F3 := evalm( &*`(transpose(phi3), F)); F4 := cvalm(C &*`(transpose(phi4), F)); FŠ := evalm(C &*`(transpose(phiŠ5), F)); F6 := evalm(C &**(transpose(phi6), F)); F7 := evalm( &**(transpose(phi7), F)); M1 := evalm( &**(evalm( &**(transpose(phil), M)), phil)); M2 := evalm( &**(evalm( &* (transpose(phi2), M)), phi2)); M3 := evalm( &**(evalm( &**(transpose(phi3), M)), phi3)); M4 := evalm(C &*`(evalm( &*`(transpose(phi4), M)), phi4)); M5 := evalm( &*`(evalmC &*`(transpose(phi5), M)), phi5)); M6 := evalmC &*`(evalm(&*`(transpose(phi6), M)), phi6)); M7 := evalm(C &*`(evalmC &*`*(transpose(phi7), M)), phi7)); Mc := diag(MI, M2, M3, M4, MS, M6, M7); btam := evalm( &**(evalm( &** (inverse(Mc), transpose(phi))), M)); > omegal := f14(1/2); omega2 := f2%(1/2); omega3 := f3(1/2); omega4 := f44(1/2); omegaS := f54(1/2); omega6 := f6^(1/2); omega7 := f7^(1/2); `ϖ1` := f11/2)*(1-xI^2)X(1/2); `ϖ2` := f2(1/2)*(1-xI^2)(1/2); `&varpm;3` := Í3(1/2)*(1xI^2)AX1/2); `ϖ4` := ƒ4^(1/2)*(1-x1^2)A(1/2); `ϖ5` := fŠ^(1/2)*(1-x1^2)X1/2); `ϖ6` := £64(1 /2)*(1-xi42)4(1/2); `ϖ7` := f7(1/2)*(1-xi^2}X1/2); Kdl1 := (int(exp(-omega *xi*(t-x))*sin( ϖ1°*(t-x)), x = t)/ ϖ 1°; Kd2 := (Int(exp(-omega2*x1*(t-x))#*sinC ϖ2`*(t-x)), x = t))/ ϖ2`; Kd3 := (nt(exp(-omega3*x1#(t-x))*sinC &varp;3`*(t-X)), x=0 t))/ ϖ3`; Kd4 := (mnt(exp(-omega4*xi*(t-x))*sinC ϖ4`*(t-x)), x = Ð))/@ϖ4`; Kd5 := (nt(exp(omega5*xi*(t-x))*sinC &varp;5`*(t-x)), x = t))/ &varp;5`; Kđố := (nt(exp(-omega6*xi*(t-x))*sinC ϖ6`*(t-x)), x = t))/'ϖ6`; Kd7 := (nt(exp(-omega7*xi*(t-x))*sinC ϖ7-*(t-x)), x = t))/ ϖ7°; > U := matrix(7, 1, [evalm(F1*Kd1/M1), evalm(F2*Kd2/M2), evalm(F3*Kd3/M3), evalm(F4*Kd4/M4), evalm(F5*Kd5/M5), evalm(F6*Kd6/M6), evalm(F7*Kd7/M7)]); chvi := evalm( &**(phi, U)); n:=chvi; n[1,1]; n{2,1]; n[3,1]}; n{4,1]; n[Š,1]; n[6,1]; n[7,1 );n1:=n[1,1]; n2:=n[2,1]; n3:=n[3,1]; n4:=n[4,1]; n5:=n[5,1]; n6:=n[6, 1]; n7:=n[7,1]; x1:= evalm(n1); x2:= evalm(n2); x3:= evalm(n3); x4:= evalm(n4); x5:= evalm(n5); x6:= evalm(n6); x7:= evalm(n7); > cv := matrix(7, 1, [evalm(x1), evalm(x2), evalm(x3), evalm(x4), evalm(x5), evalm(x6), evalm(x7)]); > cvd := evalm( &* (btam, cv)); > zd]t := evalm(0.9275913957e-2*cv[1, 1]4+0.1758297076e-1*cv[2, 1]+0.2547323934e-1*cv[3, 1]+0.3233038676e-1*cv[4, 1]+0.4399992700e-1*cv[5, 1}+0.5193456030e-1*cv[6, 1]4+0.4082896001e-1*cv[7, 1]}); zd2t := evalm(0.6135978462e-1*cv{[], 1]+.1028618291*cv[2, 1}]+.1196071016*cv[3, 1]+.1052860590*cv[4, 1]+0.1761035331e-1*cv[5, ]-0.7968991449e-1*cv[6, 1]0.9834333638e-1*cv[7, 1]); zd3t := evalm(.1066494804*cv[1, 1]+.1268006521*cv{[2, 1]+0.5561835228e-1*cv[3, 1]0.5711979869e-1*cv[4, 1]-.2045387212*cv[Š, 1]-0.3039193058e-1*cv[6, 1]+.1411617061*cv[7, 1]); zd4t := evalm(.2388512605*cv[1, 1]+.1672813271*cv{[2., 1]-.1033762887*cv{3, 1]-.2462866098*cv[4, 1]+0.9340683038e-1 *cv{[5, 1]+.1902215049*cv[6, 1]-.1529883122*cv[7, 1]); zd5t := evalm(0.6857295399e-1*cv[1, 1]+0.1941568113e-1*cv[2, 1]0.5964453772c-I*cv[3, 1)-0.4011404927e-1*cv[4, ]+.1405031133*cv[Š, 1]-.1675674029*cv[6, 1)+0.7899372176e-1*cv[7, 1]); zđốt := evalm(.3681044923*cv[1, 1]-.1960667682*cv[2, 1]-.2644394020*cv[3, 1]+.3200219345*cv[4, 1]-.1125775664*cv[{[ŠS, 1]+0.3962455858e-1*cv[6, 1]-0.1032480275e-1*cv[7, 1]); zdt7 := evalm(.1471861869*cv[1, 1]-.2378756804*cv[2, ]+.2267615768*cv[3, 1]-.1141180848*cv[4, 1]+0.2159618110e-1*cv[5, 1]-0.4111393937e-2*cv[6, ]+0.6720486576e-3*cv[7, 1)); >vdtl := diff(zd1t, 0); vdt2 := diff(zd2t, t); vdt3 := diff(zd3t, t); vdt4 := diff(zd4t, t); vdt5 := diff(zd5t, t); vdt6 := diff(zd6t, t); vdt7 := diff(zd7t, t); zd1 := subs(t = 1, zd1t); zd2 := subs(t = 1, zd2t); 2d3 := subs(t = 1, zd3t); zd4 := subs(t = 1, zd4t); zd5 := subs(t = 1, zd5t); zd6: subs(t = 1, zd6t); zd7 := subs(t = 1, zd7t); vd1 := subs(t = 1, vdt1l); vd2 := subs(t = 1, vdt2); vd3 := subs(t = 1, vdt3); vd4 := subs(t = 1, vdt4); vd5 := subs(t = 1, vdtS); vd6 := subs(t = 1, vdt6); vd7 := subs(t = 1, vdt7); > with(linalg); tim*nhiem-gd2; > f1 := 0.8192537790e-2*E/m; f2 := 0.5246674796e-1*E/m; f3 := 1509264763*E/m; f4 := 2496215899*E/m; f5 := 3339001967*E/m; f6 := 4986868863*E/m; f7 := 7175572683*E/m; omega! := f14(1/2); omega2 := f2(1/2); omega3 := f34(1/2); omega4 := f44(1/2); omega5 := f{ŠX1/2); omeega6 := f6^(1/2); omega7 := f7(1/2); `&varp;1` := f1^(1/2)*(1-x1^2)X1/2); `ϖ2` := f2(1/2)*(1-x1^2)X1/2); `ϖ3` := f3(1/2)*(1-xi^2)X1/2); `ϖ4` := f4^(1/2)*(1-xI^2)⁄X1/2); `&varpl;5` := fŠ5N1/2)*(1-xI^2)X1/2); `ϖ6` := f6^(1/2)*(1-x1^2)A1/2); `ϖ7` := f7A(1/2)*(1-xI1^2)X1/2); phì := Matnx(7, 7, {(1, 1) = 1, (1,2)=1, 0,3) =1, 0, 4) = 1, 0, 5) = 1, (1, 6) = 1, 1, 7) = 1, (2, 1) = 1.958104493, (2, 2) = 1.731692347, (2, 3) = }.228182952, (2, 4) = 7234695373, (2, 5) = 2924826399, (2, 6) = -.550215973S, (2, 7) = -1.669488035, (3, 1) = 2.836793912, (3, 2) = 2.013601200, (3, 3) = 5387157790, (3, 4) = -.4470887281, (3, 5) = -.8985001212, (3, 6) = -.7420879341, (3, 7) = 1.591485681, (4, 1) = 3.600431147, (4, 2) = 1.772504574, (4, 3) = -.5532586922, (4, 4) = -1.065156899, (4, 5) = -.6042879953, (4, 6) = 8980674380, (4, 7) = -.8009174235, (5, 1) = 4.899994200, (5, 2) = 2964726015, (5, 3) = -1.981148883, (5, 4) = 4039721440, (5, 5) = 2.116573774, (5, 6) = -.3159228658, (5, 7) = 1515689451, (6, 1) = 5.781396964, (6, 2) = -1.341590134, (6, 3) = -.2943742827, (6, 4) = 8226827615, (6, 5) = -2.524276951, (6, 6) = 1111971461, (6, 7) = -0.2885508503e-1, (7, 1) = 6.169421463, (7, 2) = -2.246435457, (7, 3) = 1.855201080, (7, 4) = -.8977668178, (7, 5) = 1.614628266, (7, 6) = -0.3931366020e1, (7, 7) = 0.639980] 404e-2 } ); 15 > zđ] := d]; zđ2 := d2; zd3 := đ3; zd4 := d4; zdŠ := đŠ5; zd6 := d6; zd7 := d7; vd1 := v1; vd2 := v2; vd3 := v3; vd4 := v4; vd5 := v5; vd6 := v6; vd7 := v7; > z] := e\(-xi*omegal *t)*(zd1*cos( ϖ1°*t)+(vd1+xi* omega] *zd1)*sin( ϖ1°*t)/ ϖ1>); z2 := e^(xi*omega2*t)*(zd2*cosC ϖ2`#*t)+(vd2+xi*omega2*zd2)*sinC &varpl;2`*t)/ˆ&varpl;2`); z3 := e^(- xi*omega3*t)*(zd3*cos(C ϖ 3`*t)+(vd3+xi*omega3*zd3)*sin(C &varpr;3`*t)/ˆ&varpl;3`); z4 := e^(xI*omega4*t)*(zd4*#cos( ϖ4`*t)+(vd4+xi*omega4*zd4)*sinC ϖ4`*t)/&varpl;4`); 25 := e^(- xi*omega5*t)*(zd5*cos(C ϖ5`*t)+(vd5+xI*œnega5*zdŠ)*sinC ϖ5`*t)/ˆ&varpl;5`); z6 := e^(- xi*omega6*t)*(zd6*cos( ϖ6 *t)}+(vd6+x1*omega6*zd6)* sin( ϖ6 *t)/ ϖ6); 27 := e^(- xi*omega7*t)*(zd7*cos( Sϖ7-*t)+(vd7+xi*omega7*zd7)*sin( ϖ7 *t)/ ϖ7 >); > z:= matrix(7, 1, [z1, z2, z3, z4, z5, z6, z?]); > chvi2 := evalm( &** (phi, z)); > with(Optimization); tinh*chuyen*vi-x72max-xung; print( output redirected ); # input placeholder tình chuyen vị - x/2max - xung > X7 := 6.169421463*e(-0.9051263884e-1*xì#sart(E/m)#t)*(d1 *cos(0.9051263884e-] *sart(E/m)*sart(1x1^2)*1)+(1 1.048180§2*(v1+0.9051263884e-1*xi?sart(E/m)*đ1))*smm(0.90512638§4e-1 *sqrt(E/m)*sqrt(1xi^2)*Ð0/Gqrt(E/m)*sart(1-x1^2)))-2.246435457*e^(-.22905621 14*xi*sartŒ/m)*0*(đ2*cos(.22905621 14*sart(E/m)*sart(1xi42)*t)+(4.365740592*(v2+.22905621 14*xi*sqrt(E/m)*d2))*sin(.22905621 14*sqrt(E/m)*sqrt(1-xi42)*t)/(sqrt(E/m)*sqrt(] xi42)))+1.855201080*e“(-.3884925692*xi*sqrt(E/m)*t)*(d3*cos(.3884925692* sqrt(E/m)*sqrt(1 xi42)*t)+(2.57405 1808*(v3+.3884925692*xi*sqrt(E/m)*d3))*sin(.3884925692* sqrt(E/m)*sqrt(1 -xi42)*t)/(sqrt(E/m)*sqrt(1xi42)))-.8977668 78*e*(-.49962 4466*x1* sqrt(E/m)*t)*(d4*cos(.4996214466* sqrt(E/m)*sqrt(1xi42)*t)+(2.001515361*(v4+.4996214466*xi*sqrt(E/m)*d4))*sin(.4996214466* sqrt(E/m)*sqrt(1 -xi42)*t)/(sqrt(E/m)*sqrt(1xi42)))+1.614628266*e“(-.5778409787*xi* sqrt(E/m)*t)* (d5*cos(.5778409787*sqrt(E/m)*sqrt(1 xi42)*t)+(1.730579929* (v5+.5778409787*xi*sqrt(E/m)*d5))*sin(.5778409787*sqrt(E/m)*sqrt(1 -xi42)*t)/(sqrt(E/m)*sqrt(1 x142)))-0.3931366020e-1*e%(-.7061776592*xi*sqrt(E/m)*t)*(d6*cos(.7061776592* sqrt(E/m)*sqrt(1xi42)*t)+(1.4 16074251 *(v6+.7061776592*x1* sqrt(E/m)*d6))*sin(.7061776592*sqrt(E/m)*sqrt(1 -xi42)*t)/(sqrt(E/m)*sqrt(1 x142)))+0.6399801404e-2*e*(-.8470875210*xi*sqrt(E/m)*t)*(d7*cos(.8470875210*sqrt(E/m)*sart(1 xi42)*t)+(1.180515561*(v7+.8470875210*xi*sqrt(E/m)*d7))*sin(.8470875210*sqrt(E/m)*sqrt(1-x142)*t)/(sqrt(E/m)* sqrt(1 x142))); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi=O.le-1 .1,t=.1 3, dl =0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.591 3e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 5454 .634, E = 2451.3 2848.75, e = 2.71 2.71, xi = 0.2e-1 0.875e-1, t= 264 .264, dl = 0.406e-1 0.562e-1, d2 = 0.265e1 0.381le-1, d3 = 0.165e-] 0.244e-1, d4 = 0.69e-2 0.107e-1, d5 = 0.23e-2 0.4e-2, d6 = 0.14e-2 0.38e-2, d7 = 0.6e-4 0.4e2, v1 = 7705 1.0834, v2 = 4022 .60, v3 = 0.762e-1 .2608, v4 = -0.74e-1 0.75e-2, v5 = -0.622e-1 -0.19e-2, v6 = -.1238 0.67e-1, v7 = -0.22e-1 -0.3e-2, maximize); NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, e =2.7] 2.71, xi = 0.3e-] 0.7S5e-1, t = 264 .264, dl] = 0.43e-1 0.535e-1, d2 = 0.282e-] : 0.36e-1, d3 = 0.177e-1 0.23e-1, d4 = 0.75e-2 0.10e-1, d5 = 0.26e-2 0.37e-2, d6 = 0.18e-2 0.34e-2, d7 = 0.le-3 0.3e-3, v1 = 8174 1.026, v2 = 4316 .5634, v3 = 0.962e-1 .1833, v4 = -0.584e-] -0.4]e-2, v5 = -0.547e-1 -0.12e-1, v6 = -.1129 -0.752e-1, v7 = -0.18Se-1 -0.58e-2, maximize); NLPSolve(x7, m = 5749 .604, E = 2583.8 2716.25, e = 2.71 2.71, xi = 0.4e-1 0.625e-1, t = 264 .264, dl = 0.455Se-1 0.507e-1, d2 = 0.299e-1 0.338e-1, d3 = 0.188e-1 0.215e-1, d4 = 0.81e-2 0.94e-2, d5 = 0.28e-2 0.34e-2, d6 = 0.21e-2 0.29e-2, d7 = 0.le-3 0.3e-3, v1 = 8643 .9686, v2 = 4609 .5269, v3 = 1163 .1598, v4 = -0.427e-1 -0.156e-1, v5 = -0.472e-1 -0.259e-1, v6 = 1021 -0.832e-1, v7 = -0.15e-1 -0.87e-2, maximize); % > NLPSolve(x7, m = 53066 64858, E = 2358 2915, e =2.7] 2.71, xi=0.le-1 .1,t = 264 264, dl = 0.3809e-1 0.5898e- 1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 5454 634, E= 2451.3 2848.75, e = 2.71 2.71, xi = 0.2e-1 0.875e-1, t = 264 .264, dl = 0.406e-1 0.562e-1, d2 = 0.265e-1 0.381e-1, d3 = 0.165e-1 0.244e-1, d4 = 0.69e-2 0.107e-1, dS = 0.23e-2 0.4e-2, d6 = 0.14e-2 0.38e-2, d7 = 0.6e-4 0.4e-2, v] =.7705 1.0834, v2 = 4022 .60, v3 = 0.762e-1 .2608, v4 = -0.74e-1 0.75e-2, v5 = -0.622e-1 -0.19e-2, v6 = -.1238 -0.67e-1, v7 = 0.22e-1 -0.3e-2); NLPSolve(x7, m = 5601 .619, E = 2517.5 2782.5, e = 2.71 2.71, x1 = 0.3e-1 0.7S5e-1, t= 264 .264, dl | = 0.43e-1 0.535e-1, d2 = 0.282e-1 0.36e-1, d3 = 0.177e-1 0.23e-1, d4 = 0.75e-2 0.10e-1, d5 = 0.26¢e-2 0.37e-2, d6 = - 0.18e-2 0.34e-2, d7 = 0.le-3 0.3e-3, v1 = 8174 1.026, v2 = 4316 .5634, v3 = 0.962e-1 .1833, v4 = -0.584e-1 -0.41e-2, v5 = -0.547e-1 -0.12e-1, v6 = -.1129 -0.752e-1, v7 = -0.185e-1 -0.58e-2); NLPSolve(x7, m = 5749 .604, E = 2583.8 2716.25, e =2.71 2.71, xi = 0.4e-1 0.625e-1, t = 264 .264, dl = 0.45Se-1 0.507e-1, d2 = 0.299e-1 0.338e-1, d3 = 0.188e1 0.215e-1, d4 = 0.81e-2 0.94e-2, d5 = 0.28e-2 0.34e-2, d6 = 0.21e-2 0.29e-2, d7 = 0.le-3 0.3e-3, v1 = 8643 .9686, v2 = 4609 .5269, v3 = 1163 .1598, v4 = -0.427e-1 -0.1S6e-1, v5 = -0.472e-1 -0.259e-1, v6 = -.102] -0.832e-], v7 = -0.15e-1 -0.87e-2); > NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.5e-1 0.5e-1, t = 264 .264, dl = 0.4794e-1 0.4794e-1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, dS = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 .91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9]21e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1152e-]); > with(Optimization); tinh*chuyen*vi* theo*lat*cat* thoi* gian-x 72-xung; print( output redirected -); # input placeholder tinh chuyen vi theo lat cat thoi gian - x72 -xung 16 %t=0.1; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e =2.71 2.71, xi = O.le-1 .1,t=0.1 0.1, dl =0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 v5 = -0.6972e-1 0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi = O.1e-1 .1,t = 0.1 0.1, dl = 0.3809e-1 d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 v5 = -0.6972e-1 0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.Se-1 O.5e-1, t= 1, d2 = 0.3168e-] 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, dS = 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 v5 = -0.3972e-1 -0.3972e-1, v6 = -0.912]e-1 -0.9121e-1, v7 = -0.1152e-] -0.1152e-1); %t=0.2; NLPSolve(x7, m = 53066 .64858, E = 2358 2915,e =2.71 2.71, x1 =0.le-1 .1,t=0.2 0.107e-2 0.1895e-1, 0.5898e-1, 0.107e-2 0.1895e-1, 0.1 0.1, dl = 0.4794e-1 0.4794e0.307e-2 0.307e-2, d6 = 0.250e-2 .1363, v4 = -0.2708e-1 -0.2708e-1, 0.2, dl = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, vl = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 64858, E = 2358 2915,e =2.71 2.71, xi=0.le-1 .1,t=0.2 0.2, d1 = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-1] -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = 0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 26S0 2650, e = 2.71 2.71, xi = 0.5e-1 0.Se-1, t= 1, d2 = 0.3168e-] 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 v§ =-0.3972e-1 -0.3972e-l, v6 = -0.9121e-] -0.9121e-1, v7 = -0.1152e-] -0.1152e-1); %t=0.4; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e =2.71 2.71, x1 = O.Je-1 .1,t=0.4 0.2 0.2, dl = 0.4794e-1 0.4794e0.307e-2 0.307e-2, d6 = 0.250e-2 .1363, v4 = -0.2708e-] -0.2708e-1, 0.4, dl] =0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.259]e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-] -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, x1=0.le-] .1,t=0.4 0.4, dl = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 0.591 3e-1, v7 = 0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.5e-1 O.Se-1, t= 0.4 0.4, dl] = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 0.9121e-1, v7 = -0.1152e-1 -0.1152e-1); %1=0.6; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, x1 = 0.le-1 .1,t=0.6 0.6, d1 =0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 64858, E = 2358 2915, e =2.71 2.71, xi =0.le-1] .1, t= 0.6 0.6, dl] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, x1 = 0.5e-1 0.Se-1, t= 1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, dS = 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 v5 = -0.3972e-1 -0.3972e-1, v6 = -0.912]e-1 0.9121le-1, v7 = -0.1152e-1 -0.1152e-1); %t=0.8; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e =2.71 2.71, x1 =O.Je-1 .1,t=0.8 0.6 0.6, d] = 0.4794e-1 0.4794e0.307e-2 0.307e-2, d6 = 0.250e-2 .1363, v4 = -0.2708e-1] -0.2708e-1, 0.8, d] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.259]e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 64858, E = 2358 2915, e = 2.71 2.71, xi=0.le-1 .1,t=0.8 0.8, d] = 0.3809e-1 0.S898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-] -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, x1 = 0.5e-1 0.Se-1, t = 1, d2 = 0.3168e-] 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.250e-2, d7 = 0.19e-3 0.19e-3, v] = 91116 91116, v2 = 490308 .490308, v3 = 1363 v5 = -0.3972e-1 -0.3972e-1, v6 = -0.912]e-1 -0.912le-1, v7 = -0.1152e-1 -0.1152e-1); %t=1.0; NLPSolve(x7, m = 53066 64858, E = 2358 2915, ¢ =2.71 2.71, x1 =O.le-1 .1,t=1 0.8 0.8, dl = 0.4794e-1 0.4794e0.307e-2 0.307e-2, d6 = 0.250e-2 .1363, v4 = -0.2708e-1] -0.2708e-1, 1,d1 =0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591le-1, d4 = 0.628e-2 0.1136e-1, dS = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); 17 NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi = O.le-1 .1,t=1 1, d] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v] = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.5e-] 0.Se-1, t=1 1, d1 = 0.4794e-1 0.4794e-1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121le-], v7 = -0.1152e-1 -0.1152e-1); %t=1.2; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e =2.71 2.71, xi=O.le-] .1,t=1.2 1.2, d1 = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-] 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e =2.71 2.71, x1 =0.le-1 .1,t= 1.2 1.2, d1 = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591le-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, x1 = 0.5e-1 O.Se-1, t= 1.2 1.2, d1 = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, dS = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-] -0.1152e-1); %t=].4; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, x1 = 0.1e-1 .1, t= 1.4 1.4, d] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, x1=0.1le-1 .1,t=1.4 1.4, dl = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.259 le-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, vl = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.Se-1 O.S5e-1, t= 1.4 1.4, dl = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-], v7 = -0.1152e-1 -0.1152e-1); %t=1.6; NLPSolve(x7, m = 53066 64858, E = 2358 2915, e = 2.71 2.7], xi=0.le-1 .1,t=1.6 1.6, d] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.7] 2.71, xi = O.le-1 .1,t= 1.6 1.6, dl = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, x1 = 0.5e-1 0.Se-1, t= 1, d2 = 0.3168e-] 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 V5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-] -0.9121e-1, v7 = -0.1152e-] -0.1152e-1); %t=1.8; NLPSolve(x7, m = 53066 64858, E = 2358 2915,e =2.71 2.71, xi=0.le-] .1,t= 1.8 1.6 1.6, dI = 0.4794e-1 0.4794e0.307e-2 0.307e-2, d6 = 0.250e-2 .1363, v4 = -0.2708e-1 -0.2708e-1, 1.8, dl = 0.3809e-] 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-], v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 64858, E = 2358 2915,e = 2.71 2.71, xi =O.le-1 .1,t= 1.8 1.8, d1 = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, vl = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, x1 = 0.S5e-1 0.5e-1, t= 1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, dS = 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1152e-1); %t=2.0; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi=0.le-] .1, t=2.0 1.8 1.8, d] = 0.4794e-1 0.4794e0.307e-2 0.307e-2, d6 = 0.250e-2 .1363, v4 = -0.2708e-1 -0.2708e-1, 2.0, dl = 0.3809e-] 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dŠ = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1] 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-] -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi = O.Je-] .1,t=2.0 2.0, d] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = O.Se-1 0.5e-1, t= 2.0 2.0, dl = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 18 0.250e-2, d7 = 0.19e-3 0.19e-3, v1] = 91116 91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.912]e-1, v7 = -0.1152e-1 -0.1152e-1); Yot=2.2; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi=0.1le-1 .1, t= 2.2 2.2, dl] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1] .23026, v4 = -0.8964e-] 0.1895e-1, v5 = -0.6972e-1 0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = 0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi=0.le-1 .1,t=2.2 2.2, dl] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-] 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.Se-1 0.Se-1, t= 2.2 2.2, dl] = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168e-3, d3 = 0.20024e-] 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1152e-1); %t=2.4; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, x1 =O0.le-1] .1,t=2.4 2.4, d] = 0.3809e-] 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi=0.Je-1 .1,t= 2.4 2.4, dl] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.561 2e-] .23026, v4 = -0.8964e-1 0.189Se-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi =0.5e-1 0.5e-1, t= 2.4 2.4, dl = 0.4794e-1 0.4794e1, d2 = 0.3168¢e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, vl = 91116 .91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1152e-1); %ot=2.6; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi= 0.1e-1 1,t = 2.6 2.6, dl] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi = 0.1e-1 .1,t=2.6 2.6, dl = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.5e-1 O.5e-1, t = 2.6 2.6, d] = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168¢e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 .91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-] -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1152e-1); %t=2.8; NLPSolve(x7,m = 53066 6485§, E= 2358 2915,e=2.71 2.71,xi=0.le-1 1,t=2.8 2.8, dì = 0.3§09e-1 0.5898e-1, d2 =0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-] -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi=O.le-1 .1,t=2.8 2.8, d] = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, Vd = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, x1 = 0.5e-1 O.5e-1, t= 2.8 2.8, dl = 0.4794e-1 0.4794e], d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v] = 91116 .91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1152e-1); %t=3.0; NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.71 2.71, xi = 0.1e-1 .1, t= 3.0 3.0, dl = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 ().419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.561] 2e-] .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x7, m = 53066 .64858, E = 2358 2915, e = 2.7] 2.71, xi=O.le-] .1, t= 3.0 3.0, dl = 0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1] -0.12e-3); NLPSolve(x7, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.5e-] O.5e-1, t = 3.0 3.0, d] = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, d5 = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1352e-1); > with(Optmization); unh*chuyen*vi-x 2max-xung; print( output redirected °); # input placeholder tinh chuyen vi - x1 2max — xung 19 >x1:= e^(-0.9051263884c-1 *xi*sqrtŒ/m)*t)*(đ1*cos(0.9051263884e-1*sart(E/m)*sqrt(1xi42)*t)+(1 1.048 18082*(v1+0.905 1263884e-1*xi*sqrt(E/m)*d1))*sin(0.905 1263884e-1 *sqrt(E/m)*sqrt(1x142)*t)/(sqrt(E/m)*sqrt(1-x142)))+e(-.2290562114*xi*sqrt(E/m)*t)*(d2*cos(.22905621 14*sqrt(E/m)*sqrt(1x142)*t)+(4.365740592*(v2+.2290562114*xi*sqrt(E/m)*d2))*sin(.22905621 14*sqrt(E/m)*sqrt(1-xi*2)*t)/(sqrt(E/m)*sqrt(1x142)))+e(-.3884925692*xi* sqrt(E/m)*t)*(d3*cos(.3884925692*sqrt(E/m)*sqrt(1 x142)*t)+(2.57405 1808*(v3+.3884925692*xi*sqrt(E/m)*d3))*sin(.3884925692*sqri(E/m)*sqrt(1-xi42)*t)/(sqrt(E/m)*sqrt(1x142)))+e(-.4996214466*xi* sqrt(E/m)*t)*(d4*cos(.49962 14466*sqrt(E/m)*sqrt(1 - x142)*t)+(2.001515361*(v4+.4996214466*xi*sqrt(E/m)*d4))*sin(.4996214466*sqrt(E/m)*sqrt(1 -x142)*t)/(sqrt(E/m)*sqrt(1x1^2)))+e^(-.5778409787*xi*sqrtŒE/m)*t)*(d5*cos(.5778409787*sqart(E/m)*sqrt(1XI^2)*0+(1.730579929*(v5+.5778409787*xi*sart(E/m)*dŠ))*sin(.Š778409787*sart(E/m)*sart(1-x1^2)*t)/(sqart(Œ/m)*sqrt(1xi42)))+eA(-.7061776592*xi* sqrt(E/m)*t)*(d6*cos(.7061776592*sqrt(E/m)*sqrt(1 xi42)*t)+(1.41607425 1*(v6+.7061776592*xi* sqrt(E/m)*d6))*sin(.706 776592* sqrt(E/m)*sqrt(] -x142)*t)/(sqrt(E/m)*sqrt(1xI^2)))>e^X-.8470875210*xi#sqrt(E/m)*t)*(d7*cos(.8470875210*sart(E/m)*sqrt(1xi42)*t)+(1.1805 15561 *(v7+.8470875210*xi*sqrt(E/m)*d7))*sin(.8470875210*sqrt(E/m)*sqrt(1-xi42)*t)/(sqrt(E/m)*sqrt(1 - x142))); NLPSolve(x1, m = 53066 64858, E = 2358 2915, e =2.71 2.71, xi =0.le-1 1,t=.1 1, d] =0.3809e-1 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v] = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-] 0.1895e-1, v5 = -0.6972e-] -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(x1, m = 5454 .634, E = 2451.3 2848.75, e = 2.71 2.71, xi = 0.2e-1 0.875e-1, t= 75 .75, dl = 0.406e-1 0.562e-], d2 = 0.265e-1 0.381e-1, d3 = 0.16S5e-1 0.244e-1, d4 = 0.69e-2 0.107e-1, d5 = 0.23e-2 0.4e-2, d6 = 0.14e-2 0.38e-2, d7 = 0.6e-4 0.4e-2, v1] = 7705 1.0834, v2 = 4022 .60, v3 = 0.762e-1 .2608, v4 = -0.74e-1 0.75e-2, v5 = -0.622e-1 0.19e-2, v6 = -.1238 -0.67e-1, v7 = -0.22e-] -0.3e-2); NLPSolve(x1, m = 5601 .619, E = 2517.5 2782.5, e = 2.71 2.71, xi = 0.3e-1 0.75e-1, t = 75 .75, dl = 0.43e-1 0.535e-1, d2 = 0.282e-1 0.36e-1, d3 = 0.177e-1 0.23e-1, d4 = 0.75e-2 0.10e-1, d5 = 0.26e-2 0.37e-2, d6 = 0.18e-2 0.34e-2, d7 = Q.le-3 0.3e-3, vl = 8174 1.026, v2 = 4316 .5634, v3 = 0.962e-1 .1833, v4 = -0.584e-1 -0.41e-2, v5 = -0.547e-1 -0.12e1, V6 = -.1129 -0.752e-1, v7 = -0.185e-1 -0.58e-2); NLPSolve(x1, m = 5749 .604, E = 2583.8 2716.25, e = 2.71 2.71, xi = 0.4e-1 0.625e-1, t= 75 .75, dl = 0.455e-1 0.507e-1, d2 = 0.299e-1 0.338e-1, d3 = 0.188e-1 0.21 5e-1, d4 = 0.8le-2 0.94e-2, d5 = 0.28e-2 0.34e-2, d6 = 0.2]e-2 0.29e-2, d7 = 0.1le-3 0.3e-3, v1 = 8643 .9686, v2 = 4609 .5269, v3 = 1163 .1598, v4 = -0.427e-1 -0.156e-l, v5 = 0.472e-1 -0.259e-1, v6 = -.1021 -0.832e-1, v7 = -0.15e-1 -0.87e-2); % NLPSolve(x1,m = 53066 64858, E= 235§ 2915,e=2.71 2.71,xi=0.le-l 1,t=.75 75, =0.3809e-1 0.5898e-], d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.20Se-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(x1, m = 5454 .634, E = 2451.3 2848.75, e = 2.7] 2.71, xi = 0.2e-1 0.87Se-1, t= 75 .75, d1 = 0.406e-1 0.562e-1, d2 = 0.26Se-1 0.381le-1, d3 = 0.165e-] 0.244e-1, d4 = 0.69e-2 0.107e-1, dS = 0.23e-2 0.4e-2, d6 = 0.14e-2 0.38e-2, d7 = 0.6e-4 0.4e-2, v1 = 7705 1.0834, v2 = 4022 .60, v3 = 0.762e-1 .2608, v4 = -0.74e-1 0.75e-2, v5 = -0.622e-1 -0.19e-2, v6 = -.1238 -0.67e-1, v7 = 0.22e-] -0.3e-2, maximize); NLPSolve(x1, m = 5601 .619, E = 2517.5 2782.5, e = 2.71 2.71, xi = 0.3e-1 0.75e-1, t = 75 75, dl = 0.43e-1 0.535e-1, d2 = 0.282e-1 0.36e-1, d3 = 0.177e-1 0.23e-1, d4 = 0.75e-2 0.10e-1, d5 = 0.26e-2 0.37e-2, d6 = 0.18e-2 0.34e-2, d7 = 0.1e-3 0.3e-3, v1 = 8174 1.026, v2 = 4316 .5634, v3 = 0.962e-1 .1833, v4 = -0.584e-1 -0.41e2, V5 = -0.547e-] -0.12e-1, v6 = -.1129 -0.752e-1, v7 = -0.185e-1 -0.58e-2, maximize); NLPSolve(x1, m = 5749 .604, E= 2583.8 2716.25, e =2.71 2.71, x1 = 0.4e-1 0.625e-], t= 75 .75, d1 = 0.455e-1 0.507e-1, d2 = 0.299e-1 0.338e-1, d3 = 0.188e-] 0.215e-1, d4 = 0.8le-2 0.94e-2, d5 = 0.28e-2 0.34e-2, d6 = 0.21e-2 0.29e-2, d7 = 0.le-3 0.3e-3, v1 = 8643 .9686, v2 = 4609 .5269, v3 = 1163 .1598, v4 = -0.427e-1 0.156e-1, v5 = -0.472e-1 -0.259e-1, v6 = -.1021 -0.832e-1, v7 = -0.15e-1] -0.87e-2, maximize); NLPSolve(x1, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, xi = 0.Se-1 0.S5e-1, t= 75 .75, dl = 0.4794e-1 0.4794e1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-] 0.20024e-1, d4 = 0.867e-2 0.867e-2, dS = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 .91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121le-1 -0.912]e-1, v7 = -0.1152e-1 -0.1152e-1); > %tinhmomen; NLPSolve(13.946*E*x1, m = 53066 64858, E = 2358 2915,e =2.71 2.71, xi=O.le-1 1,t=.1 1, d1 =0.3809e-] 0.5898e-1, d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-] 0.2591e-1, d4 = 0.628e-2 0.1136e-1, dS = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e+4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-1 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3); NLPSolve(13.946*E*x1, m = 5454 .634, E = 2451.3 2848.75, e = 2.71 2.71, xi = 0.2e-1 0.875e-1, t= 75 .75, d1 = 0.406e-1 0.562e-1, d2 = 0.26Se-! 0.381le-1, d3 = 0.16Se-1 0.244e-1, d4 = 0.69e-2 0.107e-1, dS = 0.23e-2 0.4e-2, d6 = 0.14e-2 0.38e-2, d7 = 0.6e-4 0.4e-2, vi =.7705 1.0834, v2 = 4022 .60, v3 = 0.762e-1 .2608, v4 = -0.74e-1 0.75e-2, v5 = -0.622e-1 -0.19e-2, v6 = -.1238 0.67e-1, v7 = -0.22e-1 -0.3e-2); NLPSolve(13.946*E*x1, m = 560] .619, E = 2517.5 2782.5, e = 2.71 2.71, xi = 0.3e-1 0).75e-1, t= 75 .75, dl = 0.43e-1 0.535e-1, d2 = 0.282e-1 0.36e-1, d3 = 0.177e-1 0.23e-1, d4 = 0.75e-2 0.10e-1, d5 = 0.26e-2 0.37e-2, d6 = 0.18e-2 0.34e-2, d7 = 0.1e-3 0.3e-3, v] = 8174 1.026, v2 = 4316 .5634, v3 = 0.962e-1 .1833, v4 = -0.584e-1 -0.41e-2, v5 = -0.547e-1 -0.12e-1, v6 = -.1129 -0.752e-1, v7 = -0.185e-1 -0.S8e-2); NLPSolve(13.946*E*x1, m = 5749 .604, E = 2583.8 2716.25, e = 2.71 2.71, xi = 0.4e-1 0.625e-1, t= 75 .75, dl = 0.45S5e-1 0.507e-1, d2 = 0.299e1 0.338e-1, d3 = 0.188e-1 0.21Se-1, d4 = 0.8le-2 0.94e-2, d5 = 0.28e-2 0.34e-2, d6 = 0.2]e-2 0.29e-2, d7 = 0.]e-3 0.3¢e3, v1 = 8643 .9686, v2 = 4609 .5269, v3 = 1163 .1598, v4 = -0.427e-1 -0.156e-1, v5 = -0.472e-1 -0.259e-1, v6 = -.1021 -0,832e-1, v7 = -0.15e-1 -0.87e-2); % NLPSolve(13.946*E*x1, m = 53066 64858§,E=2358 2915,e=2.71 2.71,xi=0.le-l 1,t=.75 75, =0.3809c-1 0.5898e-], d2 = 0.2473e-1 0.4024e-1, d3 = 0.15306e-1 0.2591e-1, d4 = 0.628e-2 0.1136e-1, d5 = 0.205e-2 0.434e-2, d6 = 0.107e-2 0.419e-2, d7 = 0.2e-4 0.47e-3, v1 = 72356 1.14083, v2 = 3728 .6365, v3 = 0.5612e-1 .23026, v4 = -0.8964e-] 0.1895e-1, v5 = -0.6972e-1 -0.1568e-1, v6 = -.13457 -0.5913e-1, v7 = -0.2548e-1 -0.12e-3, maximize); NLPSolve(13.946*E*x1, m = 5454 .634, E = 2451.3 2848.75, e = 2.71 2.71, xi = 0.2e-] 0.875e-1, t = 75 .75, dl = 0.406e-1 0.562e-1, d2 = 0.265e-1 0.381e-1, d3 = 0.165e-1 0.244e-1, d4 = 0.69e-2 0.107e-1, d5 = 0.23e-2 0.4e-2, d6 = 0.14e-2 0.38e-2, d7 = 0.6e-4 0.4e-2, v1 = 7705 1.0834, v2 = 4022 .60, v3 = 0.762e-1 .2608, v4 = -0.74e-1 0.75e-2, v5 = -0.622e-1 -0.19e-2, v6 = -.1238 -0.67e-1, v7 = -0.22e-] -0.3e-2, maximize); NLPSolve(13.946*E*x1, m = 560] 619, E= 20 2517.5 2782.5, e = 2.71 2.71, xi = 0.3e-1 0.75e-1,t=.75 75, d1 = 0.43e-1 0.535e-1, d2 = 0.282e-1 0.36e-1, d3 = 0.177e-1 0.23e-1, d4 = 0.75e-2 0.10e-1, d5 = 0.26e-2 0.37e-2, d6 = 0.18e-2 0.34e-2, d7 = 0.le-3 0.3e-3, v] = 8174 1.026, v2 = 4316 5634, v3 = 0.962e-] .1833, v4 = -0.584e-1 -0.41e-2, v5 = -0.547e-1 -0.12e-1, v6 = -.1129 -0.752e-1, v7 = -0.185e-1 -0.58e-2, maximize); NLPSolve(13.946*E*x1, m = 5749 .604, E = 2583.8 2716.25, e = 2.71 2.71, xi = 0.4e-1 0.625e-1, t= 75 .75, dl = 0.45Se-1 0.507e-], d2 = 0.299e-1 0.338e-1, d3 = 0.188e-1 0.215e-1, d4 = 0.81e-2 0.94e-2, d5 = 0.28e-2 0.34e-2, d6 = 0.21e-2 0.29e-2, d7 = 0.1e-3 0.3e-3, v1 = 8643 .9686, v2 = 4609 .5269, v3 = 1163 .1598, v4 = 0.427e-1 -0.156e-1, v5 = -0.472e-1 -0.259e-1, v6 = -.1021 -0.832e-1, v7 = -0.15e-1 -0.87e-2, maximize);NLPSolve(13.946*E*x1, m = 58962 .58962, E = 2650 2650, e = 2.71 2.71, x1 = O0.5e-1 0.5e-1, t= 75 .75, dl = 0.4794e-1 0.4794e-1, d2 = 0.3168e-1 0.3168e-1, d3 = 0.20024e-1 0.20024e-1, d4 = 0.867e-2 0.867e-2, dS = 0.307e-2 0.307e-2, d6 = 0.250e-2 0.250e-2, d7 = 0.19e-3 0.19e-3, v1 = 91116 .91116, v2 = 490308 .490308, v3 = 1363 .1363, v4 = -0.2708e-1 -0.2708e-1, v5 = -0.3972e-1 -0.3972e-1, v6 = -0.9121e-1 -0.9121e-1, v7 = -0.1152e-1 -0.1152e-1); Tải trọng động dài hạn dạng hình chữ nhật > with(Optimization); tinh*chuyen*vi-x7max -daihan; print( output redirected -); # input placeholder tinh chuyen vị - x7max — daihan > x7 := -(177.4793600*(exp(-0.9051263884e-1*xi*sqrt(E/m)*t)*sqrt(1.-1.*xi42)*cos(0.9051263884e-1 *sqrt(E/m)* sqrt(1.1.*xiA2)*t)+exp(-0.9051263884e-1 *xi* sqrt(E/m)*t)*xi*sin(0.905 1263884e-1 *sqrt(E/m)*sqrt(1.-1.*xi42)*t)-1.*sqri(1.1.*xi42)))* P/(sqrt(1-xi42)*E)+(8.205228241 *(exp(-.2290562114*xi*sqrt(E/m)*t)*sqrt(1.1.*x142)*cos(.2290562114*sqrt(E/m)*sqrt(1 -1.*x142)*t)+exp(.22905621 14*xi*sqrt(E/m)*t)*xi* sin(.22905621 14*sqrt(E/m)*sqrt(1.-1.*xiA2)*t)-1.*sqrt(1.-1.*xi42)))*P/(sqrt(1-x142)*E)(2.275837633* (ex p(-.3884925692* xi*sqrt(E/m)*t)*sqrt(1.-1.*x142)*cos(.3884925692*sqrt(E/m)*sqrt(2.-1.*xi42)*t)+exp(.3884925692*xi*sqrt(E/m)*t)* xi*sin(.3884925692* sqrt(E/m)*sqrt(1.-1.*x142)*t)-] *sqrt(1.-].*x142)))*P/(sqrt(1xi42)*E)+(.4491510769*(exp(-.49962 14466*xi* sqrt(E/m)*t)*sqrt(1.-1.*xi42)*cos(.4996214466* sqrt(E/m)*sqrt(1.1.*xi42)*t)+exp(-.49962 14466*xi*sqrt(E/m)*t)*xi*sin(.49962 14466*sqrt(E/m)*sqrt(1.-1.*xi42)*t)-1.*sqrt(1.-1.*x142)))*P/(sqrt(1 xi42)*E)-(.3199173860*(exp(-.5778409787* xi* sqrt(E/m)*t)* sqrt(1 -1.*x142)*cos(.5778409787*sqrt(E/m)*sqrt(1 1.*xi42)*t)+exp(-.5778409787* xi* sqrt(E/m)*t)*xi*sin(.5778409787 *sqrt(E/m)*sqrt(1.-].*x142)*t)-1.*sqrt(1.-1.*xi42)))*P/(sqrt(1 - xi42)*E)+(0.101616406Se-1*(exp(-.7061776592* xi*sqrt(E/m)*t)*sqrt(] -1.*xi42)*cos(.7061776592* sqrt(E/m)*sqri(1.1.*xiA2)*t)+exp(-.7061776592*xi* sqrt(E/m)*t)*xi*sin(.7061776592*sqrt(E/m)* sqrt(1 -1.*xi42)*t)-1.*sqrt(1.-1.*x142)))*P/(sqri(1xi42)*E)-(0.3179459884e-3* (exp(-.8470875210*xi*sqrt(E/m)*t)*sqrt(] -1.*xi42)*cos(.8470875210*sqrt(E/m)*sqrt(1.1.*xj42)*t)+exp(-.84708752 10*xi*sqrt(E/m)*t)*xi* sin(.8470875210*sqrt(E/m)*sqrt(1.-1.*x142)*t)-1.*sqrt(1.-1.*x142)))*P/(sqri(1xi42)*E); NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi = 0.le-1 .1, P= 22.5 27.5, t =0 1, maximize); NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi = O.Je-1 .1, P= 23.125 26.88, t = 513 513, maximize); NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi =0.1e-1 .1, P= 23.75 26.25, t= 513 .513, maximize); NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi =0.le-1 .1, P= 24.375 25.63, t= 513 .513, maximize); > NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi =0.le-1 1, P= 22.5 27.5, t= 513 .513); NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi= 0.le-1 .1, P= 23.125 26.88, t= 513 513); NLPSolve(x7, m = 53066 .64858, E = 2385 2915, xi = 0.1e-1 .1, P= 23.75 26.25, t= 513 .513); NLPSolve(x7, m = 53066 .64858, E = 2385 2915, x1 = 0.1e-] 1, P= 24.375 25.63, t= 513 513); > NLPSolve(x7, m = 58962 .58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25, t= 513 .513, maximize); > with(Optimization); tinh*chuyen* vi*theo*lat*cat*thoi* gian-x7-daihan; print( output redirected °); # input placeholder tinh chuyen vi theo lat cat thoi gian - x7 — daihan %t=0.0; NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi = 0.1le-] 1, P= 22.5 27.5,t=0 0); 1, P= 22.5 27.5, t=0 0, maximize); NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi=0.1e-1 NLPSolve(x7, m = 58962 58962, E = 2650 2650, x1 = 0.5e-1 0.5e-1, P=25 25,t =0 0, maximize); %t=0.2; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, %t=0.4; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, đ%t=0.6; m = 53066 64858, E = 2385 2915, xi=0.le-1 1, P= 22.5 27.5, t=0.2 0.2); m = 53066 64858, E = 2385 2915, xi=0.le-1 .1,P=22.5 27.5, t= 0.2 0.2, maximize); m = 58962 58962, E = 2650 2650, xi = 0.Se-] 0.5e-1, P= 25 25, t = 0.2 0.2, maximize); m = 53066 64858, E = 2385 2915, x1 = le-1 .1, P= 22.5 27.5, t= 0.4 0.4); m = 53066 64858, E = 2385 2915, x1 = O0.le-1 1, P= 22.5 27.5, t= 0.4 0.4, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.Se-1, P= 25 25,1 = 0.4 0.4, maximize); NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi = 0.1le-1 1, P= 22.5 27.5, t = 0.6 0.6); NLPSolve(x7, m = 53066 64858, E = 2385 2915, x1 = O.le-1 1, P= 22.5 27.5, t= 0.6 0.6, maximize); NLPSolve(x7, m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25, t= 0.6 0.6, maximize); #%t=0.8; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, %t=1.0; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, đt=] 2; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, %t=1.4; NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi = O.}e-] 1, P= 22.5 27.5, t= 0.8 0.8); m = 53066 64858, E = 2385 2915, xi=O0.}e-1 .1, P= 22.5 27.5, t= 0.8 0.8, maximize); m = 58962 58962, E = 2650 2650, xi = 0.Se-] 0.5e-1, P= 25 25, t= 0.8 0.8, maximize); m = 53066 64858, E = 2385 2915, xi=0.1e-1 1, P= 22.5 27.5, t= 1.0 1.0); m = 53066 64858, E = 2385 2915, xi = 0.1e-1 .1,P=22.5 27.5, t=1.0 1.0, maximize); m = 58962 58962, E = 2650 2650, xi = 0.Se-1 0.5e-1, P= 25 25,t = 1.0 1.0, maximize); m = 53066 64858, E = 2385 2915, x1 = 0.1e-1 1, P= 22.5 27.5, t= 1.2 1.2); m = 53066 64858, E = 2385 2915, xi=0.le-1 1, P= 22.5 27.5, t= 1.2 1.2, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P=25 25, t= 1.2 1.2, maximize); m = 53066 64858, E = 2385 2915, x1 = 0.1e-1 1, P= 22.5 27.5, t= 1.4 1.4); 21 NLPSolve(x7, NLPSolve(x7, %t=1.6; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, %t=].8; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, %t=2.0; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, Jot=2.2; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, %t=2.4; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, %t=2.6; NLPSolve(x7, NLPSolve(x7, NLPSolve(x7, Jort=2.8; m = 53066 64858, E = 2385 2915, x1 = 0.1e-1 » J, P= 22.5 27.5,t=1.4 1.4, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25,t=1.4 1.4, maximize); m = 53066 64858, E = 2385 2915, xi = 0.1e-1 .1,P=22.5 27.5,t= 1.6 1.6); m = 530666 64858, E = 2385 2915, xi=0.1e-1 .1,P=22.5 27.5,t= 1.6 1.6, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-], P=25 25,t=1.6 1.6, maximize); m = 53066 64858, E = 2385 2915, xi = O.Je-1 A, P=22.5 27.5, t= 1.8 1.8); m = 53066 64858, E = 2385 2915, xi=0.le-1 ,P=22.5 27.5, t=1.8 1.8, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P=25 25,t=1.8 1.8, maximize); m = 53066 64858, E = 2385 2915, xi=0.1le-1 A, P=22.5 27.5, t=2.0 2.0); m = 53066 64858, E = 2385 2915, xi =0.1le-1 A, P= 22.5 27.5, t=2.0 2.0, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25, t=2.0 2.0, maximize); m = 53066 64858, E = 2385 2915, xi = O.Je-] A, P= 22.5 27.5, t=2.2 2.2); m = 53066 64858, E = 2385 2915, xi = 0.le-1 1,P=22.5 27.5, t=2.2 2.2, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P=25 25,t=2.2 2.2, Maximize); m = 53066 64858, E = 2385 2915, xi=0.le-1 ,P=22.5 27.5, t=2.4 2.4); m = 53066 64858, E = 2385 2915, x1 = 0.1e-1 A, P= 22.5 27.5,t=2.4 2.4, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25,t=2.4 2.4,maximize); m = 53066 64858, E = 2385 2915, x1 =0.1e-1 A, P=22.5 27.5, t=2.6 2.6); m = 53066 64858, E = 2385 2915, xi = 0.le-1 A,P=22.5 27.5,t=2.6 2.6, maximize); m = 58962 58962, E = 2650 2650, x1 = 0.5e-] 0.5e-1, P= 25 25,t=2.6 2.6, maximize); NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi=0.le-1 .1,P=22.5 27.5,t=2.8 2.8); NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi=0.]e-] .1,P=22.5 27.5,t=2.8 2.8, maximize); NLPSolve(x7, m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25,1=2.8 2.8, maximize); %t=3.0; NLPSolve(x7, m = 53066 64858, E = 2385 2915, x1 = 0.1le-1 A, P=22.5 27.5, t=3.0 3.0); NLPSolve(x7, m = 53066 64858, E = 2385 2915, xi = 0.1le-1 ,P=22.5 27.5, t= 3.0 3.0, maximize); NLPSolve(x7, m = 58962 58962, E = 2650 2650, xi = 0.S5e-1 0.5e-1, P=25 25, t = 3.0 3.0, maximize); > with(Optimization); tinh*chuyen*vi-x ] max-daihan; print( output redirected °); # input placeholder tỉnh chuyen vị - x) max - daihan > x1 := -(28.76758560* (exp(-0.905 1263884e-1*xi* sqrt(E/m)*t)*sqrt(1 -1.*x142)*cos(0.9051263884e-1 *sqri(E/m)*sqrt(1.1.*xi42)*t)+exp(-0.9051263884e-1 ¥xi*sqrt(E/m)*t)*xi*sin(0.905 1263884e- *sqrt(E/m)*sqrt(1.-1.*xi42)*t)-1.*sqrt(1.1.*xi42)))*P/(sqrt(] -xi42)*E)-(3.652554635* (exp(-.2290562114*xi*sqrt(E/m)*t)*sqrt(1.1.*x1^2)*cos(.22905621 14*sart(E/m)*sqrt(1.-1.*xi^2)*t)+exp(.2290562114*xi*sqrt(E/m)*t)*xi*sin(.22905621 14*sqrt(E/m)*sqri(1.-1 *xI1^2)*9)-1.*sart(1.-1.*xI^2)))*P/Gqrt(1-x1^2)*E)(1.226733672*(exp(-.3884925692*xi*sqart(E/m)*t)*sqart(1.-1.*x1^2)*cos(.3884925692*sart(E/m)*sart(1.-1.#*x1^2)*)+exp(.3884925692* xi*sqrt(E/m)*t)* xi*sin(.3884925692* (.5002981487*(exp(-.49962 14466* sqrt(E/m)*sqrt(1.-1.*xi42)*t)-1 *sqrt(1.-1.*x142)))*P/(sqrt(1-x142)*E)- xi*sqrt(E/m)*t)*sqrt(1.-1.*xi42)*cos(.49962 14466* sqrt(E/m)*sqrt(1.-1.*xi42)*t)+exp(- 49962 14466*xi*sqrt(E/m)*t)*xi*sin(.49962 14466* sqrt(E/m)*sqrt(1.-1.*xi42)*t)-1.*sqrt(1 -1.*xi42)))*P/(sqrt( -xi42)*E)(.1981368670* (ex p(-.5778409787* xi* sqrt(E/m)*t)*sqrt(1 - *xi42)*cos(.5778409787* sqrt(E/m)*sqrt(1.-1.*xi42)*t)+exp(- 5778409787* xi* sqrt(E/m)*t)*xi*sin(.5778409787* sqrt(E/m)*sqrt(1.-1.*xi42)*t)-1 *sqrt(1.-1.*xi%2)))*P/(sqrt(1-x142)*E)(.2584760768* (exp(- 706 1776592* xi* sqrt(E/m)*t)*sqrt(1 -1.*x142)*cos(.7061776592*sqrt(E/m)* sqrt(1 -1.*xi42)*t)+exp(.1061776592* xi* sqr(E/m)*t)*x1*sin(.7061 776592*sqrt(E/m)*sqrt(] -1.*xi42)*t)-1.*sqrt(1.-1.*x142)))*P/(sqrt(1 -xi42)*E)- (0.496806023 le-1*(exp(-.8470875210*x1*sqrt(E/m)*t)*sqrt(1.-1.*x142)*cos(.8470875210*sqrt(E/m)*sqrt(1.-1.*x142)*1t)+exp(.84708752 10*xi*sqrt(E/m)*t)*xi*sin(.8470875210*sqrt(E/m)*sqrt(1.-1.*xi42)*t)-1 ¥sqrt(1 -1.*xi42)))*P/(sqrt(1-x142)*E); NLPSolve(x1, m = 53066 .64858, E = 2385 2915, x1 = 0.le-1] .1, P= 22.5 27.5, t= 1, maximize); NLPSolve(x1, m = 53066 .64858, E = 2385 2915, xi = 0.]e-] .1, P= 23.125 26.88, t = 617 .617, maximize); NLPSolve(x1, m = 53066 .64858, E = 2385 2915, xi = 0.1le-1 .1, P= 23.75 26.25, t= 617 .617, maximize); NLPSolve(x1], m = 53066 .64858, E = 2385 2915, xi = 0.1le-1 .1, P= 24.375 25.63, t = 617 .617, maximize); NLPSolve(x1, m = 53066 .64858, E = 2385 2915, x1 = 0.1e-1 .1, P= 22.5 27.5, t= 617 .617); NLPSolve(x1, m = 53066 64858, E = 2385 2915, xi = 0.1e-1 .1, P= 23.125 26.88, t = 617 .617); NLPSolve(x1, m = 53066 .64858, E = 2385 2915, xi = 0.le-] .1, P= 23.75 26.25, t = 617 .617); NLPSolve(xJ, m = 53066 .64858, E = 2385 2915, xi=0.1e-1 .1,P = 24.375 25.63, t = 617 .617); NLPSolve(x1, m = 58962 .58962, E = 2650 2650, xi = 0.Se-1 0.5e-1, P= 25 25, t = 617 617, maximize); > with(Optimization); tinh*chuyen*vi*theo*lat*cat*thon* gian-x -daihan; printC output redirected °); # input placeholder tinh chuyen vi theo lat cat thoi gian - x1 - dathan %1=0.0; NLPSolve(x1, NLPSolve(x1, NLPSolve(x1, %t=0.2; NLPSolve(x1, NLPSolve(x1, m = 53066 64858, E = 2385 2915, xi = O.le-1 .1, P= 22.5 27.5,t=0 0); m = 53066 .64858, E = 2385 2915, x1 = 0.le-1 .1, P= 22.5 27.5, t= 0, maximize); m = 58962 58962, E = 2650 2650, xi = 0.5e-} 0.S5e-1, P= 25 25, t =0 0, maximize); m = 53066 64858, E = 2385 2915, xi = 0.le-1 .1, P= 22.5 27.5, t= 0.2 0.2); 1, P= 22.5 27.5, t=0.2 0.2, maximize); m = 53066 64858, E = 2385 2915, xi = 0.1e-1 22 NLPSolve(x1,m %t=0.4; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=0.6; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=0.8; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=1.0; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=1.2:; NLPSolve(x1, m NLPSolve(x], m NLPSolve(x1, m %Mt=1.4; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x], m %t=1.6; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=1.8; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=2.0; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=2.2; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=2.4; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=2.6; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=2.8; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m %t=3.0; NLPSolve(x1, m NLPSolve(x1, m NLPSolve(x1, m = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1,P=25 25,t=0.2 0.2, maximize); = 53066 64858, E = 2385 2915, xi=0.le-] .1,P=22.5 27.5,t=0.4 0.4); = 53066 64858, E = 2385 2915, xi=0.le-1 .1,P=22.5 27.5,t=0.4 0.4, maximize); = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.Se-1, P=25 25,1=0.4 0.4, maximize); = 53066 .64858, E = 2385 2915, xi=0.le-1 .1, P= 22.5 27.5, t= 0.6 0.6); = 53066 .64858, E = 2385 2915, x1 = 0.1e-1 .1, P= 22.5 27.5, t= 0.6 0.6, maximize); = 58962 58962, E = 2650 2650, x1 = 0.5e-1 0.5e-1, P=25 25,t =0 0, maximize); = 53066 64858, E = 2385 2915, xi=0.1e-1 .1,P=22.5 27.5,t=0.8 0.8); = 53066 64858, E = 2385 2915, x1 =0.1e-1 .1,P=22.5 27.5,t=0.8 0.8, maximize); = 58962 58962, E = 2650 2650, x1 = 0.5e-1 0.5e-1,P=25 25,t=0.8 0.8, maximize); = 53066 64858, E = 2385 2915, xi=0.Je-l .1,P=22.5 27.5,t=1.0 1.0); = 53066 64858, E = 2385 2915, x1=0.le-] .1,P=22.5 27.5,t=1.0 1.0, maximize); = 58962 .58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25, t = 1.0 1.0, maximize); = 53066 .64858, E = 2385 2915, xi = 0.le-1 A,P=22.5 27.5,t=1 1.2); = 53066 64858, E = 2385 2915, xi =0.1e-1 A, P=22.5 27.5,t=1 1.2, maximize); = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1,P=25 25,1= 1.2 1.2, maximize); = 53066 64858, E = 2385 2915, xi = 0.1e-1 .1,P=22.5 27.5,t=1.4 1.4); = 53066 64858, E = 2385 2915, xi=0.le-1 1,P=22.5 27.5,t=1.4 1.4, maximize); = 58962 58962, E = 2650 2650, x1 = 0.5e-1 0.5e-1,P=25 25,t=1.4 1.4, maximize); = 53066 64858, E = 2385 2915, xi=0.1e-1 1, P=22.5 27.5,t=1.6 1.6); = 53066 64858, E = 2385 2915,xi=0.le-] 1, P=22.5 27.5,t=1.6 1.6, maximize); = 58962 98962, E = 2650 2650, xi = 0.5e-1 0.5e-1,P=25 25,t=1.6 1.6, maximize); = 53066 64858, E = 2385 2915, x1 = 0.le-1 1, P= 22.5 27.5,t= 1.8); = 53066 64858, E = 2385 2915, xi =0.le-1 .1,P=22.5 27.5,t = 1.8, maximize); = 58962 58962, E = 2650 2650, xi = 0.Se-] 0.5e-1, P= 25 25,t= 1.8 1.8, maximize); = 53066 64858, E = 2385 2915, xi = 0.1le-1 1, P= 22.5 27.5, t=2.0 2.0); = 53066 64858, E = 2385 2915, xi=0.le-1 ,P=22.5 27.5,t=2.0 2.0, maximize); = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25, t = 2.0 2.0, maximize); = 53066 64858, E = 2385 2915, x1 = 0.1le-1 1, P= 22.5 27.5, t=2.2 2.2); = 53066 64858, E = 2385 2915, xi = O.1e-1 -1, P= 22.5 27.5, t=2.2 2.2, maximize); = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25,t=2.2 2.2, maximize); = 53066 .64858, E = 2385 2915, x1=0.Je-1 .1, P= 22.5 27.5, t =2.4 2.4); = 53066 .64858, E = 2385 2915, xi =0.le-1 .1, P= 22.5 27.5, t=2.4 2.4, maximize); = 58962 58962, E = 2650 2650, xi = 0.Se-1 0.5e-1, P= 25 25,t=2.4 2.4, maximize); = 53066 64858, E = 2385 2915, xi = 0.1e-1 A, P= 22.5 27.5, t=2.6 2.6); = 53066 64858, E = 2385 2915, x1 = O.1e-1 ,P=22.5 27.5,t=2.6 2.6, maximize); = 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25,1=2.6 2.6, maximize); = 53066 64858, E = 2385 2915, xi = 0.1le-1 1, P=22.5 27.5, 1=2.8 2.8); = 53066 64858, E = 2385 2915, xi = 0.1e-1 1, P= 22.5 27.5,t=2.8 2.8, maximize); = 58962 58962, E = 2650 2650, xi = 0.S5e-1 0.5e-1, P=25 25,t=2.8 2.8, maximize); = 53066 .64858, E = 2385 2915, xi=0.Je-1 1, P= 22.5 27.5, t= 3.0 3.0); = 53066 .64858, E = 2385 2915, x1 = 0.le-1 .1, P= 22.5 27.5, t = 3.0 3.0, maximize); = 58962 .58962, E = 2650 2650, xi = 0.S5e-1 0.5e-1, P= 25 25, t = 3.0 3.0, maximize); > with(Optimization); tinh*momen*max-daihan; print( output redirected °); # input placeholder tinh momen max — dathan > X1:=-(28.76758560*(exp(-0.9051263884e-1 *xi*sartŒ/m)*t)*sart(1.-1.*xi^2)*cos(0.9051263884e-1*sartŒ/m)*sqrt(1.1.*x1^2)*t)+exp(-0.9051263884e-1 *xi*sart(E/m)*Ð)*xi*sin(0.9051263884c-1 *sart(E/m)*sart(1.-].*xi^2)*0)-1.*sart(1.1.*x1^2)))*P/(sqrt(1-xi^2)*E)-(3.652554635*%(exp(-.22905621 14*xi*sart(/m)*#t)*sqrt(1.1.*xi42)*cos(.2290562114*sqrt(E/m)*sqrt(1.-1.*xiA2)*t)+ex p(- ` 2290562114*xi*sqrt(E/m)*t)*xi*sin(.22905621 14*sqrt(E/m)*sqrt(1.-1.*xi142)*t)-1.*sqrt(1.-1.*xi42)))*P/(sqrt(1-x142)*E)(1.226733672*(exp(-.3884925692*xi*sqrt(E/m)*t)*sqrt(1.-1.*x142)*cos(.3884925692*sqrt(E/m)*sqrt(1.-1.*x142)*t)+exp(.3884925692*xi*sqrt(E/m)*t)*xi* sin(.3884925692*sqrt(E/m)*sqrt(1.-1.*xi42)*t)-] *sqrt(] -].*xi42)))*P/(sqrt(1 -xi42)*E)(.500298 487* (exp(-.49962 14 * xi*sqrt(E/m)*t)*sqrt(1.-1.*x142)*cos(.49962 14466* sqrt(E/m)*sqrt(1.-1.*xi42)*t)+exp(49962 14466*xi*sqrt(E/m)*t)*xi*sin(.49962 14466*sqrt(E/m)*sqrt(].-1.*x142)*t)-1.*sqrt(1.-1.*x142)))*P/(sqrt(1 -x142)*E)(.1981368670*(exp(-.5778409787* xi*sqrt(E/m)*t)*sqrt(1.-1.*xi42)*cos(.5778409787* sqrt(E/m)*sqrt(1.-1.*x142)*t)+ex p(.5778409787*xi*sart(E/m)*U*xi*sin(.5778409787*sqrt(E/m)*sart(1.-1.*x1^2)*t)-1.*sqrt(1.-1.*xI^2)))*P/art(1-xI^2)*E)(.2584760768*(exp(-.206 1776592*xi*sart(E/m)*t)*sqart(1.-1.*xi^2)*cos(.70617776592*sart(E/m)*sart(1.- Ì.*x1^2)*t)+exp(.?061776592*xi*sqrt(E/m)*t)*xi*sin(.70617776592*sart(E/m)*sqrt(1.-1.*xi^2)*t)-1.*sgrt(1.-1.*x1^2)))*P/(sqrt(-xi^2)*E)(0.496806023 le-1*(exp(-.8470875210*xi*sqrt(E/m)*t)*sqrt(1 -1.*x142)*cos(.8470875210* sqrt(E/m)*sqrt(1.-1.*xi42)*t)+exp(.8470875210*xi*sqrt(E/m)*t)*xi*sin(.84 708752 10*sqrt(E/m)*sqrt(J -1.*x142)*t)-1 *sqrt(1.-1.*x242)))*P/(sqrt(1-x142)*E)); 23 NLPSolve(13.946*E*x1 ,m=0.53066 0.64858, E=2385 2915,xi=0.01 0.1,P=22.5 27.5,t=0.617 0.617, maximize); NLPSolve(13.946*E*x1 1,m=0.53066 0.64858 E=2385 2915,xi=0.01 0.1,P=23.125 26.88 t=0.617 0.617, maximize); NLPSolve(13.946*E*x1 ,m=0.53066 0.64858,E=2385 2915,xi=0.01 0.1,P=23.75 26.25,t=0.617 0.617, maximize); NLPSolve(13.946*E*x1,m=0.53066 0.64858,E=2385 2915,xi=0.01 0.1,P=24.375 25.63,1=0.617 0.617, maximize); ; > NLPSolve(13.946*E*x1, m = 53066 .64858, E = 2385 2915, xi =0.le-] 1, P= 22.5 27.5, t= 617 .617); NLPSolve(13.946*E*x1, m = 53066 64858, E = 2385 2915, x1= 0.le-1 J, P= 23.125 26.88, t= 617 .617); 1, P= 23.75 26.25, t= 617 .617); NLPSolve(13.946*E*x1, m = 53066 .64858, E = 2385 2915, x1= 0.1e-1 NLPSolve(]3.946*E*x1, m = 53066 .64858, E = 2385 2915, xi==0.le-] 1, P = 24.375 25.63, t = 617 .617); > NLPSolve(13.946*E*x1, m = 58962 .58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P= 25 25, t = 617 617, maximize); > %Tinh*momen theo Jat cat thoi gian; print( output redirected ” ); # input placeholder % Tinh momen theo Jat cat thoi Jot=0.0; NLPSolve(13.946*E*x1, m = 53066 64858, E = NLPSolve(13.946*E*x1, m = 53066 64858, E = NLPSolve(13.946*E*x1, m= 58962 58962, E = %t=0.2; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x], %t=0.4; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x], m m m m m m %t=0.6; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=0.8; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=1.0; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=1.2; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=1.4; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=1.6; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1 NLPSolve(13.946*E*x1, Ft=1.8; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, m m m m m m m m m m m m gian 2385 2915, xi= O.le-1 .l,P=22.5 27.5,t =0 0); 2385 2915, xi= O.le-1 .1,P=22.5 27.5,t = 0, maximize); 2650 2650, xi = 0.5e-1 0.5e-1, P=25 25,t= 0, maximize); 53066 64858, E = 2385 2915, xi = Q.le-1 1, P=22.5 27.5,t=0.2 0.2); 53066 64858, E = 2385 2915,xi=0.l1e-l 1, P= 22.5 27.5,t=0.2 0.2, maximize); 58962 58962, E = 2650 2650, xi= 0.5e-1 0.5e-1, P= 25 25,t=0.2 0.2, maximize); 53066 64858, E = 2385 2915, xi= O.le-1 1, P= 22.5 27.5,t=0.4 0.4); 53066 64858, E = 2385 2915, x1i= O.le-1 1, P=22.5 27.5, t=0.4 0.4, maximize); 58962 58962, E = 2650 2650, xi= 0.5e-1 0.5e-1, P= 25 25,1 = 0.4 0.4, maximize); 53066 64858, E = 2385 2915, x1 = 0.le-1 1,P = 22.5 27.5,t=0.6 0.6); 53066 64858, E = 2385 2915, xi= 0.le-1 ,P=22.5 27.5,t=0.6 0.6, maximize); 38962 58962, E= 2650 2650, xi= 0.5e-1 0.S5e-1, P= 25 25, t=0 0, maximize); 33066 64858, E = 2385 2915, xi= le-1 1, P= 22.5 27.5,t=0.8 0.8); 53066 64858, E = 2385 2915, xi = 0.le-] 1, P= 22.5 27.5,t=0.8 0.8, maximize); 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P=25 25,1=0.8 0.8, maximize); 33066 64858, E = 2385 2915, xi = 0.le-1 1,P= 22.5 27.5,t=1.0 1.0); 53066 64858, E = 2385 2915, xi =0.le-l ,P=22.5 27.5,t=1.0 1.0, maximize); 58962 58962, E = 2650 2650, xi = 0.5e-1 0.5e-1, P=25 25,t=1.0 1.0, maximize); 53066 64858, E = 2385 2915, xi =0.le-]l wL, P= 22.5 27.5,t=1.2 1.2); 53066 64858, E = 2385 2915, xi =0.le-]l 1, P=22.5 27.5,t=1.2 1.2, maximize); 58962 58962, E = 2650 2650, xi =0.5e-] .0.5e-1,P=25 25,t=1.2 1.2, maximize); m 53066 64858, E = 2385 2915, xi = 0.le-1 .1,P=22.5 27.5,t=14 1.4); 53066 64858, E = 2385 2915, xi= O.Je-1 .1,P=22.5 27.5,t= 1.4 1.4, maximize); 58962 58962, E = 2650 2650, xi = 0.Se-1 0.5e-1,P=25 25,t=1.4 1.4, maximize); m ,m m 53066 64858, E = 2385 2915, xi = 0.1le-1 .1,P=22.5 27.5,t=1.6 1.6); 53066 64858, E = 2385 2915, xi = 0.1e-1 1, P=22.5 27.5,t=1.6 1.6, maximize); 1.6, maximize); 58962 58962, E = 2650 2650, x1 = 0.5e-1 0.5e-1, P=25 25,t= m m m= 53066 64858, E = 2385 2915, xi = 0.le-1 .1,P=22.5 27.5,t=1.8 m= 53066 64858, E = 2385 2915, xi = O.le-1 .l,P=22.5 27.5,t=1.8 NLPSolve(13.946*E*x1, E = 2650 2650, xi= 0.Se-1 0.5e-1, P= 25 25, t= 1.8 1.8, maximize); %t=2.0; NLPSolve(13.946*E*x1, m = 53066 64858, E = 2385 2915, xi= 0.le-]l 1, P=22.5 27.5, t=2.0 NLPSolve(13.946*E*x1, m = 53066 64858, E = 2385 2915, x1 = 0.le-1 A, P=22.5 27.5,t=2.0 NLPSolve(13.946*E*x1, m = 58962 58962, E = 2650 2650, xi= 0.Se-1 0.5e-1,P=25 25,t=2.0 %t=2.2; NLPSolve(13.946*E*xI, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=2.4; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x], NLPSolve(13.946*E*x1, %t=2.6; NLPSolve(13.946*E*xI, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=2.8: NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, %t=3.0; NLPSolve(13.946*E*x1, NLPSolve(13.946*E*x1, NLPSolve(13.946* E*x1, END; m m m m m m m m m m m m m m m 1.8); 1.8, maximize); 2.0); 2.0, maximize); 2.0, maximize); 53066 64858, E = 2385 2915, xi= O.le-1 .1,P=22.5 27.5,t=2.2 2.2); 53066 64858, E = 2385 2915, xi= 0.le-l 1,P=22.5 27.5,t=2.2 2.2, maximize); 58962 58962, E = 2650 2650, xi= 0.Se-1 0.5e-1,P=25 25,t=2.2 2.2, maximize); 53066 64858, E = 2385 2915, xi = 0.Je-] .1,P=22.5 27.5,t=2.4 2.4); 53066 64858, E = 2385 2915, xi= O.le-1 .1,P=22.5 27.5,t=2.4 2.4, maximize); 58962 58962, E= 2650 2650, xi = 0.5e-1 0.5e-1,P=25 25,t=2.4 2.4, maximize); 53066 64858, E = 2385 2915, xi = 0.le-] .1,P=22.5 27.5,t=2.6 2.6); 3066 64858, E = 2385 2915, xi =0.le-1 1,P=22.5 27.5,t=2.6 2.6, maximize); 58962 58962, E = 2650 2650, xi = 0.Se-1 0.5e-1,P=25 25,t=2.6 2.6, maximize); 53066 64858, E = 2385 2915, xi = 0.le-1 .l,P=22.5 27.5,t=2.8 2.8); 53066 64858, E = 2385 2915, x1= O.le-] .l,P=22.5 27.5,t=2.8 2.8, maximize); 58962 58962, E = 2650 2650, xi= 0.5e-1 0.5e-1,P=25 25,t=2.8 2.8, maximize); 33066 64858, E = 2385 2915, xi= O.le-] .1,P=22.5 27.5,t=3.0 3.0); 53066 64858, E = 2385 2915, x1 =(0.le-l 1,P=22.5 27.5,t =3.0 3.0, maximize); 58962 58962, E = 2650 2650, x1 = 0.Se-1 0.5e-1, P= 25 25, t = 3.0 3.0, maximize);

Ngày đăng: 20/12/2023, 15:03

Xem thêm:

TÀI LIỆU CÙNG NGƯỜI DÙNG

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

TÀI LIỆU LIÊN QUAN

w