用ANSYS计算拱形橡胶减震垫的极限荷载

使用ANSYS模拟地震荷载的方法选用东营胜利油田CB11B平台的ANSYS模型对模态分析和动力分析中的操作方法进行介绍。

渤海CB11B平台是一座4腿导管架平台，包括上部甲板模块、导管架和桩基三部分。

甲板面标高为+9.00m，水深为10.5m。

桩腿的单向斜度10：1，入泥1.5m。

模拟地震荷载首先需要有地震的加速度数据，这里采用迁安波，迁安波为渤海的地震波，见文件eqq1.txt。

其时程图见图1。

注：该文件只有一列，即加速度值。

图1.加速度时程图第一步要把地震加速度数据输入ANSYS软件。

下拉菜单中Parameters-Array Parameters-Define/Edit-Add, 在Par中输入所定义数组名称（eqq）；输入数组选择Array；在I、H、K No.中输入数组的行数、列数、维数，所输入的行数应该与eqq1文件中的加速度数据个数相等，列数与维数在这里均为1。

下拉菜单中Parameters-Array Parameters-Read From File, 选择Array，点击OK；ParR中输入数组名称（eqq）；在File, ext, dir Read from file中浏览到地震加速度文件eqq1.txt所在的位置；Ncol Number of columns中输入1；最后一行中输入数据格式后点击OK（G10.4代表加速度数据总共十位，小数点后有四位.例如如：+1.2532，即G7.4）。

下拉菜单中Parameters-Array Parameters-Define/Edit-Add，选择数组文件名eqq后点击Edit，可以看到地震的加速度数据文件eqq1.txt已经被输入到数组eqq 中了。

点击Close关闭。

第二步要把地震加速度数据输入结构。

注意首先要把water table清空。

要以命令流的方式把地震加速度数据输入结构：FINISH/PREP7NT=500 %总计算步数DT=0.01 %时间步长，NT*DT即为总的计算时间/SOLUANTYPE,TRANS %以命令流的方式选择瞬态动力学分析TRNOPT,FULL*DO,I,1,NT,1 %循环开始TIME,I*DTKBC,0NSUB,1ALPHAD,0.1 %输入阻尼系数alphaBETAD,0.0028 %输入阻尼系数betadACEL,EQQ(I),0,0 %输入X、Y、Z向的地震加速度数据，这里只在X方向加了加速度。

基于ANSY S 中随动橡胶圈的有限元分析高喜玲(河南工程技术学校,河南焦作454000)摘要:常规型井口光杆密封装置的胶皮闸门是固定件,要实现井口密封装置的随动对中,就必须解除固定胶皮闸门的约束,设计出能够实时跟踪光杆运动的随动胶皮闸门组件,通过该部分组件的作用可使密封装置密封主体能够跟踪光杆的运动,调偏自如,随着杆的移动、振动,可在任意方向内进行动态补偿,随动体是重要部件,采用的是橡胶材料,由于操作频繁,橡胶内膜成为易损件,所以随动件的性能直接影响到该装置的寿命,一旦橡胶内膜损坏,不仅会造成停工,而且容易引发事故。

所以通过有限元分析来优选橡胶材料及结构参数,就显得非常必要。

关键词:随动;建模;应力分析;有限元中图分类号:TH11711　文献标志码:A 文章编号:100320794(2008)0520079203B ased on ANSY S with Rubber Ring of Finite E lement AnalysisG AO Xi -ling(Henan Engineering T echnical School ,Jiaozuo 454000,China )Abstract :Wellhead conventional type of rod sealing device rubber gate is fixed pieces ,in order to attain the wellhead sealing device with the right dynamic ,must lift gate fixed rubber restraints ,designed to track real -time m ovement of rod with the gates of rubber com ponents ,the com ponents through the role of sealing device will be the main seal tracking rod m ovement ,offset com fortable with the shot migration ,vibration ,in any direction within the dynamic com pensation ,with the dynamic is an im portant com ponent is used in rubber materials because users frequently ,rubber endometrial become fragile pieces.With the m oving parts of a direct im pact on the life of the device once rubber endometrial damage ,it will not only cause suspension ,and could easily trigger accidents.S o by finite element analysis and structural optimization rubber material parameters ,it is very necessary.K ey w ords :following ;m odeling ;stress analysis ;finite element0　前言由于抽油机制造、安装误差及抽油机受力变形、地基下沉等原因使光杆轴线与井口装置轴线间形成固定偏移及随机偏移,造成密封装置严重偏磨。

