《桥梁课程设计》

《《桥梁课程设计》》由会员分享，可在线阅读，更多相关《《桥梁课程设计》（14页珍藏版）》请在装配图网上搜索。

1、课程名称：桥梁工程B设计题目ansys梁桥模拟计算（三跨）院系：土木工程系专业：年级：姓名：学号：指导教师：西南交通大学峨眉校区ansys 梁桥模拟计算（三跨）1.绪论1.1 设计目的桥梁结构分析计算是非常重要的一门技能。通过本课程设计，掌 握一门通用有限元软件分析工具，能够独立对桥梁结构进行静力或动 力分析。本课程具体要求掌握通用有限元软件ANSYS, 了解其前处 理，后处理过程以及单元应用。通过此课程设计的学习，初步具有独 立进行结构分析的能力，从而了解桥梁的具体设计。1.2设计内容及要求桥梁结构建模、确定边界条件、求解、后处理以及分析结论1、了解所选用 Beam4 等单元的属性和用法；2

2、、对桥梁进行结构离散化，建立三维有限元数值模型；3、正确地对桥梁有限元模型设定边界条件；4、掌握数值分析静力或动力求解方法；5、对计算结果进行后处理，掌握基本作图软件应用；6、对计算结果进行分析，得出结论。2.有限元分析2.1简介有限元分析是用较简单的问题代替复杂问题后再求解。它将求解 域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一 个合适的（较简单的）近似解，然后推导求解这个域总的满足条件（如 结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近 似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难 以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，

3、 因而成为行之有效的工程分析手段。有限元是那些集合在一起能够表示实际连续域的离散单元。有限 元的概念早在几个世纪前就已产生并得到了应用，例如用多边形（有 限个直线单元）逼近圆来求得圆的周长，但作为一种方法而被提出， 则是最近的事。有限元法最初被称为矩阵近似方法，应用于航空器的 结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研 究的科学家的浓厚兴趣。经过短短数十年的努力，随着计算机技术的 快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几 乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效 的数值分析方法。在解偏微分方程的过程中, 主要的难点是如何构造一个方程来

4、逼近原本研究的方程, 并且该过程还需要保持数值稳定性.目前有许 多处理的方法， 他们各有利弊. 当区域改变时（就像一个边界可变的 固体）, 当需要的精确度在整个区域上变化, 或者当解缺少光滑性时, 有限元方法是在复杂区域（像汽车和输油管道）上解偏微分方程的一 个很好的选择. 例如, 在正面碰撞仿真时, 有可能在重要区域（例 如汽车的前部）增加预先设定的精确度并在车辆的末尾减少精度（如 此可以减少仿真所需消耗）; 另一个例子是模拟地球的气候模式, 预 先设定陆地部分的精确度高于广阔海洋部分的精确度是非常重要的2.2 方法步骤对于不同物理性质和数学模型的问题，有限元求解法的基本步骤 是相同的，只是

5、具体公式推导和运算求解不同。有限元求解问题的基 本步骤通常为：第一步：问题及求解域定义：根据实际问题近似确定求解域的物 理性质和几何区域。第二步：求解域离散化：将求解域近似为具有不同有限大小和形 状且彼此相连的有限个单元组成的离散域，习惯上称为有限元网络划 分。显然单元越小（网格越细）则离散域的近似程度越好，计算结果 也越精确，但计算量及误差都将增大，因此求解域的离散化是有限元 法的核心技术之一。第三步：确定状态变量及控制方法：一个具体的物理问题通常可 以用一组包含问题状态变量边界条件的微分方程式表示，为适合有限 元求解，通常将微分方程化为等价的泛函形式。第四步：单元推导：对单元构造一个适合的

6、近似解，即推导有限 单元的列式，其中包括选择合理的单元坐标系，建立单元试函数，以 某种方法给出单元各状态变量的离散关系，从而形成单元矩阵（结构 力学中称刚度阵或柔度阵）。为保证问题求解的收敛性，单元推导有许多原则要遵循。 对工 程应用而言，重要的是应注意每一种单元的解题性能与约束。例如， 单元形状应以规则为好，畸形时不仅精度低，而且有缺秩的危险，将 导致无法求解。第五步：总装求解：将单元总装形成离散域的总矩阵方程（联合 方程组），反映对近似求解域的离散域的要求，即单元函数的连续性 要满足一定的连续条件。总装是在相邻单元结点进行，状态变量及其 导数（可能的话）连续性建立在结点处。第六步：联立方程

