/prep7
!DEFINE THE ELEMENTARY PARAMETERS
*DIM,L,ARRAY,10
*DIM,DISTC,ARRAY,10
N=4
NBOX=2
L(1)=20
L(2)=30
L(3)=30
L(4)=20
NUML=NBOX+3
NUMC=(L(1)+L(2)+L(3)+L(4))+1
CITA=25/180*3.1415926
DISTC(1)=2.9584
DISTC(2)=4.0802
DISTC(3)=0.00
DISTC(4)=4.0802
DISTC(5)=2.9584
!DEFINE THE NODES OF BRIDGE
N, 1, 0, 0, 0.00
N, NUMC, NUMC-1, 0, 0.00
FILL,1,NUMC
N, 200+1,
(DISTC(1)+DISTC(2))*SIN(CITA)/COS(CITA), 0.0,-(DISTC(1)+DISTC(2))
N,
200+NUMC,NUMC-1+(DISTC(1)+DISTC(2))* SIN(CITA)/COS(CITA),0.0,-(DISTC(1)+DISTC(2)) FILL,201,200+NUMC
N, 400+1, DISTC(2)*SIN(CITA)/COS(CITA), 0.0, -DISTC(2)
N,
400+NUMC,NUMC-1+DISTC(2)*SIN(CITA)/ COS(CITA), 0.0,-DISTC(2)
FILL,401,400+NUMC
N, 600+1, -DISTC(2)*SIN(CITA)/COS(CITA), 0.0, DISTC(2)
N,
600+NUMC,NUMC-1-DISTC(4)*SIN(CITA)/C OS(CITA),0.0, DISTC(4)
FILL,601,600+NUMC
N, 800+1,
-(DISTC(4)+DISTC(5))*SIN(CITA)/COS(CITA ),0.0, DISTC(4)+DISTC(5) N,
800+NUMC,NUMC-1-(DISTC(4)+DISTC(5))* SIN(CITA)/COS(CITA),0.0, DISTC(4)+DISTC(5) FILL,801,800+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.3133,0.0871,0.0959,1.0,0.25,,
RMORE ,,0.00265,,,,,
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.6048,0.5486,1.5400,2.0,1.0,,
RMORE ,,0.036,,,,,
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,2.3625,0.9885,2.2600,2.0,1.0,,
RMORE ,,0.0956,,,,
! 定义横梁单元材料、几何参数
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.4582,0.0416,1.0,1.0,,
RMORE ,,0.00879,,,,,
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.0250,0.2770,0.0495,1.0,1.0,,
RMORE ,,0.0609,,,,,
! 生成纵向梁格单元
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,800+I,800+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,600+I,600+I+1 *ENDDO
TYPE,3
MAT,3
REAL,3
*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,1
E,1,600+1
E,600+1,800+1 TYPE,4
MAT,4
REAL,4
*DO,I,2,L(1),1
E,200+I,400+I
E,400+I,I
E,I,600+I
E,600+I,800+I *ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+1,400+L(1)+1
E,400+L(1)+1,L(1)+1
E,L(1)+1,600+L(1)+1
E,600+L(1)+1,800+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,I
E,I,600+I
E,600+I,800+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,L(1)+L(2)+1
E,L(1)+L(2)+1,600+L(1)+L(2)+1
E,600+L(1)+L(2)+1,800+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,I
E,I,600+I
E,600+I,800+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,L(1)+L(2)+L(3)+1
E,L(1)+L(2)+L(3)+1,600+L(1)+L(2)+L(3)+1
E,600+L(1)+L(2)+L(3)+1,800+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,I
E,I,600+I
E,600+I,800+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,L(1)+L(2)+L(3)+L(4) +1
E,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,800+L(1)+L(2)+L(3) +L(4)+1
/SOL
ACEL,,9.81,0
!******************* ********************* ! 定义边界条件
!**************************************** D,400+1,UX,0,,,,UY,UZ
D,600+1,UX,0,,,,UY,UZ
D,21,UY,0
D,421,UY,0
D,621,UY,0
D,51,UY,0
D,651,UY,0
D,451,UY,0
D,81,UY,0
D,681,UY,0
D,481,UY,0
D,400+101,UX,0,,,,UY,UZ
D,600+101,UX,0,,,,UY,UZ
ANTYPE,0
/STATUS,SOLU
SOLVE
/POST1
/REP,FAST
A VPRIN,0, ,
ETABLE, ,SMISC, 2
A VPRIN,0, ,
ETABLE, ,SMISC,8 A VPRIN,0, , ETABLE, ,SMISC, 4 A VPRIN,0, , ETABLE, ,SMISC, 5 A VPRIN,0, , ETABLE, ,SMISC, 6 A VPRIN,0, , ETABLE, ,SMISC, 10 A VPRIN,0, , ETABLE, ,SMISC, 11 A VPRIN,0, , ETABLE, ,SMISC, 12 FINISH。