Phân tích ứng xử phi tuyến dầm BTCT chịu tải trọng va chạm

Đ A ÏI HỌC Q UỐ C GIA TP HỒCH Í MINH TR ƯỜ NG Đ ẠI HỌC BA Ù CH K HOA K HOA K ỸTH UẬ T XÂ Y D ỰNG - - NGUYỄ N THANH TÂ N PHÂ N TÍCH Ứ NG XỬPHI TUYẾ N DẦ M BTCT CHỊU T TRỌNG VA CHẠM CHUYÊ N NGÀ NH : XÂ Y DỰNG DÂ N DỤNG & CÔ NG NGHIỆ P MÃSỐNGÀ NH : 60 58 20 LUẬ N VĂ N THẠC SĨ TP HCM, 07/2010 CÔ NG TRÌNH ĐƯC HOÀ N THÀ NH TẠI TRƯỜ NG ĐẠI HỌC BÁ CH KHOA ĐẠI HỌC QUỐ C GIA TP HỒCHÍ MINH Cán bộhướng dẫn khoa học: TS HỒHỮ U CHỈ NH Cán bộchấm nhận xét 1: Cán bộchấm nhận xét 2: Luận văn thạc sónày bảo vệtại HỘ I ĐỒ NG CHẤ M BẢ O VỆLUẬ N VĂ N THẠC SĨ TRƯỜ NG ĐẠI HỌC BÁ CH KHOA TP HỒCHÍ MINH, ngày tháng năm 2010 TRƯỜ NG ĐẠI HỌC BÁ CH KHOA PHÒ NG ĐÀ O TẠO SĐH CỘ NG HÒ A XÃHỘ I CHỦNGHĨA VIỆ T NAM Độc lập – Tựdo – Hạnh phú c Tp.HCM, ngày .tháng .năm 2010 NHIỆ M VỤLUẬ N VĂ N THẠC SĨ Họtên học viên: NGUYỄ N THANH TÂ N Phái : Nam Ngày, tháng, năm sinh: 09/02/1981 Nơi sinh : Bình Trị Thiên Chuyên ngành: Xây dựng dân dụng vàcông nghiệp MSHV : 02108494 Khóa (năm trúng tuyển): 2008 I- TÊ N ĐỀTÀ I: PHÂ N TÍCH Ứ NG XỬPHI TUYẾ N DẦ M BTCT CHỊU TẢ I TRỌNG VA CHẠM II- NHIỆ M VỤVÀNỘ I DUNG: - Nghiên cứu phân tích ứng xửphi tuyến vật liệu BTCT theo Modified Compression Field Theory (MCFT) - Xây dựng phần tửph tứgiác, tam giác đểmôphỏng ứng xửphi tuyến kết cấu BTCT chịu tải trọng dạng xung (va chạm) - Lập trình tính toán với hỗtrợ phần mềm MATLAB, khảo sát sốví dụ tính toán cóđối chứng với kết quảphân tích module SAP2000 phi tuyến (Version 14) III- NGÀ Y GIAO NHIỆ M VỤ: 22/02/2010 IV- NGÀ Y HOÀ N THÀ NH NHIỆ M VỤ: 02/07/2010 V- CÁ N BỘHƯỚ NG DẪ N: TS HỒHỮ U CHỈ NH Nội dung vàđềcương Luận văn thạc só đãđược Hội Đồng Chuyên Ngành thông qua CÁ N BỘHƯỚ NG DẪ N TS HỒHỮ U CHỈ NH TRƯỞ NG BAN QUẢ N LÝCHUYÊ N NGÀ NH L L ên, tơi xin bày t ịng bi ịng Tôi xin chân thành c TP.H ã truy h nghiên c u t b ý th Trong trình th th chun mơn kh ,h ài li ên c , TS H chân thành c h ành ã sinh t ãg ,c nh bi s chu ên khích l r l ày, c chuyên sâu v th th TS H truy cho ki Trong th th , tài li s nh theo h ên c ình b ki xin , c ãg nhi , chia s ki ên quy hồn thành khóa h h ã may m , xin c t b Nh ình h ã giúp tơi ngh nghi ên môn c M õi quan tr bày t òng bi sâu s PGS.TS Chu Qu Th t ình, nghiên c ên mơn v àng, giúp th ni r chân thành c Tp.HCM Nguy i ã M L …………………………………………………………………………… i Danh m Danh m u………………………………………………………… …ii ình Ký hi ……………………………………………………………….…iii …………………………………………………………………… vi ………….………………………………………… .1 1.1 Gi i thi 1.2 T ình hình nghiên c 1.2.1 1.2.2 1.3 M nghiên c : MƠ HÌNH PHI TUY ên th ………3 …… ………… 11 ……………………….13 2.1 Các gi 2.2 Thi 2.3 Quan h 2.3.1 Quan h 2.3.2 Quan h ………………………………………………………… 13 ình cân b – bi – bi – bi ình bêtơng……………………17 T ………………………………………… … 20 ình t 3.2 Ma tr 3.2.1 Ma tr 3.2.2 Ma tr 3.3 Ma tr 3.4 Ma tr 3.4.1 Ma tr 3.4.2 Ma tr 3.4.3 Ma tr 3.5 Vector t 3.6 Kh 3.6.2 Tích phân t quát………………………………………………………… 20 thu g …………….25 ……………………………………31 …………………35 … ………………………………………………………… 39 … 39 3.7 Thu tốn phân tích…………………………………………… …… ……44 ………………………………………….45 4.1 Khái qt………………………………………………………………………….46 4.2 Ví d 01……………………….………………………………………………….48 4.3 Ví d …….………………………………………………….64 4.4 Ví d ……………………….……………….73 4.5 K ……………………………………………… ……………… 79 À KI ………………………………………80 5.1 K 5.2 Ki h PH …………………………………………………………………… …….82 TÀI LI ………………………………………………………….…………… 80 …………………………………………………………120 DANH M B B B B B B B B So sánh mơ hình phi So sánh gi So sánh k hân tích……………………………………………… 52 K h 1……… …………….….54 Gi qu th K ………………… 69 K ii ÌNH DANH M CH Hình 1.1 Hi …3 Hình 1.2 V àc …… ……5 Hình 1.3 Seaman, 1978)………………………………………………………………….… Hình 1.4 Xun th bêtơng, Agradh Laine (1999)…………… …6 Hình 1.5 Xuyên th sàn bêtơng, Teng c …… … Hình 1.6 Mơ hình DEM, Sawamoto c …… ….8 Hình 1.7 S ên d BTCT, Magnier Donzé (1998)……………………… …8 Hình 1.8 M .9 Hình 1.9 ành h àc 10 Hình 1.10 Mơ hình cho s ……………… 11 Hính 1.11 Quan h – bi ình t …… 12 CH Hình 2.1 Hình 2.2 Hình 2.3 Hình 2.4 Hình 2.5 Hình 2.6 Hình 2.7 Hình 2.8 Mơ hình v Ph Quan h Quan h Quan h …………… … 13 àng……………………………………………… … 14 nd …… …………14 …… …… 15 ình bêtông sau n …….….……16 – bi ……… … 17 – bi êtông b …… 18 – bi êtơng…………………………… … 19 CH Hình 3.1 S … …20 Hình 3.2 Thu g hình ch ……………………….21 Hình 3.3 Thu g hình tam giác……………………………… 22 Hình 3.4 Bi ên c … 24 Hình 3.5 pra , 2001)……………………………………………….25 Hình 3.6 Ph CT b …25 Hình 3.7 Vịng trịn Morh bi ình………………………… …… ….26 Hình 3.8 Các liên h module cát …….27 Hình 3.9 Vịng trịn Morh ình……………………… ……… … 28 Hình 3.10 Các h ………………………………….…29 Hình 3.11 Ph ình ch iii Hình 3.12 Hình 3.13 Hình 3.14 Hình 3.15 Ph ình tam giác…………………………………………….………35 Gi odule cát Các giá tr ti ……………….40 Gi g… ….44 CH Hình 4.1 Chi ti Hình 4.2 HH d Hình 4.3 g ph Hình 4.4 Quan h – chuy 50 Hình 4.5 Mơ hình SAP2000 d Hình 4.6 Quan h – ích theo SAP2000 d Hình 4.7 So sánh k – SAP2000 d ……………………………52 Hình 4.8 Mơ ph 53 Hình 4.9 d Hình 4.10 , phân tích phi Hình 4.11 tích Hình 4.12 …………… .58 Hình 4.13 S …………… …………………….…………… .59 Hình 4.14 Hình 4.15 Hình 4.16 Hình 4.17 Hình 4.18 Hình 4.19 Hình 4.20 Hình 4.21 Hình 4.22 Hình 4.23 Hình 4.24 Hình 4.25 Hình 4.26 Hình 4.27 Hình 4.28 Hình 4.29 Hình 4.30 Hình 4.31 t nk Chi ti PTHH d ong ph Quan h – chuy Các hình thái phá ho Mơ ph va ch Mơ hình SAP2000 d mút d ………………………………… 66 67 0, phân tích phi 0, phân tích Chi ti Mơ ph Mơ hình SAP2000 d ………….….75 ………….….76 iv Hình 4.32 Hình 4.33 hân tích v title('DO THI CHUYEN VI - THOI GIAN'); axis([0 -12.5 12.5]); grid on; hold on; % -END -Các function % MA TRAN DO CUNG PHAN TU function [K]=Estiff(D,a,b,d) t0 = d/(12*a*b); k11=4*b^2*t0*D(1,1)+6*a*b*t0*D(1,3)+4*a^2*t0*D(3,3); k12=3*a*b*t0*D(1,2)+4*b^2*t0*D(1,3)+4*a^2*t0*D(2,3)+3*a*b*t0*D(3,3); k13=-4*b^2*t0*D(1,1)+2*a^2*t0*D(3,3); k14=3*a*b*t0*D(1,2)-4*b^2*t0*D(1,3)+2*a^2*t0*D(2,3)-3*a*b*t0*D(3,3); k15=-2*b^2*t0*D(1,1)-6*a*b*t0*D(1,3)-2*a^2*t0*D(3,3); k16=-3*a*b*t0*D(1,2)-2*b^2*t0*D(1,3)-2*a^2*t0*D(2,3)-3*a*b*t0*D(3,3); k17=2*b^2*t0*D(1,1)-4*a^2*t0*D(3,3); k18=-3*a*b*t0*D(1,2)+2*b^2*t0*D(1,3)-4*a^2*t0*D(2,3)+3*a*b*t0*D(3,3); k22=4*a^2*t0*D(2,2)+6*a*b*t0*D(2,3)+4*b^2*t0*D(3,3); k23=-3*a*b*t0*D(1,2)-4*b^2*t0*D(1,3)+2*a^2*t0*D(2,3)+3*a*b*t0*D(3,3); k24=2*a^2*t0*D(2,2)-4*b^2*t0*D(3,3); k25=k16; k26=-2*a^2*t0*D(2,2)-6*a*b*t0*D(2,3)-2*b^2*t0*D(3,3); k27=3*a*b*t0*D(1,2)+2*b^2*t0*D(1,3)-4*a^2*t0*D(2,3)-3*a*b*t0*D(3,3); k28=-4*a^2*t0*D(2,2)+2*b^2*t0*D(3,3); k33=4*b^2*t0*D(1,1)-6*a*b*t0*D(1,3)+4*a^2*t0*D(3,3); k34=-3*a*b*t0*D(1,2)+4*b^2*t0*D(1,3)+4*a^2*t0*D(2,3)-3*a*b*t0*D(3,3); k35=k17; k36=k27; k37=-2*b^2*t0*D(1,1)+6*a*b*t0*D(1,3)-2*a^2*t0*D(3,3); k38=3*a*b*t0*D(1,2)-2*b^2*t0*D(1,3)-2*a^2*t0*D(2,3)+3*a*b*t0*D(3,3); k44=4*a^2*t0*D(2,2)-6*a*b*t0*D(2,3)+4*b^2*t0*D(3,3); k45=k18; k46=k28; k47=k38; k48=-2*a^2*t0*D(2,2)+6*a*b*t0*D(2,3)-2*b^2*t0*D(3,3); k55=k11; k56=k12; k57=k13; k58=k14; k66=k22; k67=k23; k68=k24; k77=k33; k78=k34; k88=k44; K= [k11 k12 k13 k14 k15 k16 k17 k18 k12 k22 k23 k24 k25 k26 k27 k28 k13 k23 k33 k34 k35 k36 k37 k38 k14 k24 k34 k44 k45 k46 k47 k48 k15 k25 k35 k45 k55 k56 k57 k58 k16 k26 k36 k46 k56 k66 k67 k68 k17 k27 k37 k47 k57 k67 k77 k78 k18 k28 k38 k48 k58 k68 k78 k88]; % -End of Function -% HAM TAO MA TRAN CHI SO PHAN TU function [ix]=indexos(endoe,nonpe,nodofpn) k=0; for i=1:nonpe %bien chay nut s=(endoe(i)-1)*nodofpn; for j=1:nodofpn %bien chay chi so moi nut k=k+1; ix(k)=s+j; end end % -End of Function 106 % MA TRAN CUNG TONG THE function [KOS]=Tostiff(KOS,K,ix) for i=1:length(ix) %bien chay tren ma tran is=ix(i); for j=1:length(ix) %bien chay tren ma tran js=ix(j); KOS(is,js)=KOS(is,js)+K(i,j); %bien chay tren ma tran he thong end end % -End of Function -% HAM TINH MA TRAN KHOI LUONG PHAN TU function [km] = mamass(den,a,b,d) m = den*a*b*d; h = 8; km = zeros(h,h); for i =1:h for j = 1:h if i==j km(i,j) = m/4; else km(i,j) = 0; end end end % -End of Function %GHEP NOI MA TRAN KHOI LUONG TONG THE function [MOS] = ToMaMass(MOS,km,ix) for i = 1: length(ix) is = ix(i); for j = 1: length (ix) js = ix(j); % bien chay tren ma tran %bien chay tren ma tran MOS(is,js) = MOS(is,js)+km(i,j); %bien chay tren ma tran he thong end end % -End of Function -% HAM XAC XAC DINH MA TRAN CAN THEO MA TRAN KHOI LUONG VA MA TRAN DO CUNG function [C] = matrixC(MOS,KOS) % ti so can ket cau shi = 0.03; shj = 0.05; % tan so vong dao dong omega % tan so co ban omegai = 6.2; omegaj = 117.5; % tan so thoa 90% khoi luong % tinh cac he so ti le a1 = (2*shj*omegaj-2*shi*omegai)/(omegaj^2 - omegai^2); a0 = 2*shi*omegai-a1*omegai^2; % MA TRAN CAN KET CAU THEO RAYLEIGH [C] = a0*MOS + a1*KOS; % -End of Function -% HAM XAC XAC DINH MA TRAN CAN THEO MA TRAN KHOI LUONG VA MA TRAN DO CUNG function [C] = matrixC(MOS,KOS) % ti so can ket cau shi = 0.01; shj = 0.05; 107 % tan so vong dao dong omega omegai = 6.2; % tan so co ban omegaj = 117.5; % tan so thoa 90% khoi luong % tinh cac he so ti le a1 = (2*shj*omegaj-2*shi*omegai)/(omegaj^2 - omegai^2); a0 = 2*shi*omegai-a1*omegai^2; % MA TRAN CAN KET CAU THEO RAYLEIGH [C] = a0*MOS + a1*KOS; % -End of Function -% XU LY RANG BUOC MA TRAN CUNG TINH function [KOS,f]=straintKOS(KOS,f,ixres,vodof) m=length(ixres); n=length(KOS); for i=1:m c=ixres(i); for j=1:n KOS(c,j)=0; KOS(j,c)=0; end %gan lai cac so hang tren duong cheo chinh =1 KOS(c,c)=1; f(c)=vodof(i); %xoa luc tap trung ung voi chi so bi rang buoc end % End of Function - % XU LY RANG BUOC MA TRAN CUNG DONG HOC function [KD,f]= straintKD(KD,f,ixres,vodof) m=length(ixres); n=length(KD); for i=1:m c=ixres(i); for j=1:n KD(c,j)=0; KD(j,c)=0; end KD(c,c)=1; %gan lai cac so hang tren duong cheo chinh =1 %xoa luc tap trung ung voi chi so bi rang buoc f(c)=vodof(i); end % End of Function % XU LY RANG BUOC MA TRAN KHOI LUONG function [MOS,f]=straintMOS(MOS,f,ixres,vodof) m=length(ixres); n=length(MOS); for i=1:m c=ixres(i); for j=1:n MOS(c,j)=0; MOS(j,c)=0; end MOS(c,c)=1; %gan lai cac so hang tren duong cheo chinh =1 %xoa luc tap trung ung voi chi so bi rang buoc f(c)=vodof(i); end % -End of Function 108 Bài s % -%CHUONG TRINH PHAN TICH PHI TUYEN DAM BTCT CHIU TAI TRONG DONG %Hoc vien NGUYEN THANH TAN - MSHV: 02108494 - K2008 - DHBK TP.HCM %Bai so 03 - Bai toan phan tu tam giac % -clc clear format compact % -%Dac trung vat lieu %MPa hay N/mm2 fc=25; eo=-0.002; v=0.25; %He so poisson %MPa hay N/mm2 Es=195000; px1=0.02; px2=0.00; py=0.0019; fyx=400; %MPa hay N/mm2 fyy=400; %MPa hay N/mm2 %MPa hay N/mm2 Ec=2*fc/abs(eo); fcr=0.33*sqrt(fc); %MPa hay N/mm2 ecr=fcr/Ec; % -%Dac trung hinh hoc d=200; %be day ket cau (mm) anpha1=0; %goc huong cot thep doc %(rad)goc huong cot thep dung anpha2=pi/2; %goc huong dam so voi toa tong the beta=0; % -%so nut non=63; noe=91; %so phan tu %so nut moi phan tu nonpe=3; nodofpn=2; %so bac tu cua moi nut nodofos=non*nodofpn; %so bac tu cua he thong % -%Toa nut coord(1,1)=1100; coord(1,2)=0; coord(2,1)=1250; coord(2,2)=0; coord(3,1)=1400; coord(3,2)=0; coord(4,1)=1550; coord(4,2)=0; coord(5,1)=1700; coord(5,2)=0; coord(6,1)=1850; coord(6,2)=0; coord(7,1)=2000; coord(7,2)=0; coord(8,1)=0; coord(9,1)=150; coord(10,1)=300; coord(11,1)=450; coord(12,1)=600; coord(13,1)=750; coord(14,1)=900; coord(15,1)=1100; coord(16,1)=1250; coord(17,1)=1400; coord(18,1)=1550; coord(19,1)=1700; coord(20,1)=1850; coord(21,1)=2000; coord(8,2)=100; coord(9,2)=100; coord(10,2)=100; coord(11,2)=100; coord(12,2)=100; coord(13,2)=100; coord(14,2)=100; coord(15,2)=100; coord(16,2)=100; coord(17,2)=100; coord(18,2)=100; coord(19,2)=100; coord(20,2)=100; coord(21,2)=100; coord(22,1)=0; coord(23,1)=150; coord(24,1)=300; coord(25,1)=450; coord(26,1)=600; coord(27,1)=750; coord(22,2)=200; coord(23,2)=200; coord(24,2)=200; coord(25,2)=200; coord(26,2)=200; coord(27,2)=200; 109 coord(28,1)=900; coord(29,1)=1100; coord(30,1)=1250; coord(31,1)=1400; coord(32,1)=1550; coord(33,1)=1700; coord(34,1)=1850; coord(35,1)=2000; coord(28,2)=200; coord(29,2)=200; coord(30,2)=200; coord(31,2)=200; coord(32,2)=200; coord(33,2)=200; coord(34,2)=200; coord(35,2)=200; coord(36,1)=0; coord(37,1)=150; coord(38,1)=300; coord(39,1)=450; coord(40,1)=600; coord(41,1)=750; coord(42,1)=900; coord(43,1)=1100; coord(44,1)=1250; coord(45,1)=1400; coord(46,1)=1550; coord(47,1)=1700; coord(48,1)=1850; coord(49,1)=2000; coord(36,2)=400; coord(37,2)=400; coord(38,2)=400; coord(39,2)=400; coord(40,2)=400; coord(41,2)=400; coord(42,2)=400; coord(43,2)=400; coord(44,2)=400; coord(45,2)=400; coord(46,2)=400; coord(47,2)=400; coord(48,2)=400; coord(49,2)=400; coord(50,1)=0; coord(51,1)=150; coord(52,1)=300; coord(53,1)=450; coord(54,1)=600; coord(55,1)=750; coord(56,1)=900; coord(57,1)=1100; coord(58,1)=1250; coord(59,1)=1400; coord(60,1)=1550; coord(61,1)=1700; coord(62,1)=1850; coord(63,1)=2000; coord(50,2)=500; coord(51,2)=500; coord(52,2)=500; coord(53,2)=500; coord(54,2)=500; coord(55,2)=500; coord(56,2)=500; coord(57,2)=500; coord(58,2)=500; coord(59,2)=500; coord(60,2)=500; coord(61,2)=500; coord(62,2)=500; coord(63,2)=500; % -%Khai bao phan tu elem(1,1)=1; elem(2,1)=1; elem(3,1)=2; elem(4,1)=2; elem(5,1)=3; elem(6,1)=3; elem(7,1)=4; elem(8,1)=4; elem(9,1)=5; elem(10,1)=5; elem(11,1)=6; elem(12,1)=6; elem(13,1)=7; elem(14,1)=8; elem(15,1)=9; elem(16,1)=9; elem(17,1)=10; elem(18,1)=10; elem(19,1)=11; elem(20,1)=11; elem(21,1)=12; elem(22,1)=12; elem(23,1)=13; elem(24,1)=13; elem(1,2)=15; elem(2,2)=2; elem(3,2)=16; elem(4,2)=3; elem(5,2)=17; elem(6,2)=4; elem(7,2)=18; elem(8,2)=5; elem(9,2)=19; elem(10,2)=6; elem(11,2)=20; elem(12,2)=7; elem(13,2)=21; elem(14,2)=9; elem(15,2)=23; elem(16,2)=10; elem(17,2)=24; elem(18,2)=11; elem(19,2)=25; elem(20,2)=12; elem(21,2)=26; elem(22,2)=13; elem(23,2)=27; elem(24,2)=14; elem(1,3)=14; elem(2,3)=15; elem(3,3)=15; elem(4,3)=16; elem(5,3)=16; elem(6,3)=17; elem(7,3)=17; elem(8,3)=18; elem(9,3)=18; elem(10,3)=19; elem(11,3)=19; elem(12,3)=20; elem(13,3)=20; elem(14,3)=22; elem(15,3)=22; elem(16,3)=23; elem(17,3)=23; elem(18,3)=24; elem(19,3)=24; elem(20,3)=25; elem(21,3)=25; elem(22,3)=26; elem(23,3)=26; elem(24,3)=27; 110 elem(25,1)=14; elem(26,1)=14; elem(27,1)=15; elem(28,1)=15; elem(29,1)=16; elem(30,1)=16; elem(31,1)=17; elem(32,1)=17; elem(33,1)=18; elem(34,1)=18; elem(35,1)=19; elem(36,1)=19; elem(37,1)=20; elem(38,1)=20; elem(39,1)=21; elem(25,2)=28; elem(26,2)=15; elem(27,2)=29; elem(28,2)=16; elem(29,2)=30; elem(30,2)=17; elem(31,2)=31; elem(32,2)=18; elem(33,2)=32; elem(34,2)=19; elem(35,2)=33; elem(36,2)=20; elem(37,2)=34; elem(38,2)=21; elem(39,2)=35; elem(25,3)=27; elem(26,3)=28; elem(27,3)=28; elem(28,3)=29; elem(29,3)=29; elem(30,3)=30; elem(31,3)=30; elem(32,3)=31; elem(33,3)=31; elem(34,3)=32; elem(35,3)=32; elem(36,3)=33; elem(37,3)=33; elem(38,3)=34; elem(39,3)=34; elem(40,1)=22; elem(41,1)=23; elem(42,1)=23; elem(43,1)=24; elem(44,1)=24; elem(45,1)=25; elem(46,1)=25; elem(47,1)=26; elem(48,1)=26; elem(49,1)=27; elem(50,1)=27; elem(51,1)=28; elem(52,1)=28; elem(53,1)=29; elem(54,1)=29; elem(55,1)=30; elem(56,1)=30; elem(57,1)=31; elem(58,1)=31; elem(59,1)=32; elem(60,1)=32; elem(61,1)=33; elem(62,1)=33; elem(63,1)=34; elem(64,1)=34; elem(65,1)=35; elem(40,2)=23; elem(41,2)=37; elem(42,2)=24; elem(43,2)=38; elem(44,2)=25; elem(45,2)=39; elem(46,2)=26; elem(47,2)=40; elem(48,2)=27; elem(49,2)=41; elem(50,2)=28; elem(51,2)=42; elem(52,2)=29; elem(53,2)=43; elem(54,2)=30; elem(55,2)=44; elem(56,2)=31; elem(57,2)=45; elem(58,2)=32; elem(59,2)=46; elem(60,2)=33; elem(61,2)=47; elem(62,2)=34; elem(63,2)=48; elem(64,2)=35; elem(65,2)=49; elem(40,3)=36; elem(41,3)=36; elem(42,3)=37; elem(43,3)=37; elem(44,3)=38; elem(45,3)=38; elem(46,3)=39; elem(47,3)=39; elem(48,3)=40; elem(49,3)=40; elem(50,3)=41; elem(51,3)=41; elem(52,3)=42; elem(53,3)=42; elem(54,3)=43; elem(55,3)=43; elem(56,3)=44; elem(57,3)=44; elem(58,3)=45; elem(59,3)=45; elem(60,3)=46; elem(61,3)=46; elem(62,3)=47; elem(63,3)=47; elem(64,3)=48; elem(65,3)=48; elem(66,1)=36; elem(67,1)=37; elem(68,1)=37; elem(69,1)=38; elem(70,1)=38; elem(71,1)=39; elem(72,1)=39; elem(73,1)=40; elem(74,1)=40; elem(75,1)=41; elem(76,1)=41; elem(77,1)=42; elem(78,1)=42; elem(79,1)=43; elem(80,1)=43; elem(81,1)=44; elem(82,1)=44; elem(83,1)=45; elem(84,1)=45; elem(85,1)=46; elem(86,1)=46; elem(87,1)=47; elem(66,2)=37; elem(67,2)=51; elem(68,2)=38; elem(69,2)=52; elem(70,2)=39; elem(71,2)=53; elem(72,2)=40; elem(73,2)=54; elem(74,2)=41; elem(75,2)=55; elem(76,2)=42; elem(77,2)=56; elem(78,2)=43; elem(79,2)=57; elem(80,2)=44; elem(81,2)=58; elem(82,2)=45; elem(83,2)=59; elem(84,2)=46; elem(85,2)=60; elem(86,2)=47; elem(87,2)=61; elem(66,3)=50; elem(67,3)=50; elem(68,3)=51; elem(69,3)=51; elem(70,3)=52; elem(71,3)=52; elem(72,3)=53; elem(73,3)=53; elem(74,3)=54; elem(75,3)=54; elem(76,3)=55; elem(77,3)=55; elem(78,3)=56; elem(79,3)=56; elem(80,3)=57; elem(81,3)=57; elem(82,3)=58; elem(83,3)=58; elem(84,3)=59; elem(85,3)=59; elem(86,3)=60; elem(87,3)=60; 111 elem(88,1)=47; elem(89,1)=48; elem(90,1)=48; elem(91,1)=49; elem(88,2)=48; elem(89,2)=62; elem(90,2)=49; elem(91,2)=63; elem(88,3)=61; elem(89,3)=61; elem(90,3)=62; elem(91,3)=62; % -% Khai bao dieu kien bien nores=6; ixres(1)=13; vodof(13)=0; ixres(2)=20; vodof(20)=0; ixres(3)=41; vodof(41)=0; ixres(4)=69; vodof(69)=0; ixres(5)=97; vodof(97)=0; ixres(6)=125; vodof(125)=0; % -%CAC THONG SO BAI TOAN DONG %Khoi luong rieng den = 2.638*10^(-7); %kg.mm %Cac hang so NEWMARK gama = 1/2; bela = 1/4; %XAC DINH MA TRAN KHOI LUONG VA MA TRAN CAN KET CAU % Xac dinh ma tran khoi luong km = zeros(6,6); MOS=zeros(nodofos,nodofos); ix=zeros(1,nonpe*nodofpn); %ma tran cung cua ca he thong %vector chi so for ie=1:noe %bien chay phan tu endoe(1)=elem(ie,1); endoe(2)=elem(ie,2); endoe(3)=elem(ie,3); x1=coord(endoe(1),1); y1=coord(endoe(1),2); x2=coord(endoe(2),1); y2=coord(endoe(2),2); x3=coord(endoe(3),1); y3=coord(endoe(3),2); % -a1=x3-x2; a2=x1-x3; a3=x2-x1; b1=y2-y3; b2=y3-y1; b3=y1-y2; dtich=0.5*det([1 x1 y1;1 x2 y2;1 x3 y3]); ix=indexos(endoe,nonpe,nodofpn); % -%khoi luong phan tu m = dtich*d*den; % % Dinh nghia cac vetor v1 = [x2-x1, y2-y1]; v2 = [x3-x2, y3-y2]; v3 = [x1-x3,y1-y3]; %Tinh cac tich vecto vh1 = dot(v1,v2); vh2 = dot(v2,v3); 112 vh3 = dot(v3,v1); if vh1 == km(1,1)=m/4; km(2,2)=m/4; km(3,3)=m/2; km(4,4)=m/2; km(5,5)=m/4; km(6,6)=m/4; elseif vh2 == km(1,1)=m/4; km(2,2)=m/4; km(3,3)=m/4; km(4,4)=m/4; km(5,5)=m/2; km(6,6)=m/2 ; else km(1,1)=m/2; km(2,2)=m/2; km(3,3)=m/4; km(4,4)=m/4; km(5,5)=m/4; km(6,6)=m/4; end; % -km; MOS = Tomass(MOS,km,ix); end % Ma tran khoi luong tong the MOS; % % Xac dinh ma tran cung ban dau %ma tran cung cua ca he thong KOS=zeros(nodofos,nodofos); ix=zeros(1,nonpe*nodofpn); %vector chi so K=zeros(nonpe*nodofpn,nonpe*nodofpn); %ma tran cung cua tung phan tu module=zeros(4,noe); %module dan hoi % -%Chay phan tu for ie=1:noe %bien chay phan tu endoe(1)=elem(ie,1); endoe(2)=elem(ie,2); endoe(3)=elem(ie,3); x1=coord(endoe(1),1); y1=coord(endoe(1),2); x2=coord(endoe(2),1); y2=coord(endoe(2),2); x3=coord(endoe(3),1); y3=coord(endoe(3),2); % -a1=x3-x2; a2=x1-x3; a3=x2-x1; b1=y2-y3; b2=y3-y1; b3=y1-y2; dtich=0.5*det([1 x1 y1;1 x2 y2;1 x3 y3]); ix=indexos(endoe,nonpe,nodofpn); % -MTB=1/(2*dtich)*[b1 b2 b3 0; a1 a2 a3; a1 b1 a2 b2 a3 b3]; % D=Ec/(1-v^2)*[1 v 0;v 0;0 (1-v)/2]; % -K=d*dtich*MTB'*D*MTB; KOS=Tostiff(KOS,K,ix); end % Ma tran cung tong the ban dau 113 KOS; % Xac dinh ma tran can ket cau f=zeros(nodofos,1); f(126) =1; C = matrixC(MOS,KOS); % Xu ly rang buoc cho ma tran can [C,f] = straintC(C,f,ixres,vodof); %ma tran can he thong C; % PHAN TICH BAI TOAN -sbien=0; ii =1; charP = zeros(1,2000); charCV70 = zeros(1,2000); chart = zeros(1,2000); % -%Ham tai % Tai dang xung tam giac, xung dai t1 = 3; % sec , thoi gian phan tich t = 0.0001; %Buoc thoi gian phan tich dt = 0.0001; %sec % -% chuyen vi ban dau u = zeros(nodofos,1); du = zeros(nodofos,1); % van toc ban dau % gia toc ba dau ddu= zeros(nodofos,1); while t