Ansys桥梁计算

根据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桥梁计算汇编

Ansys桥梁计算汇编摘要本文将介绍 Ansys 软件在桥梁计算方面的应用。

桥梁是人类文明的重要标志之一，而现代桥梁的建造离不开计算机技术的支持。

Ansys 软件作为一种常用的工程计算软件，可以在计算桥梁结构强度、稳定性和振动等方面发挥重要作用。

本文将从桥梁计算的基础原理、Ansys 软件的应用方法、实际案例等方面进行详细介绍。

桥梁计算基础原理在桥梁计算中，常采用有限元分析（FEA）方法。

该方法将桥梁分割成许多小的有限元，再对每个元素进行计算，最终得到整个桥梁的力学特性。

桥梁的有限元模型需要包括以下几个方面：1.桥梁的几何形状、材料特性、截面性质等参数；2.桥梁外部荷载的大小和作用方向；3.桥梁支座的约束条件。

在FE 模型中，桥梁结构被分成许多小的单元。

通过对每个单元进行力学计算，可以得到整个结构的力学性质。

FEA 方法可以更准确地预测桥梁在不同荷载条件下的强度、稳定性和振动等特性。

使用 FEA 方法可以更好地辅助设计桥梁结构，提升桥梁的施工安全性和使用寿命。

Ansys 软件在桥梁计算中的应用方法Ansys 软件是工程领域中常用的通用有限元分析软件之一。

该软件具有强大的计算能力和众多的功能模块，适用于结构分析、热力学分析、流体力学分析等广泛领域的计算。

在桥梁计算中，Ansys 软件通常可以用于以下方面：•结构强度分析：Ansys 可以通过求解桥梁分析中的应力、应变和变形等参数，来评估桥梁结构的强度和稳定性。

•静态和动态分析：静态分析是针对荷载作用下改桥梁的响应，而动态分析则更加关注桥梁在振动状态下的响应。

•稳定性分析：Ansys 软件可以对桥梁结构在不同荷载情况下的稳定性问题进行分析。

•模态分析：模态分析用于确定桥梁在不同工作状态下的振动频率和模态形态。

Ansys 软件的桥梁计算模块通常包括几何建模、网格划分、材料参数设定、荷载设定、约束条件设定、分析求解等模块。

下面给出一些使用 Ansys 进行桥梁计算的简单步骤：1.准备桥梁几何模型，在 Ansys 中对桥梁进行建模；2.将桥梁模型划分为许多小的单元，确定每个单元的材料特性、截面特性等参数；3.设定荷载和约束条件，确定计算桥梁的工作状态；4.运行计算，解析桥梁在不同工况下的强度、稳定性和振动等特性。

ansys梁系计算

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

道路与桥梁电算(ANSYS作业)

课程名称：道路与桥梁电算（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。

曲线梁桥ANSYS计算命令流

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

曲线梁桥ANSYS计算命令流

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

某桥ANSYS计算

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的一个特点是笼统与通用，因此很多东西被掩盖到背后去了。

ANSYS四跨连续梁的内力计算教程

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。

ansys连续梁桥的影响线计算

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的简支板梁桥横向分布系数计算

基于ANSYS的简支板梁桥横向分布系数计算作者：李东明刘俭王飞关明仙来源：《价值工程》2012年第20期摘要：本文采用大型有限元计算软件对简支板梁桥的横向分布系数计算进行了探讨，并通过实例对此计算过程进行了说明，同时讨论了横向连接的刚度对横向分布系数的影响，计算结果表明，采用ANSYS软件计算横向分布系数具有一定的理论依据和工程实践性。

