Ansys桥梁计算汇编

ANSYS实例分析75道(含结果)【【ANSYS算例算例】】3.4.2(1)基于图形界面的桁架桥梁结构分析基于图形界面的桁架桥梁结构分析(stepbystep)下面以一个简单桁架桥梁为例，以展示有限元分析的全过程。

背景素材选自位于密执安的“OldNorthParkBridge“(1904-1988)，见图3-22。

该桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表3-6。

桥长L=32m,桥高H=5.5m。

桥身由8段桁架组成，每段长4m。

该桥梁可以通行卡车，若这里仅考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg，若取一半的模型，可以将卡车对桥梁的作用力简化为P1，P2和P3，其中P1=P3=5000N,P2=10000N，见图3-23。

图3-22位于密执安的“OldNorthParkBridge“(1904-1988)图3-23桥梁的简化平面模型(取桥梁的一半)表3-6桥梁结构中各种构件的几何性能参数构件惯性矩m4横截面积m2顶梁及侧梁桥身弦梁底梁解答解答以下为基于ANSYS 图形界面(GraphicUserInterface,GUI)的菜单操作流程。

(1)进入进入ANSYS（设定工作目录和工作文件）（设定工作目录和工作文件）程序程序→→ANSYS→→ANSYSInteractive→→Workingdirectory（设置工作目录）→Initialjobname（设置工作文件名）:TrussBridge→→Run→→OK(2)设置计算类型设置计算类型：Preferences…→→Structural→→OK(3)定义单元类型定义单元类型ANSYSMainMenu：Preprocessor→→ElementType→→Add/Edit/Delete.→→Add…→→Beam:2delastic3→→OK（返回到ElementTypes窗口）→→Close(4)定义实常数以确定梁单元的截面参数定义实常数以确定梁单元的截面参数ANSYSMainMenu:Preprocessor→→RealConstants…→→Add/Edit /Delete→→Add…→→selectType1Beam3→→OK→→RealConsta ntsSetNo.:1,AREA:2.19E-3，，Izz:3.83e-6(1号实常数用于顶梁和侧梁)→→Apply→→RealConstantsSetNo.:2,AREA:1.185E-3，，Izz:1.87E-6(2号实常数用于弦杆)→→Apply→→RealConstantsSetNo.:3,AREA:3.031E-3，，Izz:8.47E-6(3号实常数用于底梁)→→OK(backtoRealConstantswindow)→Close(theRealConstant swindow)(5)定义材料参数定义材料参数ANSYSMainMenu:Preprocessor→→MaterialProps→→MaterialMo dels→→Structural→→Linear→→Elastic→→Isotropic→→EX:2.1e11,PRXY:0.3(定义泊松比及弹性模量)→→OK→→Density(定义材料密度)→DENS:7800,→→OK→→Close（关闭材料定义窗口）(6)构造桁架桥模型构造桁架桥模型生成桥体几何模型ANSYSMainMenu：Preprocessor→→Modeling→→Create→→Keypoints→→InActive CS→→NPTKeypointnumber：：1，，X，，Y，，ZLocationinactiveCS：：0，，0→→Apply→→同样输入其余15个特征点坐标（最左端为起始点，坐标分别为(4,0),(8,0),(12,0),(16,0),(20,0),(24,0),(28,0),(32,0),(4,5.5),(8,5.5),(12 ,5.5),(16.5.5),(20,5.5),(24,5.5),(28,5.5)）→Lines→Lines→→StraightLine→→依次分别连接特征点→→OK网格划分ANSYSMainMenu:Preprocessor→→Meshing→→MeshAttributes→→PickedLines→→选择桥顶梁及侧梁→→OK→→selectREAL:1,TYPE:1→→Apply→→选择桥体弦杆→→OK→→selectREAL:2,TYPE:1→→Apply→→选择桥底梁→→OK→→selectREAL:3,TYPE:1→→OK→→ANSYSMainMen u：Preprocessor→→Meshing→→MeshTool→→位于SizeControls 下的Lines：：Set→→ElementSizeonPicked→→Pickall→→Apply→→NDIV：：1→→OK→→Mesh→→Lines→→Pickall→→OK(划分网格)(7)模型加约束模型加约束ANSYSMainMenu:Solution→→DefineLoads→→Apply→→Struct ural→→Displacement→→OnNodes→→选取桥身左端节点→→OK→→selectLab2:AllDOF(施加全部约束)→→Apply→→选取桥身右端节点→→OK→→selectLab2:UY(施加Y方向约束)→→OK(8)施加载荷施加载荷ANSYSMainMenu:Solution→→DefineLoads→→Apply→→Struct ural→→Force/Moment→→OnKeypoints→→选取底梁上卡车两侧关键点（X坐标为12及20）→→OK→→selectLab:FY，，Value:-5000→→Apply→→选取底梁上卡车中部关键点（X坐标为16）→→OK→→selectLab:FY，，Value:-10000→→OK→→ANSYSUtilityMenu：→→Select→→Everything(9)计算分析计算分析ANSYSMainMenu：Solution→→Solve→→CurrentLS→→OK(10)结果显示结果显示ANSYSMainMenu：GeneralPostproc→→PlotResults→→Deedshape→→Defshapeonly →→OK（返回到PlotResults）→→ContourPlot→→NodalSolu→→DOFSolution,Y-Componentof Displacement→→OK（显示Y方向位移UY）(见图3-24(a))定义线性单元I节点的轴力ANSYSMainMenu→GeneralPostproc→→ElementTable→→Define Table→→Add→→Lab:[bar_I],Bysequencenum:[SMISC,1]→→OK →→Close定义线性单元J节点的轴力ANSYSMainMenu→→GeneralPostproc→→ElementTable→→Def ineTable→→Add→→Lab:[bar_J],Bysequencenum:[SMISC,1]→→OK→→Close画出线性单元的受力图(见图3-24(b))ANSYSMainMenu→→GeneralPostproc→→PlotResults→→ContourPlot→→LineElemRes→→LabI:[bar_I],LabJ:[bar_J],Fact :[1]→→OK(11)退出系统退出系统ANSYSUtilityMenu：File→→Exit→→SaveEverything→→OK(a)桥梁中部最大挠度值为0.003374m(b)桥梁中部轴力最大值为25380N图3.24桁架桥挠度UY以及单元轴力计算结果【【ANSYS算例算例】】3.4.2(2)基于命令流方式的桁架桥梁结构分析基于命令流方式的桁架桥梁结构分析!%%%%%[ANSYS 算例]3.4.2(2)%%%%%begin%%%%%%!------注：命令流中的符号$，可将多行命令流写成一行------/prep7!进入前处理/PLOPTS,DATE,0!设置不显示日期和时间!=====设置单元和材料ET,1,BEAM3!定义单元类型R,1,2.19E-3,3.83e-6,,,,,!定义1号实常数用于顶梁侧梁R,2,1.185E-3,1.87e-6,0,0,0,0,!定义2号实常数用于弦杆R,3,3.031E-3,8.47E-6,0,0,0,0,!定义3号实常数用于底梁MP,EX,1,2.1E11!定义材料弹性模量MP,PRXY,1,0.30!定义材料泊松比MP,DENS,1,,7800!定义材料密度!-----定义几何关键点K,1,0,0,,$K,2,4,0,,$K,3,8,0,,$K,4,12,0,,$K,5,16,0,,$K,6,20,0,,$K,7,2 4,0,,$K,8,28,0,,$K,9,32,0,,$K,10,4,5.5,,$K,11,8,5.5,,$K,12,12,5.5,,$K,13,16,5.5,,$K,14,20,5.5,,$K,15,24,5.5,,$K,16,28,5.5,,!-----通过几何点生成桥底梁的线L,1,2$L,2,3$L,3,4$L,4,5$L,5,6$L,6,7$L,7,8$L,8,9!------生成桥顶梁和侧梁的线L,9,16$L,15,16$L,14,15$L,13,14$L,12,13$L,11,12$L,10,11$L,1,10! ------生成桥身弦杆的线L,2,10$L,3,10$L,3,11$L,4,11$L,4,12$L,4,13$L,5,13$L,6,13$L,6,14 $L,6,15$L,7,15$L,7,16$L,8,16!------选择桥顶梁和侧梁指定单元属性LSEL,S,,,9,16,1,LATT,1,1,1,,,,!-----选择桥身弦杆指定单元属性LSEL,S,,,17,29,1,LATT,1,2,1,,,,!-----选择桥底梁指定单元属性LSEL,S,,,1,8,1,LATT,1,3,1,,,,!------划分网格AllSEL,all!再恢复选择所有对象LESIZE,all,,,1,,,,,1!对所有对象进行单元划分前的分段设置LMESH,all!对所有几何线进行单元划分!=====在求解模块中，施加位移约束、外力，进行求解/soluNSEL,S,LOC,X,0!根据几何位置选择节点D,all,,,,,,ALL,,,,,!对所选择的节点施加位移约束AllSEL,all!再恢复选择所有对象NSEL,S,LOC,X,32!根据几何位置选择节点D,all,,,,,,,UY,,,,!对所选择的节点施加位移约束ALLSEL,all!再恢复选择所有对象!------基于几何关键点施加载荷FK,4,FY,-5000$FK,6,FY,-5000$FK,5,FY,-10000/replot!重画图形Allsel,all!选择所有信息(包括所有节点、单元和载荷等)solve!求解!=====进入一般的后处理模块/post1!后处理PLNSOL,U,Y,0,1.0!显示Y方向位移PLNSOL,U,X,0,1.0!显示X方向位移!------显示线单元轴力------ETABLE,bar_I,SMISC,1ETABLE,bar_J,SMISC,1PLLS,BAR_ I,BAR_J,0.5,1!画出轴力图finish!结束!%%%%%[ANSYS算例]3.4.2(2)%%%%%end%%%%%%【【ANSYS算例算例】】3.2.5(3)四杆桁架结构的有限元分析四杆桁架结构的有限元分析下面针对【典型例题】3.2.5(1)的问题，在ANSYS平台上，完成相应的力学分析。

钢桁架桥静力受力分析对一架钢桁架桥进行具体静力受力分析，分别采用GUI方式和命令流方式。

A 问题描述图6-15 钢桁架桥简图已知下承式简支钢桁架桥桥长72米，每个节段12米，桥宽10米，高16米。

设桥面板为0.3米厚的混凝土板，当车辆行驶于桥梁上面时，轴重简化为一组集中力作用于梁上，来计算梁的受力情况。

桁架杆件规格有三种，见下表：杆件截面号形状规格端斜杆 1 工字形400×400×16×16上下弦 2 工字形400×400×12×12 横向连接梁 2 工字形400×400×12×12其他腹杆 3 工字形400×300×12×12所用材料属性如下表：表6-3 材料属性参数钢材混凝土弹性模量EX 2.1×1011 3.5×1010泊松比PRXY 0.3 0.1667密度DENS 7850 2500B GUI操作方法1.创建物理环境1）过滤图形界面：GUI：Main Menu> Preferences，弹出“Preferences for GUI Filtering”对话框，选中“Structural”来对后面的分析进行菜单及相应的图形界面过滤。

