标题: ANSYS中的梁格法 [打印本页]
作者： even2004 标题:ANSYS中的梁格法 时间： 2009/04/30/08:27
/prep7
!DEFINE THE ELEMENTARY PARAMETERS
*DIM,L,ARRAY,10
*DIM,DISTC,ARRAY,10
N=5
NBOX=4
L(1)=20
L(2)=30
L(3)=40
L(4)=30
L(5)=20
R=100
NUML=NBOX+3
NUMC=(L(1)+L(2)+L(3)+L(4)+L(5))+1
DISTC(1)=1.1375
DISTC(2)=1.850
DISTC(3)=2.000
DISTC(4)=0.0
DISTC(5)=-2.00
DISTC(6)=-1.85
DISTC(7)=-1.1375
LOCAL,11,1,0,0,0,,90
CSYS,11
!DEFINE THE NODES OF BRIDGE
N, 1, R, -80.214, 0.00
N, NUMC, R, 0, 0.00
FILL,1,NUMC
N, 200+1, R+DISTC(1)+DISTC(2)+DISTC(3), -80.214, 0.00
N, 200+NUMC,R+DISTC(1)+DISTC(2)+DISTC(3), 0, 0.00
FILL,201,200+NUMC
N, 400+1, R+DISTC(2)+DISTC(3), -80.214, 0.00
N, 400+NUMC,R+DISTC(2)+DISTC(3), 0, 0.00
FILL,401,400+NUMC
N, 600+1, R+DISTC(3), -80.214, 0.00
N, 600+NUMC,R+DISTC(3), 0, 0.00
FILL,601,600+NUMC
N, 800+1, R+DISTC(5), -80.214, 0.00
N, 800+NUMC,R+DISTC(5), 0, 0.00
FILL,801,800+NUMC
N, 1000+1, R+DISTC(5)+DISTC(6), -80.214, 0.00
N, 1000+NUMC,R+DISTC(5)+DISTC(6), 0, 0.00
FILL,1001,1000+NUMC
N, 1200+1, R+DISTC(5)+DISTC(6)+DISTC(7), -80.214, 0.00
N, 1200+NUMC,R+DISTC(5)+DISTC(6)+DISTC(7), 0, 0.00
FILL,1201,1200+NUMC
! 定义纵梁单元材料、几何参数
ET,1,BEAM4
MP,EX , 1, 3.45E+10
MP,NUXY, 1, 0.2000
MP,DENS, 1, 2600
MP,ALPX, 1, 1.00E-05
R,1,0.1812,0.1126,0.007939,0.725,0.25,,
RMORE ,,0.002964,,,,,
ET,2,BEAM4
MP,EX , 2, 3.45E+10
MP,NUXY, 2, 0.2000
MP,DENS, 2, 2600
MP,ALPX, 2, 1.00E-05
R,2,1.1813,0.6270,0.1920,1.85,1.0,,
RMORE ,,0.0306,,,,,
ET,3,BEAM4
MP,EX , 3, 3.45E+10
MP,NUXY, 3, 0.2000
MP,DENS, 3, 2600
MP,ALPX, 3, 1.00E-05
R,3,1.3750,0.8081,0.2670,1.85,1.0,,
RMORE ,,0.0351,,,,,
! 定义横梁单元材料、几何参数
ET,4,BEAM4
MP,EX , 4, 3.45E+10
MP,NUXY, 4, 0.2000
MP,DENS, 4, 0
MP,ALPX, 4, 1.00E-05
R,4,0.500,0.3854,0.0417,1.0,1.0,,
RMORE ,,0.00878,,,,,
ET,5,BEAM4
MP,EX , 5, 3.45E+10
MP,NUXY, 5, 0.2000
MP,DENS, 5, 0
MP,ALPX, 5, 1.00E-05
R,5,1.250,0.5260,0.0573,1.0,1.0,,
RMORE ,,0.0798,,,,,
! 生成纵向梁格单元
TYPE,1
MAT,1
REAL,1
*DO,I,1,NUMC-1,1
E,200+I,200+I+1
*ENDDO
*DO,I,1,NUMC-1,1
E,1200+I,1200+I+1
*ENDDO
TYPE,2
MAT,2
REAL,2
*DO,I,1,NUMC-1,1
E,400+I,400+I+1
*ENDDO
*DO,I,1,NUMC-1,1
E,1000+I,1000+I+1
*ENDDO
TYPE,3
MAT,3
REAL,3
*DO,I,1,NUMC-1,1
E,600+I,600+I+1
*ENDDO
*DO,I,1,NUMC-1,1
E,800+I,800+I+1
*ENDDO
*DO,I,1,NUMC-1,1
E,I,I+1
*ENDDO
! 生成横向梁格单元
TYPE,5
MAT,5
REAL,5
E,200+1,400+1
E,400+1,600+1
E,600+1,1
E,1,800+1
E,800+1,1000+1
E,1000+1,1200+1
TYPE,4
MAT,4
REAL,4
*DO,I,2,L(1),1
E,200+I,400+I
E,400+I,600+I
E,600+I,I
E,I,800+I
E,800+I,1000+I
E,1000+I,1200+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+1,400+L(1)+1
E,400+L(1)+1,600+L(1)+1
E,600+L(1)+1,L(1)+1
E,L(1)+1,800+L(1)+1
E,800+L(1)+1,1000+L(1)+1
E,1000+L(1)+1,1200+L(1)+1
TYPE,4
MAT,4
REAL,4
*DO,I,L(1)+2,L(1)+L(2),1
E,200+I,400+I
E,400+I,600+I
E,600+I,I
E,I,800+I
E,800+I,1000+I
E,1000+I,1200+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+L(2)+1,400+L(1)+L(2)+1
E,400+L(1)+L(2)+1,600+L(1)+L(2)+1
E,600+L(1)+L(2)+1,L(1)+L(2)+1
E,L(1)+L(2)+1,800+L(1)+L(2)+1
E,800+L(1)+L(2)+1,1000+L(1)+L(2)+1
E,1000+L(1)+L(2)+1,1200+L(1)+L(2)+1
TYPE,4
MAT,4
REAL,4
*DO,I,L(1)+L(2)+2,L(1)+L(2)+L(3),1
E,200+I,400+I
E,400+I,600+I
E,600+I,I
E,I,800+I
E,800+I,1000+I
E,1000+I,1200+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+L(2)+L(3)+1,400+L(1)+L(2)+L(3)+1
E,400+L(1)+L(2)+L(3)+1,600+L(1)+L(2)+L(3)+1
E,600+L(1)+L(2)+L(3)+1,L(1)+L(2)+L(3)+1
E,L(1)+L(2)+L(3)+1,800+L(1)+L(2)+L(3)+1
E,800+L(1)+L(2)+L(3)+1,1000+L(1)+L(2)+L(3)+1
E,1000+L(1)+L(2)+L(3)+1,1200+L(1)+L(2)+L(3)+1
TYPE,4
MAT,4
REAL,4
*DO,I,L(1)+L(2)+L(3)+2,L(1)+L(2)+L(3)+L(4),1
E,200+I,400+I
E,400+I,600+I
E,600+I,I
E,I,800+I
E,800+I,1000+I
E,1000+I,1200+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+L(2)+L(3)+L(4)+1,400+L(1)+L(2)+L(3)+L(4)+1
E,400+L(1)+L(2)+L(3)+L(4)+1,600+L(1)+L(2)+L(3)+L(4)+1
E,600+L(1)+L(2)+L(3)+L(4)+1,L(1)+L(2)+L(3)+L(4)+1
E,L(1)+L(2)+L(3)+L(4)+1,800+L(1)+L(2)+L(3)+L(4)+1
E,800+L(1)+L(2)+L(3)+L(4)+1,1000+L(1)+L(2)+L(3)+L(4)+1
E,1000+L(1)+L(2)+L(3)+L(4)+1,1200+L(1)+L(2)+L(3)+L(4)+1
TYPE,4
MAT,4
REAL,4
*DO,I,L(1)+L(2)+L(3)+L(4)+2,L(1)+L(2)+L(3)+L(4)+L(5),1
E,200+I,400+I
E,400+I,600+I
E,600+I,I
E,I,800+I
E,800+I,1000+I
E,1000+I,1200+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+L(2)+L(3)+L(4)+L(5)+1,400+L(1)+L(2)+L(3)+L(4)+L(5)+1
E,400+L(1)+L(2)+L(3)+L(4)+L(5)+1,600+L(1)+L(2)+L(3)+L(4)+L(5)+1
E,600+L(1)+L(2)+L(3)+L(4)+L(5)+1,L(1)+L(2)+L(3)+L(4)+L(5)+1
E,L(1)+L(2)+L(3)+L(4)+L(5)+1,800+L(1)+L(2)+L(3)+L(4)+L(5)+1
E,800+L(1)+L(2)+L(3)+L(4)+L(5)+1,1000+L(1)+L(2)+L(3)+L(4)+L(5)+1
E,1000+L(1)+L(2)+L(3)+L(4)+L(5)+1,1200+L(1)+L(2)+L(3)+L(4)+L(5)+1
/SOL
ACEL,,9.81,0
!*********************************************************************************
!
! 定义边界条件
!
!*********************************************************************************
D,600+1,UX,0,,,,UY,UZ
D,800+1,UX,0,,,,UY,UZ
D,21,UX,0,,,,UY,UZ
D,51,UX,0,,,,UY,UZ
D,91,UX,0,,,,UY,UZ
D,121,UX,0,,,,UY,UZ
D,600+141,UX,0,,,,UY,UZ
D,800+141,UX,0,,,,UY,UZ
ANTYPE,0
/STATUS,SOLU
SOLVE
/POST1
/REP,FAST
AVPRIN,0, ,
ETABLE, ,SMISC, 2
!*
AVPRIN,0, ,
ETABLE, ,SMISC,8
!*
AVPRIN,0, ,
ETABLE, ,SMISC, 4
!*
AVPRIN,0, ,
ETABLE, ,SMISC, 5
!*
AVPRIN,0, ,
ETABLE, ,SMISC, 6
!*
AVPRIN,0, ,
ETABLE, ,SMISC, 10
!*
AVPRIN,0, ,
ETABLE, ,SMISC, 11
!*
AVPRIN,0, ,
ETABLE, ,SMISC, 12
!*
FINISH
阳光精品论坛 - Powered By BMForum v5.0