!****************************************************************************** ***********************! case2:无偏载(以跨径布置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=3 !跨数L(1)=30 !第一跨L(2)=40 !第二跨L(3)=30 !第三跨e1=1.25 !1#墩处内支座到中心线的间距e2=1.25 !1#墩处外支座到中心线的间距e3=0 !2#墩处的支座偏心距（正的表示外偏）e4=0 !3#墩处的支座偏心距e5=1.25 !4#墩处内支座到中心线的间距e6=1.25 !4#墩处外支座到中心线的间距R=10000 !曲线桥半径H0=1.0 !梁底到截面形心处的高度M=16146 !mass21单元质量J=27246.38 !mass21单元转动惯量!*************************************LL=0.0*DO,I,1,NLL=LL+L(I)CITA(I)=L(I)/R/3.1415925*180*ENDDOCITA0=LL/R/3.1415925*180LOCAL,11,1,0,0,0,0,270,0CSYS,11!DEFINE THE NODES OF BRIDGEN, 1, R, -179.9999, H0N, L(1)+1, R, -179.9999+CITA(1), H0N, 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, H0FILL,1,L(1)+1FILL,L(1)+1,L(1)+L(2)+1FILL,L(1)+L(2)+1,L(1)+L(2)+L(3)+1!DEFINE THE REFERENCE POINT OF MAIN DECKN,2000,0,0,0!***************************************************************** !DEFINE THE NODES OF RIGID BEAM OF MAIN GIRDERN, 2001, R-e1, -179.9999, H0N, 2002, R+e2, -179.9999, H0N, 2003, R+e3, -179.9999+CITA(1), H0N, 2004, R+e4, -179.9999+CITA(1)+CITA(2), H0N, 2005, R-e5, -179.9999+CITA0, H0N, 2006, R+e6, -179.9999+CITA0, H0NROTA T,ALL!*****************************************************! 定义主梁单元材料、几何参数! 混凝土标号为C40!**************梁段截面几何特性*******************ET,1,BEAM4MP,EX , 1, 3.25E+10MP,NUXY, 1, 0.2000MP,DENS, 1, 3950.50MP,ALPX, 1, 1.00E-05R,1,4.39129,2.4428,20.7089,8.5,2.0,, !跨中梁断面RMORE ,,5.0388,,,,,!**************刚臂截面几何特性*******************! 增设抗扭或偏心支座时采用ET,3,BEAM4MP,EX , 3, 3.25E+10MP,NUXY, 3, 0.2000MP,DENS, 3, 0.0MP,ALPX, 3, 1.0E-08R,3,100,1000,1000,0.5,0.5,,RMORE ,,1000,,,,,!****************生成主梁单元**********************TYPE,1MA T,1REAL,1*DO,I,1,LL,1E,I,I+1,2000*ENDDO!*********************************************************** ! 生成主梁刚臂单元!*********************************************************** TYPE,3MA T,3REAL,3E,2001,1E,1,2002*IF,e3,GT,0.0,THENE,2003,L(1)+1*ELSENSEL,U, , ,2003*ENDIF*IF,e4,GT,0.0,THENE,2004,L(1)+L(2)+1*ELSENSEL,U, , ,2004*ENDIFE,2005,LL+1E,LL+1,2006!**************************************************! 支墩处的横隔梁采用Mass21单元来模拟!**************************************************ET,4,MASS21R,100, M ,M ,M ,J ,,,TYPE,4REAL,100E,1E,L(1)+1E,L(1)+L(2)+1E,LL+1!****施加约束条件**********D,2001,UX,,,,,UY,UZD,2002,UZ*IF,e3,GT,0.0,THEND,2003,,,,,,,UZ*ELSED,L(1)+1,,,,,,,UZ*ENDIF*IF,e4,GT,0.0,THEND,2004,,,,,,,UZ*ELSED,L(1)+L(2)+1,,,,,,,UZ*ENDIFD,2005,UX,,,,,,UZD,2006,UZ!********************************************************** ! 对桥梁桥面各节点逐点加载!********************************************************** /SOLUANTYPE,0 ! 静力分析ACEL,0,0,0K=1*DO,I,1,LL+1,1TIME,KFDELE,ALL,ALLF,I,FZ,-1000SOLVEK=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/POST1SET,I!*****************竖向位移***************DO,J,1,LL+1,1*GET,DISP1(J,I),NODE,J,U,Y*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!*****************支座反力************GET,REACTION(1,I),NODE, 2001, RF,FZ*GET,REACTION(2,I),NODE, 2002, RF,FZ*IF,e3,GT,0.0,THEN*GET,REACTION(3,I),NODE, 2003, RF,FZ*ELSE*GET,REACTION(3,I),NODE, L(1)+1, RF,FZ*ENDIF*IF,e4,GT,0.0,THEN*GET,REACTION(4,I),NODE, 2004, RF,FZ*ELSE*GET,REACTION(4,I),NODE, L(1)+L(2)+1, RF,FZ *ENDIF*GET,REACTION(5,I),NODE, 2005, RF,FZ*GET,REACTION(6,I),NODE, 2006, RF,FZKK(I)=I*ENDDO!**************************公路I级汽车活载影响线加载效应计算*******************************! 车道宽度假定为3.75m! 车道布置：横向NN=2车道，不考虑横向折减系数;跨径小于150m，不需要考虑纵向折减。

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

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