有限元模型修正中若干重要问题

有限元分析可能会出现什么错误

图一应变能随单元尺寸变化图有限元分析可能会出现什么错误？当有限元分析扩散向可能没有正式数字化程序培训的设计者的时候，专业人员必须问“最适当的方法是否被采用？这些方法是否产生了精确的结果？” 在今天的设计领域，有限元方法被广泛的应用，其中包括各种各样的通用商业软件和适合某专业领域的专业软件。这些方法日益增长的被用，在帮助确定好的新的设计的同时改良了设计性能和成本。考虑到有限元方法在各个设计领域起着重要的作用，专业人员需要问他们自己“他们的设计程序是否是可获得的最适合的技术？这些方法是否会产生精确的结果？”。这些问题是非常重要的，因为越来越多的设计人员不见得受过数字化程序培训，但是他们却在他们的工作中应用有限元方法。右图是应变能量随单元尺寸减小的变化示意图。从图中可以看出应变能会随着单元尺寸的减小而收敛。当这些有限元方法被向越来越多且越来越广的群体广泛的应用的时候，用户必须问有限元分析会出现什么样的错误。本文目的不是在广义上解决这个问题，更恰当的说，我们必须集中焦点于有限元方法的可靠性和准确应用方面。为了便于说明，我们考虑线弹性问题，假设有限元的代数方程精确地被求解。对于复杂的分析，考虑这些条件的同时，还有一些额外的要求也有必要得到。数学模型首先，设计人员应该记住有限元方法是为了求解数学模型，这数学模型是实际物理问题的理想化结果。数学模型是建立在考虑几何、材料特性、加载条件和位移边界条件等假设的基础上的。数学模型的指导方程是考虑到边界条件的偏微分方程。这些方程不能用封闭的解析方式求解，因此，设计人员要借助有限元方法获得一个数值解。例如，考虑一个几何和载荷为轴对称的阀套。在这种条件下，考虑轴对称分析条件是合理的。分析的数学模型可以通过指定几何尺寸、支撑条件、材料常数和加载条件来获得。虽然通常情况下设计人员不能用封闭方程的数学模型的精确解，但是这个数学模型的精确解是存在的，且是唯一的。高精度的精确解的近似解可以用有限元方法求得。为了充分理解这些观察到的，必须要有收敛概念。这里E 表示应变能的数学模型的精确解（未知的），Eh 表示对应于单元尺寸h 的应变能的有限元解。那么收敛表示为：上述图一中的示意图表明了收敛是如何达到。当单元网格趋于精确时（这就意味着单元尺寸h 的减少），应变能Eh 将趋向于数值E 。E 和Eh 之间的误差值的减少速度视解答的题目，也依靠采用单元类型和网格质量。明显地，网格细化过程中高阶单元减少误差率比低阶单元快。可靠性问题有限元方法的可靠性是指，在提出很好的数学模型求解时，有限元程序有两个特性。第一，在任何材料特性、位移边界条件和加载条件下，当单元尺寸h 趋向0时，有限元的结

有限元分析报告样本

《有限元分析》报告基本要求： 1. 以个人为单位完成有限元分析计算，并将计算结果上交；（不允许出现相同的分析模型，如相同两人均为不及格） 2. 以个人为单位撰写计算分析报告； 3. 按下列模板格式完成分析报告； 4. 计算结果要求提交电子版，报告要求提交电子版和纸质版。（以上文字在报告中可删除）《有限元分析》报告一、问题描述（要求：应结合图对问题进行详细描述，同时应清楚阐述所研究问题的受力状况和约束情况。图应清楚、明晰，且有必要的尺寸数据。）一个平面刚架右端固定，在左端施加一个y 方向的-3000N 的力P1，中间施加一个Y 方向的-1000N 的力P2，试以静力来分析，求解各接点的位移。已知组成刚架的各梁除梁长外，其余的几何特性相同。横截面积：A=0.0072 m2 横截高度：H=0.42m 惯性矩：I=0.0021028m4ｘ弹性模量： E=2.06ｘ10n/ m2/ 泊松比：u=0.3 二、数学模型（要求：针对问题描述给出相应的数学模型，应包含示意图，示意图中应有必要的尺寸数据；如进行了简化等处理，此处还应给出文字说明。）（此图仅为例题）

三、有限元建模（具体步骤以自己实际分析过程为主，需截图操作过程）用ANSYS 分析平面刚架 1．设定分析模块选择菜单路径：MainMenu—preference 弹出“PRreferences for GUI Filtering”对话框，如图示，在对话框中选取：Structural”,单击[OK]按钮，完成选择。 2．选择单元类型并定义单元的实常数（1）新建单元类型并定（2）定义单元的实常数在”Real Constants for BEAM3”对话框的AREA中输入“0。0072”在IZZ 中输入“0。0002108”，在HEIGHT中输入“0.42”。其他的3个常数不定义。单击[OK]按钮，完成选择 3．定义材料属性在”Define Material Model Behavier”对话框的”Material Models Available”中，依次双击“Structural→Linear→Elastic→Isotropic”如图

有限元模型如何查错

