Thông tin tài liệu
BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ QUỐC PHÒNG HỌC VIỆN KỸ THUẬT QUÂN SỰ Lê Xuân Thùy 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 Lê Xuân Thùy 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 Người hướng dẫn khoa học: GS.TS Nguyễn Thái Chung HÀ NỘI - 2022 i LỜI CAM ĐOAN Tôi Lê Xuân Thùy, xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết luận án trung thực chưa cơng bố cơng trình TÁC GIẢ Lê Xuân Thùy 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 Nguyễn Thái Chung 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ả hồ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.NGND Hoàng Xuân Lượng cố GS.TS Vũ Đình Lợi – 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 Nguyễn Văn Lệ - Đạ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ả 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 .5 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 .5 1.1.1 Tải trọng sóng xung kích .5 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 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 tố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 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 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 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 e Véc tơ lực bề mặt phần tử, f Véc tơ lực tập trung phần tử, f sbe Véc tơ lực nút phần tử vỏ có gân gia cường e b su e c 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 kk, Kuu L e e Kuu N Ma trận độ cứng tuyến tính uốn, Ma trận độ cứng phi tuyến uốn, QuQ K e K K s b e e , uu N Ma trận độ cứng cắt, Ma trận khối lượng phần tử 3D cong, phần tử vỏ Hệ số kể đến tính chất phi tuyến hệ, Ma trận độ cứng phần tử vỏ có gân gia cường hệ K e sb tọa độ phần tử, Ma trận độ cứng phần tử vỏ có gân gia cường K sb e hệ tọa độ tổng thể, K e pill K Ma trận độ cứng phần tử liên kết đàn hồi, Ma trận độ cứng tổng thể hệ (chưa khử biên), Ma trận độ cứng tổng thể hệ (sau khử biên), K K tt (i1 ) Ma trận độ cứng tiếp tuyến hiệu quả, M M Ma trận khối lượng phần tử 3D cong, phần tử vỏ s e, b e Ma trận khối lượng phần tử vỏ có gân gia cường M e sb hệ tọa độ phần tử, M e sb 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, 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 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 *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 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 *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 /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 *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)) 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 Nt=2 !So dot SXK tac dung
Ngày đăng: 06/12/2022, 07:15
Xem thêm:
