资源描述：
《ansys连续梁桥的影响线计算》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、FINI/CLEARSPAN=3!定义连续梁跨数L=20!定义单跨跨径/FILNAME,MODEL/TITLE,The%SPAN%*%L%MContinuousBridgeAnalysis/NOPR/PREP7*AFUN,DEG/VIEW,1,0,0,1!!!!初始化结束!!!!********************************!!!!********************************!!!!定义相关参数SEGMENT=5!每米单元个数NMAX=SEGMENT*L*SPAN+1!计算最大节点号X0=0.0Y0=0.0!定义原点
2、SEC_H=L/20SEC_W=L/10!计算截面高、宽!!!!定义相关参数定义结束!!!!********************************!!!!********************************!!!!建立有限元模型ET,1,BEAM3MP,EX,1,3.5E7MP,PRXY,1,0.167!定义单元和材料属性R,1,SEC_W*SEC_H,SEC_H**3*SEC_W/12,SEC_H!定义梁截面N,1,X0,Y0N,NMAX,X0+L*SPAN,Y0FILL,1,NMAX!定义节点TYPE,1MAT,1REAL,1E
3、SYS,0!定义单元属性*DO,I,1,NMAX-1E,I,I+1*ENDDO!生成单元*DO,J,1,NMAX,SEGMENT*LD,J,UY*IF,J,EQ,1,THEND,J,UX*ENDIF*ENDDO!对节点施加约束ALLS/PBC,ALL,,1EPLOTSAVE!存储模型!!!!模型建立结束!!!!********************************!!!!********************************!!!!模型求解/SOLU*DO,I,1,NMAXF,I,FY,-1*IF,I,GT,1,THENFDELE,I
4、-1,FY*ENDIF!对节点施加单元荷载LSWRITE,I!存储载荷步*ENDDOANTYPE,0OUTRES,ALL,LAST,OUTPR,BASIC,ALL,TIME,1AUTOTS,-1NSUBST,1,,,1KBC,0!定义有关求解参数ALLS,ALLLSSOLVE,1,NMAX!求解所有载荷步!!!!模型求解结束!!!!********************************!!!!********************************!!!!进入后处理器计算影响线/POST1*DIM,NODE_X,ARRAY,NMAX!定
5、义节点X坐标存储数组*DIM,N_QY,TABLE,NMAX,NMAX!定义QY影响线存储表*DIM,N_MZ,TABLE,NMAX,NMAX!定义MZ影响线存储表*DIM,QY_EXTR,ARRAY,4,NMAX!定义QY影响线极值存储数组*DIM,MZ_EXTR,ARRAY,4,NMAX!定义MZ影响线极值存储数组*DIM,TRANS_V,ARRAY,NMAX-1!定义中间传递数据向量ESEL,ALL*DO,I,1,NMAXN_QY(0,I)=IN_MZ(0,I)=I*ENDDO!给N_QY和N_MZ表的0行赋值（节点号）*VGET,NODE_X,N
6、ODE,ALL,LOC,X!给NODE_X赋值（节点X坐标）*ABSET,'CreatingEfficacyData......',BAR!初始化状态条*DO,I,1,NMAX*ABCHECK,NINT(100*I/NMAX)!更新状态条SET,IETABLE,QYI,SMISC,2ETABLE,QYJ,SMISC,8ETABLE,MZI,SMISC,6ETABLE,MZJ,SMISC,12!建立Beam3的单元表数据N_QY(I,0)=NODE_X(I,1)N_MZ(I,0)=NODE_X(I,1)!给N_QY和N_MZ的0列赋值（节点X坐标）*VGE
7、T,TRANS_V,ELEM,ALL,ETAB,QYI!从读取各单元I端剪力QYI赋予TRANS_V*MFUN,N_QY(I,1),TRAN,TRANS_V!将TRANS_V中的值赋予N_QY的第一列*GET,N_QY(I,NMAX),ELEM,NMAX-1,SMISC,8!读取最后一个单元的J端剪力QY*VGET,TRANS_V,ELEM,ALL,ETAB,MZI!从读取各单元I端弯矩MZI赋予TRANS_V*MFUN,N_MZ(I,1),TRAN,TRANS_V!将TRANS_V中的值赋予N_MZ的第一列*GET,N_MZ(I,NMAX),ELEM,
8、NMAX-1,SMISC,12!读取最后一个单元的J端弯矩MZ*ENDDO*AB