2）定义工作标题：GUI：Utility Menu> File> Change Title，在弹出的对话框中输入“Truss Bridge StaticAnalysis”，单击“OK”。

如图6-16（a）。

指定工作名：GUI：Utility Menu> File> Change Jobname，弹出一个对话框，在“Enter new Name”后面输入“Structural”，“New log and error files”选择yes，单击“OK”。

如图6-16（b）。

图6-16（a）定义工作标题图6-16（b）指定工作名3）定义单元类型和选项：GUI：Main Menu> Preprocessor> Element Type> Add/Edit/Delete，弹出“Element Types”单元类型对话框，单击“Add”按钮，弹出“Library of Element Types”单元类型库对话框。

基于ansys的钢桁架桥的分析和计算姓名: 马彦学院：建筑与环境专业：工程力学学号：1043055033指导老师：朱哲明2013/6/151.问题简述钢桁架桥简图如下，尺寸如图，单元长12m，高16m。

设桥面板为0.3m厚的混凝土板。

杆件截面号形状规格端斜杆 1 工字梁400*400*16*16上下弦 2 工字梁400*400*12*12横向连接梁 2 工字梁400*400*12*12其他腹杆 3 工字梁400*300*12*12参数钢材混凝土EX 2.1x1011 3.5x1010PRXY 0.3 0.1667DENS 7850 25002.材料实常数3.半横架桥模型镜面对称，生成整体模型3.施加约束及受力4.计算及分析结果◆整体位移云图◆结点总位移矢量图◆单元第一主应力云图◆单元第二主应力云图◆单元第三主应力云图◆节点位移结果PRINT U NODAL SOLUTION PER NODE***** POST1 NODAL DEGREE OF FREEDOM LISTING *****LOAD STEP= 1 SUBSTEP= 1 TIME= 1.0000 LOAD CASE= 0THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATESYSTEMNODE UX UY UZ USUM1 0.18808E-02-0.20919E-01 0.70316E-03 0.21015E-012 0.11411E-02-0.21354E-01 0.59772E-03 0.21393E-013 0.14813E-02-0.20809E-01 0.11202E-02 0.20892E-014 0.15919E-02-0.20373E-01 0.11392E-02 0.20467E-015 0.22549E-02-0.18918E-01 0.10528E-02 0.19081E-016 0.23458E-02-0.18310E-01 0.10055E-02 0.18487E-017 -0.10050E-02-0.18459E-01-0.38731E-02 0.18887E-018 -0.11376E-02-0.19066E-01-0.38598E-02 0.19486E-019 0.24977E-02-0.12074E-01 0.72603E-03 0.12351E-0110 0.29237E-02-0.11079E-01 0.68719E-03 0.11479E-0111 -0.35033E-02-0.10438E-01-0.84626E-02 0.13887E-0112 -0.38537E-02-0.10965E-01-0.84226E-02 0.14353E-0113 0.27521E-02 0.0000 0.0000 0.27521E-0214 0.34768E-02 0.0000 0.0000 0.34768E-0215 0.82671E-03-0.17947E-01 0.14911E-03 0.17967E-0116 0.67748E-03-0.19250E-01 0.10648E-03 0.19262E-0117 0.42077E-02-0.19398E-01 0.59595E-02 0.20725E-0118 0.40812E-02-0.18095E-01 0.59727E-02 0.19488E-0119 0.40101E-03-0.10784E-01 0.34385E-04 0.10791E-0120 0.34470E-03-0.12307E-01 0.25523E-06 0.12312E-0121 0.69212E-02-0.11199E-01 0.10204E-01 0.16656E-0122 0.65820E-02-0.10142E-01 0.10244E-01 0.15847E-0123 0.0000 0.0000 0.0000 0.000024 0.0000 0.0000 0.0000 0.0000MAXIMUM ABSOLUTE VALUESNODE 21 2 22 2VALUE 0.69212E-02-0.21354E-01 0.10244E-01 0.21393E-01◆单元受力结果PRINT ELEMENT TABLE ITEMS PER ELEMENT***** POST1 ELEMENT TABLE LISTING *****STAT CURRENT CURRENTELEM ZHOU_I ZHOU_J1 -49659. 7936.32 -42695. -3502.73 -9873.9 -28642.4 9567.9 -51440.5 -15016. 23374.6 -22120. -5510.47 -26981. -11385.8 -33355. 18549.9 -17656. -15556.10 -16095. -16301.11 -16203. -16943.12 -12683. -20132.13 4836.6 5157.114 -17901. -18351.15 -2331.6 23001.16 -18331. -20015.17 -6067.9 50464.18 -19568. -26493.19 -5052.8 51411.20 -26836. -34142.21 -23626. -29919.22 -32522. -21349.23 -35649. -25215.24 -699.47 1061.525 690.13 -1048.326 5802.4 -1462.327 -9677.8 5182.928 16212. -4765.129 -4310.8 3979.130 -25.038 0.000031 -9.3064 0.000032 23.898 0.000033 -3569.2 -42609.34 8110.9 -49823.35 -5544.6 -22051.36 -11343. -27005.37 18453. -33238.38 -28592. -9977.139 -51593. 9648.540 23614. -15193.41 -16998. -16116.***** POST1 ELEMENT TABLE LISTING *****STAT CURRENT CURRENTELEM ZHOU_I ZHOU_J42 -20120. -12682.43 -15489. -17761.44 -16350. -16082.45 5157.1 4836.646 -18351. -17901.47 -2225.2 22850.48 -18463. -19869.49 -6087.5 50530.50 -19228. -26843.51 -5332.4 51796.52 -21374. -32473.53 -25205. -35655.54 -34114. -26894.55 -29953. -23607.56 -1061.5 699.4757 1048.3 -690.1358 5171.8 -9672.159 -1448.6 5796.560 3928.8 -4269.361 -4732.8 16215.62 -20.844 0.000063 -5.2944 0.000064 36.585 0.0000MINIMUM VALUESELEM 39 4VALUE -51593. -51440.MAXIMUM VALUESELEM 40 51VALUE 23614. 51796.5.命令流文件/FILNAM,Structural/TITLE,Truss Bridge Static Analysis/COM,Structural/prep7et,1,beam4et,2,shell63sectype,1,beam,i,,0 !定义工字型截面secoffset,cent !截面至心不偏移secdata,0.4,0.4,0.4,0.016,0.016,0.016,0,0,0,0 !定义工字型截面参数sectype,2,beam,i,,0secoffset,centsecdata,0.4,0.4,0.4,0.012,0.012,0.012,0,0,0,0sectype,3,beam,i,,0secoffset,centsecdata,0.3,0.3,0.4,0.012,0.012,0.012,0,0,0,0r,1,0.0187,0.00017,0.00054,0.4,0.4,0, !定义单元实常数r,2,0.0141,0.128e-3,0.415e-3,0.4,0.4,,r,3,0.0117,0.541e-4,0.324e-3,0.3,0.4,,r,4,0.3,,,,,,MP,EX,1,2.1E11MP,PRXY,1,0.3MP,DENS,1,7850MP,EX,2,3.5E10MP,PRXY,2,0.1667MP,DENS,2,2500N,,0,0,-5,,,, !创建节点，复制结点NGEN,4,4,ALL,,,12,,,1,NGEN,2,1,ALL,,,,,10,1,NGEN,2,1,2,10,4,,16,,1,NGEN,2,1,3,11,4,,,-10,1,TYPE,1MAT,1REAL,1ESYS,0 !单元坐标系SECNUM,1TSHAP,LINEE,11,14 !建立单元TYPE,1 MAT,1 REAL,1 ESYS,0 SECNUM,2 TSHAP,LINE E,2,6E,6,10E,10,14 E,1,5E,5,9E,9,13E,3,7E,7,11E,4,8E,8,12E,1,2E,3,4E,5,6E,7,8E,9,10E,13,14 TYPE,1 MAT,1 REAL,1 ESYS,0 SECNUM,3 TSHAP,LINE E,3,6E,6,11E,4,5E,5,12E,2,3E,1,4E,6,7E,5,8E,10,11 E,9,12TYPE,2MAT,2REAL,1 ESYS,0TSHAP,QUADE,1,2,6,5E,5,6,10,9E,9,10,14,13NSYM,X,14,ALLESYM,,14,ALLNUMMRG,ALL,,,,LOW NUMCMP,ALL FINISH/SOLNSEL,S,,,23,24D,ALL,,,,,,UX,UY,UZ,,, NSEL,S,,,13,14D,ALL,,,,,,UY,UZ,,, NSEL,S,,,1,2F,ALL,FY,-100000 ALLSEL,ALL ACEL,0,10,0, ANTYPE,0SOLVEFINISH/POST1PLDISP,2PLNSOL,U,SUM,0,1PLVECT,U,,,,VECT,NODE,ON,0ETABLE,zhou_i,SMISC,1ETABLE,zhou_j,SMISC,7ETABLE,zhou_i,SMISC,2ETABLE,zhou_j,SMISC,8ETABLE,zhou_i,SMISC,6ETABLE,zhou_j,SMISC,12PRETAB,ZHOU_I,ZHOU_J,JIAN_I,JIAN_J,WAN_I,WAN_J PLLS,ZHOU_I,ZHOU_J,1,0PRNSOL,U,COMPFINISH/EXIT。

2 梁系问题和杆系相比，梁系结构杆件的弯曲效应就需要考虑了。

这时，结点上的广义位移为线位移的基础上增加了角位移。

在Ａｎｓｙｓ中二维梁单元和三维梁单元的分别是Ｂｅａｍ３和Ｂｅａｍ４，它们结点上的位移分别是３个和６个。

建立梁系结构的方法和求解步骤与前面的杆系结构完全相同，唯一不同的是它们单元特性和单元内力属性有所改变。

这里给出二维超静梁和三维精简后塔机算例来说明梁系结构的有限元计算。

　下面几个实例的代码均在Ａｎｓｙｓ５．６的ＥＤ版中调试通过。

　２．１　Ｅｘ３－３５Ｓ．ｔｘｔ 直接建模求解超静定梁（结点少）　如图所示的两跨超静定梁，整个梁上作用有线形分布的荷载，梁的中点位置上有一个集中力作用。

梁的几何尺寸和其上的荷载和边界条件如图所示，试利用Ａｎｓｙｓ计算该梁的变形和内力。

这里我们采用了两种不同的单元划分方案，分别计算出了变形和内力，并给出了内力图。

从内力图的比较可以看出，不同单元划分对内力图的形状有明显影响。

所以在利用直线连接绘制内力图时需要比较多的节点和单元才可以得出比较准确的内力图。

　第一个例子对两跨超静定梁采用了最少的结点（３个结点）和单元（２个单元）个数来说明梁系结构的有限元计算过程。

由于采用了太少的结点和单元，计算结果画出的的剪力图和弯矩图是不准确的，请读者注意。

　　上图是设置了三个结点，两个单元的模型，完全固定左侧结点，右侧结点不能移动，但可以转动。

了各单元上作用有线性分布的荷载，在中间结点上有向下的集中力。

