Bộ chương trình tính toán kết cấu tấm làm từ vật liệu đồng nhất đẳng hướng chịu tác dụng của hệ sóng xung kích được lập trình trên nền ngôn ngữ lập trình của phần mềm ANSYS. Với công cụ lập trình này, việc giải quyết bài toán tính kết cấu tấm trở nên ngắn gọn và dễ dàng (so với sử dụng công cụ lập trình Matlab). Điều kiện liên kết của tấm có thể thay đổi dễ dàng bằng cách thay đổi giá trị của biến Bien. Kết quả tính cho độ tin cậy cao và đa dạng (có thể xuất mọi kết quả theo yêu cầu của bài toán cơ học). Qua đây, tác giả giới thiệu phương pháp lập trình cho dạng tải trọng động trong ngôn ngữ lập trình ANSYS. Liên hệ thêm để được hướng dẫn chi tiết: thuylxmtagmail.com
Đề bài: Tính toán kết cấu đồng đẳng hướng chịu tác dụng tải trọng sóng xung kích Kích thước tấm: - Chiều dài: m - Chiều rộng: m - Bề dày: 0,05 m Vật liệu tấm: thép Liên kết: Hai cạnh ngắn ngàm Tải trọng: Hai lớp sóng xung kích tác dụng liên tiếp nhau, tác dụng lên bề mặt Thời gian tác dụng lớp sóng 0,05s Siêu áp mặt sóng: Pmax = 3.105 N/m2 Mô hình Kết yêu cầu: Chuyển vị ứng suất điểm Chuyển vị điểm Ứng suất Xicmax Ứng suất Xicmay Chương trình tính: /PREP7 smrt,off /TITLE, ANALYSIS OF PLATE UNDER BLAST LOADING !/TITLE, WRITED BY LE XUAN THUY ! JULY 2013 KEYW,PR_SET,1 KEYW,PR_STRUC,1 /GO /VIEW,1,1,1,1 /VUP,1,Z /REPLOT ! THONG SO HINH HOC TAM L=8 !m - LENGTH OF PALTE B=3 !m - WIDTH OF PALTE TH=0.05 !m - THICKNESS OF PLATE ! CHON DIEU KIEN BIEN CHO TAM !BIEN=1 !4CANH NGAM !BIEN=2 !2CANH Y=CONST NGAM,2CANH X=CONST TU DO BIEN=3 !2CANH X=CONST NGAM,2CANH Y=CONST TU DO !BIEN=4 !4CANH CHIU LIEN KET BAN LE !BIEN=5 !2CANH Y=CONST BAN LE,2CANH X=CONST TU DO !BIEN=6 !2CANH X=CONST BAN LE,2CANH Y=CONST TU DO ! ET,1,SHELL63,,,,,,2 R,1,TH MP,EX,1,2.1E11 !N/m2 MP,NUXY,1,0.3 !MP,ALPX,1,0.1E-5 MP,DENS,1,7.850E3 !Kg/m3 K,1,0,0,0 K,2,0,B,0 K,3,L,B,0 K,4,L,0,0 A,1,2,3,4 ESIZE,0.25 AMESH,ALL MESHING ! DIEU KIEN BIEN *IF,BIEN,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,BIEN,EQ,2,THEN NSEL,S,LOC,Y,0 D,ALL,ALL NSEL,S,LOC,Y,B D,ALL,ALL *ELSEIF,BIEN,EQ,3,THEN NSEL,S,LOC,X,0 D,ALL,ALL NSEL,S,LOC,X,L D,ALL,ALL *ELSEIF,BIEN,EQ,4,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,BIEN,EQ,5,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,BIEN,EQ,6,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 ! MODAL ANALYSIS /SOLU NSEL,ALL antype,2 modopt,lanb,2 !mxpand,2 solve *get,f1,mode,1,freq *get,f2,mode,2,freq FINISH ! APPLY BLAST LOADING !/OUTPUT,SXK /SOLU ANTYPE,TRANS !NLGEOM,ON /UNITS,USER TRNOPT,FULL !Cu the hoa lua chon phan tich qua ! Dat tai txy=0.05 dampbetad=2*txy/(40*(f1+f2)) dampalphad=dampbetad*40*f1*f2 alphad,dampalphad !He so can khoi luong betad,dampbetad !He so can cung ! !TIMINT,OFF ! Tat/mo hieu ung qua TIME,1E-9 SOLVE ! pmax1=3E3 ! pmax pmax2=3E5 ! pmax teta1=0.05 ! teta2=0.05 ! epsilon=0.025 ! deltat=0.0025 ! Buoc thoi gian ! Bat dau lap -OUTRES,ALL,ALL !Dieu khien viec ghi du lieu database TIMINT,ON *DO,I,1,200,1 TIME,I*deltat NSUBST,2 *IF,(I*deltat),LT,epsilon,THEN SFA,1,,PRES,pmax1*(1-((I-1)*deltat/teta1)) *ELSEIF,(I*deltat),GE,epsilon,AND,(I*deltat),LE,teta1,THEN SFA,1,,PRES,(pmax1*(1-((I-1)*deltat/teta1))+ pmax2*(1-((I-41)*deltat/teta2))) *ELSEIF,(I*deltat),GT,teta1,AND,(I*deltat),LE,teta2,THEN SFA,1,,PRES,pmax2*(1-((I-41)*deltat/teta2)) *ELSE SFA,1,,PRES,0 *ENDIF solve *ENDDO /POST26 NSOL,2,258,U,Z, UZ_2 STORE,MERGE XVAR,1 PLVAR,2, FORCE,TOTAL SHELL,TOP EPLOT !* !* ESOL,3,187,258 ,S,X,SX_3 STORE,MERGE FORCE,TOTAL SHELL,TOP !* !* ESOL,4,187,258 ,S,Y,SY_4 STORE,MERGE