第23卷第3期海南大学学报自然科学版V ol.23N o.3 2005年9月NATURA L SCIENCE JOURNA L OF H AINAN UNIVERSIT Y Sep.2005 文章编号:1004-1729(2005)03-0238-04用ANSYS计算拱形橡胶减震垫的极限荷载李正中1,陈玉骥2(1.怀化铁路总公司张家界工务段,湖南张家界427001;2.中南大学铁道校区土建学院,湖南长沙410075)摘　要:通过实验确定了橡胶材料的M ooney2Rivlin常数,再以拱形橡胶减震垫为对象,建立了用于计算的空间有限元模型.考虑材料非线性和几何非线性,用工程软件ANSY S得到了拱形橡胶减震垫的极限荷载以及与之相对应的拱形橡胶减震垫的压缩量.关键词:拱形橡胶减震垫;超弹性材料;极限荷载中图分类号:O344.5 文献标识码:A在铁路线路上,钢轨和混凝土轨枕之间设置了减震垫,目的在于减少火车在运行时对混凝土轨枕面的冲击,改善行车的舒适度.目前,主要应用的减震垫有板式橡胶垫和圆柱式橡胶垫,其各有优缺点.板式橡胶减震垫的承载力较大,但是减震效果欠佳;而圆柱式橡胶减震垫,因为采用空心的构造,具有较好的减震效果,但由于构造上的原因,在圆柱与上部块体的结合处易产生应力集中现象,因而极限荷载偏低,疲劳寿命也较短.拱形橡胶减震垫是一种新型减震产品(见图1),与板式橡胶减震垫和圆柱式橡胶减震垫相比,拱形橡胶减震垫具有独特的优点.由于拱形橡胶减震垫采用半球形空腔的构造,减少了局部应力集中所带来的危害,从而提高了极限荷载和使用寿命,而且空心结构在工作时所产生的空气压力可以起到气体弹簧的作用,有助于减少振动的传递和冲击.可见,拱形橡胶减震垫是一种较好的减震装置.对于这种新型减震垫的力学性能研究,还未见公开报道.为了考察拱形橡胶减震垫的强度储备和变形情况,本文采用工程软件ANSY S对其极限荷载进行了分析,以下分4个方面进行讨论.1　拱形橡胶减震垫有限元模型的建立1.1　有关橡胶材料的超弹性理论　拱形橡胶垫的橡胶是主要的受力材料.橡胶不同于一般的弹性材料,它属于超弹性材料,即它的应力应变关系不是线性规律,也不同于应力应变关系曲线的切线斜率逐渐下降的塑性材料.而是随着应力的增长,应变的增长速度逐渐减小的超弹性材料.橡胶具有超弹性,在较小的外力作用下就能显示出高度的变形能力,由此决定了橡胶在变形过程中会呈现较强的几何非线性和材料非线性[1].对于橡胶超弹性性质的描述有不同的方式[2～4],其中,有从统计理论入手分析的,也有从分子之间的连接力入手分析的.随着计算机技术的发展以及有限元理论的广泛应用,一种被称之为唯象理论的概念逐渐引入了有关橡胶材料的理论计算中,此理论的基本思想是直接从橡胶材料的大变形现象出发,不考虑具体的分子结构,推导出可以完全或比较近似表达橡胶材料的应力应变关系的一种数学格式,从而同经典的弹性理论一样,不用考虑分子和结构等细节.这种方收稿日期:2004-11-26作者简介:李正中(1967-),男,四川广元人,怀化铁路总公司张家界工务段助理工程师.法由M ooney (1940)创始,而后Rivlin 将其推广到了一种更通用更合理的形式.图1　拱形橡胶减震垫的尺寸(单位:mm )Rivlin (1956)假设橡胶材料是不可压缩的,在小应变状态下是各向同性的,而各向同性不可压缩材料的应变能函数W 可以用下式表达W =∑∞i =0,j =0C ij (I 1-3)i (I 2-3)j ,(1)式中I 1和I 2为应变不变量.(1)式中每一项都是(I 1-3)和(I 2-3)的幂函数,此处之所以选用(I 1-3)和(I 2-3)而不用I 1和I 2是为了在零应变时(即I 1=I 3=3),W 将自动消失.1.2　拱形橡胶减震垫的材料常数的确定　M ooney 2Rivlin 常数是一种适用于不可压缩橡胶类材料的材料常数,可使用2、5、9等常数来描述.在有关手册中通常查不到超弹材料的M ooney 2Rivlin常数,但可通过对材料的标准试样的实验所得的实测数据来确定这些常数[5].为得到橡胶材料的物理特性,选取了部分样品做单向压缩实验,测得材料的应力应变关系,并将所测得的数据通过工程软件ANSY S 自带的超弹性分析程序进行材料特性的拟合计算[2],得出应变能函数所需要的常数,从而确定所用橡胶材料的材料属性.按上述办法得到的由M ooney 2Rivlin 模型所描述的应力应变关系如图2所示.　图2　在ANSY S 中定义的应力应变曲线 图3　1Π4拱形橡胶减震垫有限元模型1.3　拱形橡胶垫有限元建模　由于拱形橡胶减震垫具有对称性,故采用1Π4拱形橡胶减震垫建立有限元模型.建模时,采用H yper158单元,材料属性由以上实验所得的数据以及在ANSY S 中计算所得的M ooney 2Rivlin 超弹性模型确定.拱形橡胶减震垫的网格划分有限元模型如图3所示.实际的橡胶减震垫是支撑在混凝土轨枕上的,但为了简化计算,不考虑橡胶减震垫与轨枕的接触摩932　第3期 李正中等:用ANSY S 计算拱形橡胶减震垫的极限荷载擦,故在模型底面只约束竖向位移,而在对称面上则约束在对称面法线方向上的位移.2　荷载施加形式的确定作用于拱形橡胶减震垫的荷载有顶面铁轨所传递的列车轴重荷载以及半球内部的气体压力荷载.轴重荷载采用标准特种荷载,即火车的轴重为25t ,每个拱形橡胶减震垫上的荷载为25Π2=12.5t.半球内部的气体压力荷载是由顶面荷载的作用使半球体积发生变化所引起,它与半球体积的变化有关,故事先无法知道该气体压力的值.为了得到半球内气体的荷载,需进行必要的试算.为此,先做如下假设:1)在施加荷载的过程中,半球内气体不会漏掉;2)半球内的气体变化为等温变化;3)气体的压力随着轴重的增加而线性增加;4)半球内气体的初始压力为标准大气压.试算时,首先忽略半球内气体的存在,对橡胶减震垫施加顶面的列车轴重荷载,加载采用分级加载方式,每级加载10%(即1.25t ),求解过程考虑橡胶材料的材料非线性和由于大变形引起的几何非线性.然后根据橡胶减震垫在12.5t 荷载作用下半球的变形情况,计算半球在变形后的体积,通过与变形前的体积做比较,得到此时的半球内气体的压力.计算表明,所有半球的变化形状都很接近,几乎都是椭球形,故为了简化计算,假设所有半球变化后的尺寸都一样,并取平均尺寸进行计算.根据假设2计算出在顶面荷载12.5t 作用下,有如下关系:P 0V 0=P 1V 1,(2)式中P 0为标准大气压,取为1.01×105Pa ,P 1为变形后的半球内的气体压力,V 0、V 1分别为半球变形之前、后的体积.由于半球在变形后成为一个椭球形的空腔,所以V 0=43πR 3,V 1=43πab 2.(3)式中R 为半球在变形之前的半径;a ,b 分别为半球在变形之后椭球长半轴和短半轴的长度.由(2)、(3)式得P 1=V 0V 1P 0=R 3ab 2P 0,(4)通过上述公式,得到受压变形后半球内的气体压力P 1.在拱形橡胶减震垫未受到荷载时,整个橡胶减震垫边界(包括半球内)都受到一个标准大气压作用,该标准大气压为一平衡荷载,故这时的气压不计入荷载,即以后加载时所有气体压力荷载均应以标准大气压为起点开始计算,所以在变形后半球内的气体压力应该用以下气体压力作为荷载P =P 1-P 0.(5)由于该气体压力的数值是在没有考虑气体压力的情况下得到的,与实际情况的气体压力有所不同.因此,应该将气体压力按假定3施加于半球表面上,再进行求解计算,求得的结果再与上一次的计算结果比较,依此循环计算,直到相邻两次计算结果相差较小时方可认为这时的压力为顶面荷载12.5t 荷载时的气体压力.经过反复试算,最终得到在12.5t 荷载时的平均气体压力为0.35个大气压.3　拱形橡胶减震垫的极限荷载为简化计算,先引入以下假设:042海南大学学报自然科学版 2005年　1)计算时,只考虑到材料的超弹性性质,不考虑到其它的因素,如橡胶材料的疲劳以及其工作环境对材料的影响.图4　顶面位移随荷载变化2)半球内气体压力在加载过程中随顶面荷载同比例逐渐加大,到正常使用荷载时,假设气体压力达到最大值,然后随着荷载的增加,认为气体出现了泄漏现象,气体压力逐渐减小,当达到极限荷载时,气体正好全部漏完,即此时半球内无气体压力.采用上一节的有限元模型,考虑材料非线性,以及橡胶材料受力过程中产生的大变形,通过迭代计算,得到拱形橡胶垫的极限荷载.具体计算结果见图4.从以上计算结果可以看出,拱形橡胶减震垫顶面的位移(压缩量)随着荷载的增加而增加,在荷载均匀增加的情况下,压缩量的增加率是先快后慢,表现出超弹性材料的性质.该拱形橡胶减震垫的极限荷载为37.5t ,相对应的顶面位移为6.67m m.4　结　论1)拱形橡胶减震垫的橡胶材料属于超弹性材料,在计算过程中必须考虑材料的非线性特征,同时还要考虑超弹性材料的大变形能力;2)拱形橡胶减震垫顶面的位移随着荷载的增加逐渐增加,但增量则逐渐减小;3)该拱形橡胶减震垫的极限荷载为37.5t ,安全系数37.5Π12.5=3,顶面的极限位移大小为6.67mm ,约为其厚度的一半.参考文献:[1]史守峡.平面应力不可压缩橡胶薄片的非线性有限元分析[J ].哈尔滨工程大学学报,1998,19(3):11-15.[2]特雷劳尔L G R.橡胶弹性物理学[M].王梦蛟,等译.北京:化学工业出版社,1982.[3]弗雷克利P K,佩恩A R.橡胶在工程中应用的理论与实践[M].杜承泽,等译.北京:化学工业出版社,1985.[4]朱敏.橡胶化学与物理[M].北京:化学工业出版社,1984.[5]易太连,翁雪涛,朱石坚.不可压缩橡胶体的静态性能分析[J ].海军工程大学学报,2002,1(14):76-80.C alculation on the U ltim ate Load of the Arch Rubber Shock P adby E ngineering Softw are ANSYSLI Zheng 2zhong 1,CHE N Y u 2ji 2(1.Zhangjiajie Engineering Workgroup ,Huaihua Railway Parent 2company ,Zhangjiajie 427001,China ;2.Institute of Civil Engineering ,Railway Campus of Mid 2South University ,Changsha 410075,China )Abstract :This paper determined the M ooney 2Rivlin constants for rubber materials by test ,and then taking the arch rubber shock pad as an object ,established a spatial finite element m odel used for calculation.When their material non 2linearity and geometrical non 2linearity were taken into account ,the ultimate load and the corresponding com pressive load of the arch rubber shock pad could be obtained by engineering s oft 2ware ANSY S.K ey w ords :arch rubber shock pad ;super 2elastic materials ;ultimate load 142　第3期 李正中等:用ANSY S 计算拱形橡胶减震垫的极限荷载。