　经过计算，利用单元表定义技术，利用单元上的画直线方法得到了该梁剪力图（如上图）和弯矩图（如下图）。

　显示这两个内力图的趋势是不对的，但在结点上的内力值是正确的。

　下面是所有命令：　ＦＩＮＩＳＨ　！退出以前模块　／ＣＬＥＡＲ　／ＦＩＬＮＡＭＥ，　ＥＸ３－３５Ｓ ／ＵＮＩＴＳ，　ＣＧＳ ！　单位采用ＣＧＳ　／ＴＩＴＬＥ，　ＥＸ３－３５Ｓ．ｔｘｔ，　３　ＮＯＤＥＳ　ＳＩＭＵＬＡＴＩＯＮ　／ＰＲＥＰ７　！进入前处理模块ＰＲＥＰ７ Ｎ，　１，　０，　０　！定义各个结点　Ｎ，　２，　１５，　０ Ｎ，　３，　３０，　０ ＮＰＬＯＴ ！　图形中不显示结点号码　ＮＰＬＯＴ，　１ ！　图形中显示结点号码　ＮＬＩＳＴ ！　在直角坐标系统下，列出结点资料　ＤＳＹＳ，　１ ！　列出资料时，转换至圆柱坐标　ＮＬＩＳＴ ！　在圆柱坐标系统下，列出结点资料　ＥＴ，　１，　ＢＥＡＭ３ ！定义第一类单元为平面梁单元ＢＥＡＭ３　ＭＰ，　ＥＸ，　１，　２０７ｅ５　！定义第一类材料的弹性模量ＥＸ　Ｒ，　１，　１，　０．０２０８３３，　０．５　！定义单元的第一类实常数：Ａｒｅａ，Ｉｎｅｒｔｉａ，Ｈｅｉｇｈｔ　Ｅ，　１，　２ ！定义各个单元两端结点的连接关系 Ｅ，　２，　３　ＥＰＬＯＴ　！用图形显示单元　ＥＬＩＳＴ　！列出单元信息　ＦＩＮＩＳＨ　！退出前处理模块　／ＳＯＬＵ　！进入求解模块ＳＯＬＵＴＩＯＮ ＡＮＴＹＰＥ，　ＳＴＡＴＩＣ　！申明求解类型是静力分析　ＯＵＴＰＲ，　ＢＡＳＩＣ，　ＡＬＬ ！　在输出结果中，列出元素的结果　Ｄ，　１，　ＵＸ，　０，　，　，　，　，　ＵＹ，　ＲＯＴＺ ！定义１号结点约束所有自由度ＵＸ，ＵＹ和ＲＯＴＺ　Ｄ，　３，　ＵＸ，　０，　，　，　，　，ＵＹ　！约束３号结点的ＵＸ和ＵＹ自由度　ＳＦＢＥＡＭ，　１，　１，　ＰＲＥＳ，　０，　３００　！定义各个单元上的线性分布荷载　ＳＦＢＥＡＭ，　２，　１，　ＰＲＥＳ，　３００，　６００　Ｆ，　２，　ＦＹ，　－１０００ ！定义在２号结点上的沿－Ｙ方向大小为１０００的集中力　ＳＯＬＶＥ　！开始求解 ＦＩＮＩＳＨ ！退出求解模块 ／ＰＯＳＴ１　！进入后处理模块ＰＯＳＴ１ ＰＲＤＩＳＰ ！　显示变形数据　ＥＴＡＢＬＥ，ＩＭＯＭＥＮＴ，　ＳＭＩＳＣ，　６　！　建立元素结果表，元素Ｉ点力矩　ＥＴＡＢＬＥ，ＪＭＯＭＥＮＴ，　ＳＭＩＳＣ，　１２　！　建立元素结果表，元素Ｊ点力矩　ＥＴＡＢＬＥ，　ＩＳＨＥＡＲ，　ＳＭＩＳＣ，　２　！　建立元素结果表，元素Ｉ点剪力　ＥＴＡＢＬＥ，　ＪＳＨＥＡＲ，　ＳＭＩＳＣ，　８　！　建立元素结果表，元素Ｊ点剪力　ＰＲＥＴＡＢ ！　显示单元表资料　／ＴＩＴＬＥ，　ＳＨＥＡＲ　ＦＯＲＣＥ　ＤＩＡＧＲＡＭ　！定义图形窗口标题　ＰＬＬＳ，　ＩＳＨＥＡＲ，　ＪＳＨＥＡＲ ！　结构剪力分布图　／ＴＩＴＬＥ，　ＢＥＮＤＩＮＧ　ＭＯＭＥＮＴ　ＤＩＡＧＲＡＭ　！定义图形窗口标题　ＰＬＬＳ，　ＩＭＯＭＥＮＴ，　ＪＭＯＭＥＮＴ ！　结构弯矩分布图　ＦＩＮＩＳＨ ！退出后处理模块ＰＯＳＴ１　２．２　Ｅｘ３－３５Ｃ．ｔｘｔ 直接建模求解超静定梁（结点较多）　从上面例子看出，要得到比较符合实际情况的内力图，需要设置比较多的结点和单元。

课程名称：道路与桥梁电算（ANSYS作业）题目要求：1.选择合适模型参数及尺寸2.建模求解过程借助ANSYS有限元软件完成题目作业格式要求：1.对所分析的结构做简要说明，说明计算参数取值2.建模过程所需菜单须详细说明3.结果分析中重点对变形、应力、应变规律进行阐述，必要时可借助相关图形进行说明4.作业可手写也可打印。

题目：路面结构内力计算20m1：1.537m尺寸说明：具体建模尺寸以题目尺寸要求为准（1）.材料参数=0.18m-0.2m，E=1500MPa-1800MPa，u=0.3面层：h1=0.3-0.4m，E=1200MPa-1500MPa，u=0.2基层：h1=0.15m-0.2m，E=800-1200MPa，u=0.2底基层：h1土基：E=50MPa，u=0.3（2）路基顶面宽度取16-18米，路基填筑高度6-8米，边坡1：1.2-1.5要求：建模时必须明确说明具体尺寸。

（3）建模要求沥青层网格划分使用map（映射）型网格划分形式（4）荷载形式荷载大小0.7MPa，宽度0.3m（下面是解答，先确定参数，再解答）解：一．结构：面层 h=0.18m，E=1.6E11Pa，u=0.3 基层 h=0.35m，E=1.2E11Pa，u=0.2 底基层h=0.18m，E=1.2E11Pa，u=0.2 土基 E=0.05E11Pa，u=0.3 宽度为18米,土层底层宽度为37.5米,边坡1:1.3，施加荷载0.7Mpa，宽度0.3m。

第一步，建立模型1、选择主菜单Preferences，弹出对话框中点击 Structure，选取结构分析，点OK关闭对话框。

2、定义单元类型。

依次点击主菜单中的Preferences ， ElementType，add/Edit/Delete，点击add按钮，弹出单元选择对话框，单元类型选择Solid组的Quad 4node 42，单击OK按钮，Close关闭对话框。

3.定义材料特性，按照Preprocessor，Material Props，Material Models依次点击，弹出材料属性对话框，对话框右侧选择Structure，双击后选择Linear ， Elastics，Isotropics，双击之后弹出对话框，定义材料的E和u ，然后关闭，再选择Metarial按钮选择New Material,依次建立四个路基层的E，u。

!****************************************************************************** ***********************! 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，不需要考虑纵向折减。

A桥主拱空间有限元静力计算1．模型建立根据A桥初步设计图纸以及横撑ANSYS优化计算结果，建立空间有限元模型。

鉴于A桥为有推力钢管混凝土拱桥，本次计算不考虑下部结构对上部结构的影响，即采用拱脚固接约束的方式进行建模。

全桥采用BEAM188梁单元来模拟主拱、桥面系钢梁、横撑以及其他钢管构件；用SHELL63壳单元来模拟桥面板；用LINK8杆单元来模拟吊杆。

用COMBIN14弹簧元来模拟立柱与桥面横梁的支撑约束。

其中采用双梁共节点方式来模拟钢管与混凝土的组合结构；采用面单元加偏置梁单元来模拟钢梁和混凝土面板的空间受力特性。

用约束节点自由度来模拟简支桥面系。

考虑有限元模型构件体积与实际构造物体积之间的差别，采用不同的换算容重来定义材料特性。

钢材弹性模量Es采用2.06E11Pa，波松比采用0.333，热膨胀系数采用1.2E-6。

主拱混凝土弹性模量（考虑收缩徐变效应）采用Ec=Es/15。

2．计算方法采用线性叠加原理结合主拱主要控制点纵向轴力影响线进行工况计算；采用容许应力法进行应力、反力计算。

本次计算主要考虑结构恒载、汽车活载、人群活载、整体升温、整体降温。

按照《公路桥涵设计通用规范》(JTG D60-2004)进行荷载取值。

3．应力计算结果根据相关规范，采用影响线布载计算各工况的主拱应力，再进行工况组合，对主拱八分点控制应力进行检算。

本次计算根据0/8点、1/8点、2/8点、3/8点、4/8点的影响线，拟订了9种加载方式，得到26种计算工况，最后得到46种工况组合计算结果。

具体如下：3．1影响线计算用ANSYS瞬态分析计算主拱轴力点纵向影响线，结果如下：0/8点轴力影响线1/8点轴力影响线2/8点轴力影响线3/8点轴力影响线4/8点轴力影响线3．2 布载方式的确定根据以上主拱影响线确定如下9种布载方式：4#加载方式5#加载方式8#加载方式3．3 工况确定3．4 组合确定3．5 主拱钢管、混凝土应力计算结果：4．主拱位移计算结果5．吊杆内力结果6．支座反力结果。

＊＊＊桥梁仿真单元类型 (1)一、建议选用的单元类型 (1)二、常见桥梁连接部位 (2)三、桥梁基础的处理方式 (2)＊＊＊桥梁常见模型处理 (2)一、桥梁中常用的模型可以用相应的单元 (2)二、桥梁建模要综合运用各种合适的单元 (3)三、选用合适的分析方法 (3)施加预应力的方式 (3)一、预应力的模拟方式 (3)二、建立预应力的模型 (5)＊＊＊土弹簧的模拟 (5)＊＊＊混凝土的模拟 (5)工况组合 (6)一、典型的荷载工况步骤 (6)二、存储组合后的荷载工况 (6)风荷载的确定 (7)地震波的输入 (7)初应力荷载 (8)Ansys可采用两种方法来实现铰接： (8)AUTOCAD模型输入 (9)用ANSYS作桥梁计算十三（其他文件网格划分） (12)(一)时间选项 (13)(二)子步数和时间步大小 (13)(三)自动时间步长 (14)(四)阶跃或递增载荷 (14)关于阶跃载荷和逐渐递增载荷的说明： (14)一、用于动态和瞬态分析的命令 (14)二、非线性选项 (14)三、输出控制 (15)重新启动一个分析 (16)一、重启动条件 (16)二、一般重启动的步骤 (17)三、边界条件重建 (17)在Ansys单元库中，有近200种单元类型，在本章中将讨论一些在桥梁工程中常用到的单元，包括一些单元的输人参数，如单元名称、节点、自由度、实常数、材料特性、表面荷载、体荷载、专用特性、关键选项KEYOPl等。