7、组求解和结果解释：有限元法最终导致联立方 程组。联立方程组的求解可用直接法、迭代法和随机法。求解结果是 单元结点处状态变量的近似值。对于计算结果的质量，将通过与设计 准则提供的允许值比较来评价并确定是否需要重复计算。简言之，有限元分析可分成三个阶段，前置处理、计算求解和后 置处理。前置处理是建立有限元模型，完成单元网格划分；后置处理 则是采集处理分析结果，使用户能简便提取信息，了解计算结果。3.ANSYS 模拟3.1ANSYS 历史ANSYS 公司成立于1970年，总部位于美国宾西法尼亚的匹兹堡 致力于CAE技术的研究和发展，并通过CAE技术帮助企业优化设计流 程ANSYS软件的创始人是美国匹

8、斯堡大学力学系教授、著名有限元 权威John Swanson博士。ANSYS灵活、开放的解决方案为概念设计 到最终测试的设计全过程提供了有效的CAE协同环境，使客户可以在 设计的各个阶段大规模采用CAE技术，最大程度地发挥CAE对设计流 程的贡献，从而大幅度地缩短研发流程、降低研发费用、提高设计质 量。3.2ANSYS 软件介绍ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大 型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的 美国ANSYS开发。它能与多数CAD软件接口，实现数据的共享和交 换，如 Pro/Engineer, NASTRAN, Alogor, IDEAS

9、, Auto CAD 等，是 现代产品设计中的高级CAE工具之一。软件主要包括三个部分：前处理模块，分析计算模块和后处理模 块。前处理模块:它提供了一个强大的实体建模及网格划分工具，用 户可以方便地构造有限元模型；ANSYS软件提供的分析类型如下：1结构静力分析2结构动力分析3结构非线性分析4动力学分析5热分析6电磁场分析7流体动力学分析8声场分析9压电分析后处理模块：它可将计算结果以彩色等值线显示、梯度显示、矢 量显示、粒子流迹显示、立体切片显示、透明及半透明显示（可看到 结构内部）等图形方式显示出来，也可将计算结果以图表、曲线形式 显示或输出。软件提供了100种以上的单元类型，用来模拟工程

10、中的各种结构 和材料。该软件有多种不同版本，可以运行在从个人机到大型机的多 种计算机设备上，如 PC，SGI，HP，SUN，DEC，IBM，CRAY 等。3.3结构分析流程图3.4 命令介绍3.5单元的介绍SHELL43 单元描述：SHELL43适合模拟线性、弯曲及适当厚度的壳体结构。单元中每 个节点具有六个自由度：沿x、y和z方向的平动自由度以及绕x、y 和z轴的转动自由度。平面内两个方向的形状必变都是线性的。对 于平面外的运动，用张量组的混合内插法(a mixed interpolation of tensorial components)。单元具有塑性、蠕变、应力刚化、大变形和大应变的特

11、性。关于 此单元更详细的性能可参见ANSYS, Inc. Theory Reference中的 SHELL43。如果是薄壳或者塑性和蠕变不需考虑，弹性的四边性壳单 元(SHELL63)就可以了。如果遇到收敛困难或者需要考虑大应变时， 可选择SHELL181单元。当然，对于非线性结构分析我们推荐选择 SHELL181 单元。图43.1 SHEll43单元几何图示其中：XIJ =没有定义单元坐标系时的X轴乂 =定义了单元坐标系时的X轴SHELL43 单元的输入数据：图43-1给出了此单元的几何形状、节点位置和坐标系设置。单 元由四个节点、四个壳厚度以及正交各向异性的材料特性确定。在三 角形、棱形和

12、四面体单元一章中已经提到，当把节点K和节点L定 义为同一个节点时就形成了三角形单元。正交各向异性材料的方向与单元坐标系的方向一致。单元坐标系 的方向已经在坐标系一章中描述过。单元X轴可以从X轴向Y轴旋 转一个角度THETA。单元需要有一个有效的厚度。随着在每个角节点处输入的厚度值 的不同，假定厚度在单元面积上平滑变化。如果单元厚度不变，只输 入TK (I)就可以了。如果厚度不是常数，必须分别输入四个节点的 厚度值。单元名义上绕Z轴的平面内的转动刚度由KEY0PT(3) (=0或1) 确定。另外一个真实的转动刚度(Allman转动)相应在由KEYOPT(3) =2来定义。这样的话，实常数ZSTI