有限元模型如何查错作者：PAUL KUROWSKI 在建立有限元模型的过程中很容易出错，如果你知道如何查错，修正这些错误将会变得很简单有限元分析的第一步就是建立被分析对象的数学模型，这要求我们思索建模的理论基础如弹性理论，板的Reissner理论，塑性变形理论等，和考虑问题的其它信息如几何描述、材料特性，约束和荷载等等。分析的目的就是由这些条件，计算得到精确解u_EX并同时得到位移u_EX的应力函数 F(u_EX)如Von Mises应力等。应力函数F (u_EX) 仅仅依赖于数学模型的定义，而与求解该数学问题的数值近似计算方法无关；同时应力函数F(u_EX)也不依赖于网格划分、网格类型和单元尺寸。函数F(u_EX)与模型实体物理性质之间的差异，被称为“模型错误”。下一步就是使用有限元方法去找到精确解u_EX的近似值u_FE。这个过程包括选择网格划分和构件类型，如对二维板用八节点（矩形）单元，依此类推。网格划分＆单元定义被称为有限元的离散化。离散化产生的误差可以被定义为：大部分的分析应该把这个误差控制在10％以内。同时由于建立模型和模型的离散化一定会产生这个误差，正确运用有限元分析就包括对这两类误差进行评估和控制。有限元分析结果中的名义误差＆真实误差是有区别的，最好能够加以区别：名义误差可以比建模误差和离散误差的总和小，二者可能反号而相互抵消。结果的好坏取决于模型是否反应实际（模型误差的大小）和有限元软件在转化过程中的精度控制（离散

化误差的大小）。 WHAT IS MODELING ERROR? 何为模型误差？假设要分析一个支架，我们首先考虑到的问题应该包括：我们想得到什么结果？是最大应力还是最大变形？是固有频率、弯曲刚度、还是温度分布？支架是否处于弹性变形阶段？极限荷载形式有几种情况？如何模拟支撑条件等等。有了一个明确的目标和对我们使用的理论自身局限性的把握，分析者就可以建立模型了。有时这个模型与CAD模型是相似的，但相当多的情况是，为了简化网格的划分，我们有必要修改模型的拓扑描述。部分建模的过程包括以下一些问题：用壳单元模拟薄壁墙体，对对称性、反对称性或两者的运用，是否考虑细部及忽略不重要的特征等。比如，选用壳单元而不用实体单元意味着我们考虑到数学模型和相应的有限元软件的运作方式而作出了一个重要的决定。当(研究对象的)拓扑描述已经比较理想后，我们还需要对材料属性(选择线弹性、弹塑性或其他)、荷载及支撑条件进行理想的简化。我们认为这些简化精确反应了所需模型的重要数据，而建模当中的一些重要决策有时并未过多的考虑这些(方面)。简化了的模型经常是概念错误的，一个检验模型是否不合理的方法是其解析解对应的应变能是否无穷大或趋近于零；另一个方法是对应于数学模型的我们感兴趣的数据在结果没有得到体现。很多分析者认为一个有效的网格生成器可以生成高质量的网格并降低模型误差，其实不尽然，模型是在网格划分前假定的，因此，最合理的网格划分也无法修正一个简化不合理的数学模型。 A SYSTEMATIC APPROACH 一个系统的方法确保模型误差较小的唯一方式是把所需研究的数据放在对模型假设不敏感之处。类似地，通过把所需研究的数据放在对离散不敏感之处（不敏感的表现是：结果对更细的网格划分或更大的p值并不发生明显的改变），以减少离散误差。举个例子：比如说我们对一块简支板沿着边缘方向的剪力感兴趣，那么经典的克西霍夫板模型(Kirchhoff’s plate)是不可用的，可以通过一个Reissner模型或一个全3D的弹性模型轻而易举地检验出来。一个关于板弯曲的Reissner模型假设所有平面内位移沿厚度方向呈线性变化、剪应变沿厚度方向保持不变。若采用更厚的板的话会迫使人们去置疑简支的意义、同时会置疑是否可以给出一

GARTEUR 有限元模型修正与确认研究