＊＊＊关于单元选择问题这是一个大问题，方方面面很多，主要是掌握有限元的理论知识。

首先当然是由问题类型选择不同单元，二维还是三维，梁，板壳，体，细梁，粗梁，薄壳，厚壳，膜等等，再定义你的材料：各向同性或各向异性，混凝土的各项’参数，粘弹性等等。

接下来是单元的划分与网格、精度与求解时间的要求等选择，要对各种单元的专有特性有个大概了解。

使用Ansys，还要了解Ansys的一个特点是笼统与通用，因此很多东西被掩盖到背后去了。

本文档仅需在参数设定部分输入结构尺寸数据，即可完成球铰模型。

!--------------------------------------------------------------!本命令流用四面体自动划分网格计算，体积比较大，!无撑脚!桩基5*5，下承台厚度5m，弹性模量*1!切出桩轴位置!修正上转盘和球铰偏心，墩身不偏心!加钢束，先拉直的，以后再完善出带弯。

单根钢绞线直径5mm，7根成一股，139mm2。

一孔12根或15根。

!--------------------------------------------------------------FINISH/CLEAR,START/COM, Structural*afun,deg !默认单位转为角度!!0.参数设定，定义变量-------------------------------------------------------------/prep7QJR1=2.25 !球铰与下转盘接触面半径QJR2=2.60 !球铰半径QJH1=0.05 !球铰下部高度QJH2=0.85 !球铰上部高度SZPR1=6.25 !上转盘半径SZPH1=1.0 !上转盘下部圆盘部分高度SZPH2=1.8 !上转盘上部矩形柱的高度SZPA1=14 !上转盘上部矩形柱的顺桥边长SZPB1=18.5 !上转盘上部矩形柱的横桥边长XZPH1=5.0 !下转盘的高度XZPA1=19.2 !下转盘上部矩形柱的顺桥边长XZPB1=23.2 !下转盘上部矩形柱的横桥边长ZL1= 9.5 !桩长-----------ZR1= 0.9 !桩半径ZDA1=4.0 !桩顺桥向间距ZDB1=4.7 !桩横桥向间距DSA1=6.0 !墩身的顺桥边长DSA2=2.0 !墩身的顺桥单肢边长DSB1=16 !墩身的横桥边长DSH1=6 !墩身的高度------------NA=5 !桩顺桥向排数NB=5 !桩横桥向排数PXZ=0 !偏心距离---------------CJR1=5.25 !撑脚中心距CJR2=0.51 !撑脚外径CJR3=CJR2-0.025 !撑脚内径CJN1=4 !撑脚对数DSM=40 !墩身材料C40SZPSM=45 !上转盘上部材料C45SZPXM=50 !上转盘下部材料C50QJM=50 !球铰材料C50XZPM=40 !下转盘材料C45ZM=35 !桩基材料C35!!预应力筋，直筋emst=1.95e8 !预应力筋弹性模量kPaas=139e-6 !预应力筋面积m2tf=180 !预应力筋张拉力kNdenss=7.921 !预应力筋密度t/m3GSDA1=2 !预应力筋顺桥向间距m ------------GSDB1=2 !预应力筋横桥向间距m --------------GSNA=3 !预应力筋顺桥向排数（钢筋横桥向布置）------------ GSNB=3 !预应力筋桥横向排数（钢筋顺桥向布置）------------ GSH1=0.2 !预应力筋桥距下承台地面距离m ------GSREAL=7 !预应力筋实常数编号GSM=7 !预应力筋材料编号!自动计算的参数*afun,deg !默认单位转为角度ZDA2=(XZPA1-(NA-1)*ZDA1)/2 !桩中心与下转盘边顺桥向距ZDB2=(XZPB1-(NB-1)*ZDB1)/2 !桩中心与下转盘边横桥向距ZX=-ZL1-XZPH1 !桩底标高ZS=ZX+ZL1 !桩顶标高XZPX=ZS !下转盘底标高XZPS=XZPX+XZPH1 !下转盘顶标高QJX=XZPS !球铰底标高QJS=QJX+QJH1+QJH2 !球铰顶标高SZPX=QJS !上转盘底标高SZPS=SZPX+SZPH1+SZPH2!上转盘顶标高DSX=SZPS !墩身底标高DSS=DSX+DSH1 !墩身顶标高sita=asin(CJR2/CJR1)+0.5 !撑脚与其对称轴夹角!荷载FF=220000 !最大悬臂主梁的重量，单位kN-------MM=0 !偏心弯矩，单位kN/m!材料参数ET,1,45ET,2,65ET,3,63ET,4,link8MP,EX,2,2.1e8 !钢MP,PRXY,2,0.3MP,DENS,2,7.850MP,EX,35,3.3e7 !C35 MP,PRXY,35,0.2MP,DENS,35,2.650MP,EX,40,3.4e7 !C40 MP,PRXY,40,0.2MP,DENS,40,2.650MP,EX,45,3.45e7 !C45 MP,PRXY,45,0.2MP,DENS,45,2.650MP,EX,50,3.55e7 !C50 MP,PRXY,50,0.2MP,DENS,50,2.650MP,EX,55,3.70e7 !C55MP,PRXY,55,0.2MP,DENS,55,2.650MP,EX,GSM,emst !钢束MP,PRXY,GSM,prxyMP,DENS,GSM,denssR,GSREAL,AS,TF/(EMST*AS)*1.0271ALLSEL/VIEW,1,1,1,1/ANG,1/REP,FAST/REPLO!!1.建模--------------------------------------------------------------------------------------------------- !桩基*DO,J,1,NB,1*DO,I,1,NA,1K,,-(NA-1)*ZDA1/2+(I-1)*ZDA1,-ZL1-XZPH1,-(NB-1)*ZDB1/2+(J-1)*ZDB1*ENDDO*ENDDOKGEN,2,ALL,,,,ZL1*DO,I,1,NA*NB,1L,I,I+NA*NB !*ENDDOWPOFF,,-ZL1-XZPH1wpro,,-90.000000, !旋转工作平面*DO,J,1,NB,1*DO,I,1,NA,1!CYL4,-(NA-1)*ZDA1/2+(I-1)*ZDA1,-(NB-1)*ZDB1/2+(J-1)*ZDB1,ZR1,,,,ZL1 !生成桩CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长CYL4,-(NA-1)*ZDA1/2+(I-1)*ZDA1,-(NB-1)*ZDB1/2+(J-1)*ZDB1,ZR1,,,,ZL1+XZPH1 !生成桩CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长*ENDDO*ENDDO!下转盘ALLSEL*get,kmax,kp,,num,maxK,,-XZPA1/2,XZPX,-XZPB1/2K,, XZPA1/2,XZPX,-XZPB1/2K,,-XZPA1/2,XZPX, XZPB1/2K,, XZPA1/2,XZPX, XZPB1/2K,,-XZPA1/2,XZPS,-XZPB1/2K,, XZPA1/2,XZPS,-XZPB1/2K,,-XZPA1/2,XZPS, XZPB1/2K,, XZPA1/2,XZPS, XZPB1/2V,kmax+1,kmax+2,kmax+4,kmax+3,kmax+5,kmax+6,kmax+8,kmax+7 !球铰及上转盘下部建模*get,kmax,kp,,num,maxK,, , QJX !球铰关键点K,, QJR1, QJXK,, , QJX+QJH1K,, QJR2, QJX+QJH1K,, , QJSK,, QJR2, QJSK,,SZPR1, SZPX !上转盘下部关键点K,, , SZPX+SZPH1K,,SZPR1, SZPX+SZPH1A,kmax+1,kmax+2,kmax+4,kmax+3 !球铰下部A,kmax+3,kmax+4,kmax+6,kmax+5 !球铰上部A,kmax+5,kmax+6,kmax+7,kmax+9,kmax+8 !上转盘下部*get,Amax,AREA,,num,maxVROTAT,AMAX,AMAX-1,AMAX-2,,,,KMAX+1,KMAX+3,, !面旋转成体!上转盘上部矩形建模ALLSEL*get,kmax,kp,,num,maxK,,-SZPA1/2,SZPX+SZPH1,-SZPB1/2K,, SZPA1/2,SZPX+SZPH1,-SZPB1/2K,,-SZPA1/2,SZPX+SZPH1, SZPB1/2K,, SZPA1/2,SZPX+SZPH1, SZPB1/2K,,-SZPA1/2,SZPS,-SZPB1/2K,, SZPA1/2,SZPS,-SZPB1/2K,,-SZPA1/2,SZPS, SZPB1/2K,, SZPA1/2,SZPS, SZPB1/2V,kmax+1,kmax+2,kmax+4,kmax+3,kmax+5,kmax+6,kmax+8,kmax+7 !墩身建模ALLSEL*get,kmax,kp,,num,maxK,,-DSA1/2,DSX,-DSB1/2K,, DSA1/2,DSX,-DSB1/2K,,-DSA1/2,DSX, DSB1/2K,, DSA1/2,DSX, DSB1/2K,,-DSA1/2,DSS,-DSB1/2K,, DSA1/2,DSS,-DSB1/2K,,-DSA1/2,DSS, DSB1/2K,, DSA1/2,DSS, DSB1/2V,kmax+1,kmax+2,kmax+4,kmax+3,kmax+5,kmax+6,kmax+8,kmax+7 !!生成撑脚,命令有错，无法自动划分网格!CSYS,0!系统坐标系置为当前坐标系!WPCSYS,-1 !工作平面复原!WPOFF,,-XZPH1!wpro,,-90.000000, !旋转工作平面!*afun,deg !默认单位转为角度!*DO,I,1,CJN1*2,1!XXX1=CJR1*COS((I-1)*180/CJN1+SITA)!YYY1=CJR1*SIN((I-1)*180/CJN1+SITA)!XXX2=CJR1*COS((I-1)*180/CJN1-SITA)!YYY2=CJR1*SIN((I-1)*180/CJN1-SITA)!CYL4,XXX1,YYY1,CJR2,,CJR3,,SZPS-XZPX !生成撑脚CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长!CYL4,XXX2,YYY2,CJR2,,CJR3,,SZPS-XZPX !生成撑脚CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长!*ENDDO!!!移动偏心CSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原vsel,s,loc,y,0.1-XZPH1/2,QJH1+QJH2+SZPH1+SZPH2 !偏心移动CM,SHAGNBU-V,VOLUVGEN, ,SHAGNBU-V, , , , ,PXZ, , ,1ALLSELvplot!!切出双肢墩ALLSCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原wpro,,,-90.000000, !旋转工作平面WPROTA,THXY,THYZ,THZX WPOFF,,,DSA1/2-DSA2 !WPOFF,XOFF,YOFF,ZOFFVSBW,ALLWPOFF,,,-(DSA1-2*DSA2) !WPOFF,XOFF,YOFF,ZOFF VSBW,ALLCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原VSEL,S,LOC,Y,DSX,DSSVSEL,R,LOC,X,-(DSA1/2-DSA2),DSA1/2-DSA2ALLSEL,BELOW,VOLUVPLOTCM,VD,VOLUALLSELVDELE,VDALLSELVPLOTVPTN,ALL !体分割!!预应力筋allsCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原!顺桥向钢筋CSYS,0!系统坐标系置为当前坐标系*DO,I,1,GSNB,1WPCSYS,-1 !工作平面复原WPOFFS,,,-(GSNB-1)*GSDB1/2+(I-1)*GSDB1 VSBW,ALL*ENDDO!横桥向钢筋*DO,I,1,GSNA,1WPCSYS,-1 !工作平面复原WPROTA,,,90WPOFFS,,,-(GSNA-1)*GSDA1/2+(I-1)*GSDA1 VSBW,ALL*ENDDOWPCSYS,-1 !工作平面复原WPOFFW,,XZPX+GSH1 WPROTA,,90VSBW,ALL!!工作面切割实体!切出特征面ALLSELCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原VSBW,ALLWPOFF,,,PXZ !切出偏心VSBW,ALLWPROTA,,,-90VSBW,ALL!横桥向CSYS,0!系统坐标系置为当前坐标系*DO,I,1,NB,1WPCSYS,-1 !工作平面复原WPOFFS,,,-(NB-1)*ZDB1/2+(I-1)*ZDB1 VSBW,ALL!WPOFFS,,,-ZDB2!VSBW,ALL!WPOFFS,,,2*ZDB2!VSBW,ALL*ENDDO!顺桥向*DO,I,1,NA,1WPCSYS,-1 !工作平面复原WPROTA,,,90WPOFFS,,,-(NA-1)*ZDA1/2+(I-1)*ZDA1 VSBW,ALL!WPOFFS,,,-ZDA2!VSBW,ALL!WPOFFS,,,2*ZDA2!VSBW,ALL*ENDDOALLSELVPLOTCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原VSEL,S,LOC,Y,ZX,ZSCM,Z-V,VOLUVSEL,S,LOC,Y,XZPX,XZPSCM,XZP-V,VOLUVSEL,S,LOC,Y,XZPX,XZPSCM,XZP-V,VOLUVSEL,S,LOC,Y,QJX,QJSCM,QJ-V,VOLUVSEL,S,LOC,Y,SZPX,SZPSCM,SZP-V,VOLUVSEL,S,LOC,Y,SZPX,SZPX+SZPH1CM,SZPX-V,VOLUVSEL,S,LOC,Y,SZPS-SZPH2,SZPSCM,SZPS-V,VOLUVSEL,S,LOC,Y,DSX,DSSCM,DS-V,VOLU!顺桥向钢筋LSEL,S,LOC,Z,-(GSNB-1)*GSDB1/2*DO,I,1,GSNB,1LSEL,A,LOC,Z,-(GSNB-1)*GSDB1/2+(I-1)*GSDB1*ENDDOLSEL,R,LOC,Y,XZPX+GSH1-0.01,XZPX+GSH1+0.01 CM,GSS-L,LINE!横桥向钢筋LSEL,S,LOC,X,-(GSNA-1)*GSDA1/2*DO,I,1,GSNA,1LSEL,A,LOC,X,-(GSNA-1)*GSDA1/2+(I-1)*GSDA1*ENDDOLSEL,R,LOC,Y,XZPX+GSH1-0.01,XZPX+GSH1+0.01 CM,GSH-L,LINEALLSEL!!桩单元!CMSEL,S,Z-V!CMPLOT!ALLSEL,BELOW,VOLU!VATT, 1, ,1, 0!LESIZE,ALL,0.5, , , , , , ,0!VSWEEP,Z-V!!下转盘单元!ALLS!CMSEL,S,XZP-V!CMPLOT!ALLSEL,BELOW,VOLU!VATT, 1, ,1, 0!LESIZE,ALL,0.25, , , , , , ,0!VSWEEP,XZP-V!VMESH,XZP-V!!生成单元ALLSELvplot!!划分混凝土网格!ALLSELVATT, 2, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS CM,LL,LINELESIZE,LL,0.5, , , , , , ,0LSEL,U,LOC,Y,ZX,SZPX+0.01CM,LL1,LINELESIZE,LL1,1, , , , , , ,0!墩身C40CMSEL,S,DS-VVATT, DSM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS!上转盘上部C45CMSEL,S,SZPS-VVATT, SZPSM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS !上转盘下部CMSEL,S,SZPX-VVATT, SZPXM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS !球铰CMSEL,S,QJ-VVATT, QJM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS!下转盘CMSEL,S,XZP-VVATT, XZPM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS !桩C35CMSEL,S,Z-VVATT,ZM, , 1, 0 !C35 VATT,MAT,REAL,TYPE,ESYSALLSELMSHAPE,1,3DMSHKEY,0VMESH,ALL!预应力钢筋CMSEL,S,GSS-LCMSEL,A,GSH-LLATT,GSM,GSREAL,4LMESH,ALL!加约束NSEL,S,LOC,Y,ZXCM,Z-D,NODED,Z-D,ALL!!求解!墩顶节点组NSEL,S,LOC,Y,DSSnplotCM,DSS-D,NODENSEL,R,LOC,Z,-DSB1/2-0.1+PXZ,-DSB1/2+0.1+PXZ CM,DSS1-D,NODE*GET,NNUM1,NODE,0,COUNT !获得节点组的节点数量CMSEL,S,DSS-DNSEL,R,LOC,Z,DSB1/2-0.1+PXZ,DSB1/2+0.1+PXZ CM,DSS2-D,NODE*GET,NNUM2,NODE,0,COUNT !获得节点组的节点数量!/ESHAPE,1!ALLSELVPLOT/VIEW,1,1,1,1/ANG,1/REP,FAST/REPLO SAVE,'qiujiao','db' !。

