搜档网
当前位置:搜档网 › 有限元非线性计算特点

有限元非线性计算特点

有限元非线性计算特点
有限元非线性计算特点

有限元非线性计算特点

文章通过几个典型的工程计算模型,分析比较有限元线性与非线性计算结果,阐释了有限元非线性计算的特点及优点。

标签:工程计算;线性;非线性

1 引言

有限元单元法已成为强有力的数值解法来解决工程中遇到的大量问题,其应用范围从固体到流体,从静力到动力,从力学问题到非力学问题,有限元的线性分析已被广泛采用。但对于许多航空工程中遇到的问题,如进气道等,仅仅采用线性求解是不真实的,而采用非线性计算将更符号实际情况。本文借助MSC/NASTRAN有限元分析程序,对于典型的工程计算模型分析比较线性与非线性计算结果,从而给出非线性计算相对于线性计算的优点及特点。

2 有限元非线性计算的特点及优点

为了明确有限元非线性计算结果与线性计算结果的差异,更好的展现有限元非线性计算的特点,本节将借助于有限元分析软件MSC/NASTRAN,对一受外载的矩形薄板根据不同的边界条件,进行非线性及线性静力分析,通过分析比较计算结果,说明有限元非线性静力计算中的一些特点。

2.1 非线性与线性计算结果随载荷的变化

首先,给出薄板尺寸、载荷。

模型尺寸:薄板尺寸为500×500×1.5mm。

载荷:受法向气动压力(pressure),气动压力由小到大变化依次为0.01MPa、0.02MPa、0.04MPa、0.08MPa、0.16MPa。

取薄板中央节点位移、应力及薄板边缘中部节点位移,比较线性计算结果和非线性计算结果。在分别进行有限元线性及非线性分析后,给出位移、应力及支反力结果随载荷的变化曲线。图1、图3、图5分别为采用限元线性计算得到的参考点的位移、应力及支反力变化曲线;图2、图4、图6分别为采用有限元非线性计算得到的参考点的位移、应力及支反力变化曲线。

由圖1、3、5可见,采用线性静力分析后,参考点位移、应力、支反力均随载荷增加而线性增大,位移、应力、支反力与载荷呈明显的线性关系,这是线性静力分析的特点。对于本例,可以预言,在其它条件不变的情况下,计算出一套载荷下的结果,就可以按照线性关系求出压力载荷下的位移、应力及支反力结果。

有限元非线性计算特点

有限元非线性计算特点 文章通过几个典型的工程计算模型,分析比较有限元线性与非线性计算结果,阐释了有限元非线性计算的特点及优点。 标签:工程计算;线性;非线性 1 引言 有限元单元法已成为强有力的数值解法来解决工程中遇到的大量问题,其应用范围从固体到流体,从静力到动力,从力学问题到非力学问题,有限元的线性分析已被广泛采用。但对于许多航空工程中遇到的问题,如进气道等,仅仅采用线性求解是不真实的,而采用非线性计算将更符号实际情况。本文借助MSC/NASTRAN有限元分析程序,对于典型的工程计算模型分析比较线性与非线性计算结果,从而给出非线性计算相对于线性计算的优点及特点。 2 有限元非线性计算的特点及优点 为了明确有限元非线性计算结果与线性计算结果的差异,更好的展现有限元非线性计算的特点,本节将借助于有限元分析软件MSC/NASTRAN,对一受外载的矩形薄板根据不同的边界条件,进行非线性及线性静力分析,通过分析比较计算结果,说明有限元非线性静力计算中的一些特点。 2.1 非线性与线性计算结果随载荷的变化 首先,给出薄板尺寸、载荷。 模型尺寸:薄板尺寸为500×500×1.5mm。 载荷:受法向气动压力(pressure),气动压力由小到大变化依次为0.01MPa、0.02MPa、0.04MPa、0.08MPa、0.16MPa。 取薄板中央节点位移、应力及薄板边缘中部节点位移,比较线性计算结果和非线性计算结果。在分别进行有限元线性及非线性分析后,给出位移、应力及支反力结果随载荷的变化曲线。图1、图3、图5分别为采用限元线性计算得到的参考点的位移、应力及支反力变化曲线;图2、图4、图6分别为采用有限元非线性计算得到的参考点的位移、应力及支反力变化曲线。 由圖1、3、5可见,采用线性静力分析后,参考点位移、应力、支反力均随载荷增加而线性增大,位移、应力、支反力与载荷呈明显的线性关系,这是线性静力分析的特点。对于本例,可以预言,在其它条件不变的情况下,计算出一套载荷下的结果,就可以按照线性关系求出压力载荷下的位移、应力及支反力结果。

matlab有限元分析实例

MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件,用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人,控制系统等领域。 MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂(矩阵实验室),软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。 MATLAB有限元分析与应用:

《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书,作者是卡坦,译者是韩来彬。 内容简介: 《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用,书中的每个示例都以交互的方式求解,使读者很容易就能把MATLAB用于有限分析和应用。另外,《MATLAB有限元分析与应用》还提供了大量免费资源。 《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂,循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题,具有很高的工程应用价值。

汽车成功案例

汽车成功案例 安全性问题 竞争优势 全球汽车工业对汽车安全性越来越重视,与安全强制法规相关的试验也在大量增加。目前碰撞安全问题在碰撞前、碰撞中和碰撞后阶段同时展开研究。在碰撞前阶段利用主动避撞系统;在碰撞中阶段利用车身结构、气囊展开、安全带张紧等措施减小伤害;在碰撞后阶段,主要关心油箱是否破裂以防止爆炸或起火。MSC.Software虚拟产品开发设计能够对每一个阶段进行设计研究。 碰撞前阶段 避免碰撞发生当然是车辆交通中最有效的降低伤亡的方法。而车辆的行为,例如车辆打滑、侧翻、或者车轮遇到冰路面将会发生何种状况等等可以利用虚拟样机来预测。在ADAMS/Car中结合多刚体和控制的仿真可以模拟从主动悬架到ABS制动器等系统的试验来增加主动安全性。通过同步调整机械、控制系统对车辆进行优化,可以大大缩短设计周期。 碰撞中阶段 一旦碰撞不可避免,气囊展开和座椅安全带的预张紧就成为减小伤害的关键因素,虚拟产品开发能够对这些系统进行优化。气囊展开可以利用SimOffice中的MSC Dytran,安全带约束系统的力可以利用多体仿真分析软件。在样车建造和法规试验之前进行虚拟试验可以大大地降低开发费用。法规试验中车辆各种性能可以用SimOffice中提供的有限元方法来进行精确地预测和研究。

碰撞仿真流程通常需要大量人力,管理仿真产生的海量数据也是一个挑战。模型组装、质量检查、定义工况、报告准备等方面如果引入流程自动化和数据管理则可以节省大量的人力。MSC.Software是领先的流程管理和自动化工具供应商,其产品MSC SOFY 和MSC SimManager都提供了汽车碰撞流程自动化的环境。将工作流程确定下来并进行客户化配置后,软件工具可以自动地生成代码来指导用户完成工作流程。例如,德国宝马(BMW)公司利用MSC SimManager建立碰撞仿真自动化流程,管理海量仿真数据,并且可以和供应商合作,使供应商可以上载各自相关的部件。 LSTC公司的领先的碰撞求解器LS-Dyna可以通过MSC Nastran(Sol700)的标准格式来调用。因此,适撞性和显著非线性问题都可以采用和NVH部门同样的模型,这样通过不同部门的协作可以节省大量的时间和费用。 碰撞后阶段 避免碰撞后起火取决于供油系统的完整性,该项安全要求 已在美国安全法规FMVSS301中有明确规定。车辆碰撞 后的燃油泄漏必须避免,MSC.Dytran采用拉格朗日和欧 拉技术,可以模拟碰撞中和碰撞后油箱的液固作用、结构 大变形、结构接触等问题。 MSC.SimManager也可以集成到碰撞后开发流程中,一 级供应商TI汽车公司采用MSC.SimManager管理油箱 开发过程中的冲击、压力真空、跌落、下陷等试验。 车辆动力学问题 矛盾 汽车工业需要在开发过程中减少时间和费用,同时推出创 新的产品。当前比较通用的策略是利用通用的开发平台、 共享部件开发众多系列车型。这就导致出现两个相互矛盾 的目标:一个是新系统的开发,另一个是通过共用平台和 零部件减少系统的变型。借助于虚拟产品开发可以有效地 满足这两个目标。

有限元复习题库

有限元复习 一、选择题(每题1分,共10分) 二、判断题(每空1分,共10分) 三、填空题(每空1分,共10分) 三、简答题(共44分)共6题 四、综述题(共26分)两题 一.基本概念 1. 平面应力/平面应变问题;空间问题/轴对称问题;杆梁问题;线 性与非线性问题 平面应力问题 (1) 均匀薄板(2)载荷平行于板面且沿厚度方向均匀分布 在六个应力分量中,只需要研究剩下的平行于XOY 平面的三个应力分量,即x y xy yx σσττ=、、 (000z zx xz zy yz σττττ=====,,)。 一般0z σ=,z ε并不一定等于零,但可由x σ及y σ求得,在分析问题时不必 考虑。于是只需要考虑 x y xy εεγ、、三个应变分量即可。 平面应变问题

