ansys连续梁桥的影响线计算.docVIP

FINI /CLEAR SPAN=3 !定义连续梁跨数 L=20 !定义单跨跨径 /FILNAME,MODEL /TITLE,The %SPAN%*%L%M Continuous Bridge Analysis /NOPR /PREP7 *AFUN,DEG /VIEW,1,0,0,1 !!!! 初始化结束 !!!!******************************** !!!!******************************** !!!!定义相关参数 SEGMENT=5 !每米单元个数 NMAX=SEGMENT*L*SPAN+1 !计算最大节点号 X0=0.0 Y0=0.0 !定义原点 SEC_H=L/20 SEC_W=L/10 !计算截面高、宽 !!!!定义相关参数定义结束 !!!!******************************** !!!!******************************** !!!!建立有限元模型 ET,1,BEAM3 MP,EX,1,3.5E7 MP,PRXY,1,0.167 !定义单元和材料属性 R,1,SEC_W*SEC_H,SEC_H**3*SEC_W/12,SEC_H !定义梁截面 N,1,X0,Y0 N,NMAX,X0+L*SPAN,Y0 FILL,1,NMAX !定义节点 TYPE,1 MAT,1 REAL,1 ESYS,0 !定义单元属性 *DO,I,1,NMAX-1 E,I,I+1 *ENDDO !生成单元 *DO,J,1,NMAX,SEGMENT*L D,J,UY *IF,J,EQ,1,THEN D,J,UX *ENDIF *ENDDO !对节点施加约束 ALLS /PBC,ALL,,1 EPLOT SAVE !存储模型 !!!!模型建立结束 !!!!******************************** !!!!******************************** !!!!模型求解 /SOLU *DO,I,1,NMAX F,I,FY,-1 *IF,I,GT,1,THEN FDELE,I-1,FY *ENDIF !对节点施加单元荷载 LSWRITE,I !存储载荷步 *ENDDO ANTYPE,0 OUTRES,ALL,LAST, OUTPR,BASIC,ALL, TIME,1 AUTOTS,-1 NSUBST,1, , ,1 KBC,0 !定义有关求解参数 ALLS,ALL LSSOLVE,1,NMAX !求解所有载荷步 !!!!模型求解结束 !!!!******************************** !!!!******************************** !!!!进入后处理器计算影响线 /POST1 *DIM,NODE_X,ARRAY,NMAX !定义节点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,NMAX N_QY(0,I)=I N_MZ(0,I)=I *ENDDO !给N_QY和N_MZ表的0行赋值（节点号） *VGET,NODE_X,NODE,ALL,LOC,X !给NODE_X赋值（节点X坐标） *ABSET,Creating Efficacy Data ......,BAR !初始化状态条 *DO,I,1,NMAX *ABCHECK,NINT(100*I/NMAX) !更新状态条 SET,I ETABLE,QYI,SMISC,2 ETABLE,QYJ,SMISC,8 ETABLE,MZI,SMISC,6 ETABLE,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坐 标） *VGET,TRANS_V,ELEM,ALL,ETAB,QYI !从读取各单元I端剪力QYI赋予TRANS_V *MFUN,N_QY(I,1),TRAN,TRANS_V !将TRANS