ansys土木经典算例ANSYS是一种常用的有限元软件，广泛应用于土木工程领域。

下面列举了一些经典的ANSYS土木算例。

1. 悬臂梁的静力分析：通过给定的悬臂梁的几何尺寸和加载条件，计算梁的应力、位移和变形情况，以评估梁的承载能力和结构安全性。

2. 桁架结构的模态分析：对给定的桁架结构进行模态分析，确定结构的固有频率和振型，以预测结构的动力响应和抗震性能。

3. 地基基础的强度分析：通过模拟地基基础的加载情况，计算地基的应力分布和变形情况，以评估地基的稳定性和承载能力。

4. 混凝土梁的弯曲分析：对给定的混凝土梁进行弯曲分析，确定梁的应力分布和变形情况，以评估梁的承载能力和结构安全性。

5. 土体的固结分析：通过模拟土体的加载和固结过程，计算土体的应力、应变和固结沉降情况，以评估土体的稳定性和变形特性。

6. 钢结构的塑性极限分析：对给定的钢结构进行塑性极限分析，确定结构的塑性极限强度和塑性变形分布，以评估结构的安全性和可靠性。

7. 水坝的渗流分析：通过模拟水坝的渗流过程，计算水坝的渗流速度、渗流压力和渗流路径，以评估水坝的渗漏风险和抗渗性能。

8. 地下管道的地震响应分析：对给定的地下管道进行地震响应分析，确定管道的动力响应和变形情况，以评估管道的抗震性能和结构安全性。

9. 地铁隧道的地表沉降分析：通过模拟地铁隧道的开挖和施工过程，计算地表沉降和沉降分布，以评估隧道施工对周围建筑物和地下管线的影响。

10. 高层建筑的风载分析：对给定的高层建筑进行风载分析，确定建筑物的风压分布和结构的响应，以评估建筑物的抗风性能和结构安全性。

以上是一些经典的ANSYS土木算例，涵盖了土木工程中的静力分析、动力分析、地基基础分析、固结分析、渗流分析、地震响应分析、地表沉降分析和风载分析等多个方面，可以帮助工程师们更好地理解和评估土木结构的性能和安全性。

本文档仅需在参数设定部分输入结构尺寸数据，即可完成球铰模型。

!--------------------------------------------------------------!本命令流用四面体自动划分网格计算，体积比较大，!无撑脚!桩基5*5，下承台厚度5m，弹性模量*1!切出桩轴位置!修正上转盘和球铰偏心，墩身不偏心!加钢束，先拉直的，以后再完善出带弯。

单根钢绞线直径5mm，7根成一股，139mm2。

一孔12根或15根。