(1) 纵向很长,且横截面沿纵向不变。(2)载荷平行于横截面且沿纵向 均匀分布 z yz zx εγγ===只剩下三个应变分量x y xy εεγ、、。也只需要考虑x y xy σστ、、三个应力分量即可 轴对称问题 物体的几何形状、约束情况及所受外力都对称于空间的某一根轴。 轴对称单元的特点(与平面三角形单元的区别):轴对称单元为圆环体,单元与单元间为节圆相连接;节点力与节点载荷是施加于节圆上的均布力;单元边界是一回转面;应变不是常量。 在轴对称问题中,周向应变分量θε是与r 有关。 板壳问题 一个方向的尺寸比另外两个方向尺寸小很多,且能承受弯矩的结构称为板壳结构,并把平分板壳结构上下表面的面称为中面。如果中面是平面或平面组成的折平面,则称为平板;反之,中面为曲面的称为壳。 杆梁问题 杆梁结构是指长度远大于其横断面尺寸的构件组成的系统。在结构力学中常将承受轴力或扭矩的杆件称为杆,而将承受横向力和弯矩的杆件称为梁。 平面(应力应变)问题与板壳问题的区别与联系 平面应力问题是指很薄的等厚度薄板,只在板边上受有平行于板面并且不沿厚度变化的面力,同时,体力也平行于板面并且不沿厚度变化。而平面应变问题是指很长的柱形体,在柱面上受有平行于横截面并且不沿长度变化的面力,同时

有限元分析系统的发展现状与展望外文翻译

Finite element analysis system development present situation and forecast Along with modern science and technology development, the people unceasingly are making the faster transportation vehicle, the large-scale building, the greater span bridge, the high efficiency power set and the preciser mechanical device. All these request engineer to be able precisely to forecast in the design stage the product and the project technical performance, needs to be static, technical parameter and so on dynamic strength to the structure as well as temperature field, flow field, electromagnetic field and transfusion carries on the analysis computation. For example analysis computation high-rise construction and great span bridge when earthquake receives the influence, has a look whether can have the destructive accident; The analysis calculates the nuclear reactor the temperature field, the determination heat transfer and the cooling system are whether reasonable; Analyzes in the new leaf blade the hydrodynamics parameter, enhances its operating efficiency. The sell may sum up as the solution physics question control partial differential equations often is not impossible. In recent years the finite element analysis which develops in the computer technology and under the numerical analysis method support(FEA, Finite Element Analysis) the side principle for solves these complex project analysis estimation problems to provide the effective way. Our country in " 95 " Plan period vigorously promotes the CAD technology, mechanical profession large and middle scalene terries CAD popular rate from " 85 " End 20% enhances that present 70%.With enterprise application of CAD, engineering and technical personnel has gradually get rid drawing board, and will join the main energy how to optimize the design, engineering and improving the quality of products, computer-aided engineering analysis (CAE. Computer Aided Engineering) method and software will be the key technical elements . ln engineering practice, finite element analysis software and CAD system integration design standards should be a qualitative leap, mainly in the following aspects : The increase design function, reduces the design cost; Reduces design and the analysis cycle period; Increase product and project reliability; Uses the optimized design, reduces the material the consumption or the cost;

非线性有限元分析

轨道结构的非线性有限元分析 姜建华 练松良 摘 要 实际轨道结构受载时的力学行为,属于典型的非线性力学问题。钢轨垫层刚度、钢轨抗扭刚度和扣件扣压力的大小是影响轨距扩大的主要因素。根据非线性有限元接触理论,建立了能准确反映扣件、钢轨与垫层的拧紧接触,以及受载车轮与钢轨侧向滑动接触的力学计算模型;并研究计算了不同扣件压力下,由于受载车轮与钢轨侧向滑动接触引起的轨距扩大问题。 关键词 轮轨关系,扣件压力,非线性弹性力学,有限元分析 1 引言 实际工程中常见的非线性问题一般可以归纳为三类:材料非线性、几何非线性以及边界条件非线性。材料非线性问题是由于材料的非线性本构关系所引起的,例如材料的弹塑性变形,材料的屈服和硬化等;几何非线性问题是由于结构的位移或变形相当大,以至必须按照变形后的几何位置来建立平衡方程;边界条件非线性问题是指边界条件随位移变化所引起的非线性问题。通常情况下,我们所遇到的非线性问题多数是上述三类非线性问题的组合[1,2]。 实际轨道结构受载时的力学行为,属于典型的非线性力学问题。比如基于轮轨接触的材料非线性、几何非线性及边界条件非线性问题,以及扣件、钢轨、垫层三者间相互作用时所表现的边界条件非线性行为等。所以,机车车辆在轨道结构上行驶时引起的力学现象是相当复杂的。以往在研究轨道各部分应力应变分布规律时,通常采用连续弹性基础梁理论或连续点支承,偶尔简单考虑扣件的作用和弹性垫层的使用。不管用哪一种支承方式建立模型,都由于这样那样的假设而带有一定程度的近似性。所以,如何利用现代力学理论的最新成果以及日益发展的计算机技术,根据轨道结构的具体情况,建立更为完整更为准确的轨道结构计算模型,为轨道设计部门提供更加可靠的设计依据或研究思路,已十分必要。 本文提出了用非线性有限元理论研究轮轨系统和轨道结构的思路。作为算例之一,本文将根据非线性有限元理论,建立能准确反映扣件、钢轨与垫层的拧紧接触,以及受载车轮与钢轨侧向滑动接触的力学计算模型。 2 轨道结构的有限元接触模型 对于非线性问题,不管是材料非线性、几何非线性,还是边界条件非线性,总是最终归结为求解一组非线性平衡方程及其控制方程。例如用位移作为未知数进行有限元分析时,最后可得到一组平衡方程及其控制方程为 : 图1 轮轨系统的对称性模型简图 [K(u)]{u}={R}(1) (u)= (u)(2)其中:{u}为节点位移列阵;{R}为节点载荷列阵; [K(u)]为总体刚度矩阵; (u)为边界条件。它们 36 姜建华:同济大学工程力学系,副教授、博士,上海200092