Abstract: This paper uses a large-scale finite element calculation software for calculating transverse distribution coefficient of plate girder bridge, and illustrates the calculation process by an example, and discusses the influence of rigidity of cross connection on transverse distribution coefficient. The calculation result shows that the use of ANSYS software for calculating transverse distribution coefficient has certain theoretical basis and the engineering practice.关键词： ANSYS；横向分布系数Key words: ANSYS；transverse distribution coefficient中图分类号：U44 文献标识码：A 文章编号：1006-4311（2012）20-0078-020 引言桥梁的计算内容很多，其中最重要的或者说最核心的就是桥梁的横向分布系数的计算或者说荷载在各主梁间如何分配的计算。

用ANSYS对钢筋混凝土梁进行计算模拟

一、用钢筋混凝土简支梁的数值模拟为实例，对ANSYS的使用方法进行说明钢筋混凝土简支梁，尺寸为长2000mm，宽150mm，高300mm。

混凝土采用C30，钢筋全部采用HRB335，跨中集中荷载P作用于一刚性垫板上，垫板尺寸为长150mm，宽100mm。

建立分离式有限元模型，混凝土采用SOLID65单元，钢筋采用LINK8单元，不考虑钢筋和混凝土之间的粘结滑移。

创建分离式模型时，将几何实体以钢筋位置切开，划分网格时将实体的边线定义为钢筋即可。

加载点以均布荷载近似代替钢垫板，支座处则采用线约束和点约束相结合。

单元尺寸以50mm左右为宜。

二、命令流!钢筋混凝土简支梁数值分析!分离式模型FINISH/CLEAR/PREP7!1.定义单元与材料属性ET,1,SOLID65,,,,,,,1ET,2,LINK8MP,EX,1,13585 !混凝土材料的初始弹模以及泊松比MP,PRXY,1,0.2FC=14.3 !混凝土单轴抗压强度和单轴抗拉强度FT=1.43TB,CONCR,1TBDA TA,,0.5,0.95,FT,-1 !定义混凝土材料及相关参数，关闭压碎TB,MISO,1,,11 !定义混凝土应力应变曲线，用MISO模型TBPT,,0.0002,FC*0.19TBPT,,0.0004,FC*0.36TBPT,,0.0006,FC*0.51TBPT,,0.0008,FC*0.64TBPT,,0.0010,FC*0.75TBPT,,0.0012,FC*0.84TBPT,,0.0014,FC*0.91TBPT,,0.0016,FC*0.96TBPT,,0.0018,FC*0.99TBPT,,0.002,FCTBPT,,0.0033,FC*0.85MP,EX,2,2.0E5 !钢筋材料的初始弹模以及泊松比MP,PRXY,2,0.3TB,BISO,2TBDA TA,,300,0 !钢筋的应力应变关系，用BISO模型PI=ACOS(-1)R,1R,2,0.25*PI*22*22R,3,0.25*PI*10*10TBPLOT,MISO,1 !混凝土材料的数据表绘图TBPLOT,BISO,2 !钢筋材料的数据表绘图!2.创建几何模型BLC4,,,150,300,2000*DO,I,1,19 !切出箍筋位置WPOFF,,,100VSBW,ALL*ENDDOWPCSYS,-1WPOFF,,,950 !切出拟加载面VSBW,ALLWPOFF,,,100VSBW,ALLWPCSYS,-1WPROTA,,-90WPOFF,,,30VSBW,ALLWPOFF,,,240VSBW,ALLWPCSYS,-1WPOFF,30WPROTA,,,90VSBW,ALLWPOFF,,,45VSBW,ALLWPOFF,,,45VSBW,ALLWPCSYS,-1!3.划分钢筋网格ELEMSIZ=50 !网格尺寸变量，设置为50mm LSEL,S,LOC,X,30LSEL,R,LOC,Y,30LA TT,2,2,2LESIZE,ALL,ELEMSIZLMESH,ALLLSEL,S,LOC,X,75 LSEL,R,LOC,Y,30LA TT,2,2,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,LOC,X,120 LSEL,R,LOC,Y,30LA TT,2,2,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,LOC,X,30 LSEL,R,LOC,Y,270LA TT,2,3,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,LOC,X,120 LSEL,R,LOC,Y,270LA TT,2,3,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,TAN1,Z LSEL,R,LOC,Y,30,270 LSEL,R,LOC,X,30,120 LSEL,U,LOC,X,75 LSEL,U,LOC,Z,0 LSEL,U,LOC,Z,2000 LSEL,U,LOC,Z,950 LSEL,U,LOC,Z,1050 LA TT,2,3,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,ALL!4.划分混凝土网格V A TT,1,1,1 MSHKEY,1ESIZE,ELEMSIZ VMESH,ALLALLSEL,ALL!5.施加荷载和约束LSEL,S,LOC,Y,0LSEL,R,LOC,Z,100DL,ALL,,UYLSEL,S,LOC,Y,0LSEL,R,LOC,Z,1900DL,ALL,,UYDK,KP(0,0,100),UX,,,,UZDK,KP(0,0,1900),UXP0=180000Q0=P0/150/100ASEL,S,LOC,Z,950,1050ASEL,R,LOC,Y,300SFA,ALL,1,PRES,Q0ALLSEL,ALLFINISH!6.求解控制设置与求解/SOLUANTYPE,0NSUBST,60OUTRES,ALL,ALLAUTOS,ONNEQIT,20CNVTOL,U,,0.015ALLSEL,ALLSOLVEFINISH!7.进入POST1查看结果/POST1SET,LASTPLDISP,1 !设定最后荷载步，查看变形ESEL,S,TYPE,,2ETABLE,SAXL,LS,1PLLS,SAXL,SAXL !绘制钢筋应力图ESEL,S,TYPE,,1/DEVICE,VECTOR,ONPLCRACK,1,1 !绘制裂缝和压碎图三、计算结果图1 混凝土材料的数据表绘图图2 钢筋材料的数据表绘图图3 钢筋的模拟图4 混凝土梁的模拟图5 梁在荷载作用下Y方向上的位移图图5 梁在荷载作用下Z方向上的应力图。

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计算例题

