孔挤压强化残余应力场的三维有限元模拟和实验研究
有限元分析 均布荷载作用下深梁的变形和应力
有 限 元 分 析 上 级 报 告 学院： 专业： 姓名： 班级： 学号：
均布荷载作用下深梁的变形和应力 两端简支，长度l=5m，高度h=1m的深梁，在均布荷载q =5000N/m作用下发生平面弯曲（如图4.1所示）。已知弹性模量为30Gpa，泊松比为0.3，试利用平面应力单元PLANE82，确定跨中的最大挠度，和上下边缘的最大拉压应力。 4.1 均布荷载作用下深梁计算模型 1．理论解 具有两个简支支座支承的简支梁，它的变形和应力分布在理论上是没有解析表达式。 在一般的弹性力学教科书中，只有将两边支座简化为等效力的条件，即在两个支座的侧表面上作用有均匀分布的剪力情况，才可以得到理论解答。 (1) 设定应力函数。 获得这种情况下的解答的主要思路是：按照应力解法，考虑到应力分量关于该梁中心 位置（x=2.5，y=0.5）有对称和反对称关系。可以首先假定一个应力函数为： Φ = A(y - 0.5)5+ B(x - 2.5)2 (y -0.5)3 +C(y -0.5)3+ D(x- 2.5)2+ E(x -2.5)2 (y - 0.5) (4.1) 依据这个应力函数，可以获得各个应力分量，按照上表面受均布压力作用简支梁的上 下表面和左右侧表面的应力边界条件，确定出应力函数（4.1）中的各个待定系数A，B，C，D和E。 按照应力求解平面应力问题方法，应力函数应该满足双调和函数： ?2?2Φ = 0 (4.2) 将（4.1）应力函数代入上式后，得到： 24 B( y - 0.5) +120A(y - 0.5) = 0 (4.3) 即： B = -5A (4.4) (2)确定应力分量。 应力函数与应力分量之间的关系为： (3) 利用梁的上下表面边界条件确定积分常数。 上表面受均布压力作用简支梁的上表面（y=h=1m）的应力边界条件:
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有限元分析与应用》由简单到复杂，循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题，具有很高的工程应用价值。
有限元法的基本思想及计算 步骤
有限元法的基本思想及计算步骤 有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散化。这些单元仅在顶角处相互联接，称这些联接点为结点。离散化的组合体与真实弹性体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单元之间只能通过结点来传递内力。通过结点来传递的内力称为结点力，作用在结点上的荷载称为结点荷载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改进，近似解最终将收敛于精确解。 用有限元法求解问题的计算步骤比较繁多，其中最主要的计算步骤为： 1）连续体离散化。首先，应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。 2）单元分析。所谓单元分析，就是建立各个单元的结点位移和结点力之间的关系式。现以三角形单元为例说明单元分析的过程。如图1所示，三角形有三个结点i，j，m。在平面问题中每个结点有两个位移分量u，v和两个结点力分量F x，F y。三个结点共六个结点位移分量可用列
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疲劳计算，进一部验证产品的设计寿命和可靠性。
solidworks进行有限元分析的一般步骤
1.软件形式： ㈠. SolidWorks的内置形式： ◆COSMOSXpress——只有对一些具有简单载荷和支撑类型的零件的静态分析。 ㈡. SolidWorks的插件形式： ◆COSMOSWorks Designer——对零件或装配体的静态分析。 ◆COSMOSWorks Professional——对零件或装配体的静态、热传导、扭曲、频率、掉落测试、优化、疲劳分析。 ◆COSMOSWorks Advanced Professional——在COSMOSWorks Professional的所有功能上增加了非线性和高级动力学分析。 ㈢. 单独发行形式： ◆COSMOS DesignSTAR——功能与COSMOSWorks Advanced Professional相同。 2.使用FEA的一般步骤： FEA=Finite Element Analysis——是一种工程数值分析工具，但不是唯一的数值分析工具！其它的数值分析工具还有：有限差分法、边界元法、有限体积法… ①建立数学模型——有时，需要修改CAD几何模型以满足网格划分的需要， （即从CAD几何体→FEA几何体），共有下列三法： ▲特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理（注：如果选中了“使用中面的壳网格”做为“网格类型”，COSMOSWorks会自动地创建曲面几何体）。▲清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。如模型中的细长面、多重实体、移动实体及其它质量问题会造成网格划分的困难甚至无法划分网格—这时我们可以使用CAD质量检查工具（即SW菜单: Tools→Check…）来检验问题所在，另外含有非常短的边或面、小的特征也必须清除掉（小特征是指其特征尺寸相对于整个模型尺寸非常小！但如果分析的目的是找出圆角附近的应力分布，那么此时非常小的内部圆角应该被保留）。 ②建立有限元模型——即FEA的预处理部分，包括五个步骤： ▲选择网格种类及定义分析类型（共有静态、热传导、频率…等八种类别）——这时将产生一个FEA算例，左侧浏览器中之算例名称之后的括号里是配置名称； ▲添加材料属性: 材料属性通常从材料库中选择，它不并考虑缺陷和表面条件等因素，与几何模型相比，它有更多的不确定性。 ◇右键单击“实体文件夹”并选择“应用材料到所有”——所有零部件将被赋予相同的材料属性。 ◇右键单击“实体文件夹”下的某个具体零件文件夹并选择“应用材料到所有实体”——某个零件的所有实体(多实体)将被赋予指定的材料属性。 ◇右键单击“实体文件夹”下具体零件的某个“Body”并选择“应用材料到实体”——只有
有限元分析中的一些问题
有限元分析的一些基本考虑-—-—-单元形状对于计算精度的影响 笔者发现,在分析复杂问题时,我们所可能出现的错误，竟然是一些很根本的错误,这些根本错误是由于对有限元的基本理论理解不清晰而造成的。 鉴于这个原因，笔者决定对一些基本问题(例如单元形状问题,单元大小问题,应力集中问题等)展开调查,从而形成了一系列文章，本篇文章是这些系列文章中的第一篇. 本篇文章先考虑有限元分析中的第一个基本问题：单元形状问题。 我们知道,单元形状对于有限元分析的结果精度有着重要影响,而对单元形状的衡量又有着诸多指标,为便于探讨,这里首先只讨论第一个最基本的指标：长宽比(四边形单元的最长尺度与最短尺度之比),而且仅考虑平面单元的长宽比对于计算精度的影响。 为此,我们给出一个成熟的算例。该算例是一根悬臂梁,在其端面施加竖直向下的抛物线分布载荷,我们现在考察用不同尺度的单元划分该梁时,对于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的单元进行划分的话，那么我们的结果可以说是完全错误的. 下面按照上表绘制出一张图,该图从形象的角度表达了上表的含义.
ANSYS 有限元分析基本流程
第一章实体建模 第一节基本知识 建模在ANSYS系统中包括广义与狭义两层含义，广义模型包括实体模型和在载荷与边界条件下的有限元模型，狭义则仅仅指建立的实体模型与有限元模型。建模的最终目的是获得正确的有限元网格模型，保证网格具有合理的单元形状，单元大小密度分布合理，以便施加边界条件和载荷，保证变形后仍具有合理的单元形状，场量分布描述清晰等。 一、实体造型简介 1．建立实体模型的两种途径 ①利用ANSYS自带的实体建模功能创建实体建模： ②利用ANSYS与其他软件接口导入其他二维或三维软件所建立的实体模型。 2．实体建模的三种方式 (1)自底向上的实体建模 由建立最低图元对象的点到最高图元对象的体，即先定义实体各顶点的关键点，再通过关键点连成线，然后由线组合成面，最后由面组合成体。 (2)自顶向下的实体建模 直接建立最高图元对象，其对应的较低图元面、线和关键点同时被创建。 (3)混合法自底向上和自顶向下的实体建模 可根据个人习惯采用混合法建模，但应该考虑要获得什么样的有限元模型，即在网格划分时采用自由网格划分或映射网格划分。自由网格划分时，实体模型的建立比较1e单，只要所有的面或体能接合成一体就可以：映射网格划分时，平面结构一定要四边形或三边形的面相接而成。 二、ANSYS的坐标系 ANSYS为用户提供了以下几种坐标系，每种都有其特定的用途。 ①全局坐标系与局部坐标系：用于定位几何对象(如节点、关键点等)的空间位置。 ②显示坐标系：定义了列出或显示几何对象的系统。 ③节点坐标系：定义每个节点的自由度方向和节点结果数据的方向。 ④单元坐标系：确定材料特性主轴和单元结果数据的方向。 1．全局坐标系 全局坐标系和局部坐标系是用来定位几何体。在默认状态下，建模操作时使用的坐标系是全局坐标系即笛卡尔坐标系。总体坐标系是一个绝对的参考系。ANSYS提供了4种全局坐标系：笛卡尔坐标系、柱坐标系、球坐标系、Y-柱坐标系。4种全局坐标系有相同的原点，且遵循右手定则，它们的坐标系识别号分别为：0是笛卡尔坐标系(cartesian)，1是柱坐标系 (Cyliadrical)，2是球坐标系(Spherical)，5是Y-柱坐标系(Y-aylindrical)，如图2-1所示。
基于有限元ANSYS的压力容器应力分析报告
压力容器分析报告
目录 1 设计分析依据 (1) 1.1 设计参数 (1) 1.2 计算及评定条件 (1) 1.3 材料性能参数 (1) 2 结构有限元分析 (2) 2.1 理论基础 (2) 2.2 有限元模型 (2) 2.3 划分网格 (3) 2.4 边界条件 (5) 3 应力分析及评定 (5) 3.1 应力分析 (5) 3.2 应力强度校核 (6) 4 分析结论 (8) 4.1 上封头接头外侧 (9) 4.2 上封头接头内侧 (11) 4.3 上封头壁厚 (13) 4.4 筒体上 (15) 4.5 筒体左 (17) 4.6 下封头接着外侧 (19) 4.7 下封头壁厚 (21)
1 设计分析依据 （1）压力容器安全技术监察规程 （2）JB4732-1995 《钢制压力容器-分析设计标准》-2005确认版 1.1 设计参数 表1 设备基本设计参数 正常设计压力MPa 7.2 最高工作压力MPa 6.3 设计温度℃0~55 工作温度℃5~55 工作介质压缩空气46#汽轮机油 焊接系数φ 1.0 腐蚀裕度mm 2.0 容积㎡ 4.0 容积类别第二类 计算厚度mm 筒体29.36 封头29.03 1.2 计算及评定条件 （1）静强度计算条件 表2 设备载荷参数 设计载荷工况工作载荷工况 设计压力7.2MPa 工作压力6.3MPa 设计温度55℃工作温度5~55℃ 注：在计算包括二次应力强度的组合应力强度时，应选用工作载荷进行计算，本报告中分别选用设计载荷进行计算，故采用设计载荷进行强度分析结果是偏安全的。 1.3 材料性能参数 材料性能参数见表3，其中弹性模型取自JB4732-95表G-5，泊松比根据JB4732-95的公式（5-1）计算得到，设计应力强度分别根据JB4732-95的表6-2、表6-4、表6-6确定。 表3 材料性能参数性能
有限元分析案例
有限元分析案例 图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
基于有限元的钢板弹簧应力分析
基于有限元的钢板弹簧应力分析 蒋阳 西华大学交通与汽车工程学院 摘 要：本文讨论了利用ANSYS 软件对钢板弹簧进行映射网格划分，并在两簧片的接触区域生成ANSYS 软件所提供的接触单元，建立起多片钢板弹簧的有限元模型。分析了施加预负荷和工作负荷时，板簧应力值显著增长的部位，从而预测板簧产生断裂的部位，可为改进设计提供指导作用。关键词：板簧、仿真、模态 1引言 钢板弹簧具有结构简单，制造、维修方便，除了作为弹性元件外，还可兼起导向和传递侧向、纵向力和力矩的作用，其片间的接触、摩擦在弹簧振动时还将起到阻尼的作用，是重要的高负荷安全部件，目前在商用车上仍被广泛采用[1]。 传统的钢板弹簧设计方法分为：三角形板计算法，板端接触法，共同曲率法[2]。上述三种计算方法对实际工作中的钢板弹簧进行简化，并不能反应实际工作中存在的复杂的非线性状态以及接触情况。 本文利用有限元分析软件ANSYS，对十片钢板弹簧的装配过程和工作过程进行计算分析的基础上， 求得在预负荷和工作负荷作用下的应力与位移等响应情况，为实际钢板弹簧的设计中确定参数提供了依据。 2计算模型的建立 某车型的板簧总成的三维实体结构见图1。板簧建模时，考虑板簧总成对称性，同时为了方便建模，取其1/4为研究对象。通过单元solid45划分网格之后，得到12390个单元和12348个节点。 图1 板簧的三维模型图 3钢板弹簧的材料属性和网格划分 钢板弹簧的材料为60CrMnBa，弹性模量为206GPa，泊松比为0.26。 solid45单元用于构造三维实体结构。单元通过8个节点来定义,每个节点有3个沿着xyz方向平移的自由度。该六面体单元有8个节点，每个节点具有X、Y、Z 三个方向的平动自由度，可以进行塑性分析、蠕变分析、膨胀分析、应力硬化分析、大变形分析和大应变能力[1]。 用单元solid45划分网格之后，得到12390个单元和12348个节点。 4钢板弹簧的片间接触单元的建立 钢板弹簧总成的片与片之间，接触与否事先未知，而且接触后存在着滑移，所以在片与片的节点间建立接触单元，模拟片间的作用力。在ANSYS 中，接触单元是覆盖在分析模型接触面之上的一层单元，主要通过识别接触对与生成接触单元，设置单元关键字和实常数来创建。利用接触单元可以跟踪接触位置、保证接触协调性防止接触表面相互穿透；并在接触表面之间传递接触应力(正压力和摩擦力)。本文选用的接触单元分别是TARGE170 （ 三维目标单元） 与CONTA173（三维8节点面与面接触单元）[2]。本文单元的实常数采用默认值。摩擦采用库仑模型，钢板弹簧之间的摩 101
板结构有限元分析实例详解
板结构有限元分析实例详解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对话框
基于ANSYS的齿轮应力有限元分析
本科毕业设计 论文题目：基于ansys的齿轮应力有限元分析 学生姓名： 所在院系：机电学院 所学专业：机电技术教育 导师姓名： 完成时间：
摘要 本文主要分析了在ansys中齿轮参数化建模的过程。通过修改参数文件中的齿轮相关参数，利用APDL语言在ANSYＳ软件中自动建立齿轮的渐开线。再利用图形界面操作模式，通过一系列的镜像、旋转等命令，生成两个相互啮合的大小齿轮。运用有限元分析软件ANSYS对齿轮齿根应力和齿轮接触应力进行分析计算,得出两个大小齿轮的接触应力分布云图。通过与理论分析结果的比较,验证了ANSYS在齿轮计算中的有效性和准确性。 关键词 :ANSYS,APDL，有限元分析,渐开线，接触应力。
Modeling and Finite Element Analysis of Involute Spur Gear Based on ANSYS Abstract We have mainly analyzed spur gear parametrization modelling process in the ansys software. using the APDL language through revises the gear related parameter in the parameter document,we establishes gear's involute automatically in the ANSYS software.Then, using the graphical interface operator schema, through a series of orders ,mirror images, revolving and so on, we produce the big and small gear which two mesh mutually. Carring on the stress analysis of the gear by using the finite element analysis software-- ANSYS, we obtain two big and small gear's contact stress distribution cloud charts. through with the theoretical analysis result's comparison,we explain ANSYS in the gear computation validity and the accuracy. Keywords: ANSYS; APDL;finite element analysis;involute line;contact stress
基于ANSYS的有限元分析
有限元大作业 基于ansys的有限元分析 班级： 学号： 姓名： 指导老师： 完成日期：
ANSYS软件是美国ANSYS公司研制的大型通用有限元分析（FEA）软件，是世界范围内增长最快的计算机辅助工程（CAE）软件，能与多数计算机辅助设计（CAD，computer Aided design）软件接口，实现数据的共享和交换，如Creo,NASTRAN, Alogor, I－DEAS, AutoCAD 等。是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。在核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电等领域有着广泛的应用。ANSYS功能强大，操作简单方便，现在已成为国际最流行的有限元分析软件，在历年的FEA评比中都名列第一。目前，中国100多所理工院校采用ANSYS软件进行有限元分析或者作为标准教学软件。 2D Bracket 问题描述: We will model the bracket as a solid 8 node plane stress element. 1.Geometry: The thickness of the bracket is 3.125 mm 2.Material: steel with modulus of elasticity E=200 GPa. 3.Boundary conditions: The bracket is fixed at its left edge. 4.Loading: The bracket is loaded uniformly along its top surface. The load is 2625 N/m. 5.Objective: a.Plot deformed shape b.Determine the principal stress and the von Mises stress. (Use the stress plots to determine these) c.Remodel the bracket without the fillet at the corner or change the fillet radius to 0.012 and 0.006m, and see how d.principal stress and von Mises stress chang e.
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]对话框。
有限元法分析过程
有限元法分析过程 有限元法分析过程大体可分为：前处理、分析、后处理三大步骤。 对实际的连续体经过离散化后就建立了有限元分析模型，这一过程是有限元的前处理过程。在这一阶段，要构造计算对象的几何模型，要划分有限元网格，要生成有限元分析的输入数据，这一步是有限元分析的关键。 有限元分析过程主要包括：单元分析、整体分析、载荷移置、引入约束、求解约束方程等过程。这一过程是有限元分析的核心部分，有限元理论主要体现在这一过程中。 有限元法包括三类：有限元位移法、有限元力法、有限元混合法。 在有限元位移法中，选节点位移作为基本未知量； 在有限元力法中，选节点力作为未知量； 在有限元混合法中，选一部分基本未知量为节点位移，另一部分基本未知量为节点力。 有限元位移法计算过程的系统性、规律性强，特别适宜于编程求解。一般除板壳问题的有限元应用一定量的混合法外，其余全部采用有限元位移法。因此，一般不做特别声明，有限元法指的是有限元位移法。 有限元分析的后处理主要包括对计算结果的加工处理、编辑组织和图形表示三个方面。它可以把有限元分析得到的数据，进一步转换为设计人员直接需要的信息，如应力分布状态、结构变形状态等，并且绘成直观的图形，从而帮助设计人员迅速的评价和校核设计方案。 附：FELAC 2.0软件简介 FELAC 2.0采用自定义的有限元语言作为脚本代码语言，它可以使用户以一种类似于数学公式书写和推导的方式，非常自然和简单的表达待解问题的微分方程表达式和算法表达式，并由生成器解释产生完整的并行有限元计算C程序。 FELAC 2.0的目标是通过输入微分方程表达式和算法之后，就可以得到所有有限元计算的程序代码，包含串行程序和并行程序。该系统采用一种语言(有限元语言)和四种技术(对象技术、组件技术、公式库技术生成器技术)开发而成。并且基于FELAC 1.0的用户界面，新版本扩充了工作目录中右键编译功能、命令终端输入功能，并且丰
有限元分析中的一些问题
有限元分析的一些基本考虑-----单元形状对于计算精度的影响 笔者发现，在分析复杂问题时，我们所可能出现的错误，竟然是一些很根本的错误，这些根本错误是由于对有限元的基本理论理解不清晰而造成的。 鉴于这个原因，笔者决定对一些基本问题（例如单元形状问题，单元大小问题，应力集中问题等）展开调查，从而形成了一系列文章，本篇文章是这些系列文章中的第一篇。 本篇文章先考虑有限元分析中的第一个基本问题：单元形状问题。 我们知道，单元形状对于有限元分析的结果精度有着重要影响，而对单元形状的衡量又有着诸多指标，为便于探讨，这里首先只讨论第一个最基本的指标：长宽比（四边形单元的最长尺度与最短尺度之比），而且仅考虑平面单元的长宽比对于计算精度的影响。 为此，我们给出一个成熟的算例。该算例是一根悬臂梁，在其端面施加竖直向下的抛物线分布载荷，我们现在考察用不同尺度的单元划分该梁时，对于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，则误差越大。
ansys有限元分析作业经典案例
有 限 元 分 析 作 业 作业名称 输气管道有限元建模分析 姓 名 陈腾飞 学 号 3070611062 班 级 07机制（2）班 宁波理工学院
题目描述： 输气管道的有限元建模与分析 计算分析模型如图1所示 承受内压：1.0e8 Pa R1=0.3 R2=0.5 管道材料参数：弹性模量E=200Gpa；泊松比v=0.26。 图1受均匀内压的输气管道计算分析模型（截面图） 题目分析： 由于管道沿长度方向的尺寸远远大于管道的直径，在计算过程中忽略管道的断面效应，认为在其方向上无应变产生。然后根据结构的对称性，只要分析其中1/4即可。此外，需注意分析过程中的单位统一。 操作步骤 1.定义工作文件名和工作标题 1.定义工作文件名。执行Utility Menu-File→Chang Jobname-3070611062，单击OK按钮。 2.定义工作标题。执行Utility Menu-File→Change Tile-chentengfei3070611062，单击OK按钮。 3.更改目录。执行Utility Menu-File→change the working directory –D/chen 2.定义单元类型和材料属性 1.设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK
2.选择单元类型。执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →apply Add/Edit/Delete →Add →select Solid Brick 8node 185 →OK Options…→select K3: Plane strain →OK→Close如图2所示，选择OK接受单元类型并关闭对话框。 图2 3.设置材料属性。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic，在EX框中输入2e11，在PRXY框中输入0.26，如图3所示，选择OK并关闭对话框。 图3 3.创建几何模型 1. 选择ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3) →OK
matlab有限元分析实例
1.物理现象：这个对工程师来说是直观的物理现象和物理量，温 度多少度，载荷是多大等等。通常来说，用户界面中呈现的、用户对工程问题进行设置时输入的都是此类信息。 2.数学方程：将物理现象翻译成相应的数学方程，例如流体对应 的是NS方程，传热对应的是传热方程等等；大部分描述这些现象的方程在空间上都是偏微分方程，偶尔也有ODE（如粒子轨迹、化学反应等）。在这个层面，软件把物理现象“翻译” 为以解析式表示的数学模型。 3.数值模型：在定义了数学模型，并执行了网格剖分后，商业软 件会将数学模型离散化，利用有限元方法、边界元法、有限差分法、不连续伽辽金法等方法生成数值模型。软件会组装并计算方程组雅可比矩阵，并利用求解器求解方程组。这个层面的计算通常是隐藏在后台的，用户只能通过一些求解器的参数来干预求解。 有限元是一种数值求解偏微分方程的方法。 基本过程大致是设置形函数，离散，形成求解矩阵，数值解矩阵，后处理之类的。 MATLAB要把这些过程均自己实现，不过在数值求解矩阵时可以调用已有函数。可以理解为MATLAB是一个通用的计算器，当然它的功能远不止如此。
而ANSYS之类的叫做通用有限元软件，针对不同行业已经将上述过程封装，前后处理也比较漂亮，甚至不太了解有限元理论的人也能算些简单的东西，当然结果可靠性又另说了。 比较两者，ANSYS之类的用起来容易得多，但灵活性不如MATLAB。MATLAB用起来很困难，也有人做了一些模块，但大多数只能解决一些相对简单的问题。 对于大多数工程问题，以及某些领域的物理问题，一般都用通用有限元软件，这些软件还能添加一些函数块，用以解决一些需要额外设置的东西。但是对于非常特殊的问题，以及一般性方程的有限元解，那只能用MATLAB或C，Fortran之类的了。
有限元分析的一般过程
一、结构的离散化 将结构或弹性体人为地划分成由有限个单元，并通过有限个节点相互连接的离散系统。 这一步要解决以下几个方面的问题: 1、选择一个适当的参考系，既要考虑到工程设计习惯，又要照顾到建立模型的方便。 2、根据结构的特点，选择不同类型的单元。对复合结构可能同时用到多种类型的单元，此时还需要考虑不同类型单元的连接处理等问题。 3、根据计算分析的精度、周期及费用等方面的要求，合理确定单元的尺寸和阶次。 4、根据工程需要，确定分析类型和计算工况。要考虑参数区间及确定最危险工况等问题。 5、根据结构的实际支撑情况及受载状态，确定各工况的边界约束和有效计算载荷。 二、选择位移插值函数 1、位移插值函数的要求 在有限元法中通常选择多项式函数作为单元位移插值函数，并利用节点处的位移连续性条件，将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。 位移插值函数需要满足相容（协调）条件，采用多项式形式的位移插值函数，这一条件始终可以满足。 但近年来有人提出了一些新的位移插值函数，如：三角函数、样条函数及双曲函数等，此时需要检查是否满足相容条件。 2、位移插值函数的收敛性（完备性）要求： 1）位移插值函数必须包含常应变状态。 2）位移插值函数必须包含刚体位移。 3、复杂单元形函数的构造 对于高阶复杂单元，利用节点处的位移连续性条件求解形函数，实际上是不可行的。因此在实际应用中更多的情况下是利用形函数的性质来构造形函数。 形函数的性质： 1）相关节点处的值为 1，不相关节点处的值为 0。 2）形函数之和恒等于 1。 1、建立数学模型（特征消隐，理想化，清除）（（即从CAD 几何体→FEA 几何体），共 有下列三法：▲ 特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲ 理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理▲ 清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。） 2、建立有限元模型：（选择网格种类及定义分析类型；添加材料属性；施加约束；定义载 荷；网格划分） 3、求解有限元模型：再在此基础上计算应变和应力等其它物理量；在热分析中，FEA 首先 计算的是网格中每个节点的温度（标量），再在此基础上计算温度梯度和热流等其它物理量. 一般如果模型可划分网格，那么它就可以求解，但如果没有定义材料或载荷，则求解会终止。 4、结果分析：材料线性假设、小变形假设、静态载荷假设等等。
Ansys有限元分析实例
课程论文 (2015-2016学年第一学期) 有限元理论在软件中的应用与刚度矩阵的求解 学生：张贺
有限元分析案例：打点喷枪模组（用于手机平板电脑等电子元件粘接），该产品主要是使用压缩空气推动模组内的顶针作高频上下往复运动，从而将高粘度的胶水从喷嘴中打出（喷嘴尺寸￠0.007”）。顶针是这个产品中的核心零件，设计使用材料是：AISI 4140 最高工作频率是160HZ（一个周期中3ms开3ms关），压缩空气压力3-8bar, 直接作用在顶针活塞面上，用Ansys仿真模拟分析零件的强度是否符合要求。 1. 零件外形设计图： 2. 简化模型特征后在Ansys14.0 中完成有限元几何模型创建： 3. 选择有限元实体单元并设定，单元类型是SOILD185，由于几何建模时使用的长度单位是mm, Ansys采用单位是长度：mm 压强：Mpa 密度：Ton/M3。根据题目中的材料特性设置该计算模型使用的材料属性：杨氏模量 2.1E5；泊松比：0.29； 4. 几何模型进行切割分成可以进行六面体网格划分的规则几何形状后对各个实体进行六面体网格划分，网格结果：
5. 依据使用工况条件要求对有限元单元元素施加约束和作用载荷： 说明：约束在顶针底端球面位移全约束； 分别模拟当滑块顶断面分别以8Bar，5Bar，4Bar和3Bar时分析顶针的内应力分布，根据计算结果确定该产品允许最大工作压力范围。 6. 分析结果及讨论： 当压缩空气压力是8Bar时：
当压缩空气压力是5Bar时： 当压缩空气压力是4Bar时：
结论： 通过比较在不同压力载荷下最大内应力的变化发现，顶针工作在8Bar时最大应力达到250Mpa，考虑到零件是在160HZ高频率在做往返运动，疲劳寿命要求50百万次以上，因此采用允许其最大工作压力在5Mpa，此时内应力为156Mpa，按线性累积损伤理论[3 ]进行疲劳寿命L-N疲劳计算，进一部验证产品的设计寿命和可靠性。