!--------------------------------------------------------------FINISH/CLEAR,START/COM, Structural*afun,deg !默认单位转为角度!!0.参数设定，定义变量-------------------------------------------------------------/prep7QJR1=2.25 !球铰与下转盘接触面半径QJR2=2.60 !球铰半径QJH1=0.05 !球铰下部高度QJH2=0.85 !球铰上部高度SZPR1=6.25 !上转盘半径SZPH1=1.0 !上转盘下部圆盘部分高度SZPH2=1.8 !上转盘上部矩形柱的高度SZPA1=14 !上转盘上部矩形柱的顺桥边长SZPB1=18.5 !上转盘上部矩形柱的横桥边长XZPH1=5.0 !下转盘的高度XZPA1=19.2 !下转盘上部矩形柱的顺桥边长XZPB1=23.2 !下转盘上部矩形柱的横桥边长ZL1= 9.5 !桩长-----------ZR1= 0.9 !桩半径ZDA1=4.0 !桩顺桥向间距ZDB1=4.7 !桩横桥向间距DSA1=6.0 !墩身的顺桥边长DSA2=2.0 !墩身的顺桥单肢边长DSB1=16 !墩身的横桥边长DSH1=6 !墩身的高度------------NA=5 !桩顺桥向排数NB=5 !桩横桥向排数PXZ=0 !偏心距离---------------CJR1=5.25 !撑脚中心距CJR2=0.51 !撑脚外径CJR3=CJR2-0.025 !撑脚内径CJN1=4 !撑脚对数DSM=40 !墩身材料C40SZPSM=45 !上转盘上部材料C45SZPXM=50 !上转盘下部材料C50QJM=50 !球铰材料C50XZPM=40 !下转盘材料C45ZM=35 !桩基材料C35!!预应力筋，直筋emst=1.95e8 !预应力筋弹性模量kPaas=139e-6 !预应力筋面积m2tf=180 !预应力筋张拉力kNdenss=7.921 !预应力筋密度t/m3GSDA1=2 !预应力筋顺桥向间距m ------------GSDB1=2 !预应力筋横桥向间距m --------------GSNA=3 !预应力筋顺桥向排数（钢筋横桥向布置）------------ GSNB=3 !预应力筋桥横向排数（钢筋顺桥向布置）------------ GSH1=0.2 !预应力筋桥距下承台地面距离m ------GSREAL=7 !预应力筋实常数编号GSM=7 !预应力筋材料编号!自动计算的参数*afun,deg !默认单位转为角度ZDA2=(XZPA1-(NA-1)*ZDA1)/2 !桩中心与下转盘边顺桥向距ZDB2=(XZPB1-(NB-1)*ZDB1)/2 !桩中心与下转盘边横桥向距ZX=-ZL1-XZPH1 !桩底标高ZS=ZX+ZL1 !桩顶标高XZPX=ZS !下转盘底标高XZPS=XZPX+XZPH1 !下转盘顶标高QJX=XZPS !球铰底标高QJS=QJX+QJH1+QJH2 !球铰顶标高SZPX=QJS !上转盘底标高SZPS=SZPX+SZPH1+SZPH2!上转盘顶标高DSX=SZPS !墩身底标高DSS=DSX+DSH1 !墩身顶标高sita=asin(CJR2/CJR1)+0.5 !撑脚与其对称轴夹角!荷载FF=220000 !最大悬臂主梁的重量，单位kN-------MM=0 !偏心弯矩，单位kN/m!材料参数ET,1,45ET,2,65ET,3,63ET,4,link8MP,EX,2,2.1e8 !钢MP,PRXY,2,0.3MP,DENS,2,7.850MP,EX,35,3.3e7 !C35 MP,PRXY,35,0.2MP,DENS,35,2.650MP,EX,40,3.4e7 !C40 MP,PRXY,40,0.2MP,DENS,40,2.650MP,EX,45,3.45e7 !C45 MP,PRXY,45,0.2MP,DENS,45,2.650MP,EX,50,3.55e7 !C50 MP,PRXY,50,0.2MP,DENS,50,2.650MP,EX,55,3.70e7 !C55MP,PRXY,55,0.2MP,DENS,55,2.650MP,EX,GSM,emst !钢束MP,PRXY,GSM,prxyMP,DENS,GSM,denssR,GSREAL,AS,TF/(EMST*AS)*1.0271ALLSEL/VIEW,1,1,1,1/ANG,1/REP,FAST/REPLO!!1.建模--------------------------------------------------------------------------------------------------- !桩基*DO,J,1,NB,1*DO,I,1,NA,1K,,-(NA-1)*ZDA1/2+(I-1)*ZDA1,-ZL1-XZPH1,-(NB-1)*ZDB1/2+(J-1)*ZDB1*ENDDO*ENDDOKGEN,2,ALL,,,,ZL1*DO,I,1,NA*NB,1L,I,I+NA*NB !*ENDDOWPOFF,,-ZL1-XZPH1wpro,,-90.000000, !旋转工作平面*DO,J,1,NB,1*DO,I,1,NA,1!CYL4,-(NA-1)*ZDA1/2+(I-1)*ZDA1,-(NB-1)*ZDB1/2+(J-1)*ZDB1,ZR1,,,,ZL1 !生成桩CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长CYL4,-(NA-1)*ZDA1/2+(I-1)*ZDA1,-(NB-1)*ZDB1/2+(J-1)*ZDB1,ZR1,,,,ZL1+XZPH1 !生成桩CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长*ENDDO*ENDDO!下转盘ALLSEL*get,kmax,kp,,num,maxK,,-XZPA1/2,XZPX,-XZPB1/2K,, XZPA1/2,XZPX,-XZPB1/2K,,-XZPA1/2,XZPX, XZPB1/2K,, XZPA1/2,XZPX, XZPB1/2K,,-XZPA1/2,XZPS,-XZPB1/2K,, XZPA1/2,XZPS,-XZPB1/2K,,-XZPA1/2,XZPS, XZPB1/2K,, XZPA1/2,XZPS, XZPB1/2V,kmax+1,kmax+2,kmax+4,kmax+3,kmax+5,kmax+6,kmax+8,kmax+7 !球铰及上转盘下部建模*get,kmax,kp,,num,maxK,, , QJX !球铰关键点K,, QJR1, QJXK,, , QJX+QJH1K,, QJR2, QJX+QJH1K,, , QJSK,, QJR2, QJSK,,SZPR1, SZPX !上转盘下部关键点K,, , SZPX+SZPH1K,,SZPR1, SZPX+SZPH1A,kmax+1,kmax+2,kmax+4,kmax+3 !球铰下部A,kmax+3,kmax+4,kmax+6,kmax+5 !球铰上部A,kmax+5,kmax+6,kmax+7,kmax+9,kmax+8 !上转盘下部*get,Amax,AREA,,num,maxVROTAT,AMAX,AMAX-1,AMAX-2,,,,KMAX+1,KMAX+3,, !面旋转成体!上转盘上部矩形建模ALLSEL*get,kmax,kp,,num,maxK,,-SZPA1/2,SZPX+SZPH1,-SZPB1/2K,, SZPA1/2,SZPX+SZPH1,-SZPB1/2K,,-SZPA1/2,SZPX+SZPH1, SZPB1/2K,, SZPA1/2,SZPX+SZPH1, SZPB1/2K,,-SZPA1/2,SZPS,-SZPB1/2K,, SZPA1/2,SZPS,-SZPB1/2K,,-SZPA1/2,SZPS, SZPB1/2K,, SZPA1/2,SZPS, SZPB1/2V,kmax+1,kmax+2,kmax+4,kmax+3,kmax+5,kmax+6,kmax+8,kmax+7 !墩身建模ALLSEL*get,kmax,kp,,num,maxK,,-DSA1/2,DSX,-DSB1/2K,, DSA1/2,DSX,-DSB1/2K,,-DSA1/2,DSX, DSB1/2K,, DSA1/2,DSX, DSB1/2K,,-DSA1/2,DSS,-DSB1/2K,, DSA1/2,DSS,-DSB1/2K,,-DSA1/2,DSS, DSB1/2K,, DSA1/2,DSS, DSB1/2V,kmax+1,kmax+2,kmax+4,kmax+3,kmax+5,kmax+6,kmax+8,kmax+7 !!生成撑脚,命令有错，无法自动划分网格!CSYS,0!系统坐标系置为当前坐标系!WPCSYS,-1 !工作平面复原!WPOFF,,-XZPH1!wpro,,-90.000000, !旋转工作平面!*afun,deg !默认单位转为角度!*DO,I,1,CJN1*2,1!XXX1=CJR1*COS((I-1)*180/CJN1+SITA)!YYY1=CJR1*SIN((I-1)*180/CJN1+SITA)!XXX2=CJR1*COS((I-1)*180/CJN1-SITA)!YYY2=CJR1*SIN((I-1)*180/CJN1-SITA)!CYL4,XXX1,YYY1,CJR2,,CJR3,,SZPS-XZPX !生成撑脚CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长!CYL4,XXX2,YYY2,CJR2,,CJR3,,SZPS-XZPX !生成撑脚CYL4,X坐标,y坐标,半径1,起始角度,半径2,终止角度,桩长!*ENDDO!!!移动偏心CSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原vsel,s,loc,y,0.1-XZPH1/2,QJH1+QJH2+SZPH1+SZPH2 !偏心移动CM,SHAGNBU-V,VOLUVGEN, ,SHAGNBU-V, , , , ,PXZ, , ,1ALLSELvplot!!切出双肢墩ALLSCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原wpro,,,-90.000000, !旋转工作平面WPROTA,THXY,THYZ,THZX WPOFF,,,DSA1/2-DSA2 !WPOFF,XOFF,YOFF,ZOFFVSBW,ALLWPOFF,,,-(DSA1-2*DSA2) !WPOFF,XOFF,YOFF,ZOFF VSBW,ALLCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原VSEL,S,LOC,Y,DSX,DSSVSEL,R,LOC,X,-(DSA1/2-DSA2),DSA1/2-DSA2ALLSEL,BELOW,VOLUVPLOTCM,VD,VOLUALLSELVDELE,VDALLSELVPLOTVPTN,ALL !体分割!!预应力筋allsCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原!顺桥向钢筋CSYS,0!系统坐标系置为当前坐标系*DO,I,1,GSNB,1WPCSYS,-1 !工作平面复原WPOFFS,,,-(GSNB-1)*GSDB1/2+(I-1)*GSDB1 VSBW,ALL*ENDDO!横桥向钢筋*DO,I,1,GSNA,1WPCSYS,-1 !工作平面复原WPROTA,,,90WPOFFS,,,-(GSNA-1)*GSDA1/2+(I-1)*GSDA1 VSBW,ALL*ENDDOWPCSYS,-1 !工作平面复原WPOFFW,,XZPX+GSH1 WPROTA,,90VSBW,ALL!!工作面切割实体!切出特征面ALLSELCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原VSBW,ALLWPOFF,,,PXZ !切出偏心VSBW,ALLWPROTA,,,-90VSBW,ALL!横桥向CSYS,0!系统坐标系置为当前坐标系*DO,I,1,NB,1WPCSYS,-1 !工作平面复原WPOFFS,,,-(NB-1)*ZDB1/2+(I-1)*ZDB1 VSBW,ALL!WPOFFS,,,-ZDB2!VSBW,ALL!WPOFFS,,,2*ZDB2!VSBW,ALL*ENDDO!顺桥向*DO,I,1,NA,1WPCSYS,-1 !工作平面复原WPROTA,,,90WPOFFS,,,-(NA-1)*ZDA1/2+(I-1)*ZDA1 VSBW,ALL!WPOFFS,,,-ZDA2!VSBW,ALL!WPOFFS,,,2*ZDA2!VSBW,ALL*ENDDOALLSELVPLOTCSYS,0!系统坐标系置为当前坐标系WPCSYS,-1 !工作平面复原VSEL,S,LOC,Y,ZX,ZSCM,Z-V,VOLUVSEL,S,LOC,Y,XZPX,XZPSCM,XZP-V,VOLUVSEL,S,LOC,Y,XZPX,XZPSCM,XZP-V,VOLUVSEL,S,LOC,Y,QJX,QJSCM,QJ-V,VOLUVSEL,S,LOC,Y,SZPX,SZPSCM,SZP-V,VOLUVSEL,S,LOC,Y,SZPX,SZPX+SZPH1CM,SZPX-V,VOLUVSEL,S,LOC,Y,SZPS-SZPH2,SZPSCM,SZPS-V,VOLUVSEL,S,LOC,Y,DSX,DSSCM,DS-V,VOLU!顺桥向钢筋LSEL,S,LOC,Z,-(GSNB-1)*GSDB1/2*DO,I,1,GSNB,1LSEL,A,LOC,Z,-(GSNB-1)*GSDB1/2+(I-1)*GSDB1*ENDDOLSEL,R,LOC,Y,XZPX+GSH1-0.01,XZPX+GSH1+0.01 CM,GSS-L,LINE!横桥向钢筋LSEL,S,LOC,X,-(GSNA-1)*GSDA1/2*DO,I,1,GSNA,1LSEL,A,LOC,X,-(GSNA-1)*GSDA1/2+(I-1)*GSDA1*ENDDOLSEL,R,LOC,Y,XZPX+GSH1-0.01,XZPX+GSH1+0.01 CM,GSH-L,LINEALLSEL!!桩单元!CMSEL,S,Z-V!CMPLOT!ALLSEL,BELOW,VOLU!VATT, 1, ,1, 0!LESIZE,ALL,0.5, , , , , , ,0!VSWEEP,Z-V!!下转盘单元!ALLS!CMSEL,S,XZP-V!CMPLOT!ALLSEL,BELOW,VOLU!VATT, 1, ,1, 0!LESIZE,ALL,0.25, , , , , , ,0!VSWEEP,XZP-V!VMESH,XZP-V!!生成单元ALLSELvplot!!划分混凝土网格!ALLSELVATT, 2, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS CM,LL,LINELESIZE,LL,0.5, , , , , , ,0LSEL,U,LOC,Y,ZX,SZPX+0.01CM,LL1,LINELESIZE,LL1,1, , , , , , ,0!墩身C40CMSEL,S,DS-VVATT, DSM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS!上转盘上部C45CMSEL,S,SZPS-VVATT, SZPSM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS !上转盘下部CMSEL,S,SZPX-VVATT, SZPXM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS !球铰CMSEL,S,QJ-VVATT, QJM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS!下转盘CMSEL,S,XZP-VVATT, XZPM, , 1, 0 !C40 VATT,MAT,REAL,TYPE,ESYS !桩C35CMSEL,S,Z-VVATT,ZM, , 1, 0 !C35 VATT,MAT,REAL,TYPE,ESYSALLSELMSHAPE,1,3DMSHKEY,0VMESH,ALL!预应力钢筋CMSEL,S,GSS-LCMSEL,A,GSH-LLATT,GSM,GSREAL,4LMESH,ALL!加约束NSEL,S,LOC,Y,ZXCM,Z-D,NODED,Z-D,ALL!!求解!墩顶节点组NSEL,S,LOC,Y,DSSnplotCM,DSS-D,NODENSEL,R,LOC,Z,-DSB1/2-0.1+PXZ,-DSB1/2+0.1+PXZ CM,DSS1-D,NODE*GET,NNUM1,NODE,0,COUNT !获得节点组的节点数量CMSEL,S,DSS-DNSEL,R,LOC,Z,DSB1/2-0.1+PXZ,DSB1/2+0.1+PXZ CM,DSS2-D,NODE*GET,NNUM2,NODE,0,COUNT !获得节点组的节点数量!/ESHAPE,1!ALLSELVPLOT/VIEW,1,1,1,1/ANG,1/REP,FAST/REPLO SAVE,'qiujiao','db' !。

