标题: 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