13、F1和ZSTIF2就被用来控制 Allman转动理论中的两个伪零能量模态。ZSTIF1和ZSTIF2的默认 值分别为1.0E-6和1.0E-3。ADMSUA是每单位面积的质量。单元荷载在节点和单元荷载一章中已经讲过。压力被当做单元面 上的表面荷载输入，如图43-1的圆圈数字所示。正向表示指向单元 表面。边界压力按单位长度的力来输入。温度可当成作用在角点位置 （1-8）的单元体力输入，如图43-1所示。第一个角温度T1默认值 是TUNIF,如果其它的温度都不指定，它们的值也自动按T1取值。 如果只输入了 T1和T2，T1值会赋给T1、T2、T3和T4,而输入的 T2值会赋给T5、T6、T7和T8

14、。对于任何其它的输入方案，未指定 温度时默认都是TUNIF。SHELL43输入总结.给出了此单元的输入参数总结。一般的单元 输入描述请参考单元输入一章。xrj = Element x-axis if ESYS is not supplied.x = Element x-axis if ESYS is supplied.1.1.1. BEAM188BEAM188适用于分析细长的梁。元素是基于Timoshenko梁理 论的。具有扭切变形效果。BEAM188是一个二节点的三维线性梁。BEAM188在每个节点上 有6或7个自由度，（自由度）数目的变化是由KEY0PT（ 1）来控制的。 当KEYOPT（

15、1） =0时（默认），每节点有6个自由度。分别是沿x,y,z 的位移及绕其的转动。当KEYOPT（1） = 1时,会添加第七个自由度 （翘曲量）。此元素能很好的应用于线性（分析），大偏转，大应力的非线性 （分析）。BEAM188包含应力刚度，在默认情况下，在某些分析中由NLGEOM 来打开。在进行弯曲（flexural），侧向弯曲（lateral）,和扭 转稳定性（torsional stability）分析时，应力刚度应该是被打开 的。BEAM188 能够采用 SECTYPE, SECDATA, SECOFFSET, SECWRITE, 和SECREAD来定义任何截面（形状）。.弹性（ela

16、sticity）,蠕变 （creep）,和塑性（plasticity）模型都是允许的（不考虑次截 面形状）。图1. BEAM188_ 3-D线性有限应变梁1.1.2. LINK10LINK10单元独一无二的双线性刚度矩阵特性使其成为一个轴向 仅受拉或仅受压杆单元。使用只受拉选项时，如果单元受压，刚度就 消失，以此来模拟缆索的松弛或链条的松弛。这一特性对于将整个钢 缆用一个单元来模拟的钢缆静力问题非常有用。当需要松弛单元的性 能，而不是关心松弛单元的运动时，它也可用于动力分析（带有惯性 或阻尼效应）。如果分析的目的时研究单元的运动（没有松弛单元），那么应该 使用类似于LINK10的不能松弛的单元

17、，比如：LINK8或PIPE59。对 于最终收敛结果为绷紧状态的结构，如果迭代过程中可能出现松弛状 态，那么这种静力收敛问题也不能使用LINK10单元。这时候应该采 用其它单元或者采用“缓慢动力”技术。LINK10单元在每个节点上有三个自由度：沿节点坐标系X、Y、Z 方向的平动，不管是仅受拉（缆）选项，还是仅受压（裂口）选项， 本单元都不包括弯曲刚度。本单元具有应力刚化、大变形功能。该单元的几何，节点位置以及坐标系见图1,单元通过两个节点、 横截面、初始应变或间隙以及各项同性材料特性来定义。单元的X轴 是沿着节点I到节点J的单元长度方向。单元的初始应变（ISTRN）由A/L给出，这里A是单元长

18、度L（由节点I和J的位置来定义的）和零应变长度Lo之间的差值。对 于缆选项，负的应变值表示其处于松弛状态。对于裂口选项，正的应 变值表示其处于裂开状态。这里裂口的值必须作为每单位长度的值输 入。I专d 0QVd! 24 . Q. L AE 1-07KU,_12flX.14E4-Q7. 13 SE4 07. 272E+07.-LdQE+Cn.4E4-O*7b711Al-20 44 07.3 4LE+a7.-4i?7E4-Q7three fpan ucintmii吕 grider5. 附件5 命令流/prep7/title,three span continus grider k,1,0,0k,2

