Phân tích động lực học vỏ trụ có gân gia cường trên liên kết đàn hồi chịu tác dụng của hệ sóng xung kích

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHÒNG HỌC VIỆN KỸ THUẬT QUÂN SỰ PHÂN TÍCH ĐỘNG LỰC HỌC VỎ TRỤ CÓ GÂN GIA CƯỜNG TRÊN LIÊN KẾT ĐÀN HỒI CHỊU TÁC DỤNG CỦA HỆ SĨNG XUNG KÍCH LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI - 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHÒNG HỌC VIỆN KỸ THUẬT QUÂN SỰ PHÂN TÍCH ĐỘNG LỰC HỌC VỎ TRỤ CÓ GÂN GIA CƯỜNG TRÊN LIÊN KẾT ĐÀN HỒI CHỊU TÁC DỤNG CỦA HỆ SĨNG XUNG KÍCH Chun ngành: Cơ kỹ thuật Mã ngành: 9.52.01.01 LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI - 2022 i LỜI CAM ĐOAN Tơi , xin cam đoan cơng trình nghiên cứu riêng Các số liệu, kết luận án trung thực chưa cơng bố cơng trình ii LỜI CẢM ƠN Tác giả luận án xin bày tỏ lòng biết ơn chân thành GS.TS tận tình hướng dẫn, giúp đỡ cho nhiều dẫn khoa học có giá trị giúp cho tác giả hoàn thành luận án Tác giả trân trọng động viên, khuyến khích kiến thức khoa học chuyên môn mà thầy hướng dẫn chia sẻ cho tác giả nhiều năm qua, giúp cho tác giả nâng cao lực khoa học, phương pháp nghiên cứu lòng yêu nghề Tác giả trân trọng cảm ơn lãnh đạo Học viện Kỹ thuật quân sự, tập thể Bộ mơn Cơ học vật rắn, Phịng thí nghiệm Sức bền vật liệu, Khoa Cơ khí, phịng Sau đại học - Học viện Kỹ thuật quân tạo điều kiện thuận lợi, hợp tác trình nghiên cứu Tác giả xin trân trọng cảm ơn GS.TS.N – Học viện Kỹ thuật Quân sự, GS.TSKH.NGND Đào Huy Bích - Đại học Quốc gia Hà Nội, GS.TS - Đại học Thủy Lợi cung cấp cho tác giả nhiều tài liệu quý hiếm, kiến thức khoa học đại nhiều lời khuyên bổ ích, dẫn khoa học có giá trị để tác giả hồn thành luận án Cuối cùng, tác giả xin bày tỏ lòng biết ơn người thân gia đình thơng cảm, động viên chia sẻ khó khăn với tác giả suốt thời gian nghiên cứu, hoàn thành luận án Tác giả iii MỤC LỤC Trang Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục ký hiệu, chữ viết tắt .vii Danh mục bảng .xii Danh mục hình vẽ, đồ thị xiv Mở đầu CHƯƠNG TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Tải trọng sóng xung kích kết cấu vỏ chịu tác dụng tải trọng sóng xung kích 1.1.1 Tải trọng sóng xung kích 1.1.2 Một số kết nghiên cứu kết cấu chịu tác dụng tải trọng sóng xung kích 11 1.2 Tổng quan tính tốn kết cấu tấm, vỏ 18 1.2.1 Phương pháp giải tích 18 1.2.2 Phương pháp số 21 1.3 Kết đạt từ cơng trình cơng bố 24 1.4 Các vấn đề cần tiếp tục nghiên cứu 25 1.5 Các vấn đề luận án tập trung giải 25 1.6 Kết luận rút từ tổng quan 26 CHƯƠNG PHÂN TÍCH PHI TUYẾN ĐỘNG LỰC HỌC VỎ TRỤ CÓ GÂN GIA CƯỜNG TRÊN LIÊN KẾT ĐÀN HỒI CHỊU TÁC DỤNG CỦA HỆ SĨNG XUNG KÍCH DO NỔ 28 2.1 Đặt vấn đề 28 2.2 Đặt toán, giả thiết 28 iv 2.3 Ứng xử phi tuyến phần tử vỏ 29 2.3.1 Quan hệ biến dạng chuyển vị 29 2.3.2 Quan hệ ứng suất biến dạng 34 2.3.3 Các thành phần nội lực 34 2.3.4 Phương trình dao động phần tử vỏ 35 2.4 Ứng xử phi tuyến phần tử gân gia cường 41 2.4.1 Quan hệ ứng xử học gân dọc theo trục Ox 41 2.4.2 Quan hệ ứng xử học gân theo phương trục Oy 48 2.4.3 Quan hệ tương thích chuyển vị gân vỏ 50 2.5 Thiết lập phương trình mơ tả dao động phi tuyến phần tử vỏ có gân gia cường 51 2.5.1 Phương trình dao động phần tử vỏ có gân gia cường 51 2.5.2 Phương trình mơ tả dao động phần tử hệ tọa độ tổng thể 52 2.6 Phương trình dao động phi tuyến vỏ có gân gia cường liên kết đàn hồi 54 2.6.1 Xây dựng ma trận tổng thể kết cấu từ ma trận phần tử 54 2.6.2 Phương trình tổng thể 55 2.6.3 Véc tơ tải trọng hệ sóng xung kích tác dụng 56 2.6.4 Điều kiện biên 57 2.6.5 Phương trình mơ tả dao động phi tuyến vỏ có gân gia cường liên kết đàn hồi 58 2.7 Thuật toán giải phương trình dao động phi tuyến vỏ có gân gia cường liên kết đàn hồi 58 2.7.1 Bài toán dao động tự 58 2.7.2 Bài toán dao động cưỡng 59 2.8 Giới thiệu kiểm tra độ tin cậy chương trình tính 61 2.8.1 Giới thiệu chương trình tính 61 v 2.8.2 Kiểm tra độ tin cậy chương trình 66 2.9 Kết luận chương 69 CHƯƠNG KHẢO SÁT SỐ VÀ THẢO LUẬN 70 3.1 Đặt vấn đề 70 3.2 Tính tốn số 70 3.3 Ảnh hưởng số yếu tố đến đáp ứng động vỏ trụ thoải có gân gia cường chịu tác dụng hệ sóng xung kích 74 3.3.1 Ảnh hưởng loại phần tử mô vỏ 74 3.3.2 Ảnh hưởng cách bố trí gân 76 3.3.3 Ảnh hưởng kích thước gân (tỷ số hg/bg) 78 3.3.4 Ảnh hưởng cường độ tải trọng pm 80 3.3.5 Ảnh hưởng số lượng tải trọng 82 3.3.6 Ảnh hưởng thời gian chênh lớp sóng nổ 84 3.3.7 Ảnh hưởng độ cứng gối tựa đàn hồi 87 3.3.8 Ảnh hưởng điều kiện biên 89 3.3.9 Ảnh hưởng lỗ vỏ 91 3.3.10 Ảnh hưởng bán kính cong vỏ R 93 3.3.11 Ảnh hưởng chiều dày vỏ 96 3.4 Kết luận chương 98 CHƯƠNG NGHIÊN CỨU ĐÁP ỨNG ĐỘNG CỦA VỎ TRỤ CÓ GÂN GIA CƯỜNG CHỊU TÁC DỤNG CỦA SÓNG XUNG KÍCH BẰNG THỰC NGHIỆM 101 4.1 Mục đích thí nghiệm 101 4.2 Địa điểm, kết cấu thí nghiệm 101 4.2.1 Địa điểm thí nghiệm 101 4.2.2 Kết cấu thí nghiệm 102 4.3 Cơ sở lý thuyết thực thí nghiệm 104 vi 4.4 Thí nghiệm xác định tính vật liệu làm vỏ, gân độ cứng kéo (nén) lò xo đàn hồi 105 4.5 Thí nghiệm trường xác định đáp ứng động lực học vỏ 107 4.5.1 Thiết bị thí nghiệm 107 4.5.2 Sơ đồ bố trí thí nghiệm 109 4.5.3 Thực hành thí nghiệm đánh giá kết 110 4.6 Kết luận chương 117 Kết luận kiến nghị 118 Danh mục công trình tác giả 121 Tài liệu tham khảo 123 Phụ lục Mã nguồn chương trình vii DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT Danh mục ký hiệu 1.1 Các ký hiệu chữ La tinh Ae Diện tích phần tử vỏ, bg Bề rộng gân, cij Cosin phương góc hợp trục toạ độ tự nhiên xi trục tọa độ tổng thể Xj,  C Ma trận cản tổng thể hệ (chưa khử biên),  C Ma trận cản tổng thể hệ (sau khử biên), E Mô đun đàn hồi vật liệu, f  Véc tơ lực khối phần tử, f  Véc tơ lực bề mặt phần tử, f  Véc tơ lực tập trung phần tử,  f sb Véc tơ lực nút phần tử vỏ có gân gia cường e b e su e c e hệ tọa độ tổng thể, f0 Độ vồng vỏ, G Mô đun đàn hồi trượt vật liệu, h Chiều dày vỏ, He Hàm nội suy Hamilton, hg Chiều cao gân, [Ik] Ma trận đơn vị kích thước kk,  K uu   L   K uu   N e e Ma trận độ cứng tuyến tính uốn, Ma trận độ cứng phi tuyến uốn, viii  K QuQ    e e Ma trận độ cứng cắt, e  K  b ,  K  s Ma trận khối lượng phần tử 3D cong, phần tử vỏ  uu N Hệ số kể đến tính chất phi tuyến hệ, e  K  sb Ma trận độ cứng phần tử vỏ có gân gia cường hệ tọa độ phần tử,  K sbe Ma trận độ cứng phần tử vỏ có gân gia cường hệ tọa độ tổng thể, e  K  pill Ma trận độ cứng phần tử liên kết đàn hồi, K Ma trận độ cứng tổng thể hệ (chưa khử biên), K Ma trận độ cứng tổng thể hệ (sau khử biên),  K t t (i1) Ma trận độ cứng tiếp tuyến hiệu quả, e e  M  b ,  M  Ma trận khối lượng phần tử 3D cong, phần tử vỏ s e  M  sb Ma trận khối lượng phần tử vỏ có gân gia cường hệ tọa độ phần tử,  M sbe Ma trận khối lượng phần tử vỏ có gân gia cường hệ tọa độ tổng thể, M Ma trận khối lượng tổng thể (chưa khử biên), M Ma trận khối lượng tổng thể (sau khử biên), M  x, y, t  Véc tơ mô men uốn xoắn, Mx, My, Mxy Mô men uốn mô men xoắn phân bố theo chiều dài, N  Ma trận hàm dạng, MAT,1 LMESH,_Y1 ! -TYPE,1 REAL,1 MAT,1 AMESH,ALL 11 Chương trình Shell_curve1_without_stif CSYS,5 K,1,R1,(-Phi/2-90) K,2,R1,(Phi/2-90) K,3,R1,(-Phi/2-90),L1 K,4,R1,(Phi/2-90),L1 ESIZE,L1/14 A,1,3,4,2 TYPE,1 REAL,1 MAT,1 AMESH,ALL MESHING CSYS,0 12 Chương trình Shell_curve1_with_phistif CSYS,5 ! CYLINDRICAL CO-ORDINATE SYSTEM K,1,R1,phi/(NPHI-1)-90-phi/2,0 K,2,R1,-90-phi/2,0 K,3,R1,-90-phi/2,L/(Nlen-1) K,4,R1,phi/(NPHI-1)-90-phi/2,L/(Nlen-1) A,1,2,3,4 ! NDIV_PHI=18 !khoang chia NDIV_Y=2 !* FLST,5,2,4,ORDE,2 FITEM,5,1 FITEM,5,3 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,NDIV_PHI, , , , ,1 FLST,5,2,4,ORDE,2 FITEM,5,2 FITEM,5,4 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,NDIV_Y, , , , ,1 ! -TYPE,1 REAL,1 MAT,1 AMESH,1 vo FLST,3,1,5,ORDE,1 FITEM,3,1 AGEN,(Nlen-1),P51X, , , ,,L/(Nlen-1),,0 ! -NUMMRG,ALL, , , ,LOW NUMCMP,ALL ! chia phan tu gan FLST,5,20,4,ORDE,20 FITEM,5,1 *DO,k,2,20,1 FITEM,5,3*(k-1) *ENDDO CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y TYPE,3 SECNUM,1 MAT,1 LMESH, _Y1 13 Chương trình Shell_curve1_with_lenstif CSYS,5 ! CYLINDRICAL CO-ORDINATE SYSTEM K,1,R1,Phi/(NPHI-1)-90-Phi/2,0 K,2,R1,-90-Phi/2,0 K,3,R1,-90-Phi/2,L/(Nlen-1) K,4,R1,Phi/(NPHI-1)-90-Phi/2,L/(Nlen-1) A,1,2,3,4 ! NDIV_PHI=2 !khoang chia NDIV_Y=38 !* FLST,5,2,4,ORDE,2 FITEM,5,1 FITEM,5,3 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,NDIV_PHI, , , , ,1 ! -Thiet lap mo hinh khoi co ban -FLST,5,2,4,ORDE,2 FITEM,5,2 FITEM,5,4 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,NDIV_Y, , , , ,1 ! -! -chia phan tu vo TYPE,1 REAL,1 MAT,1 AMESH,1 ! Tao khoi phan phan vo !FLST,3,1,5,ORDE,1 !FITEM,3,1 !AGEN,(Nlen-1),P51X, , , ,,L/(Nlen-1),,0 ! FLST,3,1,5,ORDE,1 FITEM,3,1 AGEN,NPHI-1,P51X, , ,,Phi/(NPHI-1), , ,0 ! -NUMMRG,ALL, , , ,LOW NUMCMP,ALL ! chia phan tu gan FLST,5,10,4,ORDE,10 FITEM,5,2 FITEM,5,4 FITEM,5,7 FITEM,5,10 FITEM,5,13 FITEM,5,16 FITEM,5,19 FITEM,5,22 FITEM,5,25 FITEM,5,28 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y TYPE,3 SECNUM,1 MAT,1 LMESH,_Y1 14 Chương trình Shell_curve1_with_lpstif CSYS,5 ! CYLINDRICAL CO-ORDINATE SYSTEM K,1,R1,phi/(NPHI-1)-90-phi/2,0 K,2,R1,-90-phi/2,0 K,3,R1,-90-phi/2,L/(Nlen-1) K,4,R1,phi/(NPHI-1)-90-phi/2,L/(Nlen-1) A,1,2,3,4 ! NDIV_PHI=2 !khoang chia NDIV_Y=2 !* FLST,5,2,4,ORDE,2 FITEM,5,1 FITEM,5,3 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,NDIV_PHI, , , , ,1 FLST,5,2,4,ORDE,2 FITEM,5,2 FITEM,5,4 10 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,NDIV_Y, , , , ,1 TYPE,1 REAL,1 MAT,1 AMESH,1 vo -FLST,3,1,5,ORDE,1 FITEM,3,1 AGEN,(Nlen-1),P51X, , , ,,L/(Nlen-1),,0 ! FLST,3,19,5,ORDE,2 FITEM,3,1 FITEM,3,-19 AGEN,NPHI-1,P51X, , ,,phi/(NPHI-1), , ,0 ! -NUMMRG,ALL, , , ,LOW NUMCMP,ALL ! chia phan tu gan TYPE,3 SECNUM,1 MAT,1 LMESH,ALL 15 Chương trình BoundaryModel_Plate MULTIPRO,'START',9 *CSET,1,3,B1,'SFSF',1.0 *CSET,4,6,B2,'FSFS',0 *CSET,7,9,B3,'SSSS',0 *CSET,10,12,B4,'CFCF',0 *CSET,13,15,B5,'FCFC',0 *CSET,16,18,B6,'CCCC',0 *CSET,19,21,B7,'EEEE',0 *CSET,22,24,B8,'EFEF',0 *CSET,25,27,B9,'FEFE',0 *CSET,61,62,'ENTER NUMBER INTO THE ','BOUNDARY TYPE' *CSET,63,64,'X = 0; X = B1; Y = 0; Y = L1','' MULTIPRO,'END' *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF 16 Chương trình BoundaryModel_curve1 MULTIPRO,'START',10 *CSET,1,3,B1,'SFSF',0 *CSET,4,6,B2,'FSFS',0 *CSET,7,9,B3,'SSSS',0 *CSET,10,12,B4,'CFCF',0 *CSET,13,15,B5,'FCFC',0 *CSET,16,18,B6,'CCCC',0 *CSET,19,21,B7,'EEEE',0 *CSET,22,24,B8,'EFEF',0 *CSET,25,27,B9,'FEFE',0 *CSET,28,30,B10,'SPR_leng',1 11 *CSET,31,33,B11,'SPR_phi',0 *CSET,34,36,B12,'SPR_lp',0 *CSET,61,62,'ENTER NUMBER INTO THE ','BOUNDARY TYPE' *CSET,63,64,'X = -Phi/2; X = Phi/2; Y = 0; Y = L1','' MULTIPRO,'END' *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF 17 Chương trình Boundaries_Plate *IF,B1,EQ,1,THEN NSEL,S,LOC,X,0 D,ALL,ALL NSEL,S,LOC,X,L D,ALL,ALL NSEL,S,LOC,Y,0 D,ALL,ALL NSEL,S,LOC,Y,B D,ALL,ALL *ELSEIF,B2,EQ,1,THEN NSEL,S,LOC,Y,0 D,ALL,ALL NSEL,S,LOC,Y,B D,ALL,ALL *ELSEIF,B3,EQ,1,THEN NSEL,S,LOC,X,0 D,ALL,ALL NSEL,S,LOC,X,L D,ALL,ALL *ELSEIF,B4,EQ,1,THEN NSEL,S,LOC,X,0 D,ALL,UX D,ALL,UY D,ALL,UZ NSEL,S,LOC,X,L D,ALL,UX D,ALL,UY D,ALL,UZ NSEL,S,LOC,Y,0 D,ALL,UX D,ALL,UY D,ALL,UZ NSEL,S,LOC,Y,B D,ALL,UX D,ALL,UY D,ALL,UZ *ELSEIF,B5,EQ,1,THEN NSEL,S,LOC,Y,0 D,ALL,UX D,ALL,UY D,ALL,UZ NSEL,S,LOC,Y,B D,ALL,UX D,ALL,UY D,ALL,UZ *ELSEIF,B6,EQ,1,THEN 12 NSEL,S,LOC,X,0 D,ALL,UX D,ALL,UY D,ALL,UZ NSEL,S,LOC,X,L D,ALL,UX D,ALL,UY D,ALL,UZ *ENDIF ALLSEL FINISH 18 Chương trình Boundaries_curve1 !CSYS,5 !NSEL,S,LOC,Y,-90-PHI/2 !NSEL,R,LOC,X,R1 !CM,_N1,NODE !NSEL,S,LOC,z,0 !CSYS,0 !NGEN,2,3000,_N1,,,,,-L3 *DIM,T,ARRAY,20,1,,,, *SET,T(1,1,1),2154 *SET,T(2,1,1),2142 *SET,T(3,1,1),2130 *SET,T(4,1,1),2118 *SET,T(5,1,1),2106 *SET,T(6,1,1),2094 *SET,T(7,1,1),2082 *SET,T(8,1,1),2070 *SET,T(9,1,1),2058 *SET,T(10,1,1),2046 *SET,T(11,1,1),2034 *SET,T(12,1,1),2022 *SET,T(13,1,1),2010 *SET,T(14,1,1),1998 *SET,T(15,1,1),1986 *SET,T(16,1,1),1974 *SET,T(17,1,1),1962 *SET,T(18,1,1),1950 *SET,T(19,1,1),1938 *SET,T(20,1,1),1934 ! *DIM,S,ARRAY,20,1,, , , *SET,S(1,1,1),2 *SET,S(2,1,1),6 *SET,S(3,1,1),22 *SET,S(4,1,1),38 *SET,S(5,1,1),54 *SET,S(6,1,1),70 *SET,S(7,1,1),86 *SET,S(8,1,1),102 *SET,S(9,1,1),118 *SET,S(10,1,1),134 *SET,S(11,1,1),150 *SET,S(12,1,1),166 *SET,S(13,1,1),182 *SET,S(14,1,1),198 13 *SET,S(15,1,1),214 *SET,S(16,1,1),230 *SET,S(17,1,1),246 *SET,S(18,1,1),262 *SET,S(19,1,1),278 *SET,S(20,1,1),294 ! -CSYS,0 TYPE,2 REAL,2 *DO,k,1,20,1 NGEN,2,3000,T(K,1,1),,,,,-L3 NGEN,2,3000,S(K,1,1),,,,,-L3 E,T(K,1,1),T(K,1,1)+3000 E,S(K,1,1),S(K,1,1)+3000 D,T(K,1,1)+3000,UX,0 D,T(K,1,1)+3000,UY,0 D,T(K,1,1)+3000,UZ,0 D,S(K,1,1)+3000,UX,0 D,S(K,1,1)+3000,UY,0 D,S(K,1,1)+3000,UZ,0 *ENDDO CSYS,5 NSEL,S,LOC,Y,-90-PHI/2 D,ALL,UX,0 ! NSEL,S,LOC,Y,-90+PHI/2 D,ALL,UX,0 ! NSEL,S,LOC,z,0 NSEL,R,LOC,X,R1 D,ALL,Uy,0 NSEL,S,LOC,z,L1 NSEL,R,LOC,X,R1 D,ALL,Uy,0 NSEL,ALL ALLSEL EPLOT FINISH 19 Chương trình Modal_Analysis /SOLU ANTYPE,MODAL MODOPT,LANB,8 SOLVE FINISH 20 Chương trình Trans_Analysis /SOLU ANTYPE,MODAL MODOPT,LANB,2 SOLVE *GET,F1,MODE,1,FREQ *GET,F2,MODE,2,FREQ FINISH 14 /OUTPUT,THUYHVKT MULTIPRO,'START',5 *CSET,1,3,T1,'P0 = constant at center',0 *CSET,4,6,T2,'pr = constant',0 *CSET,7,9,T3,'P0 = Psin(omegat) at center',1 *CSET,10,12,T4,'pr = Prsin(omegat)',0 *CSET,13,15,T5,'Blast loading',0 *CSET,61,62,'SELECT THE TYPE OF LOAD EFFECT','' *CSET,63,64,'ENTER NUMBER INTO THE ','CORRESPONDING BOX' MULTIPRO,'END' !* *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF 21 Chương trình Analysis_Trans_P0 MULTIPRO,'START',1 *CSET,1,3,P0,'Load amplitude[N]:',100 *CSET,61,62,'DEFINES THE LOAD PARAMETER','' *CSET,63,64,'Pt = P0','at center' MULTIPRO,'END' *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF /SOLU ANTYPE,0 OUTPR,BASIC,LAST NLGEOM,ON NSEL,S,Loc,X,B1/2 NSEL,R,Loc,y,L1/2 CM,_Z1,NODE /UNITS,USER F,_Z1,FZ,P0 SOLVE 22 Chương trình Analysis_Trans_pr MULTIPRO,'START',1 *CSET,1,3,pt0,'Load amplitude[N/m2]:',100 *CSET,61,62,'DEFINES THE LOAD PARAMETER','' *CSET,63,64,'Pt = Pt0','at area' MULTIPRO,'END' *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF /SOLU ANTYPE,0 OUTPR,BASIC,LAST NLGEOM,ON ASEL,S,Loc,Z,0 CM,_Z1,AREA /UNITS,USER SFA,_Z1,,PRES,pt0 SOLVE 23 Chương trình Analysis_Trans_P0sin MULTIPRO,'START',2 *CSET,1,3,P0,'Load amplitude[N]:',100 *CSET,4,6,Omega,'Omega = 2.pi.f',10 15 *CSET,61,62,'DEFINES THE LOAD PARAMETER','' *CSET,63,64,'Pt = P0sin(Omega.t)','at center' MULTIPRO,'END' *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF /SOLU ANTYPE,TRANS !NLGEOM,OFF NLGEOM,ON /UNITS,USER TRNOPT,FULL GXY=0.05 DAMPBETAD=2*GXY/(40*(F1+F2)) DAMPALPHAD=DAMPBETAD*40*F1*F2 ALPHAD,DAMPALPHAD BETAD,DAMPBETAD TIMINT,OFF TIME,1E-9 SOLVE NSEL,S,Loc,X,B1/2 NSEL,R,Loc,y,L1/2 CM,_Z1,NODE OUTRES,ALL,ALL TIMINT,ON *DO,I,1,50,1 TIME,I*deltat NSUBST,2 F,_Z1,FZ,-P0*SIN(Omega*I*deltat) solve *ENDDO 24 Chương trình Analysis_Trans_prsin MULTIPRO,'START',2 *CSET,1,3,P0,'Load amplitude[N]:',100 *CSET,4,6,Omega,'Omega = 2.pi.f',10 *CSET,61,62,'DEFINES THE LOAD PARAMETER','' *CSET,63,64,'Pt = P0sin(Omega.t)','at center' MULTIPRO,'END' *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF /SOLU ANTYPE,TRANS !NLGEOM,OFF NLGEOM,ON /UNITS,USER TRNOPT,FULL GXY=0.05 DAMPBETAD=2*GXY/(40*(F1+F2)) 16 DAMPALPHAD=DAMPBETAD*40*F1*F2 ALPHAD,DAMPALPHAD BETAD,DAMPBETAD TIMINT,OFF TIME,1E-9 SOLVE ASEL,S,Loc,Z,0 CM,_Z1,AREA OUTRES,ALL,ALL TIMINT,ON *DO,I,1,50,1 TIME,I*deltat NSUBST,2 SFA,_Z1,,PRES,-Pt0*SIN(Omega*I*deltat) solve *ENDDO 25 Chương trình Analysis_Trans_blast MULTIPRO,'START',2 *CSET,1,3,P0,'Load amplitude[N]:',100 *CSET,4,6,Omega,'Omega = 2.pi.f',10 *CSET,61,62,'DEFINES THE LOAD PARAMETER','' *CSET,63,64,'Pt = P0sin(Omega.t)','at center' MULTIPRO,'END' *IF,_BUTTON,EQ,1,THEN /EOF *ENDIF /SOLU ANTYPE,TRANS !NLGEOM,OFF NLGEOM,ON /UNITS,USER TRNOPT,FULL GXY=0.05 DAMPBETAD=2*GXY/(40*(F1+F2)) DAMPALPHAD=DAMPBETAD*40*F1*F2 ALPHAD,DAMPALPHAD BETAD,DAMPBETAD TIMINT,OFF TIME,1E-9 SOLVE !CSYS,5 ASEL,S,LOC,X,R1 CM,_Z1,AREA !ASEL,ALL Teta1=0.012 Teta2=0.012 17 Nt=2 !So dot SXK tac dung