课程名称：桥梁工程B设计题目：ansys梁桥模拟计算（三跨）院系：土木工程系专业：年级：姓名：学号：指导教师：西南交通大学峨眉校区年月日ansys梁桥模拟计算（三跨）1.绪论1.1设计目的桥梁结构分析计算是非常重要的一门技能。

通过本课程设计，掌握一门通用有限元软件分析工具，能够独立对桥梁结构进行静力或动力分析。

本课程具体要求掌握通用有限元软件ANSYS, 了解其前处理，后处理过程以及单元应用。

通过此课程设计的学习，初步具有独立进行结构分析的能力，从而了解桥梁的具体设计。

1.2设计内容及要求桥梁结构建模、确定边界条件、求解、后处理以及分析结论1、了解所选用Beam4等单元的属性和用法；2、对桥梁进行结构离散化，建立三维有限元数值模型；3、正确地对桥梁有限元模型设定边界条件；4、掌握数值分析静力或动力求解方法；5、对计算结果进行后处理，掌握基本作图软件应用；6、对计算结果进行分析，得出结论。

2.有限元分析2.1简介有限元分析是用较简单的问题代替复杂问题后再求解。

它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。

这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。

由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。

有限元是那些集合在一起能够表示实际连续域的离散单元。

有限元的概念早在几个世纪前就已产生并得到了应用，例如用多边形（有限个直线单元）逼近圆来求得圆的周长，但作为一种方法而被提出，则是最近的事。

有限元法最初被称为矩阵近似方法，应用于航空器的结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。

经过短短数十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。

【ANSYS 算例】3.4.2(1) 基于图形界面的桁架桥梁结构分析(step by step)下面以一个简单桁架桥梁为例，以展示有限元分析的全过程。

背景素材选自位于密执安的"Old North Park Bridge" (1904 - 1988)，见图3-22。

该桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表3-6。

桥长L=32m,桥高H=5.5m 。

桥身由8段桁架组成，每段长4m 。

该桥梁可以通行卡车，若这里仅考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg ，若取一半的模型，可以将卡车对桥梁的作用力简化为P 1 ，P 2和P 3 ，其中P 1= P 3=5000 N, P 2=10000N ，见图3-23。

图3-22位于密执安的"Old North Park Bridge" (1904 - 1988)图3-23 桥梁的简化平面模型(取桥梁的一半)表3-6 桥梁结构中各种构件的几何性能参数构件 惯性矩m 4 横截面积m 2顶梁及侧梁(Beam1) 643.8310m -⨯322.1910m -⨯ 桥身弦梁(Beam2) 61.8710-⨯31.18510-⨯ 底梁(Beam3)68.4710-⨯ 33.03110-⨯解答 以下为基于ANSYS 图形界面(Graphic User Interface , GUI)的菜单操作流程。

(1) 进入ANSYS （设定工作目录和工作文件）程序 → ANSYS → ANSYS Interactive → Working directory （设置工作目录）→ Initial jobname （设置工作文件名）:TrussBridge → Run → OK(2) 设置计算类型ANSYS Main Menu ：Preferences … → Structural → OK(3) 定义单元类型ANSYS Main Menu：Preprocessor → Element Type → Add/Edit/Delete... → Add…→ Beam: 2d elastic 3 → OK（返回到Element Types窗口）→ Close(4) 定义实常数以确定梁单元的截面参数ANSYS Main Menu: Preprocessor → Real Constants…→ Add/Edit/Delete → Add…→ select Type 1 Beam 3 → OK → input Real Constants Set No. : 1 , AREA: 2.19E-3，Izz: 3.83e-6(1号实常数用于顶梁和侧梁) → Apply → input Real Constants Set No. : 2 , AREA: 1.185E-3，Izz: 1.87E-6 (2号实常数用于弦杆) → Apply → input Real Constants Set No. : 3, AREA: 3.031E-3，Izz: 8.47E-6 (3号实常数用于底梁) → OK (back to Real Constants window) →Close(the Real Constants window) (5) 定义材料参数ANSYS Main Menu: Preprocessor → Material Props → Material Models → Structural → Linear → Elastic → Isotropic → input EX: 2.1e11, PRXY: 0.3(定义泊松比及弹性模量) → OK → Density (定义材料密度) →input DENS: 7800, → OK → Close（关闭材料定义窗口）(6) 构造桁架桥模型生成桥体几何模型ANSYS Main Menu：Preprocessor → Modeling → Create → Keypoints → In Active CS → NPT Keypoint number：1，X，Y，Z Location in active CS：0，0 → Apply →同样输入其余15个特征点坐标（最左端为起始点，坐标分别为 (4,0), (8,0), (12,0), (16,0), (20,0), (24,0), (28,0), (32,0), (4,5.5), (8,5.5), (12,5.5), (16.5.5), (20,5.5), (24,5.5), (28,5.5)）→Lines →Lines →Straight Line →依次分别连接特征点→ OK网格划分ANSYS Main Menu: Preprocessor → Meshing → Mesh Attributes → Picked Lines →选择桥顶梁及侧梁→ OK → select REAL: 1, TYPE: 1 → Apply →选择桥体弦杆→ OK → select REAL: 2, TYPE: 1 → Apply →选择桥底梁→ OK → select REAL: 3, TYPE:1 → OK →ANSYS Main Menu：Preprocessor → Meshing → MeshTool →位于Size Controls下的Lines：Set → Element Size on Picked → Pick all → Apply → NDIV：1 → OK → Mesh → Lines → Pick all → OK (划分网格) (7) 模型加约束ANSYS Main Menu: Solution → Define Loads → Apply → Structural→ Displacement → On Nodes →选取桥身左端节点→ OK → select Lab2: All DOF(施加全部约束)→ Apply →选取桥身右端节点→ OK → select Lab2: UY(施加Y方向约束)→ OK(8) 施加载荷ANSYS Main Menu: Solution → Define Loads → Apply → Structural → Force/Moment → On Keypoints →选取底梁上卡车两侧关键点（X坐标为12及20）→ OK → select Lab: FY，Value: -5000 → Apply →选取底梁上卡车中部关键点（X坐标为16）→ OK → select Lab: FY，Value: -10000 → OK →ANSYS Utility Menu：→ Select → Everything(9) 计算分析ANSYS Main Menu：Solution → Solve → Current LS → OK(10) 结果显示ANSYS Main Menu：General Postproc → Plot Results → Deformed shape → Def shape only →OK（返回到Plot Results）→ Contour Plot → Nodal Solu → DOF Solution, Y-Component of Displacement → OK（显示Y方向位移UY）(见图3-24(a))定义线性单元I节点的轴力ANSYS Main Menu →General Postproc →Element Table → Define Table →Add →Lab: [bar_I], By sequence num: [SMISC,1] → OK → Close定义线性单元J节点的轴力ANSYS Main Menu →General Postproc →Element Table → Define Table →Add →Lab: [bar_J], By sequence num: [SMISC,1] → OK → Close画出线性单元的受力图(见图3-24(b))ANSYS Main Menu → General Postproc → Plot Results → Contour Plot → Line Elem Res →LabI: [ bar_I], LabJ: [ bar_J], Fact: [1] → OK(11) 退出系统ANSYS Utility Menu：File → Exit → Save Everything → OK(a)桥梁中部最大挠度值为0.003 374m (b)桥梁中部轴力最大值为25 380N图3.24 桁架桥挠度UY以及单元轴力计算结果【ANSYS算例】3.4.2(2) 基于命令流方式的桁架桥梁结构分析!%%%%% [ANSYS算例]3.4.2(2) %%%%% begin %%%%%%!------注：命令流中的符号$，可将多行命令流写成一行------/prep7 !进入前处理/PLOPTS,DATE,0 !设置不显示日期和时间!=====设置单元和材料ET,1,BEAM3 !定义单元类型R,1,2.19E-3,3.83e-6, , , , , !定义1号实常数用于顶梁侧梁R,2,1.185E-3,1.87e-6,0,0,0,0, !定义2号实常数用于弦杆R,3,3.031E-3,8.47E-6,0,0,0,0, !定义3号实常数用于底梁MP,EX,1,2.1E11 !定义材料弹性模量MP,PRXY,1,0.30 !定义材料泊松比MP,DENS,1,,7800 !定义材料密度!-----定义几何关键点K,1,0,0,, $ K,2,4,0,, $ K,3,8,0,, $K,4,12,0,, $K,5,16,0,, $K,6,20,0,, $K,7,24,0,, $K,8,28,0,, $K,9,32,0,, $K,10,4,5.5,, $K,11,8,5.5,, $K,12,12,5.5,, $K,13,16,5.5,, $K,14,20,5.5,, $K,15,24,5.5,, $K,16,28,5.5,,!-----通过几何点生成桥底梁的线L,1,2 $L,2,3 $L,3,4 $L,4,5 $L,5,6 $L,6,7 $L,7,8 $L,8,9!------生成桥顶梁和侧梁的线L,9,16 $L,15,16 $L,14,15 $L,13,14 $L,12,13 $L,11,12 $L,10,11 $L,1,10!------生成桥身弦杆的线L,2,10 $L,3,10 $L,3,11 $L,4,11 $L,4,12 $L,4,13 $L,5,13 $L,6,13 $L,6,14 $L,6,15 $L,7,15 $L,7,16 $L,8,16!------选择桥顶梁和侧梁指定单元属性LSEL,S,,,9,16,1,LATT,1,1,1,,,,!-----选择桥身弦杆指定单元属性LSEL,S,,,17,29,1,LATT,1,2,1,,,,!-----选择桥底梁指定单元属性LSEL,S,,,1,8,1,LATT,1,3,1,,,,!------划分网格AllSEL,all !再恢复选择所有对象LESIZE,all,,,1,,,,,1 !对所有对象进行单元划分前的分段设置LMESH,all !对所有几何线进行单元划分!=====在求解模块中，施加位移约束、外力，进行求解/soluNSEL,S,LOC,X,0 !根据几何位置选择节点D,all,,,,,,ALL,,,,, !对所选择的节点施加位移约束AllSEL,all !再恢复选择所有对象NSEL,S,LOC,X,32 !根据几何位置选择节点D,all,,,,,,,UY,,,, !对所选择的节点施加位移约束ALLSEL,all !再恢复选择所有对象!------基于几何关键点施加载荷FK,4,FY,-5000 $FK,6,FY,-5000 $FK,5,FY,-10000/replot !重画图形Allsel,all !选择所有信息(包括所有节点、单元和载荷等)solve !求解!=====进入一般的后处理模块/post1 !后处理PLNSOL, U,Y, 0,1.0 !显示Y方向位移PLNSOL, U,X, 0,1.0 !显示X方向位移!------显示线单元轴力------ETABLE,bar_I,SMISC, 1ETABLE,bar_J,SMISC, 1PLLS,BAR_I,BAR_J,0.5,1 !画出轴力图finish !结束!%%%%% [ANSYS算例]3.4.2(2) %%%%% end %%%%%%【ANSYS算例】3.2.5(3) 四杆桁架结构的有限元分析下面针对【典型例题】3.2.5(1)的问题，在ANSYS平台上，完成相应的力学分析。