19、,-2.1,0k,3,-2.6,-0.125 k,4,-2.8,-0.125 k,5,-3,-0.125 k,6,-3.4857,-0.1036 k,7,-3.9714,-0.0821 k,8,-4.4571,-0.0607 k,9,-4.9429,-0.0393 k,10,-5.4286,-0.0179 k,11,-5.9143,0.0036 k,12,-6.4,0.025 k,13,-2.800,-1.85 k,14,0.0000,-1.85 kgen,9,1,12,1,0,0,49/8,100 kgen,2,1,12,1,0,0,50,900 kgen,9,901,912,1,0,0,3

20、4.5/8,100 kgen,2,901,912,1,0,0,35,900c1=0.000843399 c2=0.001701323 *dim,x1,array,8 x1(1)=49/8 x1(2)=2*49/8 x1(3)=3*49/8 x1(4)=4*49/8 x1(5)=5*49/8 x1(6)=6*49/8 x1(7)=7*49/8 x1(8)=8*49/8 *dim,x2,array,8 x2(1)=34.5/8 x2(2)=2*34.5/8 x2(3)=3*34.5/8 x2(4)=4*34.5/8 x2(5)=5*34.5/8 x2(6)=6*34.5/8 x2(7)=7*34.

21、5/8 x2(8)=8*34.5/8*dim,yb,array,8yb(1)=c1*x1(1)*x1(1) yb(2)=c1*x1(2)*x1(2) yb(3)=c1*x1(3)*x1(3) yb(4)=c1*x1(4)*x1(4) yb(5)=c1*x1(5)*x1(5) yb(6)=c1*x1(6)*x1(6) yb(7)=c1*x1(7)*x1(7) yb(8)=c1*x1(8)*x1(8) *dim,ym,array,8 ym(1)=-c2*(x2(1)-34.5)*2+2.025 ym(2)=-c2*(x2(2)-34.5)*2+2.025 ym(3)=-c2*(x2(3)-34.5

22、)*2+2.025 ym(4)=-c2*(x2(4)-34.5)*2+2.025 ym(5)=-c2*(x2(5)-34.5)*2+2.025 ym(6)=-c2*(x2(6)-34.5)*2+2.025 ym(7)=-c2*(x2(7)-34.5)*2+2.025 ym(8)=-c2*(x2(8)-34.5)*2+2.025 kgen,2,13,14,1,0,-yb(1),x1(1),100 kgen,2,13,14,1,0,-yb(2),x1(2),200 kgen,2,13,14,1,0,-yb(3),x1(3),300 kgen,2,13,14,1,0,-yb(4),x1(4),400

23、 kgen,2,13,14,1,0,-yb(5),x1(5),500 kgen,2,13,14,1,0,-yb(6),x1(6),600 kgen,2,13,14,1,0,-yb(7),x1(7),700 kgen,2,13,14,1,0,-yb(8),x1(8),800 kgen,2,813,814,1,0,0,1,100 kgen,2,913,914,1,0,ym(1),x2(1),100 kgen,2,913,914,1,0,ym(2),x2(2),200 kgen,2,913,914,1,0,ym(3),x2(3),300 kgen,2,913,914,1,0,ym(4),x2(4),

24、400 kgen,2,913,914,1,0,ym(5),x2(5),500 kgen,2,913,914,1,0,ym(6),x2(6),600 kgen,2,913,914,1,0,ym(7),x2(7),700 kgen,2,913,914,1,0,ym(8),x2(8),800 kgen,2,1,14,1,0,0,-0.5,50 et,1,shell63 mp,ex,1,3.5e10 mp,dens,1,2500 mp,prxy,0.1667 r,1,0.25*do,i,0,16,1 a,1+i*100,1+(i+1)*100,2+(i+1)*100,2+i*100 *enddoa,1

25、,51,52,2 aatt,1,1,1 asel,none r,2,0.375 *do,i,0,16,1 a,2+i*100,2+(i+1)*100,3+(i+1)*100,3+i*100 *enddo a,2,52,53,3 aatt,1,2,1 asel,none r,3,0.375 *do,i,0,16,1 a,3+i*100,3+(i+1)*100,4+(i+1)*100,4+i*100 a,4+i*100,4+(i+1)*100,5+(i+1)*100,5+i*100 *enddoa,3,53,54,4 a,4,54,55,5 aatt,1,3,1 asel,none*dim,hd,