02地震分析算例(ANSYS)土木工程中除了常见的静力分析以外，动力分析，特别是结构在地震荷载作用下的受力分析，也是土木工程中经常遇到的问题。

结构的地震分析根据现行抗震规范要求，一般分为以下两类：基于结构自振特性的地震反应谱分析和基于特定地震波的地震时程分析。

本算例将以一个4质点的弹簧－质点体系来说明如何使用有限元软件进行地震分析。

更复杂结构的分析其基本过程也与之类似。

关键知识点：(a)模态分析(b)谱分析(c)地震反应谱输入(d)地震时程输入(e)时程动力分析(1)在ANSYS窗口顶部静态菜单，进入Parameters菜单，选择Scalar Parameters选项，在输入窗口中填入DAMPRATIO=0.02，即所有振型的阻尼比为2%(2)ANSYS主菜单Preprocessor->Element type->Add/Edit/Delete，添加Beam 188单元(3)在Element Types窗口中，选择Beam 188单元，选择Options，进入Beam 188的选项窗口，将第7个和第8个选项，Stress/Strain (Sect Points) K7, Stress/Strain (Sect Nods) K8，从None改为Max and Min Only。

即要求Beam 188单元输出积分点和节点上的最大、最小应力和应变(4)在Element Types窗口中，继续添加Mass 21集中质量单元(5)下面输入材料参数，进入ANSYS主菜单Preprocessor->Material Props-> MaterialModels菜单，在Material Model Number 1中添加Structural-> Linear-> Elastic->Isotropic属性，输入材料的弹性模量EX和泊松比PRXY分别为210E9和0.3。