E,2002,601 E,601,2001 E,2001,801 E,400+L(1)+1,2003 E,2003,800+L(1)+1 E,400+L(1)+L(2)+1,2004 E,2004,800+L(1)+L(2)+1 E,400+L(1)+L(2)+L(3)+1,2005 E,2005,800+L(1)+L(2)+L(3)+1 E,400+L(1)+L(2)+L(3)+L(4)+1,2006 E,2006,800+L(1)+L(2)+L(3)+L(4)+1 E,401+LL,2008 E,2008,601+LL E,601+LL,2007 E,2007,801+LL !*********************************** *********** ! 生成无横隔板处梁格横梁单元 !*********************************** ************ TYPE,2 MAT,2 REAL,4 *DO,I,1,Nbeam-1,1 *DO,J,2,L(1),1 E,200*(I-1)+J,200*I+J *ENDDO *ENDDO *DO,I,1,Nbeam-1,1 *DO,J,L(1)+2,L(1)+L(2),1 E,200*(I-1)+J,200*I+J *ENDDO *ENDDO *DO,I,1,Nbeam-1,1 *DO,J,L(1)+L(2)+2,L(1)+L(2)+L(3)+2,1 E,200*(I-1)+J,200*I+J *ENDDO *ENDDO

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梁桥模拟计算..

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

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

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

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

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

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

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

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

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

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

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

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

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

桥梁计算（常用的计算方法)在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)。

基于ANSYS的桥梁计算模块开发