Ansys有限元分析实例[教学]

Ansys有限元分析实例[教学] 有限元分析案例:打点喷枪模组(用于手机平板电脑等电子元件粘接),该产品主要是使用压缩空气推动模组内的顶针作高频上下往复运动,从而将高粘度的胶水从喷嘴中打出(喷嘴尺寸,0.007”)。顶针是这个产品中的核心零件,设计使用材料是:AISI 4140 最高工作频率是160HZ(一个周期中3ms开3ms关),压缩空气压力3-8bar, 直接作用在顶针活塞面上,用Ansys仿真模拟分析零件的强度是否符合要求。 1. 零件外形设计图:

2. 简化模型特征后在Ansys14.0 中完成有限元几何模型创建:

3. 选择有限元实体单元并设定,单元类型是SOILD185,由于几何建模时使用的长度单位是mm, Ansys采用单位是长度:mm 压强: 3Mpa 密度:Ton/M。根据题目中的材料特性设置该计算模型使用的材料属性:杨氏模量 2.1E5; 泊松比:0.29; 4. 几何模型进行切割分成可以进行六面体网格划分的规则几何形状后对各个实体进行六面体网格划分,网格结果: 5. 依据使用工况条件要求对有限元单元元素施加约束和作用载荷:

说明: 约束在顶针底端球面位移全约束; 分别模拟当滑块顶断面分别以8Bar,5Bar,4Bar和3Bar时分析顶针的内应力分布,根据计算结果确定该产品允许最大工作压力范围。 6. 分析结果及讨论: 当压缩空气压力是8Bar时: 当压缩空气压力是5Bar时:

当压缩空气压力是4Bar时: 结论: 通过比较在不同压力载荷下最大内应力的变化发现,顶针工作在8Bar时最大应力达到250Mpa,考虑到零件是在160HZ高频率在做往返运动,疲劳寿命要求50百万次以上,因此采用允许其最大工作压力在5Mpa,此时内应力为156Mpa,按线性累积损伤理论[3 ]进行疲劳寿命L-N疲劳计算,进一部验证产品的设计寿命和可靠性。

有限元仿真技术的发展及其应用