(6)继续给Material Model Number 1添加Density属性，输入密度为7800。

2.1 载荷概述有限元分析的主要目的是检查结构或构件对一定载荷条件的响应。

因此，在分析中指定合适的载荷条件是关键的一步。

在ANSYS程序中，可以用各种方式对模型加载，而且借助于载荷步选项，可以控制在求解中载荷如何使用。

2.2 什么是载荷在ANSYS术语中，载荷（loads）包括边界条件和外部或内部作用力函数，如图2-1所示。

不同学科中的载荷实例为：结构分析：位移，力，压力，温度（热应变），重力热分析：温度，热流速率，对流，内部热生成，无限表面磁场分析：磁势，磁通量，磁场段，源流密度，无限表面电场分析：电势（电压），电流，电荷，电荷密度，无限表面流体分析：速度，压力图2-1 “载荷”包括边界条件以及其它类型的载荷载荷分为六类：DOF约束，力（集中载荷），表面载荷，体积载荷、惯性力及耦合场载荷。

·DOF constraint（DOF约束）将用一已知值给定某个自由度。

例如，在结构分析中约束被指定为位移和对称边界条件；在热力分析中指定为温度和热通量平行的边界条件。

·Force（力）为施加于模型节点的集中载荷。

例如，在结构分析中被指定为力和力矩；在热力分析中为热流速率；在磁场分析中为电流段。