当对模型的工况全部进行设定之后，在系统目录下，用户可以观察到如图３所示的９个荷载步文件，依次对应于上述的９种基本工况。
图２公路－Ｉ级主菜单
图３荷载步文件
打开ＡｌｌＳｏｌｖｅ或ＰａｒｔＳｏｌｖｅ求解对话框，既可进行全部工况或部分工况的求解。
２抗震计算模块
期挠度预测探讨［Ｊ］．公路交通科技，２００７，２４（０１）：９１－９３．
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! （上接第１１４页）
ቤተ መጻሕፍቲ ባይዱ３结语
本文介绍了基于有限元软件ＡＮＳＹＳ二次开发的针对桥梁结构的计算模块，这些模块考虑了桥梁结构和规范的特点，既便是不熟悉ＡＮＳＹＳ的普通设计人员也可以很方便的利用这些求解模块进行桥梁结构的设计计算工作。鉴于时间和水平有限，这些模块的通用性和功能尚需进一步完善和扩展。
４结语
本文以将军河汉江大桥的施工监控为背景，介绍了
公路交通科技应用技术版
大跨度预应力混凝土连续刚构桥施工监控的基本概念和注意事项，及现今监控有待进一步讨论的问题，对做监控此类桥梁的监控有一定参考价值，将军河汉江大桥已于２００６年上半年通车，整桥线形良好，受到建设方的好评。
（６）对施工质量评定的一个误区是单纯强调合龙精度。虽然合龙精度是施工控制的重要目标，仅靠最后几个节段标高的调整而取得的高精度合龙必然会影响桥梁目标线形的实现，这也偏离了监控的内涵。除合龙精度外，成桥后还要考察实际线形和理想线形的偏差，比如说各节段节点和理想线形的偏差。由于刚成桥后，桥面的标高考虑了半活载和３年的收缩徐变，其与设计图纸上的标高有很大差别，此时的理想线形的评定有必要进一步研究。

用ANSYS软件计算桥梁结构的温度应力

第31卷　第3期2005年6月四川建筑科学研究Sichuan Building Science 收稿日期:2004205210作者简介:段　凯(1979-),男,湖北武汉人,工学硕士,主要从事桥梁结构及有限元数值计算方面的研究。

用ANS YS 软件计算桥梁结构的温度应力段　凯1,杨新华1,杨文兵1(华中科技大学土木工程与力学学院,湖北武汉　430074)摘　要:温度应力是混凝土桥梁开裂的主要因素,曾造成多座预应力混凝土桥梁结构严重损害,所以在进行桥梁结构设计时,必须考虑温度应力的影响。

本文在分析桥梁结构温度应力基本特点和ANSYS 软件特性的基础上,利用ANSYS 软件及其提供的二次开发工具开发了一个计算模块,实现了桥梁结构温度应力的求解。

关键词:温度应力;ANSYS 软件;桥梁结构中图分类号:TU311141 文献标识码:A 文章编号:1008-1933(2005)03-0055-051　概　述暴露在自然环境中的混凝土结构,受到周围环境气温以及日照等因素影响,外表面温度可能发生急变(升高或降低)。

由于混凝土材料的导热系数小(一般仅为黑色金属的几十分之一),混凝土内部的温度变化非常缓慢,从而产生明显的滞后现象,并且在混凝土结构内部形成较大的温度梯度。

当由此产生的温度变形被结构的内、外约束阻碍时,会产生相当大的温差应力。

在混凝土桥梁结构中,温度应力有时甚至比活载产生的应力还要大,不少预应力混凝土桥梁因此发生严重裂损。

随着大跨度预应力混凝土箱形桥梁结构的发展,温度应力对桥梁结构的影响和危害越来越大,因此在桥梁结构设计中,必须考虑温度应力的影响。

目前,桥梁结构温度应力的计算基本上采用力等效原理,由于过于简化,在处理复杂温度场时面临很大困难;计算使用的桥梁专用程序大多基于二维有限元理论,难以考虑梁的空间效应。

近十几年来,随着计算机技术的日益发展和有限元法的广泛应用,出现了一些大型通用的有限元分析程序,但是这些程序应用于桥梁结构温度应力计算有一定的局限性,操作过程复杂。

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：荷载-跨中挠度曲线与超筋梁的试验荷载-跨中挠度曲线很相似，在荷载达到极限情况下，没有出现屈服平台，而是突然跌落。