有限元仿真技术的发展及其应用 许荣昌 孙会朝(技术研发中心) 摘 要:介绍了目前常用的大型有限元分析软件的现状与发展,对其各自的优势进行了分析,简述了有限元软件在冶金生产过程中的主要应用领域及其发展趋势,对仿真技术在莱钢的应用进行了展望。 关键词:有限元仿真 冶金生产 发展趋势 0 前言 自主创新,方法先行,创新方法是自主创新的根本之源,同时,随着市场竞争的日益激烈,冶金企业的产品设计、工艺优化也由经验试错型向精益研发方向发展,而有限元仿真技术正是这种重要的创新方法。近年来随着计算机运行速度的不断提高,有限元分析在工程设计和分析中得到了越来越广泛的应用,比如,有限元分析在冶金、航空航天、汽车、土木建筑、电子电器、国防军工、船舶、铁道、石化、能源、科学研究等各个领域正在发挥着重要的作用,主要表现在以下几个方面:增加产品和工程的可靠性;在产品的设计阶段发现潜在的问题;经过分析计算,采用优化设计方案,降低原材料成本;缩短产品研发时间;模拟试验方案,减少试验次数,从而减少试验成本。与传统设计相比,利用仿真技术,可以变经验设计为科学设计、变实测手段为仿真手段、变规范标准为分析标准、变传统分析技术为现代的计算机仿真分析技术,从而提高产品质量、缩短新产品开发周期、降低产品整体成本、增强产品系统可靠性,也就是增强创新能力、应变能力和竞争力(如图1、2) 。 图1 传统创新产品(工艺优化)设计过程为大循环 作者简介:许荣昌(1971-),男,1994年毕业于武汉钢铁学院钢铁冶金专业,博士,高级工程师。主要从事钢铁工艺技术研究工 作。 图2 现代CA E 创新产品(工艺优化)设计过程为小循环 1 主要有限元分析软件简介 目前,根据市场需求相继出现了各种类型的应用软件,其中NASTRAN 、ADI N A 、ANSYS 、 ABAQUS 、MARC 、MAGSOFT 、COS MOS 等功能强大的CAE 软件应用广泛,为实际工程中解决复杂的理论计算提供了非常有力的工具。但是,各种软件均有各自的优势,其应用领域也不尽相同。本文将就有限元的应用范围及当今国际国内C AE 软件的发展趋势做具体的阐述,并对与冶金企业生产过程密切相关的主要有限元软件ANSYS 、AB AQUS 、MARC 的应用领域进行分析。 M SC So ft w are 公司创建于1963年,总部设在美国洛杉矶,M SC M arc 是M SC Soft w are 公司于1999年收购的MARC 公司的产品。MARC 公司始创于1967年,是全球首家非线性有限元软件公司。经过三十余年的发展,MARC 软件得到学术界和工业界的大力推崇和广泛应用,建立了它在全球非线性有限元软件行业的领导者地位。随着M arc 软件功能的不断扩展,软件的应用领域也从开发初期的核电行业迅速扩展到航空、航天、汽车、造船、铁 道、石油化工、能源、电子元件、机械制造、材料工程、土木建筑、医疗器材、冶金工艺和家用电器等,成为许多知名公司和研究机构研发新产品和新技术的重要工具。在航空业M SC N astran 软件被美国联邦航空管理局(F AA )认证为领取飞行器适 13

(完整word版)有限元分析软件的比较

有限元分析软件的比较(购买必看)-转贴 随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能,需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响,看看是否会发生破坏性事故;分析计算核反应堆的温度场,确定传热和冷却系统是否合理;分析涡轮机叶片内的流体动力学参数,以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式,这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析(FEA,Finite Element A nalysis)方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中,有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃,主要表现在以下几个方面: 增加设计功能,减少设计成本; 缩短设计和分析的循环周期; 增加产品和工程的可靠性; 采用优化设计,降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案,减少试验时间和经费; 进行机械事故分析,查找事故原因。 在大力推广CAD技术的今天,从自行车到航天飞机,所有的设计制造都离不开有限元分析计算,FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局(NASA)在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本,是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在,世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件,主要有德国的ASKA、英国的PA FEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 以下对一些常用的软件进行一些比较分析: 1. LSTC公司的LS-DYNA系列软件

非线性有限元方法及实例分析

非线性有限元方法及实例分析 梁军 河海大学水利水电工程学院,南京(210098) 摘 要:对在地下工程稳定性分析中常用的非线性方程组的求解方法进行研究,讨论了非线性计算的迭代收敛准则,并利用非线性有限元方法分析了一个钢棒单轴拉伸的实例。 关键词:非线性有限元,方程组求解,实例分析 1引 言 有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题,其应用范围从固体到流体,从静力到动力,从力学问题到非力学问题。有限元的线性分析已经设计工具被广泛采用。但对于绝大多数水利工程中遇到的实际问题如地下洞室等,将其作为非线性问题加以考虑更符合实际情况。根据产生非线性的原因,非线性问题主要有3种类型[1]: 1.材料非线性问题(简称材料非线性或物理非线性) 2.几何非线性问题 3.接触非线性问题(简称接触非线性或边界非线性) 2 非线性方程组的求解 在非线性力学中,无论是哪一类非线性问题,经过有限元离散后,它们都归结为求解一个非线性代数方程组[2]: ()()()00 021212211=… …==n n n n δδδψδδδψδδδψΛΛΛ (1.1) 其中n δδδ,,,21Λ是未知量,n ψψψ,,,21Λ是n δδδ,,,21Λ的非线性函数,引用矢量记 号 []T n δδδδΛ21= (1.2) []T n ψψψψΛ21= (1.3) 上述方程组(1.1)可表示为 ()0=δψ (1.4) 可以将它改写为 ()()()0=?≡?≡R K R F δδδδψ (1.5) 其中()δK 是一个的矩阵,其元素 是矢量的函数,n n ×ij k R 为已知矢量。在位移有限 元中,δ代表未知的结点位移,()δF 是等效结点力,R 为等效结点荷载,方程()0=δψ表示结点平衡方程。 在线弹性有限元中,线性方程组

有限元分析案例