（ii ）纵向钢筋：PIPE20 （iii ）横向箍筋：PIPE202.2 材料性质（i ）、混凝土材料表5-4 混凝土材料的输入参数一览表[16~19]·单轴受压应力-应变曲线（εσ-曲线）在ANSYS ○R程序分析中，需要给出混凝土单轴受压下的应力应变曲线。

在本算例中，混凝土单轴受压下的应力应变采用Sargin 和Saenz 模型[17,18]：221⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-+=c c s c c E E E εεεεεσ （5-30）式中取4'4')108.0028.1(c c c f f -=ε；断面图配筋图断面图配筋图断面图配筋图RCBEAM-01RCBEAM-02RCBEAM-03图5-12 各梁FEM模型断面图（a）单元网格图（b）钢筋单元划分图图5-13 算例（一）的FEM模型图2.4 模型求解在ANSYS○R程序中，对于非线性分析，求解步的设置很关键，对计算是否收敛关系很大，对于混凝土非线性有限元分析，在计算时间容许的情况下，较多的求解子步（Substeps）或较小的荷载步和一个非常大的最大子步数更容易导致收敛[2]。

在本算例中，设置了100个子步。

最终本算例收敛成功，在CPU为P41.6G、内存为256MB的微机上计算，耗时约为8小时。

2.5 计算结果及分析2.5.1 荷载—位移曲线图5-14为ANSYS○R程序所得到的各梁的荷载-跨中挠度曲线，从图中可以看出：（i）、梁RCBEAM-01：曲线形状能基本反映钢筋混凝土适筋梁剪切破坏的受力特点，而且荷载-跨中挠度曲线与钢筋混凝土梁的弯剪破坏形态非常类似，即当跨中弯矩最大截面的纵筋屈服后，由于裂缝的开展，压区混凝土的面积逐渐减小，在荷载几乎不增加的情况下，压区混凝土所受的正应力和剪应力还在不断增加，当应力达到混凝土强度极限时，剪切破坏发生，荷载突然降低。

（ii）、梁RCBEAM-02：荷载-跨中挠度曲线与超筋梁的试验荷载-跨中挠度曲线很相似，在荷载达到极限情况下，没有出现屈服平台，而是突然跌落。

ANSYS四跨连续梁的内力计算四跨连续梁模型图如下所示，各个杆件抗弯刚度EI相同,利用平面梁单元分析它的变形和内力1.结构力学分析利用结构力学方法可以求出这个连续梁的剪力图和弯矩图如下这里只给出了梁的弯曲刚度相同条件,没有指定梁截面的几何参数和材料的力学性质.从结构力学分析的条件上看,这些条件对于确定梁的内力已经足够，但是对于梁的变形分析和应力计算，还需要补充材料的力学参数和截面几何参数。

所以以下分析中，假定梁的截面面积位0。

3m2，抗弯惯性矩为0.003m4,截面高度为0.1m；材料的弹性模量为1000kN/m2，泊松比为0。

3。

补充这些参数对于梁的内力没有影响,但是对于梁的变形和应力是有影响的。

2.用节点和单元的直接建模求解按照前面模型示意图布置节点和单元，在图示坐标系里定位节点的坐标和单元连接信息，以及荷载作用情况和位移约束.由于第二跨中间有两个集中力,所以在集中力位置设置两个节点。

这样，就可以将这两个集中力直接处理成节点荷载。

对于平面梁单元的节点只需输入平面上的两个坐标值，所以这里只输入节点的x坐标和y坐标.（1）指定为结构分析运行主菜单中preference偏好设定命令,然后在对话框中,指定分析模块为structural结构分析，然后单击ok按钮（2）新建单元类型运行主菜单preprocessor—element type-add/edit/delete命令，接着在对话框中单击add按钮新建单元类型（3）定义单元类型先选择单元为beam,接着选2d elastic 3，然后单击ok按钮确定，完成单元类型的选择（4）关闭单元类型的对话框回到单元类型对话框,已经新建了beam3的单元，单击对话框close按钮关闭对话框（5）定义实力常量运行主菜单preprocessor-real constants—add/edit/delete命令，接着在对话框中单击add按钮新建实力常量接着选择定义单元beam3的实力常量，选择后单击ok按钮，然后输入该单元的截面积为0。

FINI/CLEARSPAN=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.0Y0=0.0 !定义原点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,1ESYS,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-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 !定义节点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,NODE,ALL,LOC,X !给NODE_X赋值（节点X坐标）*ABSET,'Creating Efficacy Data ......',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坐标）*VGET,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,NMAX-1,SMISC,12 !读取最后一个单元的J端弯矩MZ *ENDDO*ABFINISH !终止状态条（*ABSET *ABCHECK *ABFINISH须联用）TRANS_V(0)=!!!!QY和MZ的影响线建立结束!!!!********************************!!!!********************************!!!!计算影响线的最大、最小值并决定其位置*DIM,COLUMN_V,ARRAY,NMAX !定义中间临时数组*ABSET,'Geting Extremum ......',BAR !初始化状态条*DO,I,1,NMAX*ABCHECK,NINT(100*I/NMAX) !更新状态条*VCOL,1*MFUN,COLUMN_V(1),COPY,N_QY(1,I) !从N_QY中读取第1列数据给COLUMN_V *VSCFUN,INDEX,LMAX,COLUMN_V !从COLUMN_V中提取最大值所处位置序号QY_EXTR(1,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值*VSCFUN,QY_EXTR(2,I),MAX,COLUMN_V !从COLUMN_V表中提取QY最大值*VSCFUN,INDEX,LMIN,COLUMN_V !从COLUMN_V中提取最小值所处位置序号QY_EXTR(3,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值*VSCFUN,QY_EXTR(4,I),MIN,COLUMN_V !从COLUMN_V表中提取QY最小值*VCOL,1*MFUN,COLUMN_V(1),COPY,N_QY(1,I) !从N_MZ中读取第1列数据给COLUMN_V *VSCFUN,INDEX,LMAX,COLUMN_V !从COLUMN_V中提取最大值所处位置序号MZ_EXTR(1,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值*VSCFUN,MZ_EXTR(2,I),MAX,COLUMN_V !从COLUMN_V表中提取MZ最大值*VSCFUN,INDEX,LMIN,COLUMN_V !从COLUMN_V中提取最小值所处位置序号MZ_EXTR(3,I)=NODE_X(INDEX) !根据序号从NODE_X中提取坐标值*VSCFUN,MZ_EXTR(4,I),MIN,COLUMN_V !从COLUMN_V表中提取MZ最小值*ENDDO*ABFINISH !终止状态条COLUMN_V(0)=!!!!计算影响线的最大、最小值结束!!!!********************************/eof !结束该宏的执行，以下可以选择执行/AXLAB,X,Span_X/XRANGE,0,SPAN*L/GROPT,DIVX,12/GROPT,CURL,1/AXLAB,Y,VALUE*VPLOT,N_QY(1,0),N_QY(1,1), 151*VPLOT,N_MZ(1,0),N_MZ(1,1),51,101,151。

桥梁计算（常用的计算方法)在Ansys单元库中，有近200种单元类型,在本章中将讨论一些在桥梁工程中常用到的单元,包括一些单元的输人参数，如单元名称、节点、自由度、实常数、材料特性、表面荷载、体荷载、专用特性、关键选项KEYOPl等。

*＊＊关于单元选择问题这是一个大问题，方方面面很多，主要是掌握有限元的理论知识.首先当然是由问题类型选择不同单元，二维还是三维，梁，板壳,体，细梁，粗梁，薄壳，厚壳，膜等等，再定义你的材料：各向同性或各向异性,混凝土的各项’参数，粘弹性等等.接下来是单元的划分与网格、精度与求解时间的要求等选择，要对各种单元的专有特性有个大概了解。

使用Ansys，还要了解Ansys的一个特点是笼统与通用，因此很多东西被掩盖到背后去了。

比如单元类型，在Solid里面看到十几种选择，Solid45, Solidl85，Solid95等，看来区别只是节点数目上。

但是实际上每种类型里还有Keyopt分成多种类型,比如最常用的线性单元Solid45,其Keyopt(1）:in●cludeorexclude extradisplacement shapes，就分为非协调元和协调元，Keyopt (2）：fullintegration。

rreducedintegration其实又是两种不同的单元，这样不同组合一下这个Solid45实际上是包含了6种不同单元,各有各的不同特点和用处。

因此使用Ansys要注意各单元的Keyopt选项.不同的选项会产生不同的结果。

·举例来说：对线性元例如Solid45，要想把弯曲问题计算得比较精确,必须要采用非协调模式.采用完全积分会产生剪切锁死，减缩积分又会产生零能模式(ZEM）,非协调的线性元可以达到很高的精度，并且计算量比高阶刷、很多，在变形较大时,用Enhanced Strain比非协调位移模式(Enhaced Displacement）更好(Solidl85)。