26、array,9 hd(9)=0.2hd(8)=0.2+(3.4/7)*0.3/3.4 hd(7)=0.2+2*(3.4/7)*0.3/3.4 hd(6)=0.2+3*(3.4/7)*0.3/3.4 hd(5)=0.2+4*(3.4/7)*0.3/3.4 hd(4)=0.2+5*(3.4/7)*0.3/3.4 hd(3)=0.2+6*(3.4/7)*0.3/3.4 hd(2)=0.2+7*(3.4/7)*0.3/3.4 hd(1)=0.5*do,i,1,8,1 r,30+i,hd(i+1),hd(i+1),hd(i),hd(i) *enddo*do,i,0,16,1 *do,k,5,11,1 a

27、,k+i*100,k+(i+1)*100,k+1+(i+1)*100,k+1+i*100 a,k,k+50,k+1+50,k+1aatt,1,k+26,1 asel,none *enddo *enddo asel,none r,4,0.5 *do,i,0,16,1 a,4+i*100,4+(i+1)*100,13+(i+1)*100,13+i*100*enddo a,4,54,63,13 aatt,1,4,1 asel,none c3=6.2474E-05 c4=0.000126024 *dim,h1,array,8 h1(1)=0.25+c3*x1(1)*x1(1) h1(2)=0.25+c

28、3*x1(2)*x1(2) h1(3)=0.25+c3*x1(3)*x1(3)h1(4)=0.25+c3*x1(4)*x1(4)h1(5)=0.25+c3*x1(5)*x1(5)h1(6)=0.25+c3*x1(6)*x1(6) h1(7)=0.25+c3*x1(7)*x1(7) h1(8)=0.25+c3*x1(8)*x1(8) *dim,h2,array,9h2(1)=0.25+c4*x2(1)*x2(1) h2(2)=0.25+c4*x2(2)*x2(2) h2(3)=0.25+c4*x2(3)*x2(3) h2(4)=0.25+c4*x2(4)*x2(4) h2(5)=0.25+c4*

29、x2(5)*x2(5) h2(6)=0.25+c4*x2(6)*x2(6) h2(7)=0.25+c4*x2(7)*x2(7) h2(8)=0.25+c4*x2(8)*x2(8) h2(9)=0.4*do,i,1,8,1 r,10+i,h1(i) a,13+(i-1)*100,13+i*100,14+i*100,14+(i-1)*100 aatt,1,10+i,1asel,none *enddo *do,i,1,9,1 r,20+i,h2(i) a,13+(i+7)*100,13+(i+8)*100,14+(i+8)*100,14+(i+7)*100 a,13,63,64,14aatt,1,2

30、0+i,1 asel,none *enddo r,50,2 a,1,51,64,14 a,801,901,914,814 aatt,1,50,1 asel,none allsel esize,0.4 mshape,0 mshkey,1 amesh,all nsel,all nsym,x,10000,1,80000,1 esym,10000,1,10000,1 cskp,12,0,1714,1713,1701 nsym,z,100000,1,100000,1 esym,100000,1,1000000,1 allsel, nummrg,all csys,0nsel,s,loc,y,-1.85 n

31、sel,r,loc,x,-3,-2.6nsel,r,loc,z,-0.5,0 d,all,ux,uy,roty,rotz allsel nsel,s,loc,y,-1.85 nsel,r,loc,x,2.6,3 nsel,r,loc,z,-0.5,0 d,all,ux,uy,roty,rotz allsel nsel,s,loc,y,-1.85 nsel,r,loc,x,-3,-2.6 nsel,r,loc,z,169,169.5 d,all,ux,uy,roty,rotz allsel nsel,s,loc,y,-1.85 nsel,r,loc,x,2.6,3 nsel,r,loc,z,16

32、9,169.5 d,all,ux,uy,roty,rotz allsel nsel,s,loc,y,-3.875 nsel,r,loc,x,2.6,3 nsel,r,loc,z,49,50 d,all,ux,uy,uz,roty,rotz allsel nsel,s,loc,y,-3.875 nsel,r,loc,x,-3,-2.6 nsel,r,loc,z,49,50 d,all,ux,uy,uz,roty,rotz allselnsel,s,loc,y,-3.875 nsel,r,loc,x,-3,-2.6 nsel,r,loc,z,119,120 d,all,ux,uy,uz,roty,rotz allselnsel,s,loc,y,-3.875 nsel,r,loc,x,2.6,3 nsel,r,loc,z,119,120 d,all,ux,uy,uz,roty,rotz allsel /solu allsel acel,9.8 solve /post1 plnsol,s,1 nsel,s,loc,z,10 esln,s,all plnsol,s,1 finish save