有限元分析案例 图1 钢铸件及其砂模的横截面尺寸 砂模的热物理性能如下表所示: 铸钢的热物理性能如下表所示: 一、初始条件:铸钢的温度为2875o F,砂模的温度为80o F;砂模外边界的对流边界条件:对流系数0.014Btu/hr.in2.o F,空气温度80o F;求3个小时后铸钢及砂模的温度分布。 二、菜单操作: 1.Utility Menu>File>Change Title, 输入Casting Solidification; 2.定义单元类型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete, Add, Quad 4node 55; 3.定义砂模热性能:Main Menu>Preprocessor>Material Props>Isotropic,默认材料编号1, 在Density(DENS)框中输入0.054,在Thermal conductivity (KXX)框中输入0.025,在S pecific heat(C)框中输入0.28; 4.定义铸钢热性能温度表:Main Menu>Preprocessor>Material Props>-Temp Dependent->Temp Table,输入T1=0,T2=2643, T3=2750, T4=2875; 5.定义铸钢热性能:Main Menu>Preprocessor>Material Props>-Temp Dependent ->Prop Table, 选择Th Conductivity,选择KXX, 输入材料编号2,输入C1=1.44, C2=1.54, C3=1.22, C4=1.22,选择Apply,选择Enthalpy,输入C1=0, C2=128.1, C3=163.8, C4=174.2; 6.创建关键点:Main Menu>Preprocessor>-Modeling->Create>Keypoints>In Active

对有限元方法的认识

我对有限元方法的认识 1有限元法概念 有限元方法(The Finite Element Method, FEM)是计算机问世以后迅速发展起来的一种分析方法。每一种自然现象的背后都有相应的物理规律,对物理规律的描述可以借助相关的定理或定律表现为各种形式的方程(代数、微分、或积分)。这些方程通常称为控制方程(Governing equation)。 针对实际的工程问题推导这些方程并不十分困难,然而,要获得问题的解析的数学解却很困难。人们多采用数值方法给出近似的满足工程精度要求的解答。 有限元方法就是一种应用十分广泛的数值分析方法。 有限元方法是处理连续介质问题的一种普遍方法,离散化是有限元方法的基础。 这种思想自古有之:古代人们在计算圆的周长或面积时就采用了离散化的逼近方法:即采用内接多边形和外切多边形从两个不同的方向近似描述圆的周长或面积,当多边形的边数逐步增加时近似值将从这两个方向逼近真解。 近年来随着计算机技术的普及和计算速度的不断提高,有限元分析在工程设计和分析中得到了越来越广泛的重视,已经成为解决复杂的工程分析计算问题的有效途径,现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算,其在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器、国防军工、船舶、铁道、石化、能源、科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃。 国际上早在 60 年代初就开始投入大量的人力和物力开发有限元分析程序。“有限单元”是由Clough R W于1960年首次提出的。但真正的有限元分析软件是诞生于 70 年代初期,随着计算机运算速度的提高,内、外存容量的扩大和图形设备的发展,以及软件技术的进步,发展成为有限元分析与设计软件,但初期其前后处理的能力还是比较弱的,特别是后处理能力更弱。

板结构有限元分析实例详解

板结构有限元分析实例详解1:带孔平板结构静力分析本节介绍带孔平板结构静力分析问题,同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析 有孔的矩形平板,左侧边缘固定,长400mm,宽200 mm,厚度为10 mm,圆孔在板的正中心,半径为40 mm,左侧全约束,右侧边缘均布应力1MPa,如图8.7所示。求板的变形、位移及应力变化情况。(材料的材料属性为:弹性模量为300000 MPa,剪切模量为0.31。) 图8.7 带孔的矩形平板 由于小孔处边缘不规则,本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名 启动ANSYS Mechanical APDL Product Launcher窗口,如图8.8所示。在License下 拉选框中选择ANSYS Multiphysics产品,在Working Directory输入栏中输入工作目 录:C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1,在Job Name一栏中输入工作文件名:Chapter8-1。以上参数设置完毕后,单 击Run按钮运行ANSYS。

图8.8 ANSYS设置窗口菜单 可以先在目标文件位置建立工作目录,然后单击Browse按钮选择工作目录;也 可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性 选择Main Menu>Preferences命令,出现Preferences for GUI Filtering对话框, 如图8.9所示,在Individual discipline(s) to show in the GUI中勾选Structural,过滤掉ANSYS GUI菜单中与结构分析无关的选项,单击OK按钮关闭该对话框。 图8.9 Preferences for GUI Filtering对话框

有限元分析71831

