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
Trang 1Đề bài: Tính toán kết cấu tấm đồng nhất đẳng hướng chịu tác dụng của tải trọng sóng xung
kích
Kích thước tấm:
- Chiều dài: 8 m
- Chiều rộng: 3 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 đều lên bề mặt tấm Thời gian tác dụng của mỗi lớp sóng là 0,05s Siêu áp mặt sóng: Pmax = 3.105 N/m2
Mô hình tấm
Kết quả yêu cầu: Chuyển vị và ứng suất tại điểm giữa tấm
Chuyển vị tại điểm giữa tấm
Trang 2Ứ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 -
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
Trang 3!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
Trang 4D,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 do
! -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 do cung
! -
!TIMINT,OFF ! Tat/mo hieu ung qua do
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 ra 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))
Trang 5*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