收稿日期:2003207207;修订日期:2004203225 基金项目:教育部博士学科点专项基金(20010227012)资助项目文章编号:100026893(2004)0420372204 GARTEUR 有限元模型修正与确认研究费庆国,张令弥,郭勤涛 (南京航空航天大学振动工程研究所,江苏南京 210016) Case Study of FE Model Updating and Validation via an Air craft Model Structur e FEI Qing 2guo,Z HANG Ling 2mi,G UO Qin 2tao (Institu te o f Vi brati on Engi neering,Nanjing University of Aeronau tics and Astro nautics,Nanjing 210016,China)摘要:待修正参数的选择以及修正后模型的质量评估是有限元模型修正的两个重要问题。以欧洲学术界广泛采用的GA RTEUR 飞机模型为例,利用基于灵敏度分析的模型修正方法,通过仿真算例研究参数选择对模型修正质量的影响,并以试验数据为目标值对有限元模型进行修正与确认。为全面评估模型的修正质量,引入三级标准对修正后有限元模型进行确认。关键词:固体力学;模型确认;有限元法;模型修正;参数选择中图分类号:O 248121 文献标识码:A Abstr act:Parameter selection and quali ty validation are of g reat i mpo rtance in fini te element model updating.This paper presents so me results which demonstrate the relationship betw een parameter selection and updated model .s quality throu gh si mulation cases.Three q uali ty levels w ith corresponding validation criteria are emplo yed with an emphasis o n updated mod 2el .s predictio n ability.Results of updating based on exper i mental modal test data are sho w n as an application example.A n aircraft test structure,GA RTEUR,which is g enerally utilized in Europe,is employed in bo th the si mulation case and the exper i mental case.Sensi tivity 2based model updating appro ach is applied. Key wor ds:solid mechanics;model validation;finite element method;model updating;parameter selection 在航空工程中,准确的有限元模型对于动态响应预测以及动态设计至关重要。建模过程中的不确定因素,如离散化误差、材料物理参数的不确定性、边界条件的近似等,导致有限元模型必然存在误差。设计规范规定,有限元模型必须通过振动模态试验或者地面共振试验来检验[1]。近30年来,有限元模型修正技术得到了长足的发展[2~6]。根据修正对象的不同可将修正方法分为矩阵型方法和设计参数型方法。后者物理意义明确,更具工程应用价值。本文采用基于灵敏度分析的设计参数型修正方法。基于灵敏度分析的设计参数型修正方法主要包括待修正设计参数选择,灵敏度分析,参数修正以及模型确认等环节。待修正设计参数的选择是模型修正的起始环节。通常,候选参数是有限元模型存在不确定性因素的参数。近20年虽然发展了很多种参数选择或者误差定位的策略与算法,工程应用中仍然难以准确无遗漏地确定误差参数。因此,有必要讨论参数选择对模型修正质量的影响。模型确认是模型修正的检验环节。在当前的研究及工程应用中,通常只要求修正后模型的计算结果能够复现修正过程中利用的试验数据。事实上,为全面评估模型的质量,模型的复现能力与预测能力应予以同等重视[7]。本研究引入了三级质量标准对修正后的有限元模型进行确认。本文采用G ARTE UR 飞机模型为研究对象,通过仿真算例来研究参数选择对模型修正质量的影响,并给出了利用振动模态测试结果对G AR 2TEUR 飞机模型的有限元模型进行修正与确认的结果。 1 模型修正方法与模型确认准则 (1)模型修正方法模型修正可归结为以下的优化问题[8] Min p +R(p )+2 2,R(p )=f E -f A (p )s.t V L [p [V U (1) 其中:p 代表设计参数;f E ,f A 是结构动态特性试验与分析结果;R 代表残差;V L ,V U 是设计参数的下、上限。令设计参数的初始值为p 0,动态特性f 是设计参数的隐函数,其泰勒展开式为第25卷第4期航空学报 Vol 125N o 14 2004年 7月ACT A AERO NA U TICA E T AS TRO NA U TICA SINICA July 2004

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

板结构有限元分析实例详解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对话框

汽车车架简化模型有限元分析

汽车典型零部件简化模型有限元分析任务1：连杆简化模型的有限元分析 1. 分析任务：对图一所示的连杆的二维简化模型进行有限元分析，确定该设计是否满足结构的强度要求；若强度不够，修改设计直至最大应力减小至材料允许的范围内。在修改结构时，注意不可改变连杆小头衬套的内径和连杆大头的内径，也不可改变连杆各处的厚度和材料。 2. 分析所需数据： a．连杆采用两种材料，连杆本体用的是40Cr结构钢，左侧小头中的衬套用的是铜。 b．连杆杆身和大头的厚度为1.5mm，小头的厚度为3.0mm。注意在杆身和小头的过渡处有R2.0的过渡圆角； c．连杆结构的其它尺寸如图二所示； d．施加在大、小头内壁上的边界条件用于模拟连杆与曲轴及活塞销的连接。假定载荷分布在小头夹角为90o的内壁上，且为锥状分布；约束施加在连杆大头夹角为90o的内壁上； e．40Cr材料的弹性模量：210GPa；泊松比：0.3；屈服极限为：850MPa，设计安全系数为6；铜的弹性模量：120GPa，泊松比：0.33；屈服极限为：250MPa；设计安全系数为4。 3. 完成该分析应掌握的ANSYS技术： a．单元类型的选择；单元的尺寸控制； b．不同厚度和材料的二维实体建模； c．工作平面的灵活应用；

d．按载荷和约束的要求分割线和面； e．模型参数（材料，实常数，单元类型号等） f．粘结、合并等布尔运算操作 g．局部坐标系，旋转节点坐标系； h．线性分布载荷的施加； i．单元网格误差估计； j．Ansys 命令日志文件及其在修改设计中的应用； k．多窗口显示的功能 4. 分析报告内容的基本要求： a.对分析任务的描述；列出分析所需数据： b.利用多窗口显示的功能绘出连杆的实体模型和网格模型，在模型上能反映出连杆各部位材料、厚度的不同； c.绘图反映连杆的边界条件； d.绘出对连杆原设计进行有限元分析后得到的变形图和应力等值线图； e.图示SEPC和SERR并说明有限元分析的建模误差； f.详细说明对不符合设计要求的结构所作的设计修改；及最终符合设计要求的计算结果； g.在分析中遇到的关键问题（在实体建模、网格剖分、边界条件施加等各个步骤中出现的）及解决的办法； h.整理命令日志文件，并在每个语句后添加说明（说明该语句的功能，说明前要加！号）。注意：添加的说明（可以用中文说明）应该反映在建模中的操作步骤而不是简单的ANSYS命令定义。

有限元分析中的一些问题

有限元分析的一些基本考虑-—-—-单元形状对于计算精度的影响笔者发现,在分析复杂问题时,我们所可能出现的错误，竟然是一些很根本的错误,这些根本错误是由于对有限元的基本理论理解不清晰而造成的。鉴于这个原因，笔者决定对一些基本问题(例如单元形状问题,单元大小问题,应力集中问题等)展开调查,从而形成了一系列文章，本篇文章是这些系列文章中的第一篇. 本篇文章先考虑有限元分析中的第一个基本问题：单元形状问题。我们知道,单元形状对于有限元分析的结果精度有着重要影响,而对单元形状的衡量又有着诸多指标,为便于探讨,这里首先只讨论第一个最基本的指标：长宽比(四边形单元的最长尺度与最短尺度之比),而且仅考虑平面单元的长宽比对于计算精度的影响。为此,我们给出一个成熟的算例。该算例是一根悬臂梁,在其端面施加竖直向下的抛物线分布载荷,我们现在考察用不同尺度的单元划分该梁时,对于A点位移的影响。这五种不同的划分方式，都使用矩形单元,只不过各单元的长宽比不同。例如第一种(1)AＲ＝1．1,就是长宽比接近1；第二种（2)AＲ=1．5，就是长宽比是1。5.其它类推。第五种(5）AＲ=24，此时单元的长度是宽度的2４倍。现在我们看看按照这五种单元划分方式对于A点位移的影响，顺便我们也算出了B点的位移,结果见下表.

我们现在仔细查看一下上表，并分析其含义。我们先考虑第一行，它是第一种单元划分情况,此时每个单元的长宽比是１。１，由此我们计算出A点,B点的垂直位移，可以看到,A点的竖直位移是—1.093英寸，而B点的竖直位移是-0。34６英寸。而这两点我们都是可以用弹性力学的方式得到精确解的,其精确解分别是-1。１52以及—０。360。这样，我们可以得到此时Ａ点位移误差的百分比是[(—1.0９３)—(-1。152)]／1。1５2 =5。2%．对于其它情况,也采用类似的方式得到A点位移误差的百分比。从上表可以看出来，随着长宽比的增加,位移误差越来越大，竟然大到５6％.因此，如果我们是用长宽比为24的单元进行划分的话，那么我们的结果可以说是完全错误的. 下面按照上表绘制出一张图,该图从形象的角度表达了上表的含义.

结构有限元及其应用软件

中国海洋大学本科生课程大纲课程属性：公共基础/通识教育/学科基础/专业知识/工作技能，课程性质：必修、选修一、课程介绍 1.课程描述（中英文）：本课程是一门重要的结构计算分析课程，通过多媒体教学和上机练习，系统学习结构有限元FEM的基本原理和方法，熟悉掌握通用有限元应用软件ANSYS进行结构静力和动力分析的方法和步骤，并初步掌握使用ANSYS进行海工典型结构强度计算的方法。 Structural finite element method and its application software is an important course of structural calculation and analysis. Through multimedia teaching and computer practice, the basic principles and methods of Finite Element Method (FEM) are learned systematically. The general finite element application software ANSYS for the methods and procedures of structural static and dynamic analysis are mastered.At the same time, the strength calculation method of typical ocean engineering structures using ANSYS is preliminarily mastered. 2.设计思路：有限元方法是一种现代设计方法，应用于结构设计中，是一种具有重要经济意义和巨大潜力的先进结构设计技术。因此选择该课程作为结构设计方面的一门必修课程，主要介绍结构有限元的基本原理和方法，还选择了通用的有限元软件ANSYS进行示例分析。包括要求掌握有限元法的基本思想和基本原理、平面刚架结构的有限元法、弹

机械结构有限元分析

机械结构有限元分析有限元分析软件ANSYS在机械设计中的应用摘要：在机械设计中运用ANSYS软件进行有限元分析是今后机械设计发展的必然趋势，将有限元方法引入到机械设计课程教学中，让学生参与如何用有限元法来求解一些典型零件的应力，并将有限元结果与教材上的理论结果进行对照。这种新的教学方法可以大大提高学生的学习兴趣，增强学生对专业知识的理解和掌握，同时还可以培养学生的动手能力。在机械设计课程教学中具有很强的实用价值。关键词：机械设计有限元 Ansys 前言：机械设计课程是一门专业基础课，其中很多教学内容都涉及到如何求取零件的应力问题，比如齿轮、v带、螺栓等零件。在传统的教学过程中，都是根据零件的具体受力情况按材料力学中相应的计算公式来求解。比如，在求解齿轮的接触应力时，是把齿轮啮合转化为两圆柱体的接触，再用公式求解。这些公式本身就比较复杂，还要引入各种修正参数，因此我们在学习这些内容时普遍反映公式难记，学习起来枯燥乏味，而且很吃力。近年来有限元法在结构分析中应用越来越广泛，因此如果能将这种方法运用到机械设计课程中，求解一些典型零件的应力应变，并将分析结果和教材上的理论结果进行对比，那么无论是对于提高学生学习的热情和积极性，增强对重点、难点知识的理解程度，还是加强学生的计算机水平都是一件非常有益的事情。由于直齿圆柱齿轮的接触强度计算是机械设计课程中的一个重要内容，齿轮强度的计算也是课程中工作量最繁琐的部分。下面就以渐开线直齿圆柱齿轮的齿根弯曲疲劳强度的计算为例，探讨在机械设计课程中用ANSYS软件进行计算机辅助教学的步骤和方法，简述如何将有限元方法应用到这门课程的教学中。 1．传统的直齿圆柱齿轮齿根弯曲疲劳强度的计算传统方法把轮齿看作宽度为b的矩形截面的悬臂梁。因此齿根处为危险剖面，它可用30。切线法确定。如图l所示。作与轮齿对称中心线成30。角并与齿根过渡曲线相切的切线，通过两切点作平行与齿轮轴线的剖面，即齿根危险剖面。理论上载荷应由同时啮合的多对齿分担，但为简化计算，通常假设全部载荷作用于齿顶来进行分析，另用重合度系数E对齿根弯曲应力予以修正。由材料力学弯曲应力计算方法求得齿根最大弯曲应力为：

有限元分析中的一些问题

有限元分析的一些基本考虑-----单元形状对于计算精度的影响笔者发现，在分析复杂问题时，我们所可能出现的错误，竟然是一些很根本的错误，这些根本错误是由于对有限元的基本理论理解不清晰而造成的。鉴于这个原因，笔者决定对一些基本问题（例如单元形状问题，单元大小问题，应力集中问题等）展开调查，从而形成了一系列文章，本篇文章是这些系列文章中的第一篇。本篇文章先考虑有限元分析中的第一个基本问题：单元形状问题。我们知道，单元形状对于有限元分析的结果精度有着重要影响，而对单元形状的衡量又有着诸多指标，为便于探讨，这里首先只讨论第一个最基本的指标：长宽比（四边形单元的最长尺度与最短尺度之比），而且仅考虑平面单元的长宽比对于计算精度的影响。为此，我们给出一个成熟的算例。该算例是一根悬臂梁，在其端面施加竖直向下的抛物线分布载荷，我们现在考察用不同尺度的单元划分该梁时，对于A点位移的影响。这五种不同的划分方式，都使用矩形单元，只不过各单元的长宽比不同。例如第一种（1）AR=1.1，就是长宽比接近1；第二种（2）AR=1.5,就是长宽比是1.5.其它类推。第五种（5）AR=24，此时单元的长度是宽度的24倍。现在我们看看按照这五种单元划分方式对于A点位移的影响，顺便我们也算出了B点的位移，结果见下表。

我们现在仔细查看一下上表，并分析其含义。我们先考虑第一行，它是第一种单元划分情况，此时每个单元的长宽比是1.1，由此我们计算出A点，B点的垂直位移，可以看到，A点的竖直位移是-1.093英寸，而B点的竖直位移是-0.346英寸。而这两点我们都是可以用弹性力学的方式得到精确解的，其精确解分别是-1.152以及-0.360.这样，我们可以得到此时A点位移误差的百分比是 [(-1.093)-(-1.152)]/1.152 = 5.2%. 对于其它情况，也采用类似的方式得到A点位移误差的百分比。从上表可以看出来，随着长宽比的增加，位移误差越来越大，竟然大到56%。因此，如果我们是用长宽比为24的单元进行划分的话，那么我们的结果可以说是完全错误的。下面按照上表绘制出一张图，该图从形象的角度表达了上表的含义。由此可见，长宽比越接近于1，那么结算结果越精确，越远离1，则误差越大。

结构有限元分析的形状处理方法_杜平安

结构有限元分析的形状处理方法杜平安摘要　介绍结构形状处理的各种方法,包括类型简化、细节简化、形式变换、局部结构和利用对称性等。关键词　形状处理　有限元分析　建模 Abstract The processing method is intro-duced in the paper ,including ty pe simplifica tion 、details simplifica tio n 、fo rm tra nsfo rmatio n 、local structure a nd symm etry utiliza tion . Key words Shape processing Finite element analysis Modelling 收稿日期:1999-08-18 1　结构类型简化根据结构形状、载荷和约束条件的特点,结构类型可分为空间问题、平面问题、轴对称问题、板壳问题和杆件问题等。其中平面问题和轴对称问题的几何模型是一平面图形,在平面上划分网格比在空间内划分要容易得多,单元数量也少得多。因此将空间问题作适当近似,使其按平面问题来处理,则可使分析过程大为简化。在图1a 中,计算轮毂与轴过盈配合的接触压力时,由于辐孔尺寸较小且远离接触面,因此可以不考虑辐孔而将轮毂简化为轴对称结构。同样,在计算图1b 中螺栓与螺母螺纹面上的接触压力时,由于螺旋升角较小,也可以不考虑升角的影响,而将螺栓与螺母简化为轴对称结构。图1　结构类型简化结构 2　结构细节简化细节是结构中相对尺寸很小的局部,如倒圆、倒角、退刀槽和加工凸台等。根据网格划分特点,一条直线或曲线至少要划分一个单元边;一个平面或曲面至少要划分一个单元面;一个圆至少要用三个单元边离散,因此几何模型中的细节将限制细节处及其附近的网格大小,从而影响整个结构的网格分布和增加网格数量。图2是有无细节时自动划分出的网格,从中可以看出细节对网格划分的影响。图2　细节对网格划分的影响因此,建立几何模型时应尽量忽略一些不必要的细节。在静力分析中,高应力区域中的细节会引起应力集中,细节大小和形状对应力影响很大,这些细节不能忽略。而处于结构低应力区的细节一般可以忽略。在动力计算中,由于结构固有频率和模态振型主要取决于结构的质量分布和刚度,因此细节一般可以忽略。在热分析中,细节不会在结构中引起局部高温,这时也可以考虑较少的细节。 3　结构形式变换有些结构尽管形状不是很复杂,但划分网格却很困难。如果对结构形式作适当变换,则可使网格划分变得容易,划分出的单元更少。例如图3a 所示的带肋板,划分网格时需要用板单元和梁单元组合,且两类单元为偏心连接,自动分网难以满足这种要求。如果将带肋板变换为平板(图3b 所示),则在平板上划分网格要容易得多。由于带肋板用于焊接而成支撑箱式立柱,其特性要求主要是刚度,因此可按等刚度条件作为变换 · 26·《机械与电子》2000(1)

有限元+螺栓简化

1 概述螺栓是机载设备设计中常用的联接件之一。其具有结构简单,拆装方便,调整容易等优点,被广泛应用于航空、航天、汽车以及各种工程结构之中。在航空机载环境下，由于振动冲击的影响，设备往往产生较大的过载，对作为紧固件的螺栓带来强度高要求。螺栓是否满足强度要求，关系到机载设备的稳定性和安全性。传统力学的解析方法对螺栓进行强度校核，主要是运用力的分解和平移原理，解力学平衡方程，借助理论和经验公式，理想化和公式化。没有考虑到连接部件整体性、力的传递途径、部件的局部细节(如应力集中、应力分布)等等。通过有限元法，整体建模，局部细化，可以弥补传统力学解析的缺陷。用有限元分析软件MSC.Patran/MSC.Nastran提供的特殊单元来模拟螺栓连接，过程更方便，计算更精确，结果更可靠。因此，有限元在螺栓强度校核中的应用越来越广泛。 2 有限元模型的建立对于螺栓的模拟，有多种模拟方法，如多点约束单元法和梁元法等。多点约束单元法(MPC)即采用特殊单元RBE2来模拟螺栓连接。在螺栓连接处，设置其中一节点为从节点(Dependent)，另外一个节点为主节点(Independent)。主从节点之间位移约束关系使得从节点跟随主节点位移变化。比例因子选为1,使从节点和主节点位移变化协调一致，从而模拟实际工作状态下，螺栓对法兰的连接紧固作用。梁元法模拟即采用两节点梁单元Beam，其能承受拉伸、剪切、扭转。通过参数设置，使梁元与螺栓几何属性一致。本文分别用算例来说明这两种方法的可行性。 2.1 几何模型如图1所示组合装配体，底部约束。两圆筒连接法兰通过8颗螺栓固定。端面受联合载荷作用。

abaqus有限元建模小例子

问题一: 工字梁弯曲 1.1 问题描述: 在<<材料力学实验>>中，弯曲实验測定了工字梁弯曲应变大小及其分布，以验证弯曲正应力公式。在这里，採用ABAQUS/CAE建立试验件的有限元模型，ABAQUS/Standard模块进行分析求解，得到应力、应变分布，对比其与理论公式计算值及实验測量值的差別。弯曲实验的相关数据：材料：铝合金E=70GPa 泊松比0.3 实验装置结构简图如图所示：结构尺寸测量值：H=50（+/-0.5mm） h=46(+/-0.5mm) B=40（+/-0.5mm） b=2(+/-0.02mm) a=300(+/-1mm) F1=30N Fmax=300N N ? F100 = 1.2 ABAQUS有限元建模及分析一对象: 工字型截面铝合金梁梁的结构简图如图1所示，結构尺寸、载荷、約束根据1.1设定，L取1600mm，两端各伸出100mm。二用ABAQUS/CAE建立实验件的有限元模型,效果图如下：边界条件简化：左侧固定铰支座简化为下表面左参考点处的约束U1=U2=U3=0

右侧活动铰支座简化为下表面右参考点处的约束U1=U2=UR3=0 几何模型

有限元模型三ABAQUS有限元分析結果 ①应力云图（Z方向正应力分量）：施加载荷前 F=300N

②应变（Z方向分量）：中间竖直平面的厚度方向应变分布图： F=100N F=200N

F=300N 由上图可以看出应变沿着厚度方向呈线性比例趋势变化，与实验测得的应变值变化趋势相同。中性轴处应变均接近零值，应变与距离中性轴位移基本为正比关系。 1.3分析结果: 中间竖直截面上下边缘轴向应力数值对比：*10^-6 MPa 距中性轴距ABAQUS模拟实验测量值平均理论值 1/2H -96.182*70000 -97*70000 -6.9165=-70000*98.807 -1/2H 95.789*70000 92*70000 6.9165

ANSYS有限元分析与实体建模

第五章实体建模 5.1实体建模操作概述用直接生成的方法构造复杂的有限元模型费时费力，使用实体建模的方法就是要减轻这部分工作量。我们先简要地讨论一下使用实体建模和网格划分操作的功能是怎样加速有限元分析的建模过程。自下向上地模造有限元模型：定义有限元模型顶点的关键点是实体模型中最低级的图元。在构造实体模型时，首先定义关键点，再利用这些关键点定义较高级的实体图元（即线、面和体）。这就是所谓的自下向上的建模方法。一定要牢记的是自下向上构造的有限元模型是在当前激活的坐标系内定义的。图5-1自下向上构造模型自上向下构造有限元模型：ANSYS程序允许通过汇集线、面、体等几何体素的方法构造模型。当生成一种体素时，ANSYS程序会自动生成所有从属于该体素的较低级图元。这种一开始就从较高级的实体图元构造模型的方法就是所谓的自上向下的建模方法。用户可以根据需要自由地组合自下向上和自上向下的建模技术。注意几何体素是在工作平面内创建的，而自下向上的建模技术是在激活的坐标系上定义的。如果用户混合使用这两种技术，那么应该考虑使用CSYS，WP或CSYS，4命令强迫坐标系跟随工作平面变化。图5-2自上向下构造模型（几何体素）注意：建议不要在环坐标系中进行实体建模操作，因为会生成用户不想要的面或体。

运用布尔运算：可以使用求交、相减或其它的布尔运算雕塑实体模型。通过布尔运算用户可直接用较高级的图元生成复杂的形体。布尔运算对于通过自下向上或自上向下方法生成的图元均有效。图5-3使用布尔运算生成复杂形体。拖拉或旋转：布尔运算尽管很方便，但一般需耗费较多的计算时间。故在构造模型时，如果用拖拉或旋转的方法建模，往往可以节省计算时间，提高效率。图5-4拖拉一个面生成一个体〔VDRAG〕移动和拷贝实体模型图元：一个复杂的面或体在模型中重复出现时仅需要构造一次。之后可以移动、旋转或拷贝到所需的地方。用户会发现在方便之处生成几何体素再将其移动到所需之处，这样往往比直接改变工作平面生成所需体素更方便。图5-5拷贝一个面网格划分：实体建模的最终目的是为了划分网格以生成节点和单元。在完成了实体建模和建立了单元属性，网格划分控制之后，ANSYS程序可以轻松地生成有限元网格。考虑到要满足特定的要求，用户可以请求映射网格划分生成全部都是四边形、三角形或块单元。

有限元分析步骤

有限元建模与分析有限元分析（FEA）是一种预测结构的偏移与其它应力影响的过程，有限元建模（FEM）将这个结构分割成单元网格以形成实际结构的模型，每个单元具有简单形态（如正方形或三角形）。这样有限元程序就有了可写出在刚度矩阵结构中控制方程方面的信息。每个单元上的未知量就是在节点上的位移，这个点就是单元元的连接点。有限元程序将这些单个单元的刚度矩阵组合起来以形成整个模型的总刚度矩阵，并给予已知力和边界条件来求解该刚度矩阵以得出未知位移，从节点上位移的变化就可以计算出每个单元中的应力。有限单元由假定的应变方程式导出，有些单元可假设其应变是常量，而另外一些可采用更高阶的函数。利用给定单元的这些方程和实际几何体，则可以写出外力和节点位移之间的平衡方程。对于单元的每个节点来说，每个自由度就有一个方程，这些方程被十分便利地写成矩阵的形式以用于计算机的演算中，这个系数的矩阵就变成了一个显示出力对位移的关系的刚度矩阵：{F}=[K]、{d} 尽管求知量处于离散的自由度，内部方程仍被写成表述为连续集的应变函数。这就意味着如果选择了正确单元的话，纵然这个有限元模型有一组离散的方程，只要用有限的节点和单元也可以收敛出正确的答案。有限元模型是解决全部结构问题的完全理想的模型。这些问题包括节点的定位，单元，物理的和材料的特性，载荷和边界条件，根据分析类型的不同，如静态结构载荷，动态的或热力分析，这个模型就确定得不同。一个有限元模型常常由不止一种单元类型来建立，有限元模型是以结构的偏移来建立成数学模型，而不只是在外观上象原结构。也许某个零件用梁单元最好，而另外的零件则可能用薄壳单元最理想。对于给定的问题来讲，求解结果的准确性将取决于结构建模的好坏，负载和边界条件的确定，以及所用单元的精度。一般来讲，如模型细分更小的单元，则求解将更准确。了解你在最终的求解结果上有充分收敛的唯一确信的方法是用更细网格的单元来建立更多的模型，以检查求解结果的收敛性。新的有限元用户经常产生想象上的错误，即建立一个有限元模型的目的是建立一个看起来象这种结构的模型。有限元建模的目的是建立一个从数学意义是“相似”的模型，而不是一个外观相似的模型。一个有经验的使用者学会了怎样选择单元的正确类型，和在模型的不同区域中怎样来细分网格。一个经常忽略的错误根源是在一个模型中的负载和边界条件上进行了错误的假设。同时也很轻易地相信一个有限元模型的每个十进位的结果。以及忘掉了在负载和边界条件上粗糙的假设。如果有一个关于怎样建立边界条件模型的问题的话，宁可用你的模型以不同的方法去测试其灵敏度，而不是仅遵循一种方法，得出一种答案，

结构有限元分析

中国海洋大学本科生课程大纲一、课程介绍 1.课程描述：本课程是船舶与海洋工程专业重要的结构计算分析课程，通过多媒体教学和上机练习，系统学习结构有限元FEM的基本原理和方法，熟悉掌握通用的有限元软件ANSYS进行结构静力和动力分析的方法和步骤。 2.设计思路：有限元方法是一种现代设计方法，应用于结构设计中，是一种具有重要经济意义和巨大潜力的先进结构设计技术。因此选择该课程作为结构设计方面的一门必修课程，主要介绍结构有限元的基本原理和方法，还选择了通用的有限元软件ANSYS进行示例分析。包括要求掌握有限元法的基本思想和基本原理、平面刚架结构的有限元法、弹性力学平面问题以及结构动力学问题的有限元基本理论，并通过通用的有限元软件ANSYS了解解决相关问题的过程，同时掌握ANSYS进行结构静力和模态分析的基本步骤和方法，了解ANSYS进行结构瞬态动力分析的基本步骤。 3.课程与其他课程的关系先修课程：结构力学、弹性力学。本课程与结构力学和弹性力学相关，在掌握了结构分析的基本概念和方法之后才能很好地学习结构有限元分析。在后续课程中结构 - 3 -

有限元分析为学生在海洋平台设计课程设计及毕业设计中提供了结构分析的方法和软件工具。二、课程目标本课程的目标是学习掌握现代结构分析方法FEM，初步掌握通用的有限元软件ANSYS，为船舶与海洋工程结构设计、强度校核提供计算结果，为海洋结构动力响应分析提供建模基础。三、学习要求结构有限元分析是一门理论和实践性都很强的课程，在机房上课，人手一台计算机，强调实际ANSYS操作能力的培养。要达到以上学习任务，学生必须：（1）按时上课，上课认真听讲，积极参与结构分析典型案例分析。本课程将包含较多的课堂有限元ANSYS作业练习和课后结构有限元计算作业。（2）保质保量地按时完成课堂ANSYS作业练习和课后结构有限元计算作业，每位学生一个账号通过网络提交课堂ANSYS作业，只有在各项作业中认真练习才能够不断提高ANSYS的操作水平和结构分析的技能。四、教学内容 - 3 -

有限元分析

第2章软件介绍 2.1 建模软件 SolidWorks 是生信国际有限公司推出的机械设计软件, 是全参数化特征造型软SolidWorks有全面的零件实体建模功能,变量化的草图轮廓绘制,驱动参数改变特征的大小和位置,丰富的数据转换接口使SolidWorks可以将几乎所有的机械CAD 软件集成到现在的设计环境中来,SolidWorks支持的数据标准有：IGES,DXFDWG, ASCII, Parasolid等多种文件格式, 从而实现与有限元分析、流体分析软件的数据交换；Solidworks 环境下提供的文件格式EPRT 或EASM 格式可以将Solidworks文件转变为可执行文件,从而实现与外界的数据共享。该软件最大的特点是简单易学,Solidworks 软件从1997 年面世就受到广大工程技术人员的喜爱,它是参数化特征造型软件的新秀。 SolidWorks 软件包含零件建模、装配设计、工程图与钣金等模块，还与高级图像渲染软件PhotoWorks，高级有限元分析软件Simulation，机构运动学分析软件Motionworks，产品数据管理（PDM）软件SmarTeam，以及数控加工等软件无缝集成。SolidWorks还首创了自上而下的全相关设计，并凭借高效运行的装配设计使之成为实作技术。 2.2 有限元分析法有限单元法是在当今科学技术发展和工程分析中获得最广泛应用的数值计算方法，它是用计算机把复杂的零件形体自动分割成有限个形状简单的单元，然后逐个分析、计算这些小单元体的变形，并按一定的关系求得零件的总变形。由于它的通用性和有效性，理论基础牢靠，物理概念清晰，解决问题效率高，能为工程师在设计阶段掌握产品性能、优化产品的结构，缩短设计试验周期，使设计制造的产品具有较强的竞争力等优点，因而受到工程技术界的高度重视。伴随着计算机科学和技术的快速发展，现己成为计算机辅助工程和数值仿真的重要组成部分。 2.2.1 有限元原理有限元分析法的基本原理是用较简单的问题代替复杂问题后再求解。它的核