有限元分析 有限元法的基本思想是将结构离散化,用有限个容易分析的单元来表示复杂的对象,单元之间通过有限个节点相互连接,然后根据变形协调条件综合求解。由于单元的数目是有限的,节点的数目也是有限的,所以称为有限元法(FEM,Finite Element Method)。 有限元法是一种求解关于场问题的一系列偏微分方程的数值方法.这种类型的问题会在许多工程学科中遇到,如机械设计、声学、电磁学、岩土力学、断裂力学、流体力学等.在机械工程中,有限元分析被光分应用在结构、振动和传热问题上。 有限元法是60年代以来发展起来的新的数值计算方法,是计算机时代的产物。虽然有限元的概念早在40年代就有人提出,但由于当时计算机尚未出现,它并未受到人们的重视。随着计算机技术的发展,有限元法在各个工程领域中不断得到深入应用,现已遍及宇航工业、核工业、机电、化工、建筑、海洋等工业,是机械产品动、静、热特性分析的重要手段。早在70年代初期就有人给出结论:有限元法在产品结构设计中的应用,使机电产品设计产生革命性的变化,理论设计代替了经验类比设计。目前,有限元法仍在不断发展,理论上不断完善,各种有限元分析程序包的功能越来越强大,使用越来越方便。 大约在300年前,牛顿和莱布尼茨发明了积分法,证明了该运算具有整体对局部的可加性。虽然,积分运算与有限元技术对定义

域的划分是不同的,前者进行无限划分而后者进行有限划分,但积分运算为实现有限元技术准备好了一个理论基础。 在牛顿之后约一百年,著名数学家高斯提出了加权余值法及线性代数方程组的解法。这两项成果的前者被用来将微分方程改写为积分表达式,后者被用来求解有限元法所得出的代数方程组。在18世纪,另一位数学家拉格郎日提出泛函分析。泛函分析是将偏微分方程改写为积分表达式的另一途经。 在19世纪末及20世纪初,数学家瑞雷和里兹首先提出可对全定义域运用展开函数来表达其上的未知函数。1915年,数学家伽辽金提出了选择展开函数中形函数的伽辽金法,该方法被广泛地用于有限元。1943年,数学家库朗德第一次提出了可在定义域内分片地使用展开函数来表达其上的未知函数。这实际上就是有限元的做法。 所以,到这时为止,实现有限元技术的第二个理论基础也已确立。 20世纪50年代,飞机设计师们发现无法用传统的力学方法分析飞机的应力、应变等问题。波音公司的一个技术小组,首先将连续体的机翼离散为三角形板块的集合来进行应力分析,经过一番波折后获得前述的两个离散的成功。20世纪50年代,大型电子计算机投入了解算大型代数方程组的工作,这为实现有限元技术准备好了物质条件。1960年前后,美国的R. W. Clough教授及我国的冯康教授分别独立地在论文中提出了“有限单元”,这样的名词。此后,这样

ANSYS有限元分析实例

有限元分析 一个厚度为20mm的带孔矩形板受平面内张力,如下图所示。左边固定,右边受载荷p=20N/mm作用,求其变形情况 P 一个典型的ANSYS分析过程可分为以下6个步骤: ①定义参数 ②创建几何模型 ③划分网格 ④加载数据 ⑤求解 ⑥结果分析 1定义参数 1.1指定工程名和分析标题 (1)启动ANSYS软件,选择File→Change Jobname命令,弹出如图所示的[Change Jobname]对话框。 (2)在[Enter new jobname]文本框中输入“plane”,同时把[New log and error files]中的复选框选为Yes,单击确定 (3)选择File→Change Title菜单命令,弹出如图所示的[Change Title]对话框。 (4)在[Enter new title]文本框中输入“2D Plane Stress Bracket”,单击确定。 1.2定义单位

在ANSYS软件操作主界面的输入窗口中输入“/UNIT,SI” 1.3定义单元类型 (1)选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令,弹出如图所示[Element Types]对话框。 (2)单击[Element Types]对话框中的[Add]按钮,在弹出的如下所示[Library of Element Types]对话框。 (3)选择左边文本框中的[Solid]选项,右边文本框中的[8node 82]选项,单击确定,。 (4)返回[Element Types]对话框,如下所示 (5)单击[Options]按钮,弹出如下所示[PLANE82 element type options]对话框。

有限元分析的发展趋势