·Surface load（表面载荷）为施加于某个表面上的分布载荷。

例如，在结构分析中为压力；在热力分析中为对流和热通量。

·Body load（体积载荷）为体积的或场载荷。

例如，在结构分析中为温度和fluences；在热力分析中为热生成速率；在磁场分析中为流密度。

·Inertia loads（惯性载荷）由物体惯性引起的载荷，如重力加速度，角速度和角加速度。

主要在结构分析中使用。

·Coupled-field loads（耦合场载荷）为以上载荷的一种特殊情况，从一种分析得到的结果用作为另一分析的载荷。

例如，可施加磁场分析中计算出的磁力作为结构分析中的力载荷。

其它与载荷有关的术语的定义在下文中出现。

基于ANSYS的橡胶减震垫的设计选型李霞，周奖清（英格索兰亚太工程技术中心，上海200051）摘要：介绍了螺杆式空压机驱动机构减震垫设计选型的优化。

通过ANSYS的辅助分析，预测减震垫的最大变形，寻找最优的安装位置和最合适刚度K值减震垫，同时提到的该方法也适用于一般的减震垫的设计和选型优化。

关键词：减震垫；支座反力；重心；刚度K值；压缩量中图分类号:TB535.1文献标志码:A文章编号：1002-2333（2020）01-0116-03 Design and Selection of Rubber Cushion Based on ANSYSLI Xia,ZHOU Jiangqing(Ingersoll Rand Engineering Technology Center-Asia Pacific,Shanghai200051,China)Abstract:This paper mainly introduces the optimization of design and selection of the cushion for the drive mechanism of screw air compressor.Through the auxiliary analysis of ANSYS,the maximum deformation of the cushion is predicted,and the optimal installation position and the cushion with the most suitable stiffness K value are found.This method is also suitable for the design and selection optimization of general cushions.Keywords:cushion;support reaction force;center of gravity;stiffness K value;compression0引言螺杆式空压机在运行时，压缩机主机和主电动机本身会发生振动，传至连接的底架、机组罩壳，引起整个系统的振动和噪声。

采用ANS Y S的橡胶弹簧的有限元建模与仿真任茂文1,周长峰2(1江苏电大泗洪分校,泗洪223900;2东南大学机械工程学院,南京210096)摘要:建立AD250铰接式自卸车前悬架橡胶弹簧的参数化非线性有限元接触模型,比较橡胶弹簧实体建模与平面建模的计算精度,并分析大载荷下橡胶弹簧内部应力的分布,得到的结果与试验结果有很好的一致性。

此模型可为橡胶弹簧结构的参数优化提供理论支持。

关键词:橡胶弹簧;非线性;有限元分析中图分类号:U463133+415　文献标识码:B　文章编号:1671—3133(2008)05—0059—04Rubber spr i n gπs f i n ite elem en t si m ul a ti on m odeli n gand si m ul a ti on ba sed on ANS Y SRen Mao2wen1,Zhou Chang2feng2(1School of Sihong,J iangsu T V University,Sihong223900,J iangsu,CHN;2School ofMechanicalEngineering,Southeast University,Nanjing210096,CHN)Abstract:The AD250hinge type du mp truck fr ont sus pensi on rubber s p ringπs para metrizati on non2linear finite ele ment contact model has been established.The computati onal accuracies of the rubber s p ring entity modelling and the p lain modelling have been compared.Besides,the rubber s p ring internal stress distributi on under the big l oad has been analyzed.There are good unif or m ity bet w een the result and the experi m ental result.This model may p r ovide the theoretical support for the para meter op ti m izati on of the rubber s p ring structure.Key words:Rubber s p ring;Non2linear;Finite ele ment analysis 近年来,具有可变刚度特点的橡胶弹簧和空气弹簧在工程车辆的悬架系统中已逐渐得到应用。

如何在ANSYS中模拟橡胶隔震⽀座在ANSYS中，没有现成的隔震⽀座单元。

然⽽隔震⽀座的⼒学模型可以简化为由⽔平两⽅向的⾮线
性弹簧、粘滞阻尼器以及竖向的线性弹簧所组成，并且ANSYS提供了⼀系列的弹簧阻尼等连接单元。

这样，我们可以⽤若⼲单元相组合的⽅式来实现隔震⽀座的模拟。

竖向刚度的模拟采⽤combin14单元。

在两个⽔平⽅向采⽤combin40单元，该单元可以引⼊双线性的强化模型、粘滞阻尼的影响。

combin40⼒学原理图
隔震垫单元图
隔震垫⽔平⽅向基本参数主要有：Ku（屈服前刚度）、Kd（屈服后刚度）、Qd（屈服⼒）、阻尼⽐。

⽽combin40单元的实常数主要有：K1，K2，C，FSLIDE，GAP。

由combin40⼒学原理图，可以得到这些实常数的选取⽅法：K2=Kd，K1=Ku-Kd，FSLIDE=Qd，GAP=0，C为阻尼系数。

