BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ LÊ NGỌC TRUỜNG PHÂN TÍCH ỨNG XỬ PHI TUYẾN DẦM COMPOSITE CHỊU TÁC DỤNG TẢI TRỌNG ÐIỀU HOÀ DI ÐỘNG NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DÂN DỤNG & CƠNG NGHIỆP - 60580208 SKC005129 Tp Hồ Chí Minh, tháng 4/2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ LÊ NGỌC TRƯỜNG PHÂN TÍCH ỨNG XỬ PHI TUYẾN DẦM COMPOSITE CHỊU TÁC DỤNG TẢI TRỌNG ĐIỀU HOÀ DI ĐỘNG NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DÂN DỤNG VÀ CƠNG NGHIỆP - 60580208 Tp Hồ Chí Minh, tháng 04/2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ LÊ NGỌC TRƯỜNG PHÂN TÍCH ỨNG XỬ PHI TUYẾN DẦM COMPOSITE CHỊU TÁC DỤNG TẢI TRỌNG ĐIỀU HOÀ DI ĐỘNG NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH DÂN DỤNG VÀ CƠNG NGHIỆP - 60580208 Hướng dẫn khoa học: PGS.TS NGUYỄN TRUNG KIÊN Tp Hồ Chí Minh, tháng 04/2017 BQ GIAO DUC VA DAO TAO TRTIOIIG DAI Hoc sT'PHAM rf ruuAr THAI\H PHO-HO CHI MINH 56: 1909/QD-DHSPKTTp H6 Chi Minh, ngdy 29 thdng ndm 2016 QUYET DINH giao tld tdri lu$n vln t6t nghiQp vdr nguoi hufng ddn_Kh6a 2015 B HIEU TRITOI.{G TRUONG DAI HQC SII PHAM KY THUAT TP HO CHf MINH Cdn cti Quytit ttinh sd ll8l2000lQD-TTg ngdy 10 thSng 10 ndm 2000 cta Thri VG viQc tudng Chinh phri v0 viQc thay AOi tO chric cria Dai hqc Qu6c gia TP HO Chi Minh, t6ch Trudng Dai hgc Su ph4m K! thuft TP HO Chi Minh tryc thuQc B0 Gi6o dpc vd Ddo t4o Cdn cri Quy€t dinh sd 70l20l4lQD-TTg ngdy vd ban hanh DiAu lQ trubng Dai hqc t0ll2l20l4 cria Thri tu6ng Chinh phri Cdn cri Thdng tu s6 I112}I41TT-BGDDT ngdy 151512014 cria Ddo t4o vO viQc Ban hdnh Qui ch€ tldo tpo trinh d0 thac si;; Cdn cri vdo Bi€n b6n bdo vQ Chuy6n dA cria ngdnh trinh 86 Gi6o dpc vd K! thu4t xdy dgng c6ng & c6ng nghiQp vdo ngiy 1810812016,; X6t nhu cAu c6ng t6c vi ndng can b0; d6n dsng Xdt dA nghi cria Trucmg phong Edo t4o, Didu Giao dA tdi Lufln nrm "r" :iJrtlr"iill;, "u ngudi hu6ng ddn cao hsc 2016 cho: L0 Nggc Trwdng Kho6: 2015 -2017 B : Kj thuQt xdy dqng cAng ffinh ddn dyng & cdng nghiQp : PhAn fich hng xft phi nyAn dilm composite chlu ttic : dqtng tdi trgng didu hba di ilQng Nguoi hudng ' Di6u ddn : PGS TS Nguy1n Trung KiAn Giao cho Phong Ddo t4o qu6n ly, thUc hiqn theo ihing Qui chti ddo t4o trinh dQ th4c si ctra BQ Gi6o dpc & Ddo t4o ban hdnh Diiiu Tru&ng c6c don vi, phdng Ddo t4o, c6c Khoa quin ngdnh cao hgc vd c6c A Ong (Be) co tOn t4i Quy0t clfnh c6 hiQu lyc kC tir ngdy k!./ Noi nhfin: - BGH (AC Uie0; - Nhu'di6u2,3; - Luu: VT, SDH ;ii BQ crAo DUC & DAo rAo TRTIONG DAr HQg Sr/PrrlM xY ceNG noA xA Her cHU rnu4r THANH PHO HO CHI MINH DQc lSp xcnia vrfr xavr - Tq - Hgnh Phrfc xAc NHAN cua cAN Be HUONG DAN H9 vd t6n hgc vi0n: L0 Ngqc 'Chuy6n ngdnh: Truong Msfry: 1580822 Ki thu$t xiy dgng cdng trinh DD&CNKh6a: 20158 TOn dd tdi: PhAn tich rfrng xfr phi tuy6n dAm composite chlu tfli trgng tlidu hda di tlQng Hgc viAn ild hodn thdnh LVTN theo ilrtng yAu cdu vA nQi dung vd hinh thtbc (theo qui tlinh) crta mQt luQn vdn thgc si Gif,ng viOn hurimg dfin @j,A ghi rp hpc tAn) WfPrTla" j'i %Ӝ*,È2'Ө&9¬ Ҥ27Ҥ2 TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH PHIẾU NHẬN XÉT LUẬN VĂN THẠC SỸ (Dành cho JiảnJ viên phản biện) Tên đề tài luận văn thạc sỹ: 3KkQWtFKӭQJ[ӱSKLWX\ӃQGҫPFRPSRVLWHFKӏXWiFGөQJWҧLWUӑQJÿLӅXKzDGL ÿӝQJ MSHV: 1580822 Tên tác Jiả: /Ç1*Ӑ&75ѬӠ1* Ngành: WKXұW[k\GӵQJF{QJWUuQKGkQGөQJYjF{QJQJKLӋS Khóa: 2016 Định hướng: ӬQJGөQJ Họ tên nJười phản biện: 3*6761JX\ӉQ9 Q+LӃX Cơ quDn cônJ tác: ĈҥLKӑF.LӃQ7U~F7S+&0 Điện thoại liên hệ: 0938123299 I Ý KIẾN NHẬN XÉT Về hình thức & kết cấu luận văn: +uQKWKӭFWUuQKEj\WѭѫQJÿӕLU}UjQJV~FWtFK ӃWFҩXOXұQY QPҥFKOҥFFyWUuQKWӵ Về nội dunJ: 2.1 Nhận xét tính khoa học, rõ ràng, mạch lạc, khúc chiết luận văn &iFKGLӉQÿҥWWURQJOXұQY QFyWtQKNKRDKӑFYjNKiU}UjQJGӉWKHRG}LFiFÿӅPөFYjQӝLGXQJWUuQK bày 2.2 Nhận xét đánh giá việc sử dụng trích dẫn kết NC người khác có qui định hành pháp luật sở hữu trí tuệ /XұQY QFyVӱGөQJYjWUtFKGүQWjLOLӋXWKDPNKҧRÿ~QJTX\ÿӏQKKLӋQKjQKFӫDSKiSOXұWVӣKӳXWUtWXӋ 2.3 Nhận xét mục tiêu nghiên cứu, phương pháp nghiên cứu sử dụng LVTN 3KѭѫQJSKiSQJKLrQFӭXYjFiFKWLӃSFұQSKKӧSYӟLPөFWLrXÿӅUDFӫDÿ ӅWjL 2.4 Nhận xét Tổng quan đề tài 3KҫQ7әQJTXDQFӫDÿӅWjLQKuQFKXQJÿҥW\rXFҫX7X\QKLrQSKҫQPөFWLrXFҫQJLӟLKҥQWKHRÿ~ QJQӝL GXQJÿmÿѭӧFQJKLrQFӭXF QJQKѭWtQKPӟLFKѭDQrXU}VӵNKiFELӋWVRYӟLQKӳQJQJKLrQFӭXWUѭӟFÿk\ [HPWUDQJSKҫQPөFWLrXYjWtQKPӟLӣEҧQJ 2.5 Nhận xét đánh giá nội dung & chất lượng LVTN 1ӝLGXQJYjFKҩWOѭӧQJFӫDÿӅWjLÿҥW\rXFҫXFӫDPӝWOXұQY QWKҥFV 2.6 Nhận xét đánh giá khả ứng dụng, giá trị thực tiễn đề tài KҧQ QJӭQJGөQJYjJLiWUӏWKӵFWLӉQFӫDÿӅWjLNK{QJFDR 2.7 Luận văn cần chỉnh sửa, bổ sung nội dung (thiết sót tồn tại): 3KҫQWyPWҳWWLӃQJ$QK$EVWUDFWFҫQYLӃWOҥLFKRSKKӧSYӟLQӝLGXQJWyPWҳWWLӃQJ9L ӋWYjY QSKҥP WLӃQJ$QK - Trang 51NӃWOXұQNKLWӍVӕ/K20WKuGҫPFyÿӝY}QJQKӓWKuFyWKӇEӓTXDҧQKKѭӣQJFӫDELӃQGҥQJOӟQ GRӭQJ[ӱWX\ӃQWtQKYjSKLWX\ӃQOjQKѭQKDX LӅXQj\NK{QJSKKӧSYӟLNӃWTXҧVӕWUrQ+uQK5a 5b &iFKuQKYӁFKӭDKDLÿӗWKӏFҫQÿiQKVӕWKӭWӵDE II CÁC VẤ0Ề CẦN LÀM RÕ (Các câu hỏi giảng viên phản biện) *LҧLWKtFKELӃQGҥQJ9RQ.DUPDQNKiFELӃQGҥQJOӟQ9RQ.DUPDQӣWUDQJ17 sao? *LҧLWKtFKYuVDRӭQJ[ӱSKLWX\ӃQYjWX\ӃQWtQKFyKuQKGҥQJJLӕQJQKDX+uQK3.3 3.4)? 7iFJLҧJLҧLWKtFKYuVDRNKL/K10WKuOӡLJLҧLSKLWX\ӃQWUQJOӡLJLҧLWX\ӃQWtQK " TT Mục đánh Jiá 7tQKNKRDKӑFU}UjQJPҥFKOҥFNK~FFKLӃWWURQJOXұQY Q iQKJLiYLӋFVӱGөQJKRһFWUtFKGүQNӃWTXҧ1&FӫDQJѭӡLNKiFFy ÿӏQKKLӋQKjQKFӫDSKiSOXұWVӣKӳXWUtWXӋ 0өFWLrXQJKLrQFӭXSKѭѫQJSKiSQJKLrQFӭXVӱGөQJWURQJ/9 7әQJTXDQFӫDÿӅWjL iQKJLiYӅQӝLGXQJ FKҩWOѭӧQJFӫD/971 iQKJLiYӅNKҧQ QJӭQJGөQJJLiWUӏWKӵFWLӉQFӫDÿӅWjL 0Н± Øг0 (Giảng viên phản biện ghi rõ ý kiến “Tán thành luận văn” hay “Không tán thành luận văn”) 7iQWKjQKOXұQY Q 73+&0QJj\WKiQJQ P NJười nhận xét ê JKLU}KӑWrQ 3*6761JX\ӉQ9 Q+LӃX %Ӝ*,È2'Ө&9¬ Ҥ27Ҥ2 TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH PHIẾU NHẬN XÉT LUẬN VĂN THẠC SỸ (Dành cho JiảnJ viên phản biện) Tên đề tài luận văn thạc sỹ: 3KkQWtFKӭQJ[ӱSKLWX\ӃQGҫPFRPSRVLWHFKӏXWiFGөQJWҧLWUӑQJÿLӅXKzDGL ÿӝQJ Tên tác Jiả: /Ç1*Ӑ&75ѬӠ1* MSHV: 1580822 Ngành: WKXұW[k\GӵQJF{QJWUuQKGkQGөQJYjF{QJQJKLӋS Khóa: 2016 Định hướng: ӬQJGөQJ Họ tên nJười phản biện: 76&KkX uQK7KjQK Cơ quDn cônJ tác: KRD;k\'ӵQJ7UѭӡQJ +6ѭSKҥP WKXұW7S+&0 Điện thoại liên hệ: 0903092979 I Ý KIẾN NHẬN XÉT Về hình thức & kết cấu luận văn: +uQKWKӭFYjNӃWFҩXOXұQY QSKKӧSYӟLÿӅWjLWKHRKѭӟQJQJKLrQFӭXSKѭѫQJSKiSWtQKWRiQJӗPFiF FKѭѫQJWәQJTXDQFѫVӣOêWKX\ӃWYtGөVӕPLQKKӑDYjNӃWOXұQ Về nội dunJ: 2.1 Nhận xét tính khoa học, rõ ràng, mạch lạc, khúc chiết luận văn /XұQY QYLӃWU}UjQJPҥFKOҥFYjNK~FFKLӃW1ӝLGXQJOêWKX\ӃWÿѭӧFWKӇKLӋQEҵQJFiFNӃWTXҧVӕFyVR ViQKYӟLFiFQJKLrQFӭXNKiFKD\ELӋQOXұQNӃWTXҧFyÿѭӧF 2.2 Nhận xét đánh giá việc sử dụng trích dẫn kết NC người khác có qui định hành pháp luật sở hữu trí tuệ /XұQY QVӱGөQJNӃWTXҧKRһFWUtFKGүQQJKLrQFӭXNKiFÿ~QJTXLÿӏQKKLӋQKjQK 2.3 Nhận xét mục tiêu nghiên cứu, phương pháp nghiên cứu sử dụng LVTN 0өFWLrXQJKLrQFӭXSKѭѫQJSKiSQJKLrQFӭXVӱGөQJSKKӧSYӟLWrQÿӅWjLYjQӝLGXQJFҫQÿҥWÿ ѭӧF 2.4 Nhận xét Tổng quan đề tài ĈӅWjLJLҧLTX\ӃWEjLWRiQSKLWX\ӃQKuQKKӑFGҫPFRPSRVLWHQKLӅXOӟSFKӏXWҧLWUӑQJÿLӅXKzDGLÿӝQJG QJ OӡLJLҧL[ҩS[ӍGҥQJÿDWKӭFKRһFOѭӧQJJLiF.ӃWTXҧÿҥWÿѭӧFJӗPSKkQWtFKWҫQVӕYjFKX\ӇQYӏFӫDGҫP &iFWKDPVӕNKiFQKDXFӫDWҧLWUӑQJQKѭWҫQVӕYұQWӕFKѭӟQJVӧLÿһFWUѭQJYұWOLӋXÿmÿѭӧFNKҧRVi WYj ELӋQOXұQ 2.5 Nhận xét đánh giá nội dung & chất lượng LVTN 1ӝLGXQJYjFKҩWOѭӧQJÿӅWjLWӕW 2.6 Nhận xét đánh giá khả ứng dụng, giá trị thực tiễn đề tài ӃWTXҧÿҥWÿѭӧFWURQJÿӅWjLOjWjLOLӋXWKDPNKҧRFKRFiFQJKLrQFӭXNKiF [34] Khdeir AA, Reddy JN Free vibration of cross-ply laminated beams with arbitrary boundary conditions Int J Eng Sci 1994;32(12):1971–80 cited By (since 1996) 47 [35] Thai H-T, Vo TP Bending and free vibration beam theories.Int J Mech Sci 2012:62(1):57-66 [36] Trung-Kien Nguyen, Ngoc-Duong Nguyen, Thuc P Vo, Huu-Tai Thai Trigonometric-series solution for analysis of laminated composite beams Composite Structures 2017; 160: 142–115 [37] Jean-Marie Berthelot Composite Materials: Mechanical Behavior and Structural Analysis, Springer-Verlag Berlin Heidelberg New York in 1999 62 PHỤ LỤC Tính tần số dao động %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear all ; %% Input global P0 omega vp L N Nb hs syms lamda format long L = 20; S = input('nhap ti le L/H: h = L/S ; b = 1; E1 = 241.5 E2 = E1/40 G12 = 0.6*E2 ; G13 = G12; G23 = 0.5*E2 ; v12 = 0.25; v21 =v12*(E2/E1); Ro = 1.0 * tt1 =90*pi/180; tt2 = 0*pi/180; z1 = -h/2; z2 = +h/2; kss=5/6; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% P0=1000000; % don vi N vp=20; % don vi m/s omega=20; N=12; RL1=250; RL2=500; % Tinh toan cac hang so cung Q11 = E1/(1-v12*v21); Q12 = v12*E2/(1-v12*v21); Q13 = 0; Q14 = 0; Q15 = 0; Q16 = 0; Q22 = E2/(1-v12*v21); Q23 = 0; Q24 = 0; Q25 = 0; Q26 = 0; Q33 = 0; Q34 = 0; Q35 = 0; Q36 = 0; Q44 = G23; Q45 = 0; Q46 = 0; 63 Q55 = G13; Q56 = 0; Q66 = G12; %syms m n tt1 E1 v12 v21 E2 G12 G13 G23 Q11 Q12 Q13 Q14 Q15 Q16 Q22 Q23 Q24 Q25 Q26 Q33 Q34 Q35 Q36 Q44 Q45 Q46 Q55 Q56 Q66 Q = [Q11 Q12 Q13 Q14 Q15 Q16; Q12 Q22 Q23 Q24 Q25 Q26; Q13 Q23 Q33 Q34 Q35 Q36; Q14 Q24 Q34 Q44 Q45 Q46; Q15 Q25 Q35 Q45 Q55 Q56; Q16 Q26 Q36 Q46 Q56 Q66]; m = cos(tt1); n = sin(tt1); Qp_tt1(1,1) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; Qp_tt1(1,2) = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); Qp_tt1(2,2) =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; Qp_tt1(1,6) =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; Qp_tt1(2,6) =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; Qp_tt1(6,6) =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); Qp_tt1(4,4) =Q44*m^2+Q55*n^2; Qp_tt1(4,5) =(Q55-Q44)*n*m; Qp_tt1(5,5) =Q55*m^2+Q44*n^2; Qp_tt1; m = cos(tt2); n = sin(tt2); Qp_tt2(1,1) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; Qp_tt2(1,2) = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); Qp_tt2(2,2) =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; Qp_tt2(1,6) =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; Qp_tt2(2,6) =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; Qp_tt2(6,6) =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); Qp_tt2(4,4) =Q44*m^2+Q55*n^2; Qp_tt2(4,5) =(Q55-Q44)*n*m; Qp_tt2(5,5) =Q55*m^2+Q44*n^2; Qp_tt2; Qp = [Qp_tt1(1,1) Qp_tt2(1,1)]; % -Axx=eval(int(Qp_tt1(1,1),z,-h/2,0)+int(Qp_tt2(1,1),z,0,h/2)); Bxx=eval(int(z*Qp_tt1(1,1),z,-h/2,0)+int(z*Qp_tt2(1,1),z,0,h/2)); Dxx=eval(int(z^2*Qp_tt1(1,1),z,- h/2,0)+int(z^2*Qp_tt2(1,1),z,0,h/2)); Axz=eval(int(Qp_tt1(5,5),z,-h/2,0)+int(Qp_tt2(5,5),z,0,h/2)); Ia=eval(int(Ro,z,-h/2,0)+int(Ro,z,0,h/2)); Ib=eval(int(z*Ro,z,-h/2,0)+int(z*Ro,z,0,h/2)); Id=eval(int(z^2*Ro,z,-h/2,0)+int(z^2*Ro,z,0,h/2)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% BC='CC'; switch BC case 'SS' Nb=2; Ks=BCSS(N,L,Nb); case 'HS' Nb=3; Ks=BCHS(N,L,Nb); case 'HH' Nb=4; Ks=BCHH(N,L,Nb); case 'CF' 64 Nb=4; Ks=BCCF(N,L,Nb); case 'CS' Nb=5; Ks=BCCS(N,L,Nb); case 'CH' Nb=6; Ks=BCCH(N,L,Nb); case 'CC' Nb=8; Ks=BCCC(N,L,Nb); end Kl=LinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia, Ib, Id, Nb); M=MatrixM(N,L,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); Kn=NonLinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); J=solve(det((Kl+Ks)-lamda*M)); lamda=sort(double(sqrt(J).*L^2*sqrt(Ro/E2)./h)); lamda=sort(lamda) %Tan so dao dong khong thu nguyen cua dam Tính chuyển vị theo vận tốc khơng thứ nguyên clear all; clc global P0 omega vp L N Nb hs syms lamda N0 z format long h = 1; S = 15 %input('nhap ti le L/H: L = S*h ; b = 0.5; E1 = 241.5 * 10^6 ; E2 = E1/40 ; G12 = 0.6*E2 ; G13 = G12; G23 = 0.5*E2 ; v12 = 0.25; v21 =v12*(E2/E1); Ro = 1389 * 10^-2; % Kg/m3 => KN/m3 tt1 =90*pi/180; tt2 = 0*pi/180; z1 = -h/2; z2 = +h/2; kss=5/6; kz=0; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% P0=1000; % don vi KN omega=0; N=10; RL1=250; %so diem mau tinh toan ( RecordLength )(Linear) RL2=500; % Tinh toan cac hang so cung Q11 = E1/(1-v12*v21); Q12 = v12*E2/(1-v12*v21); Q13 = 0; Q14 = 0; Q15 = 0; Q16 = 0; Q22 = E2/(1-v12*v21); Q23 = 0; 65 Q24 Q25 Q26 Q33 Q34 Q35 Q36 Q44 Q45 Q46 Q55 Q56 Q66 %syms m n tt1 E1 v12 v21 E2 G12 G13 G23 Q11 Q12 Q13 Q14 Q15 Q16 Q22 Q23 Q24 Q25 Q = [Q11 Q12 Q13 Q14 Q15 Q16; Q12 Q13 Q14 Q15 Q16 m = cos(tt1); n = sin(tt1); Qp_tt1(1,1) Qp_tt1(1,2) Qp_tt1(2,2) Qp_tt1(1,6) Qp_tt1(2,6) Qp_tt1(6,6) Qp_tt1(4,4) Qp_tt1(4,5) Qp_tt1(5,5) Qp_tt1; m = cos(tt2); n = sin(tt2); Qp_tt2(1,1) Qp_tt2(1,2) Qp_tt2(2,2) Qp_tt2(1,6) Qp_tt2(2,6) Qp_tt2(6,6) Qp_tt2(4,4) Qp_tt2(4,5) Qp_tt2(5,5) Qp_tt2; Qp = [Qp_tt1(1,1) Qp_tt2(1,1)]; % -% Axx=eval(int(Qp_tt1(1,1),z,-h/2,0)+int(Qp_tt2(1,1),z,0,h/2)); Bxx=eval(int(z*Qp_tt1(1,1),z,-h/2,0)+int(z*Qp_tt2(1,1),z,0,h/2)); Dxx=eval(int(z^2*Qp_tt1(1,1),z,-h/2,0)+int(z^2*Qp_tt2(1,1),z,0,h/2)); Axz=eval(int(Qp_tt1(5,5),z,-h/2,0)+int(Qp_tt2(5,5),z,0,h/2)); Ia=eval(int(Ro,z,-h/2,0)+int(Ro,z,0,h/2)); Ib=eval(int(z*Ro,z,h/2,0)+int(z*Ro,z,0,h/2)); Id=eval(int(z^2*Ro,z,h/2,0)+int(z^2*Ro,z,0,h/2)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% Boundary conditions 66 Q26 Q22 Q23 Q24 Q25 Q26 BC='SS'; switch BC case 'SS' Nb=2; Ks=BCSS(N,L,Nb); case 'HS' Nb=3; Ks=BCHS(N,L,Nb); case 'HH' Nb=4; Ks=BCHH(N,L,Nb); case 'CF' Nb=4; Ks=BCCF(N,L,Nb); case 'CS' Nb=5; Ks=BCCS(N,L,Nb); case 'CH' Nb=6; Ks=BCCH(N,L,Nb); case 'CC' Nb=8; Ks=BCCC(N,L,Nb); end Kl=LinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia, Ib, Id, Nb); M=MatrixM(N,L,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); Kn=NonLinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); J=solve(det((Kl+Ks)-lamda*M)); D=sqrt(J); D=sort(D) for i=1:12 i vp=0.1*i*20.8447*L/2/pi; % beta=0.1*i dt1=L/vp/RL1; dt2=L/vp/RL2; Kl=LinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia, Ib, Id,Nb); M=MatrixM(N,L,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); Kn=NonLinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); [xl,vl,al]=Newmarkmethod_L(RL1,dt1,Kl,M,Ks); [xn,vn,an]=Newmarkmethod_NN(RL2,dt2,Kl,M,Kn,Ks); x1=max(xl(1,:)) %chuyen vi lon nhat tai giua nhip (Linear) x2=max(xn(1,:)) %chuyen vi lon nhat tai giua nhip (NonLinear) end Tính chuyển vị theo tỉ số L/h clear all ; clc %% Input global P0 omega vp L N Nb hs syms lamda N0 z format long h = 0.5; % m S = input('nhap ti le L/H: ') ; L = S*h ; %m b = 0.5; %m0 %L/h 67 E1 = 241.5 * 10^6 ; E2 = E1/40 ; G12 = 0.6*E2 ; G13 = G12; G23 = 0.5*E2 ; v12 = 0.25; v21 =v12*(E2/E1); Ro = 1389 * 10^-2; % Kg/m3 => KN/m3 tt1 =90*pi/180; tt2 = 0*pi/180; z1 = -h/2; z2 = +h/2; kss=5/6; kz=0; % (Pa) He so nen d %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% P0=1000; % don vi KN vp=20; omega=20; N=10; RL1=250; %so diem mau tinh toan ( RecordLength )(Linear) RL2=500; dt1=L/vp/RL1; dt2=L/vp/RL2; t1=0:dt1:L/vp; t2=0:dt2:L/vp; % Tinh toan cac hang so cung Q11 = E1/(1-v12*v21); Q12 = v12*E2/(1-v12*v21); Q13 = 0; Q14 = 0; Q15 = 0; Q16 = 0; Q22 = E2/(1-v12*v21); Q23 = 0; Q24 = 0; Q25 = 0; Q26 = 0; Q33 = 0; Q34 = 0; Q35 = 0; Q36 = 0; Q44 = G23; Q45 = 0; Q46 = 0; Q55 = G13; Q56 = 0; Q66 = G12; %syms m n tt1 E1 v12 v21 E2 G12 G13 G23 Q11 Q12 Q13 Q14 Q15 Q16 Q22 Q23 Q24 Q25 Q26 Q33 Q34 Q35 Q36 Q44 Q45 Q46 Q55 Q56 Q66 Q = [Q11 Q12 Q13 Q12 Q22 Q23 Q13 Q23 Q33 Q14 Q24 Q34 Q15 Q25 Q35 Q16 Q26 Q36 m = cos(tt1); n = sin(tt1); Q14 Q24 Q34 Q44 Q45 Q46 Q15 Q25 Q35 Q45 Q55 Q56 Q16; Q26; Q36; Q46; Q56; Q66]; 68 Qp_tt1(1,1) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; Qp_tt1(1,2) = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); Qp_tt1(2,2) =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; Qp_tt1(1,6) =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; Qp_tt1(2,6) =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; Qp_tt1(6,6) =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); Qp_tt1(4,4) =Q44*m^2+Q55*n^2; Qp_tt1(4,5) =(Q55-Q44)*n*m; Qp_tt1(5,5) =Q55*m^2+Q44*n^2; Qp_tt1; m = cos(tt2); n = sin(tt2); Qp_tt2(1,1) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; Qp_tt2(1,2) = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); Qp_tt2(2,2) =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; Qp_tt2(1,6) =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; Qp_tt2(2,6) =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; Qp_tt2(6,6) =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); Qp_tt2(4,4) =Q44*m^2+Q55*n^2; Qp_tt2(4,5) =(Q55-Q44)*n*m; Qp_tt2(5,5) =Q55*m^2+Q44*n^2; Qp_tt2; Qp = [Qp_tt1(1,1) Qp_tt2(1,1)]; % -Axx=eval(int(Qp_tt1(1,1),z,-h/2,0)+int(Qp_tt2(1,1),z,0,h/2)); Bxx=eval(int(z*Qp_tt1(1,1),z,-h/2,0)+int(z*Qp_tt2(1,1),z,0,h/2)); Dxx=eval(int(z^2*Qp_tt1(1,1),z,- h/2,0)+int(z^2*Qp_tt2(1,1),z,0,h/2)); Axz=eval(int(Qp_tt1(5,5),z,-h/2,0)+int(Qp_tt2(5,5),z,0,h/2)); Ia=eval(int(Ro,z,-h/2,0)+int(Ro,z,0,h/2)); Ib=eval(int(z*Ro,z,-h/2,0)+int(z*Ro,z,0,h/2)); Id=eval(int(z^2*Ro,z,-h/2,0)+int(z^2*Ro,z,0,h/2)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% BC='SS'; switch BC case 'SS' Nb=2; Ks=BCSS(N,L,Nb); case 'HS' Nb=3; Ks=BCHS(N,L,Nb); case 'HH' Nb=4; Ks=BCHH(N,L,Nb); case 'CF' Nb=4; Ks=BCCF(N,L,Nb); case 'CS' Nb=5; Ks=BCCS(N,L,Nb); case 'CH' Nb=6; Ks=BCCH(N,L,Nb); case 'CC' Nb=8; Ks=BCCC(N,L,Nb); end Kl=LinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia, Ib, Id, Nb); 69 M=MatrixM(N,L,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); Kn=NonLinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); [xl,vl,al]=Newmarkmethod_L(RL1,dt1,Kl,M,Ks); [xn,vn,an]=Newmarkmethod_NN(RL2,dt2,Kl,M,Kn,Ks); % Chuyen vi tai giua nhip cua dam theo t L/h thay doi figure(1) hold on grid on plot(t1,xl(1,:),' r','LineWidth',2); plot(t2,xn(1,:),'LineWidth',2); xlabel('Time(s)') ;ylabel('Midspan Displacement (m)'); xl=max(xl(1,:)) % Chuyen vi lon nhat tai giua nhip (linear) xn=max(xn(1,:)) % Chuyen vi lon nhat tai giua nhip (Nonlinear) Tính chuyển vị theo tần số lực kích thích clear all ; clc %% Input global P0 omega vp L N Nb hs syms lamda N0 z format long h = 1; S = 15 %input('nhap ti le L/H: L = S*h ; b = 0.5; E1 = 241.5 * 10^6 ; E2 = E1/40 ; G12 = 0.6*E2 ; G13 = G12; G23 = 0.5*E2 ; v12 = 0.25; v21 =v12*(E2/E1); Ro = 1389 * 10^-2; % Kg/m3 => KN/m3 tt1 =90*pi/180; tt2 = 0*pi/180; z1 = -h/2; z2 = +h/2; kss=5/6; kz=0; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% P0=1000; % don vi KN vp=30; % don vi m/s N=10; RL1=250; %so diem mau tinh toan ( RecordLength )(Linear) RL2=500; dt1=L/vp/RL1; dt2=L/vp/RL2; t1=0:dt1:L/vp; t2=0:dt2:L/vp; % Tinh toan cac hang so cung Q11 = E1/(1-v12*v21); Q12 = v12*E2/(1-v12*v21); Q13 = 0; Q14 = 0; Q15 = 0; Q16 = 0; Q22 = E2/(1-v12*v21); Q23 = 0; 70 Q24 = 0; Q25 = 0; Q26 = 0; Q33 = 0; Q34 = 0; Q35 = 0; Q36 = 0; Q44 = G23; Q45 = 0; Q46 = 0; Q55 = G13; Q56 = 0; Q66 = G12; %syms m n tt1 E1 v12 v21 E2 G12 G13 G23 Q11 Q12 Q13 Q14 Q15 Q16 Q22 Q23 Q24 Q25 Q26 Q33 Q34 Q35 Q36 Q44 Q45 Q46 Q55 Q56 Q66 Q = [Q11 Q12 Q13 Q14 Q15 Q16; Q12 Q22 Q23 Q24 Q25 Q26; Q13 Q23 Q33 Q34 Q35 Q36; Q14 Q24 Q34 Q44 Q45 Q46; Q15 Q25 Q35 Q45 Q55 Q56; Q16 Q26 Q36 Q46 Q56 Q66]; m = cos(tt1); n = sin(tt1); Qp_tt1(1,1) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; Qp_tt1(1,2) = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); Qp_tt1(2,2) =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; Qp_tt1(1,6) =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; Qp_tt1(2,6) =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; Qp_tt1(6,6) =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); Qp_tt1(4,4) =Q44*m^2+Q55*n^2; Qp_tt1(4,5) =(Q55-Q44)*n*m; Qp_tt1(5,5) =Q55*m^2+Q44*n^2; Qp_tt1; m = cos(tt2); n = sin(tt2); Qp_tt2(1,1) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; Qp_tt2(1,2) = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); Qp_tt2(2,2) =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; Qp_tt2(1,6) =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; Qp_tt2(2,6) =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; Qp_tt2(6,6) =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); Qp_tt2(4,4) =Q44*m^2+Q55*n^2; Qp_tt2(4,5) =(Q55-Q44)*n*m; Qp_tt2(5,5) =Q55*m^2+Q44*n^2; Qp_tt2; Qp = [Qp_tt1(1,1) Qp_tt2(1,1)]; % -Axx=eval(int(Qp_tt1(1,1),z,-h/2,0)+int(Qp_tt2(1,1),z,0,h/2)); Bxx=eval(int(z*Qp_tt1(1,1),z,-h/2,0)+int(z*Qp_tt2(1,1),z,0,h/2)); Dxx=eval(int(z^2*Qp_tt1(1,1),z,-h/2,0)+int(z^2*Qp_tt2(1,1),z,0,h/2)); Axz=eval(int(Qp_tt1(5,5),z,-h/2,0)+int(Qp_tt2(5,5),z,0,h/2)); Ia=eval(int(Ro,z,-h/2,0)+int(Ro,z,0,h/2)); Ib=eval(int(z*Ro,z,h/2,0)+int(z*Ro,z,0,h/2)); Id=eval(int(z^2*Ro,z,h/2,0)+int(z^2*Ro,z,0,h/2)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% BC='SS'; switch BC case 'SS' 71 end Nb=2; Ks=BCSS(N,L,Nb); case 'HS' Nb=3; Ks=BCHS(N,L,Nb); case 'HH' Nb=4; Ks=BCHH(N,L,Nb); case 'CF' Nb=4; Ks=BCCF(N,L,Nb); case 'CS' Nb=5; Ks=BCCS(N,L,Nb); case 'CH' Nb=6; Ks=BCCH(N,L,Nb); case 'CC' Nb=8; Ks=BCCC(N,L,Nb); Kl=LinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia, Ib, Id, Nb); M=MatrixM(N,L,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); Kn=NonLinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); J=solve(det((Kl+Ks)-lamda*M)); D=sqrt(J); D=sort(D) for i=1:6 i omega=10*i; Kl=LinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia, Ib, Id,Nb); M=MatrixM(N,L,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); Kn=NonLinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); [xl,vl,al]=Newmarkmethod_L(RL1,dt1,Kl,M,Ks); [xn,vn,an]=Newmarkmethod_NN(RL2,dt2,Kl,M,Kn,Ks); Mn=Bxx*(xn(N+2,:))+Dxx*(xn(2*N+2,:)); % Moment tai giua nhip xl=max(xl(1,:)) xn=max(xn(1,:)) Tính chuyển vị P đặt ¼ nhịp nhịp clear all; clc global P0 omega vp L N Nb hs syms lamda N0 z x hss1 hss2 hss3 format long h = 1; S = 20 %input('nhap ti le L/H: L = S*h ; b = 0.5; E1 = 241.5 * 10^6 ; E2 = E1/40 ; G12 = 0.6*E2 ; G13 = G12; G23 = 0.5*E2 ; v12 = 0.25; v21 =v12*(E2/E1); Ro = 1389 * 10^-2; % Kg/m3 => KN/m3 tt1 =90*pi/180; 72 tt2 = 0*pi/180; z1 = -h/2; z2 = +h/2; kss=5/6; % He so hieu chinh ung suat cat kz=0; % (Pa) He so nen dan hoi Winkler %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% P0=1500; % don vi KN vp=80; % don vi m/s omega=20; N=10; RL1=250; %so diem mau tinh toan ( RecordLength )(Linear) RL2=500; %Cho bai toan chuyen vi theo van toc,Omega,L/h dt1=L/vp/RL1; dt2=L/vp/RL2; % Tinh toan cac hang so cung Q11 = E1/(1-v12*v21); Q12 = v12*E2/(1-v12*v21); Q13 = 0; Q14 = 0; Q15 = 0; Q16 = 0; Q22 = E2/(1-v12*v21); Q23 = 0; Q24 = 0; Q25 = 0; Q26 = 0; Q33 = 0; Q34 = 0; Q35 = 0; Q36 = 0; Q44 = G23; Q45 = 0; Q46 = 0; Q55 = G13; Q56 = 0; Q66 = G12; %syms m n tt1 E1 v12 v21 E2 G12 G13 G23 Q11 Q12 Q13 Q14 Q15 Q16 Q22 Q23 Q24 Q25 Q26 Q33 Q34 Q35 Q36 Q44 Q45 Q46 Q55 Q56 Q66 Q = [Q11 Q12 Q13 Q14 Q15 Q12 Q22 Q23 Q24 Q25 Q13 Q23 Q33 Q34 Q35 Q14 Q24 Q34 Q44 Q45 Q15 Q25 Q35 Q45 Q55 Q16; Q26; Q36; Q46; Q56; Q16 Q26 Q36 Q46 Q56 Q66]; m = cos(tt1); n = sin(tt1); Qp_tt1(1,1) Qp_tt1(1,2) Qp_tt1(2,2) Qp_tt1(1,6) Qp_tt1(2,6) Qp_tt1(6,6) Qp_tt1(4,4) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); =Q44*m^2+Q55*n^2; Qp_tt1(4,5) =(Q55-Q44)*n*m; Qp_tt1(5,5) =Q55*m^2+Q44*n^2; Qp_tt1; m = cos(tt2); n = sin(tt2); 73 Qp_tt2(1,1) = Q11*m^4 + 2*(Q12+2*Q66)*n^2*m^2+Q22*n^4; Qp_tt2(1,2) = (Q11+Q22-4*Q66)*m^2*n^2+Q12*(m^4+n^4); Qp_tt2(2,2) =Q11*n^4+2*(Q12+2*Q66)*m^2*n^2+Q22*m^4; Qp_tt2(1,6) =(Q11-Q12-2*Q66)*n*m^3+(Q12-Q22+2*Q66)*n^3*m; Qp_tt2(2,6) =(Q11-Q12-2*Q66)*n^3*m+(Q12-Q22+2*Q66)*n*m^3; Qp_tt2(6,6) =(Q11+Q22-2*Q12-2*Q66)*n^2*m^2+Q66*(n^4+m^4); Qp_tt2(4,4) =Q44*m^2+Q55*n^2; Qp_tt2(4,5) =(Q55-Q44)*n*m; Qp_tt2(5,5) =Q55*m^2+Q44*n^2; Qp_tt2; Qp = [Qp_tt1(1,1) Qp_tt2(1,1)]; Axx=eval(int(Qp_tt1(1,1),z,-h/2,0)+int(Qp_tt2(1,1),z,0,h/2)); Bxx=eval(int(z*Qp_tt1(1,1),z,-h/2,0)+int(z*Qp_tt2(1,1),z,0,h/2)); Dxx=eval(int(z^2*Qp_tt1(1,1),z,h/2,0)+int(z^2*Qp_tt2(1,1),z,0,h/2)); Axz=eval(int(Qp_tt1(5,5),z,-h/2,0)+int(Qp_tt2(5,5),z,0,h/2)); Ia=eval(int(Ro,z,-h/2,0)+int(Ro,z,0,h/2)); Ib=eval(int(z*Ro,z,-h/2,0)+int(z*Ro,z,0,h/2)); Id=eval(int(z^2*Ro,z,-h/2,0)+int(z^2*Ro,z,0,h/2)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% BC='SS'; switch BC case 'SS' Nb=2; Ks=BCSS(N,L,Nb); case 'HS' Nb=3; Ks=BCHS(N,L,Nb); case 'HH' Nb=4; Ks=BCHH(N,L,Nb); case 'CF' Nb=4; Ks=BCCF(N,L,Nb); case 'CS' Nb=5; Ks=BCCS(N,L,Nb); case 'CH' Nb=6; Ks=BCCH(N,L,Nb); case 'CC' Nb=8; Ks=BCCC(N,L,Nb); end Kl=LinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia, Ib, Id, Nb); M=MatrixM(N,L,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); Kn=NonLinearMatrixK(N,L,kss,kz,Axx,Bxx,Dxx,Axz,Ia,Ib,Id,Nb); [xl,vl,al]=Newmarkmethod_L(RL1,dt1,Kl,M,Ks); [xn,vn,an]=Newmarkmethod_NN(RL2,dt2,Kl,M,Kn,Ks); % TINH CHUYEN VI CUA DAM KHI P O GIUA DAM a=[]; for i=N:-1:1 a=[a sym(strcat('a',num2str(i)))]; end hss1=a(end:1:1); w = sym(0); for i = a(:).' w = w*x + i; 74 end w = expand(w); xl1=subs(w,hss1,double(xl(1:N,125))); %chuyen vi P giua dam (Linear) xl2=subs(w,hss1,double(xl(1:N,62))); %chuyen vi P o 1/4 dam (Linear) xn1=subs(w,hss1,double(xn(1:N,125))); %chuyen vi P giua dam (NonLinear) xn2=subs(w,hss1,double(xn(1:N,62))); %chuyen vi P 1/4 dam (NonLinear) % HINH VE CHUYEN VI CUA DAM x=-10:2:10; xl1=subs(xl1) xn1=subs(xn1) xl2=subs(xl2) xn2=subs(xn2) %HINH VE CHUYEN VI CUA DAM KHI LUC TAI VI TRI 1/4 DAM figure(3) hold on grid on plot(x,xl2,'LineWidth',2); plot(x,xn2,'-r','LineWidth',2); xlabel('x(m)');ylabel('Displacement(m)'); %HINH VE CHUYEN VI CUA DAM KHI LUC TAI VI TRI GIUA DAM figure(6) hold on grid on plot(x,xl1,'LineWidth',2); plot(x,xn1,'-r','LineWidth',2); xlabel('x(m)');ylabel('Displacement(m)'); 75 ... Đề tài phân tích ứng xử phi tuyến dầm composite nhịp chịu tải trọng điều hịa di động Mơ hình vật lý tốn cho sau: Hình 1.8: Dầm composite chịu tác dụng tải trọng điều hòa di động P(t) 1.3 Tổng... luận văn này, ứng xử phi tuyến dầm Composite chịu tải trọng điều hòa di động phân tích dựa lý thuyết dầm Timosenko quan hệ biến dạng chuyển vị phi tuyến Von-Karman Phương trình động lực học dầm. .. dạng dầm khác Tuy vậy, ứng xử động phi tuyến dầm composite chịu tải trọng di động chưa nghiên cứu nhiều Tiếp nối với nghiên cứu ứng xử kết cấu dạng dầm làm vật liệu composite Đề tài phân tích ứng