LX Thuy_Chuyên ngành Cơ học ứng dụng, Cơ kỹ thuật: Tính toán kết cấu tấm đồng nhất đẳng hướng chịu tác dụng của hệ sóng xung kích

Đề 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


Ứ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 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))


*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