这样，⼀个隔震⽀座由三个单元所组成，combin40（X⽅向）、combin40（Y⽅
向）、combin14（Z⽅向），三单元不相交的各⾃节点约束所有⾃由度，交点处节点约束转动⾃由度（基底隔震垫），若设为层间隔震垫（如照⽚中的⼀层顶隔震），则耦合combin14两端的转动⾃由度（其中⼀端的节点与下部结构共享）。

其中：GAP－间隙⼤⼩；假如GAP等于零，接触⾯不能打开。

假如GAP是负值，有初始冲突。

建⽴⽀座上下两个节点，然后在这两个节点分间别建⽴三个弹簧单元（两个combin14与⼀
个combin40），弹簧的⽅向通过改变弹簧单元keyopt来实现，combin14修改keyopt(2)，combin40修
改keyopt(3)来指定弹簧⽅向。

在任意面施加任意方向任意变化的压力在某些特殊的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷，当然，这在一定程度上可以通过ANSYS表面效应单元实现。

如果利用ANSYS的参数化设计语言，也可以非常完美地实现此功能，下面通过一个小例子描述此方法。

!!!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来!!!finish/prep7et,500,shell63press=100e6amesh, allesla, snsla,s,1! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将! 所有节点坐标系旋转到局部坐标系下.*get,enmax,elem,,num,maxdofsel,s,fx,fy,fzfcum,add !!!将力的施加方式设置为"累加",而不是缺省的"替代"*do,i,1,enmax*if,esel,eq,1,then*get,ae,elem,i,area !此命令用单元真实面积，如用投影面积，请用下几条命令! *get,ae,elem,i,aproj,x !此命令用单元X投影面积，如用真实面积，请用上一条命令! *get,ae,elem,i,aproj,y !此命令用单元Y投影面积! *get,ae,elem,i,aproj,z !此命令用单元Z投影面积xe=centrx !单元中心X坐标(用于求解压力值)ye=centry !单元中心Y坐标(用于求解压力值)ze=centrz !单元中心Z坐标(用于求解压力值)! 下面输入压力随坐标变化的公式,本例的压力随X和Y坐标线性变化.p_e=(xe‐10)*press+(ye‐5)*pressf_tot=p_e*aeesel,s,elem,,insle,s,corner*get,nn,node,,countf_n=f_tot/nn*do,j,1,nnf,nelem(i,j),fx,f_n !压力的作用方向为X方向! f,nelem(i,j),fy,f_n !压力的作用方向为Y方向! f,nelem(i,j),fz,f_n !压力的作用方向为Z方向*enddo*endifesla,s*enddoaclear,allfcum,repl !!!将力的施加方式还原为缺省的"替代"dofsel,allallsel说明：本信息在任意面施加任意方向任意变化的压力在某些特殊的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷，当然，这在一定程度上可以通过ANSYS表面效应单元实现。

1. 荷载步中荷载的处理方式无论是线性分析或非线性分析处理方式是一样的。

①对施加在几何模型上的荷载(如 fk,sfa 等)：到当前荷载步所保留的荷载都有效。

如果前面荷载步某个自由度处有荷载，而本步又在此自由度处施加了荷载，则后面的替代前面的；如果不是在同一自由度处施加的荷载，则施加的所有荷载都在本步有效(删除除外！)。

②对施加在有限元模型上的荷载(如 f,sf,sfe,sfbeam 等)：ansys缺省的荷载处理是替代方式，可用 fcum,sfcum 命令修改，可选择三种方式：替代(repl)、累加(add)、忽略(igno)。

当采用缺省时，对于同一自由度处的荷载，后面施加的荷载替代了前面施加的荷载(或覆盖)；而对于不是同一自由度的荷载(包括集中或分布荷载)，前面的和本步的都有效。

当采用累加方式时，施加的所有荷载都在本步有效。

特别注意的是，fcum 只对在有限元模型上施加的荷载有效。

2. 线性分析的荷载步从荷载步文件(file.snn)中可以看到，本步的约束条件和荷载情况，而其处理与上述是相同的。

由于线性分析叠加原理是成立的，或者讲每步计算是以结构的初始构形为基础的，因此似乎可有两种理解。

①每个荷载步都是独立的：你可以根据你本步的约束和荷载直接求解(荷载步是可以任意求解的，例如可以直接求解第二个荷载步，而不理睬第一个荷载步:lssolve,2,2,1)，其结构对应的是你的约束和荷载情况，与前后荷载步均无关！（事实上，你本步可能施加了一点荷载，而前步的荷载继续有效，形成你本步的荷载情况）②后续荷载步是在前步的基础上计算的(形式上！)。

以荷载的施加先后出发，由于本步没有删除前面荷载步的荷载，你在本步仅仅施加了一部分荷载, 而结构效应是前后荷载共同作用的结果。