有限元分析的发展趋势 摘要:1965年“有限元”这个名词第一次出现,到今天有限元在工程上得到广泛应用,经历了三十多年的发展历史,理论和算法都已经日趋完善。有限元的核心思想是结构的离散化,就是将实际结构假想地离散为有限数目的规则单元组合体,实际结构的物理性能可以通过对离散体进行分析,得出满足工程精度的近似结果来替代对实际结构的分析,这样可以解决很多实际工程需要解决而理论分析又无法解决的复杂问题。 关键词:有限元分析结构计算结构设计 Abstract: The 1965 "finite" appeared for the first time this term, and today is widely used finite element in engineering, after more than 30 years of history, theory and algorithms have been improved. Finite element discretization of the core idea is to structure, is the actual structure of the supposed discrete combination unit for a limited number of rules, the actual structure to analyse the physical properties can be felt through a discrete body of drawn precision engineering approximation as an alternative to the analysis of actual structures, this would solve a lot of theoretical analysis and practical engineering needed to address complex problems that cannot be resolved. Key words: finite element analysis structural calculation physical design 1 有限元的发展历程 有限元法的发展历程可以分为提出(1943)、发展(1944一1960)和完善(1961-二十世纪九十年代)三个阶段。有限元法是受内外动力的综合作用而产生的。 1943年,柯朗发表的数学论文《平衡和振动问题的变分解法》和阿格瑞斯在工程学中取得的重大突破标志着有限元法的诞生。 有限元法早期(1944一1960)发展阶段中,得出了有限元法的原始代数表达形式,开始了对单元划分、单元类型选择的研究,并且在解的收敛性研究上取得了很大突破。1960年,克劳夫第一次提出了“有限元法”这个名称,标志着有限元法早期发展阶段的结束。 有限元法完善阶段(1961一二十世纪九十年代)的发展有国外和国内两条线索。在国外的发展表现为: 第一,建立了严格的数学和工程学基础;第二,应用范围扩展到了结构力学以外的领域;第三,收敛性得到了进一步研究,形成了系统的误差估计理论;第四,发展起了相应的商业软件包。 近年来随着计算机技术的普及和计算速度的不断提高,有限元分析在工程设计和分析中得到了越来越广泛的重视,已经成为解决复杂的工程分析计算问题的有效途径,现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算,其在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器,国防军工,船舶,铁道,石化,能源,科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃,主要表现在以下几个方面: 一、增加产品和工程的可靠性; 二、在产品的设计阶段发现潜在的问题 三、经过分析计算,采用优化设计方案,降低原材料成本

非线性有限元分析

非线性有限元分析 1 概述 在科学技术领域,对于许多力学问题和物理问题,人们已经得到了它们所应遵循的基本方程(常微分方程或偏微分方程)和相应的定解条件(边界条件)。但能够用解析方法求出精确解的只是少数方程性质比较简单,并且几何形状相当规则的问题。对于大多数工程实际问题,由于方程的某些特征的非线性性质,或由于求解区域的几何形状比较复杂,则不能得到解析的答案。这类问题的解决通常有两种途径。一是引入简化假设,将方程和几何边界简化为能够处理的情况,从而得到问题在简化状态下的解答。但是这种方法只是在有限的情况下是可行的,因为过多的简化可能导致误差很大甚至是错误的解答。因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。特别是五十多年来,随着电子计算机的飞速发展和广泛应用,数值分析方法已成为求解科学技术问题的主要工具。 已经发展的数值分析方法可以分为两大类。一类以有限差分法为代表,主要特点是直接求解基本方程和相应定解条件的近似解。其具体解法是将求解区域划分为网格,然后在网格的结点上用差分方程来近似微分方程,当采用较多结点时,近似解的精度可以得到改善。但是当用于求解几何形状复杂的问题时,有限差分法的精度将降低,甚至发生困难。 另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法,然后再建立近似解法并求解。如果原问题的方程具有某些特定的性质,则它的等效积分提法可以归结为某个泛函的变分,相应的近似解法实际上就是求解泛函的驻值问题。诸如里兹法,配点法,最小二乘法,伽辽金法,力矩法等都属于这一类方法。但此类方法也只能局限于几何形状规则的问题,原因在于它们都是在整个求解区域上假设近似函数,因此,对于几何形状复杂的问题,不可能建立合乎要求的近似函数。 1960年,R.W.CLOUGH发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”,这同时也标志着有限单元法(FEM)的问世。有限单元法的基本思想是将连续的求解区域离散为一组有限个,且按一定方式相互联接在一起的单元的组合体。由于单元能按不同的联结方式进行组合,且单元本身又可以有不同形状,因此可以模型化几何形状复杂的求解域。并且可以利用在每一个单元假设的近似函数来分片地表示全求解域上待求的未知场函数,从而使一个连续的无限自由度问题变成离散的有限自由度问题。 现已证明,有限单元法是基于变分原理的里兹法的另一种形式,从而使里兹法分析的所有理论基础都适用于有限单元法,确认了有限单元法是处理连续介质问题的一种普遍方法。利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的,而且事先不要求满足任何边界条件,因此可以用来处理很复杂的连续介质问题。 在短短四十余年的时间里,有限单元的分析方法已经迅速地发展为适合于使用各种类型计算机解决复杂工程问题的一种相当普及的方法。如今,有限元广泛地应用于各个学科门类,已经成为工程师和科研人员用于解决实际工程问题,进行科学研究不可或缺的有力工具。有限单元法的应用围已由弹性力学平面问题扩展到空间问题,板壳问题,由静力平衡问题扩展到稳定问题,动力问题和波动问题。分析的对象从弹性材料扩展到塑性,粘弹性,粘塑性和复合材料等,从固体

相关主题