曲线梁桥ANSYS计算命令流

!无偏载  (以跨径布置30m+40m+30m,桥宽8.5为例）

!上海城市设计研究院L1+L2+L3预应力混凝土曲线连续梁桥结构分析

!两端为抗扭支座，中间支座为点铰支座

!每次要记得修改横隔梁的参数，即Mass21单元的实常数

!*****************************************************************************************************

FINI

/CLE

/prep7

!DEFINE THE ELEMENTARY PARAMETERS

*DIM,L,ARRAY,10

*DIM,H,ARRAY,10

*DIM,CITA,ARRAY,10

!*****以下参数均可修改***************

N跨数

L第一跨

L第二跨

L第三跨

e墩处内支座到中心线的间距

e墩处外支座到中心线的间距

e墩处的支座偏心距（正的表示外偏）

e墩处的支座偏心距

e墩处内支座到中心线的间距

e墩处外支座到中心线的间距

R曲线桥半径

H梁底到截面形心处的高度

M单元质量

J单元转动惯量

!*************************************

LL=0.0

*DO,I,1,N

 LL=LL+L(I)

 CITA(I)=L(I)/R/3.1415925*180

*ENDDO

CITA0=LL/R/3.1415925*180

LOCAL,11,1,0,0,0,0,270,0

CSYS,11

!DEFINE THE NODES OF BRIDGE

N,     

N,    

N,      L(1)+L(2)+1,  R,      -179.9999+CITA(1)+CITA(2),         H0

N,      L(1)+L(2)+L(3)+1, R,  -179.9999+CITA0,                    H0

FILL,1,L(1)+1

FILL,L(1)+1,L(1)+L(2)+1

FILL,L(1)+L(2)+1,L(1)+L(2)+L(3)+1

!DEFINE THE REFERENCE POINT OF MAIN DECK

N,2000,0,0,0

!*****************************************************************

!DEFINE THE NODES OF RIGID BEAM OF MAIN GIRDER

N,  2001,   

N,  2002,   

N,  2003,   

N,  2004,  

N,  2005,   

N,  2006,   

NROTAT,ALL

!***************************************************** 

!定义主梁单元材料、几何参数 

!混凝土标号为C40

!**************梁段截面几何特性*******************

ET,1,BEAM4

MP,EX     

MP,NUXY,    

MP,DENS, 

MP,ALPX, 

R跨中梁断面

RMORE  ,,5.0388,,,,,

!**************刚臂截面几何特性*******************

!增设抗扭或偏心支座时采用

ET,3,BEAM4

MP,EX     

MP,NUXY,    

MP,DENS,    

MP,ALPX, 

R,3,100,1000,1000,0.5,0.5,,

RMORE  ,,1000,,,,,

!****************生成主梁单元**********************

TYPE,1

MAT,1

REAL,1

*DO,I,1,LL,1

 E,I,I+1,2000

*ENDDO

!***********************************************************

!生成主梁刚臂单元

!***********************************************************

TYPE,3

MAT,3

REAL,3

E,2001,1

E,1,2002

*IF,e3,GT,0.0,THEN

E,2003,L(1)+1

*ELSE

NSEL,U, , ,2003 

*ENDIF

*IF,e4,GT,0.0,THEN

E,2004,L(1)+L(2)+1

*ELSE

NSEL,U, , ,2004 

*ENDIF

E,2005,LL+1

E,LL+1,2006

!**************************************************

!支墩处的横隔梁采用Mass21单元来模拟

!**************************************************

ET,4,MASS21

TYPE,4

REAL,100

E,1

E,L(1)+1

E,L(1)+L(2)+1

E,LL+1 

!****施加约束条件**********

D,2001,UX,,,,,UY,UZ

D,2002,UZ

*IF,e3,GT,0.0,THEN

D,2003,,,,,,,UZ

*ELSE

D,L(1)+1,,,,,,,UZ

*ENDIF

*IF,e4,GT,0.0,THEN

D,2004,,,,,,,UZ

*ELSE

D,L(1)+L(2)+1,,,,,,,UZ

*ENDIF

D,2005,UX,,,,,,UZ

D,2006,UZ

!**********************************************************

!对桥梁桥面各节点逐点加载

!**********************************************************

A静力分析

ACEL,0,0,0

K=1

*DO,I,1,LL+1,1

  TIME,K

  FDELE,ALL,ALL

  F,I,FZ,-1000

  SOLVE

  K=K+1

 *ENDDO

!******定义各片纵梁的竖向挠度、弯矩、扭矩、剪力的包络图********

*DIM,AMAX,ARRAY,LL+1

*DIM,AMIN,ARRAY,LL+1

*DIM,JZMAX,ARRAY,LL+1

*DIM,JZMIN,ARRAY,LL+1

*DIM,DISP1,ARRAY,LL+1,LL+1

*DIM,MZ1,ARRAY,LL,LL+1

*DIM,MX1,ARRAY,LL,LL+1

*DIM,QY1,ARRAY,LL,LL+1

*DIM,DISP1MAX,ARRAY,LL+1

*DIM,DISP1MIN,ARRAY,LL+1

*DIM,MZ1MAX,ARRAY,LL

*DIM,MZ1MIN,ARRAY,LL

*DIM,MX1MAX,ARRAY,LL

*DIM,MX1MIN,ARRAY,LL

*DIM,QY1MAX,ARRAY,LL

*DIM,QY1MIN,ARRAY,LL

*DIM,KK,ARRAY,LL+1

*DIM,REACTION,ARRAY,6,LL+1

*DIM,REACTIONMAX,ARRAY,6

*DIM,REACTIONMIN,ARRAY,6

*DO,I,1,LL+1,1

 /POST1

 SET,I

 !*****************竖向位移**************

 *DO,J,1,LL+1,1

 *ENDDO

 !*****************竖向弯矩Mz***********

 *DO,J,1,LL,1

   *GET,MZ1(J,I),ELEM,   J, SMISC,6

 *ENDDO

!*****************扭矩Mx***********

 *DO,J,1,LL,1

   *GET,MX1(J,I),ELEM,   J, SMISC,4

 *ENDDO

 !*****************剪力Qy***********

 *DO,J,1,LL,1

   *GET,QY1(J,I),ELEM,   J, SMISC,2

 *ENDDO

 !*****************支座反力***********

*IF,e3,GT,0.0,THEN

*ELSE

*ENDIF

*IF,e4,GT,0.0,THEN

*ELSE

*ENDIF

KK(I)=I

*ENDDO

!**************************公路I级汽车活载影响线加载效应计算*******************************

!车道宽度假定为3.75m

!车道布置：横向NN=2车道，不考虑横向折减系数;跨径小于150m，不需要考虑纵向折减。

!荷载标准：PK=180+(40-5)/(50-5)*180=320kN,均布荷载标准为10.5kN/m。

!当计算剪力效应时，集中荷载标准值应乘以1.2的系数。暂不考虑冲击系数。

!*******************************************************************************************

NN=2

PK=180+(L(2)-5)/(50-5)*180

QK=10.5

!*****************主梁位移包络图*******************

*DO,I,1,LL+1,1

  AMAX(I)=0.0

  AMIN(I)=0.0

  JZMAX(I)=0.0

  JZMIN(I)=0.0

*DO,J,1,LL+1,1

  *ENDDO

*DO,J,1,LL+1,1

  *ENDDO

*ENDDO 

!*****************主梁弯矩包络图*******************

*DO,I,1,LL,1

  AMAX(I)=0.0

  AMIN(I)=0.0

  JZMAX(I)=0.0

  JZMIN(I)=0.0

*DO,J,1,LL+1,1

  *ENDDO

*DO,J,1,LL+1,1

  *ENDDO

!*****************主梁扭矩包络图*******************

*DO,I,1,LL,1

  AMAX(I)=0.0

  AMIN(I)=0.0

  JZMAX(I)=0.0

  JZMIN(I)=0.0

*DO,J,1,LL+1,1

  *ENDDO

*DO,J,1,LL+1,1

  *ENDDO

!*****************主梁剪力包络图*******************

*DO,I,1,LL,1

  AMAX(I)=0.0

  AMIN(I)=0.0

  JZMAX(I)=0.0

  JZMIN(I)=0.0

*DO,J,1,LL+1,1

  *ENDDO

*DO,J,1,LL+1,1

  *ENDDO

*ENDDO  

!*****************主梁支座反力包络值*******************

*DO,I,1,6,1

  AMAX(I)=0.0

  AMIN(I)=0.0

  JZMAX(I)=0.0

  JZMIN(I)=0.0

*DO,J,1,LL+1,1

  *ENDDO

*DO,J,1,LL+1,1

  *ENDDO

*ENDDO  

LABEL(1,2)='DISPMAX(M)'

LABEL(1,3)='DISPMIN(M)'

LABEL(1,4)='MZMAX(N.M)'

LABEL(1,5)='MZMIN(N.M)'

LABEL(1,6)='MXMAX(N.M)'

LABEL(1,7)='MXMIN(N.M)'

LABEL(1,8)='QYMAX(N)'

LABEL(1,9)='QYMIN(N)'

LABEL(1,10)='REACTIONMAX(N)'

LABEL(1,11)='REACTIONMIN(N)'

!********************************************************************

!

!将各片主梁位移及内力包络值输出到文件中

!********************************************************************

*CFOPEN,DISPMAXMIN,dat 

*VWRITE,LABEL(1,1),LABEL(1,2),LABEL(1,3)

(2X,A5,1X,A10,1X,A10)

(1X,F5.0,2X,E10.4,2X,E10.4)  

*CFCLOS 

*CFOPEN,MZMAXMIN,dat 

*VWRITE,LABEL(1,1),LABEL(1,4),LABEL(1,5)

(2X,A5,1X,A10,1X,A10)

(1X,F5.0,2X,E10.4,2X,E10.4)  

*CFCLOS 

*CFOPEN,MXMAXMIN,dat 

*VWRITE,LABEL(1,1),LABEL(1,6),LABEL(1,7)

(2X,A5,1X,A10,1X,A10)

(1X,F5.0,2X,E10.4,2X,E10.4)  

*CFCLOS 

*CFOPEN,QYMAXMIN,dat 

*VWRITE,LABEL(1,1),LABEL(1,8),LABEL(1,9)

(2X,A5,1X,A10,1X,A10)

(1X,F5.0,2X,E10.4,2X,E10.4)  

*CFCLOS 

*CFOPEN,REACTIONMAXMIN,dat 

*VWRITE,LABEL(1,1),LABEL(1,10),LABEL(1,11)

(2X,A5,1X,A10,1X,A10)

(1X,F5.0,2X,E10.4,2X,E10.4)  

*CFCLOS