不管你怎样理解，但计算结果是一样的。

(Ansys是怎样求解的，得不到证实。

是每次对每个荷载步进行求解，即[K]不变，而[P]是变化的，且[P]对应该荷载步的所有荷载向量呢？或是[P]对应一个增量呢？不用去管他，反正结果一样) 也有先生问，想在第 N 步的位移和应力的基础上，施加第 N+1 步的荷载，如何？对线性分析是没有必要的，一是线性分析的效应是可以叠加的，二是变形很小(变形大时不能采用线性分析)。

如何使用ANSYS 建模进行桥梁荷载评估摘要本文主要介绍利用ANSYS建立桥梁的结构模型，并计算出恒载、活载作用下所产生的内力及应力。

其中活载采用公路－Ⅰ级和公路－Ⅱ级的荷载标准对桥梁进行加载，并将活载作用下杆件所产生的内力及应力进行组合，将其组合值与实测应力幅值相比较，以此来作为评估荷载等级的依据。

随后对实桥的交通进行观测，在观测的基础上统计出车辆的类型、车重、轴距和过往频率，再根据统计的结果最终建立车辆模拟荷载。

关键词：钢桁架、有限元模型、车辆荷载模拟This paper mainly introduces the structure of the bridge by ANSYS established model, and calculates the constant load, under the influence of live load and the internal force produced and stress. Live load of the highway-Ⅰlevel and highway-Ⅱlevel to bridge for the standard load load, and will be under the influence of live load and stem a generated internal force and stress in combination, will the combination with measured values stress amplitude compared, as the basis for evaluating load level. Then the real bridge traffic for observation, the basis of the observation statistics a vehicle type, weight of the car, wheelbase and past frequency, again according to statistics results vehicle load establish simulation.Keywords: steel truss, finite element model, the vehicle load simulation 中图分类号：F407.9 文献标识码：A 文章编号1 概述近年来，随着我国交通运输事业的不断发展，大量低等级公路被改建扩建，服役桥梁能否继续使用已成为公路建设决策部门的一件大事。

在ANSYS里做地震分析时，需要对结构施加地震惯性荷载，地震惯性力是通过加速度的方式输入进结构的，然后与结构的质量一起形成动力计算时的惯性荷载，下面说一下在ANSYS 里施加地震惯性力的方法。

首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00 .......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流文件，并命名为acce.inp*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行，同上ACCEXYZ(0,3)=3 !08行，同上finish/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。

!荷载——结构方法计算(马蹄形断面）finish !退出当前处理程序/clear !清除以前数据,重新开始一个新的分析/COM,Structural !定义分析类型,结构分析(热分析、流体分析等)/prep7 !进入前处理器*AFUN,deg !定义角度单位为度（缺省为弧度，RAD）! 定义建模及材料参数的一些变量值*set,Py,0.13492e6 !定义垂直围岩压力大小（若有地表荷载加地表荷载值）*set,px1,0.03373e6*set,px2,0.03373e6*set,CQHD,0.45 !定义初支或二衬厚度*set,CQDYCD,0.6 !定义梁单元长度参数及弹簧单元面积（梁单元长度与弹簧单元面积相等）*set,CQETXML,31.5e9 !定义衬砌（初支或二衬）的弹性模量*set,CQUBSB,0.2 !定义衬砌的泊松比*set,WYTXKL,350e6 !定义围岩的弹性抗力系数*set,CQMD,2500 !定义衬砌的密度!定义单元类型及材料属性及单元实常数et,1,beam3 !定义1号单元为梁单元mp,ex,1,CQETXML !定义1号材料的弹性模量mp,prxy,1,CQUBSB !定义1号材料的泊松比mp,dens,1,CQMD !定义1号材料的密度R,1,CQHD,1/12*CQHD*CQHD*CQHD,CQHD !1-实常数号；第一个参数为梁截面的面积；第二个参数为梁单元的转动惯量；第三个参数为梁高。

!建立几何模型!创建关键点k,1,0,-4.42439, 0 !k-点命令；1-点号；坐标（x,y,z)k,2,2.49468,-4.19466,k,3,4.4,-3.69635,k,4,5.50615,-2.98047,k,5,6.09977,-1.97,k,6,6.41,0,k,7,0,6.41,k,8,-6.41,0,k,9,-6.09977,-1.97,k,10,-5.50615,-2.98047,k,11,-4.4,-3.69635,k,12,-2.49468,-4.19466,!二衬轴线larc,1,3,2 !三点画圆弧：起点-终点—中间点最好全按逆时针画，弯矩及轴力的方向一致larc,3,5,4larc,5,7,6larc,7,9,8larc,9,11,10larc,11,1,12!设置线单元材料属性，划分单元(初支单元)lsel,s,,,1,6 !LSEL,Type,Item,Comp,VMIN, VMAX,VINC（选择1-8单元）lATT,1,1,1 !给线单元付材料号、实常数（1只是代码，对应上面的函数值）、单元类型号LATT, MAT, REAL, TYPE, ESYSMSHKEY,1 !设置单元为映射单元（单元比较规则）lesize,1,CQDYCD !LESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE（划分线的长度，以衬砌单元厚度为lesize,2,CQDYCDlesize,3,CQDYCDlesize,4,CQDYCDlesize,5,CQDYCDlesize,6,CQDYCDlsel,alllmesh,all !划分线单元nummrg,all,,,,low !所有号数从1开始numcmp,node,eord !压缩节点号并排序!施加水平弹簧支撑单元（在梁单元的每个节点上分别产生水平方向弹簧单元）*do,i,2,22 !do循环，要对应下面一个*EDNDOa1=abs(ny(i-1)-ny(i+1))*WYTXKL !计算节点i-1和节点i+1之间y的变化量*弹性抗力kx=a1/2PSPRNG,i,TRAN,kx,1,,,, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM!弹簧单元,节点号,平移,弹性系数,弹簧水平长度产生水平弹簧*enddo*do,i,48,67,1 !do循环，要对应下面一个*EDNDOa1=abs(ny(i-1)-ny(i+1))*WYTXKL !计算节点i-1和节点i+1之间y的变化量*弹性抗力kx=a1/2PSPRNG,i,TRAN,kx,-1,,,, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM!弹簧单元,节点号,平移,弹性系数,弹簧水平长度产生水平弹簧*enddoa1=abs(ny(67)-ny(1))*WYTXKL !计算节点64上的弹性抗力kx=a1/2PSPRNG,68,TRAN,kx,-1,,,,!施加竖直弹簧支撑单元（在梁单元的每个节点上分别产生竖直方向弹簧单元）b1=abs(nx(68)-nx(2))*WYTXKL !计算节点1上的弹性抗力ky=b1/2PSPRNG,1,TRAN,ky,,-1,,,*do,i,2,17 !do循环，要对应下面一个*EDNDOb1=abs(nx(i-1)-nx(i+1))*WYTXKL !计算节点i-1和节点i+1之间x的变化量*弹性抗力ky=b1/2PSPRNG,i,TRAN,ky,,-1,,, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM!弹簧单元,节点号,平移,弹性系数,弹簧水平长度产生水平弹簧*enddo*do,i,53,67 !do循环，要对应下面一个*EDNDOb1=abs(nx(i-1)-nx(i+1))*WYTXKL !计算节点i-1和节点i+1之间x的变化量*弹性抗力ky=b1/2PSPRNG,i,TRAN,ky,,-1,,, !PSPRNG, NLOC, TYPE, K, DX, DY, DZ, ELEM!弹簧单元,节点号,平移,弹性系数,弹簧水平长度产生水平弹簧*enddob1=abs(nx(67)-nx(1))*WYTXKL !计算节点1上的弹性抗力ky=b1/2PSPRNG,68,TRAN,ky,,-1,,,!施加边界条件与围岩压力/solu !!!!!进入求解器计算!nsel,s,,,1 !选择约束的节点,仰拱中间节点!d,all,ux !施加水平方向的约束allselfcum,add,, !一定要，使荷载能叠加!竖向荷载（y方向压力，最大跨度以上节点施加竖向荷载，换算成节点荷载，在节点上施加）*do,i,18,52 !在18—48号节点施加竖向围岩压力a=-Py*0.5*abs(nx(i-1)-nx(i+1)) !计算作用在节点上的节点力f,i,fy,a*enddo!水平向荷载（X方向压力，全部节点施加水平荷载，换算成节点荷载，在节点上施加）*do,i,2,34 !在2—32号节点施加水平围岩压力a=-(px1+(px2-px1)/(ny(35)-ny(1))*(ny(35)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddo*do,i,36,67 !在36—67号节点施加水平围岩压力a=(px1+(px2-px1)/(ny(35)-ny(1))*(ny(35)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddo!设定重力加速度ACEL,0,9.8,0, !重力加速度以向下位置，不用加负号ALLSEL,ALLfinish!!!!!进入求解器计算/soluESEL,ALLNSEL,ALLsolve !求解!进入后处理器，对结果进行后处理/POST1ESEL,s,TYPE,,1 !选择单元1（隧道轴力中心线单元）NSLE,s,ALL !选择所有节点etabel,Ni,smisc,1 !定义梁i节点的轴力etabel,Nj,smisc,7 !定义梁j节点的轴力etabel,Mi,smisc,6 !定义梁i节点的弯矩etabel,Mj,smisc,12 !定义梁j节点的弯矩etable,Qi,smisc,2 !定义梁i节点的剪力etable,Qj,smisc,8 !定义梁j节点的剪力plls,ni,nj,1 !查看轴力图plls,mi,mj,-1 !查看弯矩图。