ANSYS飞机仿真方案

ANSYS航空发动机仿真方案

第三章航空发动机仿真方案 1.航空发动机行业概况

1903年12月17日，美国莱特兄弟实现了人类历史上首次有动力、载人、持续、稳定和可操作的重于空气飞行器的飞行。这使得几千年来由少数人从事的飞行探索事业在后来的百年中发展成为对世界政治、经济、军事和技术，甚至人们的生活方式都有重要影响的航空业。

因此，从狭义上说，航空发动机是航空器飞行的动力，是航空器的“心脏”；从广义上说，航空发动机是航空事业发展的推动力，对航空器的性能和研制有着决定性的影响。

1.1.世界航空发动机发展历史

航空发动机百余年的发展历史大致可分为两个时期：第一个时期从莱特兄弟的首次飞行开始到二次世界大战结束为止。这个时期内，活塞式发动机统治了40年左右的时间。第二个时期是从二次世界大战结束至今，燃气涡轮发动机取代了活塞式发动机，居航空动力的主导地位，开创了喷气时代。在燃气涡轮发动机的60多年发展历程中，大致经历了四次更新换代：

第一代是单转子亚音速喷气发动机。这一代发动机大多数在20世纪30～40年代研制，40年代末50年代初投入使用。压气机采用离心式和轴流式两种，总增压比在5左右，单管燃烧室，单级涡轮；推重比3左右。有代表性的机种有：美国的J47(TG-190)、前苏联的VK-1和法国的阿塔（Atar）发动机。

第二代是超声速涡喷发动机。这些发动机在第一代发动机的基础上有了许多创新，大都在50年代研制。主要技术特点是：双转子、进口导流叶片可调、超声速压气机、高温涡轮、推重比达到5左右。用这一代发动机装配的飞机都是超声速战斗机。代表机种有：美国的J79和前苏联的R11-300R。

第三代是超声速涡扇发动机。这一代发动机的研制始于60年代，主要技术特点是：涡扇发动机、核心机技术，2D设计、环形燃烧室、气冷涡轮、结构完整性设计、新材料、推重比8。这一代发动机的成长得益于全世界各种大型试验设备的建设、计算技术和制造技术的发展。用这一代发动机装配的飞机都是高性能超声速战斗机。代表机种有：美国的F404和F100、前苏联的AL31F和RD33、英国的RB199和法国的M88-2。

第四代是先进技术涡扇发动机。这一代从80年代中期开始发展，目前仍处于研制阶段。主要技术特点是：结构简单，抗撞击能力强，具有良好的耐久性可维护性；增加了不加力条件下的持续超声速巡航能力、采用2D喷管的有限矢量推力能力和隐身能力。第四代发动机的推重比为9～10。代表机种有：美国的

F119、前苏联的AL-41F和英国的EJ200。

航空燃气涡轮发动机在60多年的发展历史中经历了众多技术进步，如表

3-1-1所示：

表3-1-1航空燃气涡轮发动机的技术进步

目前，美、英第一、二代均已退役，第三代是现役主力机种。由于其性能先进，且还在不断改进改型，服役期比第一、二代长很多，估计将使用到2010年左右，俄、法、日、印、韩等国第二、第三代并存，以第三代为主。

第四代战斗机是美、苏冷战对抗时期开始研制的，原计划90年代中期装备部队。自苏联和华约解体后，是否还需要继续发展，在美国和西欧开展了一场大辩论。许多国会议员提出，将F15、F16经现代化改装后，就可以达到应付未来“地区冲突”的要求。在此影响下，德国曾一度退出欧洲战斗机EF2000发展计划。但辩论的结果认为：F15、F16经改装后，不能跨越“代”的鸿沟。为了满足“全球到达，全球力量”的战略目标，发展第四代战斗机是必须的。这场辩论使第四代战斗机的装备时间推迟了十年左右。第四代战斗机具有隐身、过失速机动、不加力超声速巡航、短距起降、超视距多目标攻击和装备更先进的航空电子与武器系统等许多特点，较之第三代具有全面优势。据报道，F22与F15相比，每飞行小时的维修工时降低约70%，其综合作战效能提高近10倍。但由于F22太贵，难以大量装备部队，又决定同时研制装单台F119发动机（推力增大型）的JSF联合战斗机，已有近10个国家参与该项计划。该机三种型别（常规型、短距起飞/垂直降落型和垂直起降型）的意向订货量高达3300架左右，总价达1800亿美元的订单已被马丁公司所得。美国将以联合战斗机JSF和F22这两种第四代战斗机轻重搭配跨入21世纪，从2003年起陆续装备部队，全面取代现役的第三代战斗机，成为美国、部分西方国家、甚至我国部分周边国家和地区 21世纪上半叶的主战机种。

西方各国对航空动力技术的预先研究一向给予极大重视，开展了一系列大型研究计划。如美国军方早从50年代中期就开始实施航空推进技术探索发展计划；70年代初至80年代又相继实施了先进战术战斗机发动机计划（ATFE）、先进涡轮发动机燃气发生器计划（ATEGG）和飞机推进分系统综合计划；70年代末以来，美国政府（由NASA主持）也先后实施了发动机部件改进计划、高效节能发动机计划（E3）、先进螺旋桨计划和发动机热端部件技术计划（HOST）。正是这些研究计划，为各种先进军、民用发动机提供了坚实的技术基础，才使得美国达到了当今世界的领先水平，使军、民用发动机跨上了一个又一个技术新台阶。航空动力行业已成为世界各航空强国的军事工业和国民经济的重要支柱产业之一。

美国在研制第四代F119发动机的同时，从1988年起的15年内又投入50

亿美元巨资，由军方与政府联合主持实施“综合高性能发动机技术计划”（即IHPTET计划）；英国则着手进行先进军用核心发动机第II阶段计划（ACME-II）。其共同目标是利用计算流体力学（CFD）、结构力学、燃烧、传热、新材料、新工艺、电子调节和计算机仿真等方面的最新成就，使推进系统的能力在现有基础上翻一番，预计2020年后有可能研制出第五代推重比为15～20的发动机。这意味着他们用15年左右的时间，在推重比、耗油率、成本等方面取得的技术进步，相当于过去30~40年的成就，充分表明世界航空发动机技术呈现加速发展态势。

1.2.我国航空发动机发展历史

我国航空动力行业经过几十年的建设，从无到有，由小到大，在维护修理、测绘仿制、批量生产、改进改型、民机开发等方面均取得很大成绩。新机研制和预先研究也取得可喜进展，共生产了近60000万台各型发动机，为建立一支强大的人民空军、海军和陆军航空兵，保障国家安全和促进国民经济建设做出了重大贡献。

但是，航空动力长期在测绘仿制中徘徊，走了不少弯路，与加速发展的世界先进水平相比，我国航空发动机大大落后了。我国现役军机大量使用的发动机，大都是仿制前苏联的产品或者是其改进改型，推重比为5～5.5左右，仅相当于国外早已淘汰的第二代水平；民机动力方面，目前生产的仍是仿制前苏联50～60年代的中小型涡桨发动机，干线民航机大涵道比涡扇发动机基本上还是空白；在航改燃机方面，80年代在成熟的航空发动机基础上改型了6种型号，生产了百余台，但仅占我国近800万千瓦燃机总装机容量的5%左右；从80年代初开始的高性能发动机关键技术预研，虽取得一定进展，但距工程应用还有一段较长的路要走；80年代中期开始研制的某发动机，虽已实现首飞，但距设计定型尚需时日。综合评估我国航空动力的总体技术水平，较国外相差一代半，落后约25～30年，而且这种差距还有进一步拉大的危险。我国航空发动机落后，已成为严重制约整个航空工业快速发展的“瓶颈”，这是不争的事实。造成这种局面的原因很多，“冰冻三尺，非一日之寒”，客观上航空发动机技术十分复杂，研制难度很大，花钱多，周期长，我国工业和技术基础相对薄弱，预研和型号经费投入不足；主观上对航空发动机研制的复杂性和规律性认识不足，对预研工作重视不够，技术储备少；摊子大，战线长，力量分散，型号品种多。在国际形势出现变化和动力技术快速发展的新形势下，却未能及时进行必要的结构调整和技术改造；引进仿制机种过多，总以为测仿来得快，40多年中先后测仿了来自国外的十多种型号，大都是当时该国的二、三流产品；尤为失策的是对引进技术未能很好地组织消化、吸收和创新，且引进往往挤掉或消弱了国内的新机研制；长期过份强调型号牵引，有了型号才有钱，对预先研究、打基础重视不够，预研投资强度太低，结果使基础薄弱，技术储备不足，许多关键技术没有提前突破，型号研制成了“无米之炊”，久攻不下，为了满足飞机对发动机的需求而不得不再来一轮新的引进、仿制，从而陷入“恶性循环”；在管理上，领导机构重叠，管理分散，政出多门，意见不一，缺乏稳定、权威的中长期发展规划，缺乏科学的决策程序，型号和大型预研项目的“上马”和“下马”，往往因机构调整和人事变动而出现大的变化；没有协调处理好发动机与飞机（平台）、型号与预研、全新研制与改进改型、自行研制与国外引进、工业部门和使用单位之间的关系；在工业基础方面没有安排好提前突破先进的关键材料和制造工艺。

管理上滞后和指导思想上的不适应，对发动机的发展有时甚至比技术上落后的影响更大。

事实上，早在上世纪60年代，周恩来总理就指出：我们的飞机得了“心脏病”。40多年过去了，这种状况并没有实质性的好转，国内民航干线客机动力已全部被国外所占领，而新研军机因没有自行研制的先进动力可用，不得不买装国外发动机。我军装备因动力而受制于人，形势十分严峻。根治飞机“心脏病”已刻不容缓，这是摆在我们面前严肃而又紧迫的任务。图3-1-1表明，从军用航空发动机最重要的性能指标“推重比”来看，我国航空发动机与国外发达国家对比呈现出剪刀差的落后趋势，我国航空发动机仿制生产起于50年代初，起步并不算晚，而70年代以后差距愈拉愈大，与世界航空发动机技术的加速发展态势形成巨大的反差！

值得特别指出的是，我国周边国家日本早就与美国合作生产第三代F110发动机，并参与世界一流水平的大型民用涡扇发动机的国际合作研制，目前又正在与美、英合作研制飞行速度5倍声速的HYPR-9组合循环发动机，力图在高超声速推进技术领域抢占领先地位。印度自行研制的GTX-35VS双转子涡扇发动机，推重比7.5左右，预计2002年将装在他们自行设计的LCA轻型战斗机上首飞（后改用美国的F404-F2发动机）。他们计划在未来几年内对其现有的 780架各型作战飞机进行大规模更新换代，使现役第三代飞机由目前的230架增加至560

多架，还准备引进50架装有推力矢量喷管的俄制SU-30MKT 战斗机，并计划在15年内生产150～200架。台湾地区虽未独自研制航空发动机，但在美国人帮助下，合资研制了FTE1042涡扇发动机，并成功地应用于“经国号”轻型战斗机（装备130多架）。此外继购买150架F16A/B之后，又引进60架幻影2000-5，明显地提高了装备水平，增强了空军实力。可以说，我国空军（包括海军航空兵和陆军航空兵）的装备，除了数量上仍占一定优势之外，其技术水平已经或即将被这些周边国家和地区赶上和超过。这种严峻的形势，对我国国家安全已构成严重威胁。对此我们应有高度警惕和足够的认识。

1.3.新世纪我国航空动力具有良好的发展机遇

虽然我国还没有一台自行研制的发动机投入使用，但经过半个世纪的努力，我们进行了多个型号研制，开展了推重比8和推重比10发动机这两项大型预研工作，建成了包括大型高空台在内的基本配套的试验设施，形成了较强的制造加工能力。特别是改革开放20多年来的快速发展和对外合作，为21世纪航空动力的快速发展打下了较坚实的技术和人才基础。

过去航空产品的设计主要是依赖于各种试验，使得航空产品尤其是航空发动机的研制周期长、耗资多、风险高。近年来，由于信息技术特别是计算机辅助工程仿真（CAE）和计算流体力学（CFD）技术的发展以及大量试验数据的积累，20世纪90年代以来，在西方航空发达国家引发了一场设计技术的“革命”，初步实现了从“传统设计”向依靠计算机数学模型优化计算和虚拟现实仿真“预测设计”的转变，从而大大减少了试验工作量，提高了设计的成功率，既节约了经费，又缩短了研制周期。使发动机的研制周期从过去的10～15年缩短到6～8

年甚至4～5年，试验机也从过去的40～50台减少到10台左右。

我们若能牢牢抓住这个机遇，在我国预研和型号已取得成果的基础上，再加大力度引进、消化一些国外的先进软件，建立一套中国自己的航空动力设计体系和数值仿真系统，提高设计“起跑线”，发挥“后发优势”，不再重复西方几十年所走过的老路。这是我国航空动力走出困境，以较少的投入、较快的速度缩短与国际先进水平差距的有效途径，也是实现跨越式发展的难得机遇。

2.航空发动机研制中的典型CAE问题

如前所述，计算机辅助工程分析技术已经在西方航空发达国家引发了一场设计技术的“革命”，在航空发动机领域更是如此。其最直接的效益就是：大大减少试验工作量、提高了设计成功率、节约研制经费、缩短了研制周期。

现代航空发动机设计对CAE分析的要求呈现三大趋势：一是对分析的精度要求越来越高、二是需要分析的对象越来越庞大和复杂（系统级分析）、三是对多物理场耦合分析尤其是流固耦合分析的需求越来越多。下面简单罗列一些航空发动机研制过程中的典型CAE问题以供参考，在下面几个章节中，我们会针对这些问题予以更详细的说明并提出ANSYS的解决方案。

与结构力学分析相关的典型问题有:

压气机和涡轮的叶片、轮盘等重要构件在考虑温度、气动力、高速旋转离心力等载荷作用下的强度、振动以及高低周疲劳计算；

高低压转子系统以及与支撑结构相耦合的转子系统在考虑温度、气动力、高速旋转离心力和陀螺力矩等载荷作用下的动力特性（转子动力学）和不平衡响应计算；

齿轮、轴承、传动轴、机匣、管路系统以及其它辅机系统在相应复杂载荷和约束条件下的结构强度、振动、疲劳计算；

叶片鸟撞、包容性、抗异物损伤等高度非线性结构动力计算；

……

与计算流体力学分析相关的典型问题有

单/多级风扇、单/多级压气机、单/多级涡轮、进气道、喷管等零部件气动计算；

燃油喷射和雾化、燃烧室/加力燃烧室燃烧分析；

涡轮盘、涡轮叶片、燃烧室、喷管冷却分析；

飞机/发动机匹配气动计算，以及多相流模拟发动机地面吞水试验；……

与计算电磁学分析相关的典型问题有

多/全电发动机耐高温高性能磁悬浮轴承、启动发电机等关键零部件电磁特性计算；

电气化传动附件电磁特性及EMC/EMI计算；……

需要高性能计算的多物理场耦合计算问题示例

风扇和压气机叶片的气动-结构耦合分析；

涡轮叶片、涡轮盘、燃烧室以及喷管的气动-热-结构耦合分析；

磁悬浮轴承、启动发电机电磁-热-结构耦合分析；……

3.航空发动机结构力学计算需求及ANSYS实现

计算结构力学分析是航空发动机领域最重要、也是最早并且获得了最广泛应用的CAE方向。航空燃气涡轮发动机对性能、安全性、可靠性、耐久性、适用性的要求非常高，其重点就是体现在对结构强度、振动和寿命的设计要求上。计算结构力学分析的所有应用方向在航空燃气涡轮发动机的设计中都有很强的需求，本小节中我们将对这样一些主要分析领域的特点和ANSYS解决方案予以简要说明。

3.1结构静力分析

结构静力分析是计算结构力学分析领域中的最基础应用方向，在航空燃气涡轮发动机中，几乎所有零部件都会涉及重要性、程度和精度要求不等的静力计算问题，其主要目的是获得结构件的静强度和静刚度，并为后续的其它分析要求服务。

3.1.1.线性静力分析

线性静力分析通常适用于支承及传力框架、机匣、附件设备等对保形及刚度要求较高的零部件，有时也用于对叶片、轮盘等重要构建的前期分析。比如，在发动机总体结构中，支承刚度往往是影响整机振动的关键因素，必须对支承刚度进行分析。发动机常用的弹性支承元件包括弹性环、鼠笼和拉杆等。通过有限元静态分析，求得在一定载荷作用下结构的变形，进一步得到结构的刚度。

对于CAE程序而言，线性静力分析在技术层面上早已不是难题，可能存在的问题通常仅限于对大型问题的方便建模以及有效求解（例如，复杂叶片结构的全自动六面体协调单元网格划分、发动机整机结构的刚度和模态分析等）。在对大型问题的方便建模方面，ANSYS提供诸多全球独一无二的领先技术来完美解决，比如AWE环境下与CAD的双向参数传递技术、复杂装配体自动识别技术、

AI*Environment的基于拓扑映射的网格雕塑技术、复杂结构全自动六面体网格划分（ANSYS高级网格模块），等等；在求解方面，ANSYS也提供一系列独具特色的高效求解器，比如拥有众多专利技术的迭代求解器（JCG、ICCG、PCG等），以及大量的网络分布式并行求解器（DDS、DPCG、DJCG等）。

3.1.2.弹塑性/蠕变等材料非线性静力分析

航空燃气涡轮发动机CAE分析中非常普遍地会遇到以弹塑性/蠕变为代表的材料非线性问题。从理论上说，任何结构分析都应该是非线性的，线性分析只是对其在低强度、小变形状态下的假设而已。发动机中的高速旋转部件，如叶片、盘、鼓（以弹塑性分析为主），以及高温部件，如涡轮、燃烧室、喷管等（以蠕变以及弹塑性与蠕变组合作用分析为主）均存在不可忽视的材料非线性行为，在准确的强度/寿命校核中必然要考虑到各种复杂的材料非线性因素。

材料非线性是一种非常复杂的物理现象，在工程设计和分析中，要获得满意的仿真结果，必须很好地理解和准确地描述材料行为特性，但要做到这一点是比较困难的，Lemaitre和Chaboche曾经用如下很形象的例子描述了材料特性的复

杂性：

“在室温下的一块钢材，可以认为是：

常规结构分析中的线弹性材料；

振动阻尼问题中的粘弹性材料；

在非常大的载荷作用下，是理想塑性材料；

要精确计算永久变形时，它是具有硬化特性的弹塑性材料；

对于应力松弛问题，又有弹性粘塑性特性；

当计算成型极限时，它是延展性破坏；

当计算寿命极限时，它是疲劳破坏；……”

因此，在什么情况下用什么材料模型，与分析者所关心的现象相关。ANSYS 通用结构力学分析所提供的非线性材料本构非常复杂和宽泛，计有：20种率无关弹塑性模型，125种组合蠕变模型，7种粘塑性模型，以及粘弹性、混凝土、D-P、垫片、铸铁、形状记忆合金等等，可最大限度地满足航空发动机分析中的需求。同时，ANSYS的材料模型具有很多独特特性，比如：

多线性随动硬化模型允许直接输入实验数据；

Chachoche模型可与数个非线性随动硬化模型组合以模拟材料的复杂循环

特性（如棘齿效应、循环稳定、循环硬化和软化等）；

Hill屈服准则及其与各种强化模式的组合可模拟各种各向异性材料的非线性特性（如定向结晶和单晶涡轮叶片、复合材料风扇叶片和机匣等）；

ANSYS提供线性和非线性回归算法，用户可直接输入材料的试验数据即可获得其材料参数值，该功能适用于所有超弹模型、粘弹模型和蠕变模型。

在最常用的弹塑性材料方面，ANSYS所提供的丰富模型可以综合考虑单调硬化、包辛格效应、循环硬化或软化、棘齿效应以及记忆效应（形状记忆合金）等等，已经是非常的完善了，除此之外，ANSYS还提供了诸多在航空发动机设计中会遇到的较特殊的非线性材料模型，包括：

铸铁

铸铁由于其石墨颗粒的作用使其呈现与其它金属材料显著不同的塑性特性：在拉和压两个方向上有不同的屈服强度、塑性流动特性和硬化特性；

粘弹性

ANSYS中的粘弹性模型是Maxwell模型的通用积分形式，其松弛函数由Prony 级数表示。该模型功能全面，Maxwell、Kevin和标准线性实体都是其特殊形式，全面支持亚粘弹性和大应变超粘弹性。大应变超粘弹性基于Simo建议的列式，粘弹性行为的定义分为超弹性和松弛两个部分，所有的 ANSYS超弹性材料模型都可采用粘弹性选项（PRONY）。

粘塑性和蠕变

ANSYS程序提供几个选项用于模拟材料的率相关行为，包括蠕变，蠕变选项又含一系列适用于常规蠕变分析的蠕变法则。率相关选项是一种“过应力”模式，主要用于冲击载荷问题。程序还提供Anand模型，它最初是用于模拟高温金属的成型过程，如辊压、深拉等，该模型采用一个被称为“变形阻力”的内部变量来考虑对材料非弹性流动的各向同性阻力作用，因而也能模拟材料的硬化和软化行为，现在已广泛应用于其它领域，比如分析电子封装的焊接连接等。

超弹性

超弹性材料的应用十分广泛，比如用于密封的橡胶环等。橡胶等超弹性材料的非线性是很严重的，体现在：

非常大的应变水平（可达百分之几百！）；

材料的应力和应变呈高度非线性关系；

材料近似或完全不可压；

有很强的温度相关性；

通常并不单独存在，而是与金属等其它显著不一样的材料之间有很大的相互作用。

ANSYS超弹性材料模型有很多种类，用户可根据实际材料的实验特性等来选择合适的模型。对于超弹应用而言，ANSYS程序本身从求解器、单元技术以及解算策略等方面都进行了完善的设计，具有很好的效率和效果。

垫片材料

垫片通常在装配零件之间起密封作用，是由诸如铁、铜、橡胶、复合材料等多种材料做成的一个很薄的构件。垫片的主要变形通常限制在其法线方向上，其膜（平面内）和横向剪切等方向对刚度的贡献非常小，可以忽略不计。垫片材料通常都是受压，具有很强的非线性特性，当压力撤消时，其卸载行为非常复杂。 ANSYS垫片材料允许以表格的形式直接输入试验测量的加载过程和多个卸

载过程的“压力－闭合量（垫片上下表面之间的相对位移）”曲线数据，如果无卸载数据，则材料的卸载曲线与加载曲线一致。对于复杂的垫片连接分析，ANSYS 还提供了很多高级选项，比如：允许有初始间隙、拉伸应力限值、稳定刚度等等。同时，ANSYS还专门开发了一系列“界面”单元来专门用于垫片分析，包括4节点四边形、6节点四边形、8节点六面体、16节点六面体、12节点棱形体等适合于各种二维和三维应用的单元。

ANSYS的其它非线性材料种类就不在此一一介绍了，有兴趣者可从相关资料中去了解。最后值得一提的是，ANSYS不但支持数量庞大的材料本构模型，而且提供了一个具有良好逻辑关系的树状菜单结构来让用户很方便、很准确地定义材料参数（尤其是一些复杂的组合材料本构）。这些材料本构的特性参数都可以是与温度相关，直接在软件界面上定义。

3.1.3.接触/摩擦等状态非线性静力分析

以接触为代表的状态非线性分析问题在航空燃气涡轮发动机设计中非常普遍地存在：叶片榫头与轮盘榫槽之间的接触、带冠叶片的冠间接触、传动齿轮之间的啮合接触、轴承内的接触、各种各样的装配接触等等。

模拟两个实体之间相互作用（通常还要考虑摩擦、热、电或其它形式的能量交换）的能力对于一个分析软件来说是非常关键的，事实上，接触分析能力已经成为判断非线性分析程序优劣的一个最重要标准。同时，接触分析的性能和鲁棒性、接触对的定义及其属性管理的方便性、有效的纠错工具等几个方面是考察非线性软件接触分析能力的三个指标。

早在ANSYS2.0版本，ANSYS即具有了接触分析能力，数十年来一直根据高级分析的需要而不断地采用最新的接触算法和技术，表3-3-2简要概述了ANSYS 的接触分析特性：

ANSYS最早开发的接触单元是Contac12和Contac52，分别用于模拟二维和三维状态下的点对点接触，单元的基础算法是罚函数法和弹性库仑摩擦模型，是

所有单元中最简单的。

后来开发的Contac48和Contac49是可用于通用接触问题的点对面接触单元，采用带拉格郎日增量的罚函数法以提高兼容性，其主要特点有：允许大滑动、有摩擦或无摩擦、支持热接触等等。虽然这些单元可以模拟高度非线性接触问题（比如金属成型、辊压等），但由于点对点和点对面接触单元需要人为调整罚刚度，比较难以使用，而其最大的缺点是不能很好地与高阶实体单元配合以模拟复杂曲面接触。

事实上，在现在的ANSYS程序中，12/52、26、以及48/49等接触单元的主要作用是提供与老版本的兼容性，在新版本中都推荐采用融合了最新技术和ANSYS 独特技术的新一代接触单元，下面对这些新一代单元的主要特点做简要说明。

面对面接触

接触分析功能在1997年的ANSYS5.4版本就开始了一个质的飞跃，最初，一系列面对面接触单元（169到174）提供了刚对柔面面接触能力，虽然也是基于增广拉格朗日罚函数法，但与以前有很大的不同。在缺省情况下，程序将罚刚度与多种参数相关联（包括相邻单元的尺度、附着实体的材料特性等），无需用户提供罚刚度的绝对值，但可通过一个无量纲缩放系数来修正缺省值，这对于以弯曲变形为主的接触分析是相当必要的。目前的算法已改进为可根据其附着实体单

元的应力来自动修正接触单元罚刚度。

169-174接触单元的计算结果具有很好的可视性和可用性，其输出的结果数据是基于应力的形式而非力的形式。即使针对大型问题，这些接触单元的数值算法也非常有效，无需使用诸如自适应下降等试探性的手段即可很好地收敛。

171-174接触单元在ANSYS5.5版本就扩展为支持柔对柔面面接触的通用接触单元，这些单元采用的独特的基于高斯点的接触算法具有非常强的优势：可完美地与诸如20节点六面体单元、10节点四面体单元、8节点面单元等高阶单元结合使用，这是其它接触算法无法做到的，这也是为什么很多有限元分析程序在作接触分析时都要求用户最好采用低阶单元，而这在模拟复杂结构时要想有效地获得高精度是非常困难的。当然，ANSYS的这些单元类型也提供基于节点的接触算法作为补充，以用于一些特殊情况。

采用高阶接触单元的优势是：高阶单元以二次函数的方式模拟曲面，而低阶单元只能以一个个小平面来近似模拟曲面。这一点在航空燃气涡轮发动机设计分析中是非常重要的：由于对计算精度的高要求，尤其是针对叶片、轮盘等部件，通常必须采用高阶单元进行强度分析，ANSYS唯一提供了支持这种高阶单元接触的精确算法。

基于高斯点的接触算法避免了“锐角相交时接触方向的不明确性”，这样，目标表面可以采用较粗的网格划分。同时，这种算法也避免了节点接触算法中通常会遇到的“滑过边界”问题。ANSYS的这些接触单元还提供大量的、在某些情况下非常有用的特殊接触分析能力，包括单边和单面接触、绑定接触、无分离接触、粗糙滑动接触等。同时，ANSYS的接触分析功能还扩展到模拟壳体与壳体、壳体与实体等的固连装配，即多点约束（MPC）。

ANSYS程序采用一个功能非常强大的“接触管理器”定义、修改、管理、以及后处理观察复杂构件间的复杂接触关系，其“接触向导”功能帮助用户快速地定义各接触面及其特性参数，并有很多实用工具实现诸如观察初始接触状态、接触方向、以及根据需要修改接触特性（比如反转接触法向等）。同样，ANSYS 后处理程序中也提供了很多手段来更好地观察和理解接触分析的结果。

技术的发展是无止境的，最新版本的ANSYS在新一代ANSYSWorkBenchEnvironment（AWE）的环境下，实现了自动接触定义功能。对于象发动机整体分析这样的由数百个零部件构成的复杂装配体，如果用手工的方式来一一定义接触，将是一项非常烦琐的工作，而在ANSYS中，所有的事情都由程序自动完成了。

对于耦合场分析而言，这些接触单元所具有分析功能是前所未有的：可以模拟电－磁－热－结构耦合接触问题（如点焊过程的模拟）、热接触的热传导系数可以是接触压力或温度的函数、可以定义接触面闭合时的热传导特性和接触面分开时的热对流和辐射特性、可以进行摩擦生热分析以及控制摩擦生热量在摩擦体之间的分配、与CFD功能配合模拟装配缝隙在流体压力作用下的涨开和流体泄漏……

拉格朗日乘子法

接触分析是一个很宽泛的概念，目前还没有哪一种算法能有效地涵盖所有的应用需求。前面提到的增广拉格朗日法虽然适用性很强，且能解决大量的复杂接触问题，但该算法会产生一定量的“接触穿透”（接触穿透通常都是一个非常小的数值，可忽略不计，算法本身也提供很多选项来改善该数值）。对于需要获得完美接触相容性的问题（即要求穿透量为零），最新版本的ANSYS新一代接触单元提供纯拉格朗日乘子法以增强对穿透量的约束（该算法完全避免输入罚刚度），且避免了常规拉格朗日乘子法的一些诸如收敛困难、求解量增加等缺点。这一点在航空发动机设计分析中也是非常重要的，其高精度计算要求需要获得完美的接触相容性，以此才能准确地判断诸多的接触破坏问题，比如榫头接触面的应力剥落现象等。

3.1.

4.大变形几何非线性静力分析

航空燃气涡轮发动机的叶片变形量直接关系到间隙设计的好坏，所以必须对叶片的变形进行分析，尤其是风扇叶片。现代大功率涡扇发动机由于气动性能的要求常常设计成大展弦比或掠型，叶片在离心力和气动载荷下，其变形量远远超过线性小变形有限元分析的适用范围，必须进行非线性的大变形分析。

大变形分析其实仅是几何非线性的一种，ANSYS支持所有的几何非线性分析，除了大变形以外，大转动、大位移、大应变、应力强化、旋转软化、压力载荷强化等所有的几何非线性问题在ANSYS中均有专门手段予以完善解决。这些几何非线性类型在发动机设计中都是非常有用的，比如，旋转软化这种在其它分析程序中难以见到的分析功能非常适合于叶片动力设计分析，它完整地考虑了动刚度的对叶片动态特性的影响，类似的还有应力强化等，后续我们还会详细介绍它。

3.2.结构动力分析

航空燃气涡轮发动机是高速运转的复杂机构，对动力学的分析需求历来是重中之重。而其中叶片振动分析、转子动力学分析等又是最重点关注的问题。对于叶片振动来说，据统计资料说明，叶片的损坏事故（裂纹、断裂等）绝大部分是由于振动引起的；通过试验也证明，如果叶片振动幅度降低10％，其寿命将增加10倍！由此可见振动分析、优化结构设计以降低振动的重要性。叶片的振动通常有四种类型：气流尾迹引起的强迫振动、颤振、旋转失速、以及随机振动，其中前两种振动的危害性最大。

对于实际叶片，重要的是自振频率、振型和振动应力三个因素。叶片的自振频率范围很广，从数十赫兹（如风扇叶片）到数万赫兹，但其中以低频振动最为危险。一般情况下，频率越高，振幅就越小，危险性也就越小。在叶片振动分

析中，重点是进行准确的模态分析，并以此为基础进行谐响应分析、随机振动分析、颤振分析等，以全面而准确地了解叶片的振动特性（如振动应力等）。下面就对这些分析内容作简要描述。

3.2.1.模态分析

航空燃气涡轮发动机中的模态分析通常有这样一些类型：叶片在静止和旋转预应力状态下的固有模态计算、叶片/轮盘组合体在静止和旋转预应力条件下的循环对称模态计算、整个转子系统在静止和旋转预应力状态下的固有模态计算、整个发动机系统在静止和旋转预应力状态下的固有模态计算、等等。模态分析是其它几乎所有动力分析类型的基础，其分析精度直接关系到各种后续动力分析的效果好坏，因此，需要特别高的技术要求，主要体现在以下几点：要能考虑预应力影响（应力强化）：在旋转预应力状态下，叶片振动模态和静止状态有本质差别，频率和振型都会发生变化，通常同一阶次模态的频率都会有较大幅度的提高。ANSYS通过两步走的方式可直接考虑预应力影响，即先做一个在转速和气动力作用下的静力分析，然后再在模态分析中打开预应力开关即可；

要能考虑旋转软化的影响（动刚度）：对于象风扇这样的大展弦比叶片，大变形引起的动刚度会使叶片变得更“柔软”，降低其模态特征频率。ANSYS模态分析中的“SpinSoft”选项正是用来考虑此项影响的；

要能考虑各种阻尼的影响（复模态）：任何真实结构在动力状态下都会受到阻尼的影响，其特征值体现为复数形式（虚部为特征频率，实部表征其稳定性（振幅的指数增加或降低））。ANSYS在模态分析中可以考虑多种阻尼：α阻尼（总体质量阵乘子）、β阻尼（总体刚度阵乘子）、材料相关阻尼、单元阻尼等等，在后续的响应分析中，还可以考虑频率相关阻尼（恒定阻尼比、模态阻尼比）、频率无关阻尼等等，以完整地记及各种阻尼对结构动力特性及响应的影响。值得指出的是，在航空发动机转子动力分析中，油模（薄层液体）阻尼对动力特性的影响非常显著，ANSYS提供了专门的手段来自动计算此阻尼并用于随后的动力分析中。我们在下面的“转子动力学计算”一小节中会详细谈到此点；

要能方便地进行循环对称模态分析：叶片/轮盘结构是典型的循环对称结构，利用其循环对称性进行模态分析能极大地提高解算效率。ANSYS提供了极其完善的循环对称模态分析功能（除此之外，还能进行循环对称静力分析和屈曲分析），包括专门的分网工具、定义命令和菜单等，在常规静态/预应力模态分析流程的基础上仅增加3个命令即可完成循环对称静态/预应力循环对称模态分析。同

时，ANSYS很多独特特点使循环对称模态分析更具实用性，比如，它可以不要求扇区上下截面的网格完全匹配（这在很多情况下是难以做到的）；

要具有模态综合分析能力：对于象转子系统以及整机系统这样的系统级模态分析，如果没有模态综合分析能力，在设计上哪怕仅仅只对某个部件做一定的修改，也需要整个系统重新进行模态计算。ANSYS的模态综合分析（CMS）能力可以高效地解决这个问题：各零部件独立进行模态分析，通过CMS形成整体模态，单一部件的修改仅需重新计算该部件的模态并重新综合即可。

综上所述，ANSYS为航空燃气涡轮发动机发动机的各种零部件模态分析提供了非常完善的手段，可综合考虑各种影响模态分析的因素，为准确地进行其它动力响应分析提供了可靠的保障。

3.2.2.谐响应分析

谐响应分析用以确定结构在呈正弦变化的周期性激励作用下的动力响应（变形、应力等），也就是通常所说的共振响应分析。航空燃气涡轮发动机中诸多的零部件的振动强度分析都需要使用该分析功能，比如：叶片在周期性尾流作用下的响应计算、旋转失速状态下的气流激振、转子系统不平衡响应计算、支承系统及整机系统谐响应计算，等等。

ANSYS谐响应分析提供了一套完整的技术体系满足航空发动机各种谐响应

分析的需求：

丰富的算法：ANSYS谐响应分析提供全方法、模态叠加法、凝聚法三种计算模式；

能考虑各种阻尼影响：三种计算模式均可考虑各种阻尼效应；

能考虑预应力影响：三种计算计算模式均可考虑预应力影响，这对转动部件非常重要。

3.2.3.响应谱和随机振动分析

航空燃气涡轮发动机在飞行过程中会受到诸多复杂载荷的作用，比如飞机传给发动机的复杂载荷谱、由于叶片对气流的干扰和气流振荡燃烧等因素引起的噪声激励等，这些激励载荷的频率非常宽广，激励载荷并没有一个明确的载荷－时间历程，或载荷－时间历程由于太复杂而无法一一计算，这就需要用到谱分析技术。谱分析的实质就是将一个由复杂“载荷－时间历程”所需的规模庞大的时域瞬态计算问题转化为一个相对简单的频域“响应谱”分析问题，响应谱分析所获得的是结构在载荷作用下动力响应的最大值或响应的概率值（随机振动分析）。

ANSYS响应谱分析功能同样提供了一套完整的计算技术，分析方法包括单点

响应谱、多点响应谱、PSD（功率谱密度）随机振动计算等。谱的类型可以是力、位移、速度、加速度等，且在分析中可以考虑包括材料阻尼、恒定阻尼比、模态阻尼等多种阻尼因素的影响。

3.2.

4.气弹颤振分析

气弹颤振属于“自激振动”，由气动力、结构（弹性）力、以及惯性力三者相互作用而引起。颤振时，叶片的振型及频率都大致与自由振动的情况相同（包括弯曲及扭转振动）。它与前述谐响应、响应谱等强迫振动分析的不同之处在于它不伴有任何带频率（扰频）的激振力，颤振的频率基本由叶片自身的几何尺寸及材料性质所决定。航空燃气涡轮发动机叶片的颤振通常有这样五类：亚音速失速颤振、亚音速非失速颤振、超音速失速颤振、超音速非失速颤振、堵塞颤振等，其中最重要的是亚音速失速颤振。这些颤振形式的共同特点都是在来流达到某个特定速度（颤振速度）时，流体对叶片的振动不再是起阻尼作用而是提供维持或增强叶片振动的激励力。

颤振的发生是非常危险的，ANSYS提供了一套完整的预测叶片颤振速度（颤振裕度）的技术，它用Structural结构分析模块计算叶片结构的自振特性，用LinFlow边界元法流体分析模块计算流体激振力，并用该模块提供的两种气弹颤振计算方法计算颤振速度：V-g法（也叫K方法或美国方法）和P- k法（也叫英国方法），这两种方法都是类似的，通过迭代方法求解临界颤振速度。

3.2.5.转子动力学分析

高速旋转状态下转子的振动形态是一种复合的运动状态：转子既要绕转轴本身的轴线旋转，同时又要绕两轴承中心的连线旋转，前者称为自转，后者称为公转或进动，这种复合运动的总称为涡动。在转子不平衡力驱动下，转子一般作同步正进动，当转子进动转速等于转子固有频率时，转子出现共振,相应转速就称为该转子的临界转速。

虽然转子动力学计算与转子系统在旋转预应力条件下的模态固有频率计算有相似之处，但计算转子临界转速必须要考虑旋转结构涡动时产生的陀螺效应（哥氏惯性力矩与离心惯性力矩）对转子临界转速的影响（其影响程度也随转子的具体结构而异），这是转子临界转速计算同其他非旋转结构固有频率计算的差异所在。

一般有限元程序均不具备计算转子临界转速的功能，但ANSYSStructural 通用结构力学分析模块则提供了两种临界转速计算方法供用户选择：

第一种方法是将转子系统简化为梁结构，用BEAM4或PIPE16单元建模，这两种单元都可输入其自转角速度并提供陀螺阻尼矩阵计算选项，直接考虑陀螺效应；

第二种方法不用任何的简化，直接以三维实体的方式计算转子临界转速。这主要得益于ANSYS的两项独特功能：基于部件的旋转速度施加（针对多转子系统的临界转速计算）和旋转软化（SpinSoft）计算选项。旋转软化实际上就是“动刚度法临界转速计算”，直接考虑了陀螺效应对转子刚度的影响。

对于支承刚度和阻尼，这两种方法都可用ANSYSCombin14弹簧阻尼单元来建模，直接输入刚度和阻尼即可。在支承刚度和阻尼未知的情况下，也可直接在ANSYS中予以解算：前者对支承结构（包括轴承）进行静强度（刚度）计算即可完成，对于后者，由于通常都是轴承的油膜阻尼，ANSYS专门提供了独具特色的油膜单元FLUID136、FLUID138和FLUID139，分别用于计算油模的挤压、粘性和滑动阻尼。

3.3.高度非线性冲击动力分析

航空燃气涡轮发动机在高度非线性冲击结构动力学方面的分析需求主要有三个方面：发动机鸟撞过程模拟、机匣的叶片包容性分析、以及发动机抗异物损伤能力计算。

3.3.1.航空发动机鸟撞过程仿真

据概略统计，全球每年大约发生1万次鸟撞飞机事件，国际航空联合会已把鸟害升级为“A”类航空灾难。鸟撞对航空器尤其是发动机的破坏是灾难性的：一只 0.45千克重的飞鸟与时速500公里的飞机相撞时，冲击力为8000多公斤，足以使发动机叶片或外罩等严重变形或断裂，造成灾难性后果。因此，航空发动机的抗鸟撞性能一直是设计者追求的目标，现代发动机研制时都要进行鸟撞物理试验，从最初的每只0.8公斤的单鸟撞击到现在的每只3.6公斤的多鸟撞击试验，通常要求发动机在鸟撞后仍保持足够的推力和能够继续飞行约20分钟。对于现代设计来说，由于鸟撞物理试验费用昂贵、周期长、难以提供足够的信息，因而通常都只是验证性的试验，在研究过程中，利用计算机进行有限元仿真模拟的技术得到广泛的采用，并成为主要的研究手段。

由于极高的撞击速度，飞鸟撞击发动机的过程发生在很短时间内，一般为50毫秒左右，此过程中飞鸟肌体将发生流动变形和解体而四处抛洒，发动机结构亦将产生大变形，甚至发生破坏，例如发动机风扇叶片断裂等。同时，结构的动态响应将在较长时间内持续发生，但令人感兴趣的时间段一般不超过100

毫秒。利用计算机模拟鸟撞时，分析程序必须具备的功能是：

飞鸟物理材料的描述；

飞鸟流动变形以及破碎的描述；

飞鸟与飞行器接触的描述；

飞行器结构大变形、侵彻和破坏过程的描述；

针对大型数值问题的求解能力。

当前，世界范围内对鸟撞进行分析仿真所广泛采用的工具为ANSYSLS-DYNA，LS-DYNA是著名高度非线性有限元显式求解程序，主要用于分析结构在高速撞击、爆炸等动载荷下的动态响应，同时具有强大的流体功能，可进行流体/固体耦合分析。包括Boeing、GEAircraftEngines、 Pratt&Whitney、Rolls-Royce、LockheedMartin、LoralVoughtSystems等在内的著名航空航天器及动力设备厂商，以及包括JetPropulsionLab、NASAGlennResearch、MIT、Oxford等在在内的高级研究机构，都在采用LS-DYNA进行包括鸟撞在内的发动机/飞行器结构动力及安全性分析。

飞鸟在高速撞击时

将产生强大压力，足以使金属材料发生变形和破坏。在这样的变形条件下，飞鸟的材料呈流体。LS-DYNA中的飞鸟材料采用流体动力材料，此种材料除定义一般材料性质，如密度、粘度外，附加的状态方程用于定义其流体属性，如可压缩性，

飞鸟破碎的参数等。

以前，人们在进行鸟撞问题分析或实验时主要关注结构（飞行器）的变形和响应，对飞鸟变形过程不够重视。导致这样问题可能的原因主要包括：1.早年程序中单元的变形能力不足以模拟飞鸟的流动和破碎过程；2.从鸟撞的整个物理过程来看，飞鸟对结构是唯一载荷，但载荷的大小不仅决定于飞鸟的动能，还与其流动过程以及破碎的时间密切相关。即正确描述飞鸟的流动和破碎过程对整

个分析至关重要。以前的研究对此认识有欠缺。LS-DYNA提供三种方式：采用LAGRANGE或ALE（任意欧拉-拉格朗日）；

采用EULER两种类型单元描述飞鸟的流动和破碎；

采用SPH光顺质点描述飞鸟的流动变形和破碎。

LAGRANGE或ALE的变形能力很大，足以描述与结构分离前的变形；EULER

单元可正确描述任意程度的变形；SPH算法可准确描述飞鸟的解体和碎裂。

LS-DYNA在处理飞鸟与飞行器的接触过程中亦提供三种方式：当采用LAGRANGE（或ALE）描述时，使用结构/结构接触算法；

当采用EULER描述时，使用流体/结构耦合算法；

当采用SPH描述时，使用SPH通用接触算法。

对于结构（主要是发动机叶片），可使用LS-DYNA附加破坏算法的结构材料，例如弹塑性破坏材料（金属叶片、金属外罩）或可考虑失效的叠层复合材料（如复合材料风扇叶片、复合材料外罩）等。

3.3.2.航空发动机叶片包容性分析

为了提高推进效率，现代大型航空发动机通常都由核心机和风扇系统构成，由于涵道比不断提高，风扇叶片也越来越大，通常每个风扇叶片重达20～25公斤，风扇叶尖速度可达300～600米/秒。同时，由于风扇叶片最容易受到鸟撞和异物冲击，严重时叶片发生断裂，释放出极大的动能。按照适航性要求，断裂的风扇叶片必须要包容在风扇机匣内，否则，叶片穿透机匣后一旦击中飞机油路系统或机舱，灾难性事故就难以避免。

发动机风扇叶片和风扇机匣通常由特种复合材料或钛合金等构成，要模拟这种高速冲击和侵蚀过程，要求软件必须能：

具有带失效的、能描述复合材料在高速状态下的特性的材料模式；

具有良好的变形体接触和侵蚀接触算法。

LS-DYNA具有7种带失效的复合材料模式和10种带失效的金属材料模式；53种变形体接触和3种侵蚀接触算法，综合运用这些功能，可以良好地模拟叶片与机匣的高速碰撞和破坏过程。

另外，为了增强风扇机匣的包容效果，通常会在机匣上缠绕高强度纤维织物

（如凯芙拉碳纤维等），这些质量轻、强度高、韧性好的纤维材料能有效地包容穿透了金属机匣的叶片和机匣碎片。这些纤维织物在LS-DYNA里也可进行有效的模拟，前一页右下图即为ARA公司用LS-DYNA分析的一种用于发动机包容的纤维织物结构，该研究项目由FAA资助。

3.3.3.发动机异物损伤（FOD）

航空发动机在工作过程中，由于强大气流的作用等因素，在起飞或降落过程中，往往会吸入诸如小石子、飞机轮胎碎片、螺钉螺帽等（法航协和飞机失事就是因为发动机吸入轮胎碎片造成的），高速转动的发动机部件一旦与这些异物相撞，会发生严重变形甚至破坏，与叶片包容的算法类似，LS-DYNA对发动机的FOD设计也能提供很好的仿真结果，右图即为牛津大学工程科学系用LS-DYNA 分析的一个FOD问题。

3.4.疲劳寿命分析

疲劳寿命分析在航空燃气涡轮发动机设计分析中占有非常重要的地位，实际上，在投入使用的发动机所出现的结构问题中，疲劳破坏占了相当大的一部分，发动机的叶片、轮盘、轴、传动系统等几乎所有重要构建都存在疲劳计算的问题。发动机零部件疲劳破坏可粗略地分为两大类：高周疲劳（应力疲劳）和低周疲劳（应变疲劳），前者通常发生在高振动、低应力、低温零件上，如压气机叶片、管路系统、轴及传动系统等，而后者则主要发生在高应力、高温零件上，如涡轮叶片、轮盘、燃烧室等。

材料疲劳破坏的机理相当复杂，涉及固体物理、金属学、固体力学等广泛的学术领域，具有某些边缘学科的性质。适用的疲劳寿命预测技术主要侧重于宏观领域的研究，由于影响因素众多，要想比较准确地进行寿命预测，有如下几个方面的要求：

准确的材料固有疲劳特性参数（应力－寿命曲线、应变－寿命曲线等）

准确的零件“应力－时间”或“应变－时间”历程

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疲劳计算，进一部验证产品的设计寿命和可靠性。

ANSYS复合材料仿真分析及其在航空领域的应用

ANSYS复合材料仿真分析及其在航空领域的应用 复合材料，是由两种或两种以上性质不同的材料组成。主要组分是增强材料和基体材料。复合材料不仅保持了增强材料和基体材料本身的优点，而且通过各相组分性能的互补和关联，获得优异的性能。复合材料具有比强度大、比刚度高、抗疲劳性能好、各向异性、以及材料性能可设计的特点，应用于航空领域中，可以获得显著的减重效益，并改善结构性能。目前，复合材料技术已成为影响飞机发展的关键技术之一，逐渐应用于飞机等结构的主承力构件中，西方先进战斗机上复合材料使用量已达结构总重量的25%以上。飞机结构中，复合材料最常见的结构形式有板壳、实体、夹层、杆梁等结构。板壳结构如机翼蒙皮，实体结构如结构连接件，夹层结构如某些薄翼型和楔型结构，杆梁结构如梁、肋、壁板。此外，采用缠绕工艺制造的筒身结构也可视为层合结构的一种形式。一.复合材料设计分析与有限元方法复合材料层合结构的设计，就是对铺层层数、铺层厚度及铺层角的设计。采用传统的等代设计（等刚度、等强度）、准网络设计等设计方法，复合材料的优异性能难以充分发挥。在复合材料结构分析中，已经广泛采用有限元数值仿真分析，其基本原理在本质上与各向同性材料相同，只是离散方法和本构矩阵不同。复合材 料有限元法中的离散化是双重的，包括了对结构的离散和每一铺层的离散。这样的离散可以使铺层的力学性能、铺层方向、铺层形式直接体现在刚度矩阵中。有限元分析软件，均把增强材料和基体复合在一起，讨论结构的宏观力学行为，因此可以忽略复合材料的多相性导致的微观力学行为，以每一铺层为分析单元。二.ANSYS复合材料仿真技术及其在航空领域应用复合材料具有各向异性、耦合效应、层间剪切等特殊性质，因此复合材料结构的精确仿真，已成为现代航空结构的迫切需求。许多CAE程序都可以进行复合材料的分析，但是大多程序并没有提供完备的功能，使复合材料的精确仿真难以完成。如有些程序不提供非线性分析能力，有些不提供层间剪切应力的求解能力，有些不提供考虑材料失效破坏继续计算能力等等。ANSYS作为一款著名的商业化大型通用有限元软件，广泛应用于航空航天领域，为飞机结构中的复合材料层合结构分析提供了完整精确的解决方案。1.复合材料的有限元模型建立针对飞机结构中的复合材料层合板、梁、实体以及加筋板等结构类型，ANSYS提供一种特殊的复合材料单元———层单元，以模拟各种复合材料，铺层数可达250层以上，并提供一系列技术模拟各种复杂层合结构。复合材料层单元支持非线性、振动特性、热应力、疲劳断裂等各种结构和热的分析功能和算法。2.复合材料的层合结构定义：■铺层结构：ANSYS对于每一铺层可先定义材料性质、铺层角、铺层厚度，然后通过由下到上的顺序逐层叠加组合为复合材料层合结构；也可以通过直接输入材料本构矩阵来定义复合材料性质。■板壳和梁单元截面形状：ANSYS利用截面形状工具可定义矩形、I型、槽型等各种形式；还可以定义各种函数曲线以模拟变厚度截面。3.特殊层合结构的模拟：?变厚度板壳铺层切断：将切断的某铺层厚度定义为零，即可模拟铺层切断前后的板壳实际形状。（图1上）?不同铺层板壳的节点协调：ANSYS板壳层单元的节点均可偏置到任意位置，使不同铺层数板壳的节点在中面或顶面、底面对齐。（图1下）?蜂窝/泡沫夹层结构：ANSYS通过板壳层单元来模拟夹层结构的特性，夹层面板和芯子可以是不同材料。（图2）?板－梁－实体组合结构：ANSYS将实体、板壳与梁等不同类型单元通过MPC技术相联系，各类单元的节点不需要重合并协调，便于飞机等复杂结构模型的处理。4.复合材料有限元模型的检查：复合材料结构模型建立后，可以将板壳和梁单元显示为实际形状，还可以通过图形显示和列表直观地观察铺层厚度、铺层角度和铺层组合形式，方便模型的检查及校对。（图3）5.复合材料层合结构分析ANSYS层单元支持各种静强度刚度、非线性、稳定性、疲劳断裂和振动特性等结构分析。完成分析后，可以图形显示或输出每个铺层及层间的应力和应变等结果（虽然一个单元包含许多铺层），根据这些结果可以判断结构是否失效破坏和满足设计要求。6.复合材料失效准则ANSYS已经预定义了三种复合材料破坏准则来评价复合材料结构安全性，包括最大应变/应力失效准则，蔡－吴（Tsai-Wu）准则。每种强度准则均可定义与温度相关，考虑不同温度下的材料性能。另外，用户也可自定义最多达六种的

ANSYS 有限元分析 平面薄板

《有限元基础教程》作业二：平面薄板的有限元分析 班级：机自101202班 姓名：韩晓峰 学号：201012030210 一．问题描述： P P h 1mm R1mm 10m m 10mm 条件：上图所示为一个承受拉伸的正方形板，长度和宽度均为10mm ，厚度为h 为1mm ，中心圆的半径R 为1mm 。已知材料属性为弹性模量E=1MPa ，泊松比为0.3，拉伸的均布载荷q ＝ 1N/mm 2。根据平板结构的对称性，只需分析其中的二分之一即可，简化模型如上右图所示。 二．求解过程： 1 进入ANSYS 程序 →ANSYS 10.0→ANSYS Product Launcher →File management →input job name: ZY2→Run 2设置计算类型 ANSYS Main Menu: Preferences →select Structural → OK 3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK → Options… →select K3: Plane Strs w/thk →OK →Close 4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX: 1e6, PRXY:0.3 → OK 5定义实常数以及确定平面问题的厚度 ANSYS Main Menu: Preprocessor →Real Constants …→Add/Edit/Delete →Add →Type 1→OK →Real Constant Set No.1,THK:1→OK →Close 6生成几何模型 a 生成平面方板 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Rectangle →By 2 Corners →WP X:0,WP Y:0,Width:5,Height:5→OK b 生成圆孔平面 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Circle →Solid Circle →WPX=0,WPY=0,RADIUS=1→OK b 生成带孔板 ANSYS Main Menu: Preprocessor →Modeling →Operate →Booleans → Subtract →Areas →点击area1→OK →点击area2→OK 7 网格划分 ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →(Size Controls) Global: Set →SIZE: 0.5 →OK →iMesh →Pick All → Close 8 模型施加约束

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

ansys有限元分析工程实例大作业

ansys有限元分析工程实例大作业

————————————————————————————————作者：————————————————————————————————日期：

辽宁工程技术大学 有限元软件工程实例分析 题目基于ANSYS钢桁架桥的静力分析专业班级建工研16-1班（结构工程）学号 471620445 姓名 日期 2017年4月15日

基于ANSYS钢桁架桥的静力分析 摘要：本文采用ANSYS分析程序，对下承式钢桁架桥进行了有限元建模;对桁架桥进行了静力分析，作出了桁架桥在静载下的结构变形图、位移云图、以及各个节点处的结构内力图（轴力图、弯矩图、剪切力图），找出了结构的危险截面。 关键词：ANSYS；钢桁架桥；静力分析；结构分析。 引言：随着现代交通运输的快速发展，桥梁兴建的规模在不断的扩大，尤其是现代铁路行业的快速发展更加促进了铁路桥梁的建设，一些新建的高速铁路桥梁可以达到四线甚至是六线，由于桥面和桥身的材料不同导致其受力情况变得复杂，这就需要桥梁需要有足够的承载力，足够的竖向侧向和扭转刚度，同时还应具有良好的稳定性以及较高的减震降噪性，因此对其应用计算机和求解软件快速进行力学分析了解其受力特性具有重要的意义。 1、工程简介 某一下承式简支钢桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表1，材料属性见表2。桥长32米，桥高5.5米，桥身由8段桁架组成，每个节段4米。该桥梁可以通行卡车，若只考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg，若取一半的模型，可以将卡车对桥梁的作用力简化为P1，P2，和P3，其中P1=P3=5000N，P2=10000N,见图2,钢桥的形式见图1，其结构简图见图3。

ansys大作业ANSYS电磁场分析及与ansoft仿真分析结果比较要点

期末大作业 题目：简单直流致动器 ANSYS电磁场分析及与ansoft仿真分析结果比较作者姓名：柴飞龙 学科（专业）：机械工程 学号：21225169 所在院系：机械工程学系 提交日期2013 年 1 月

1、 背景简述： ANSYS 软件是融结构、流体、电场、磁场、声场分析于一体的大型通用软件有限元分析软件，是现代产品设计中的高级CAE 工具之一。而ansoft Maxwell 软件是一款专门分析电磁场的分析软件，如传感器、调节器、电动机、变压器等。 本人在实验室做的课题涉及到电机仿真，用的较多的是ansoft 软件，因为其对电机仿真的功能更强大，电机功能模块更多，界面友好。 现就对一电磁场应用实例，用ANSYS 进行仿真分析，得到的结果与ansoft 得到的结果进行简单核对比较。 2、 问题描述： 简单直流致动器由2个实体圆柱铁芯，中间被空气隙分开的部件组成，线圈中心点处于空气隙中心。衔铁是导磁材料，导磁率为常数（即线性材料，r μ=1000），线圈是可视为均匀材料，空气区为自由空间（1=r μ），匝数为2000，线圈励磁为直流电流：2A 。模型为轴对称。 3、 ANSYS 仿真操作步骤： 第一步：Main menu>preferences

第二步：定义所有物理区的单元类型为PLANE53 Preprocessor>Element type>Add/Edit/Delete 第三步：设置单元行为 模拟模型的轴对称形状，选择Options（选项） 第四步：定义材料 Preprocessor>Material Props> ?定义空气为1号材料(MURX = 1) ?定义衔铁为2号材料(MURX = 1000) ?定义线圈为3号材料(自由空间导磁率，MURX=1)

ansys有限元分析大作业

ansys有限元分析大作业

有限元大作业 设计题目: 单车的设计及ansys有限元分析 专业班级: 姓名: 学号: 指导老师: 完成日期: 2016.11.23

单车的设计及ansys模拟分析 一、单车实体设计与建模 1、总体设计 单车的总体设计三维图如下，采用pro-e进行实体建模。 在建模时修改proe默认单位为国际主单位(米千克秒 mks) Proe》文件》属性》修改

2、车架 车架是构成单车的基体，联接着单车的其余各个部件并承受骑者的体重及单车在行驶时经受各种震动和冲击力量，因此除了强度以外还应有足够的刚度，这是为了在各种行驶条件下，使固定在车架上的各机构的相对位置应保持不变，充分发挥各部位的功能。车架分为前部和后部，前部为转向部分，后部为驱动部分，由于受力较大，所有要对后半部分进行加固。

二、单车有限元模型 1、材料的选择 单车的车身选用铝合金（6061-T6）T6标志表示经过热处理、时效。 其属性如下： 弹性模量：） .6+ 90E （2 N/m 10 泊松比：0.33 质量密度：） 3 2.70E+ N/m （2 抗剪模量：） 60E .2+ N/m （2 10 屈服强度：） .2+ （2 75E 8 N/m 2、单车模型的简化 为了方便单车的模拟分析，提高电脑的运算

效率，可对单车进行初步的简化；单车受到的力的主要由车架承受，因此必须保证车架能够有足够的强度、刚度，抗振的能力，故分析的时候主要对车架进行分析。简化后的车架如下图所示。 3、单元体的选择 单车车架为实体故定义车架的单元类型为实体单元（solid）。查资料可以知道3D实体常用结构实体单元有下表。 单元名称说明 Solid45 三维结构实体单元，单元由8个节点定义，具有塑性、蠕变、应力刚化、 大变形、大应变功能，其高阶单元是 solid95

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]对话框。

ansys有限元建模与分析实例-详细步骤

《有限元法及其应用》课程作业ANSYS应用分析 学号： 姓名： 专业：建筑与土木工程

角托架的有限元建模与分析 一 、模型介绍 本模型是关于一个角托架的简单加载，线性静态结构分析问题，托架的具体形状和尺寸如图所示。托架左上方的销孔被焊接完全固定，其右下角的销孔受到锥形压力载荷，角托架材料为Q235A 优质钢。角托架材料参数为：弹性模量366E e psi =；泊松比0.27ν= 托架图（厚度：0.5） 二、问题分析 因为角托架在Z 方向尺寸相对于其在X,Y 方向的尺寸来说很小，并且压力荷载仅作用在X,Y 平面上，因此可以认为这个分析为平面应力状态。 三、模型建立 3.1 指定工作文件名和分析标题 （1）选择菜单栏Utility Menu → 命令.系统将弹出Jobname(修改文件名)对话框,输入bracket （2）定义分析标题 GUI ：Utility Menu>Preprocess>Element Type>Add/Edit/Delete 执行命令后，弹出对话框，输入stress in a bracket 作为ANSYS 图形显示时的标题。 3.2设置计算类型 Main Menu: Preferences … →select Structural → OK 3.3定义单元类型 PLANE82 GUI ：Main Menu →Preprocessor →Element Type →Add/Edit/Delete 命令，系统将弹出Element Types 对话框。单击Add 按钮，在对话框左边的下拉列表中单击Structural Solid →Quad 8node 82，选择8节点平面单元PLANE82。单击ok ，Element Types 对话框，单击Option ，在Element behavior 后面窗口中选取Plane strs w/thk 后单击ok 完成定义单元类型。 3.4定义单元实常数 GUI ：Main Menu: Preprocessor →Real Constants →Add/Edit/Delete ，弹出定义实常数对话框，单击Add ，弹出要定义实常数单元对话框，选中PLANE82单元后，单击OK →定义单元厚度对话框，在THK 中输入0.5.

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

ANSYS经典界面与workbench结合进行仿真分析报告

联合ANSYS WORKBENCH和经典界面进行后处理 (2012-12-18 18:58:55) ▼ 标签： ansys 分类：CAE 前面几篇文章已经提到过，ANSYS WORKCENCH主要是为不大懂ANSYS命令和编程的工程师服务的，而经典界面则适用于初学者和研究人员。初学者和研究人员是完全不同的两个层次，为什么ANSYS经典界面却同时适合二者呢？ 实际上，学好ANSYS，关键并非是操作界面，而是要学好有限元。如果初学者直接从WORKBENCH来学习ANSYS，那么对于有限元就毫无收获，可以说一头雾水。而如果从经典界面进去，因为涉及到很多与有限元概念密切相关的操作，对于理解有限元很有好处。只是学到一定程度以后，需要转移到WORKBENCH中进行三维零件的分析和装配体的分析。 而当我们用到一定程度以后，发现WOKRBENCH虽然操作方便，但是的确不容易操作底层。前面的文章已经说明了如何联合二者进行仿真，以充分使用WOKRNBEHCN对于建模的方便性以及经典界面对于底层的操控性。这里再举一个例子，说明如何用WOKRBENCH进行建模，而后在经典界面中进行后处理，目的是为研究人员提供参考。 一个两边固定的梁，上面受到分布载荷作用如下图。 该分布载荷随时间而改变，其载荷的时间历程如下曲线，从0-1秒，载荷增加到1Mpa，而后保持1秒钟，接着减小到0Mpa，终止时间是3秒。 为了便于控制，这里对每个载荷步均采用自定义载荷子步的方式，划分为10个载荷子步，见下面的细节视图。

然后进行瞬态隐式动力学分析，得到该梁的位移和von mises应力。 我们现在要知道该梁上某一个应力最大的点，其应力是如何随时间而改变的。这个任务使用WOKRBENCH很难达到，但是用经典界面则轻而易举，因此我们决定使用经典界面进行后处理。 要使用经典界面后处理，只需要把WORKBENCH中生成的结果文件导入到经典界面中即可。 首先找到WORKBENCH中生成的结果文件如下图所示的路径。该文件叫file.rst， 为了方便，把file.rst拷贝到D盘的根目录下，然后启动ANSYS APDL，即经典界面。 进入经典界面后，直接进入通用后处理，并点击data file opts，以便读入结果文件。

ansys有限元分析作业

有限元分析作业 作业名称输气管道有限元建模分析 姓名邓伟 学号 p1202100706 班级：浦机械1007 题目描述： 输气管道的有限元建模与分析 计算分析模型如图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接受单元类型并关闭对话框。

ANSYS电磁兼容仿真软件解析

ANSYS 电磁兼容仿真设计软件 用途：用于电子系统电磁兼容分析，包括PCB信号完整性、电源完整性和电磁辐射协同仿真，数模混合电路的噪声分析和抑制，以及 机箱系统屏蔽效能和电磁泄漏仿真，确保系统的电磁干扰和电磁兼容性能满足要求。 一、购置理由 1 现代电子系统设计面临越来越恶劣的电磁工作环境，一方面电子系统包括了电源模块、信号处理、计算机控制、传感与机电控制、光电系统及天线与微波电路等部分，系统内部相互不发生干扰，正常工作，本身就非常困难；另一方面，在隐身、电子对抗、静放电，雷击和 电磁脉冲干扰等恶劣电磁环境下，设备还需要有足够的抗干扰能力，为电路正常工作留有足够的设计裕量。为了确保xx 系统的工作可靠性，设备必须通过相关的电磁兼容标准，如国军标 GJB151A，GJB152A。 长期以来，设备的电磁兼容设计和仿真一直缺乏必要的仿真设计 手段，只能依赖于设备后期试验测试，不仅测量成本高昂，而且，如 果EMI 测量超标，后续的查找问题和修正问题基本上依赖于经验和猜测。而解决电磁兼容问题，也只能靠经验进行猜想和诊断，采取的 措施也只能通过不断的试验进行验证，这已经成为制约我们产品进度的重要原因。 2 目前我所数字电路设计的经验和手段已经有很大改善，我们在复杂PCB布线、高速仿真方面取得了很多的成果和经验，并且已经

开始高速通道设计的预研。在相关PCB 布线工具的帮助下，将复杂 的多电源系统PCB布通，确保集成电路之间的正确连接已经基本上 没有问题。但是随着应用深入，也存在一些困难，特别在模拟数字转换、高速计算与传输PCB和系统的设计中，我们不仅要保证电路板 的正常工作，还要提高关键性的技术指标，例如数模转换电路的有效 位数、信号传输系统的速率和误码率等，此外，还要满足整个卫星电子系统的电磁兼容/电磁干扰要求，为此，我们迫切需要建立的仿真功能包括： ● 高速通道中，连接器，电缆等三维全波精确和建模仿真，这 些结构的寄生效应对于信号的传输性能有至关重要的影响； ● 有效的PCB电源完整性分析工具，对PCB 上的电源、地等 直流网络的信号质量进行仿真 ●为提高仿真精度，需要SPICE 模型，IBIS模型和S 参数模 型的混合仿真 ●需要同时进行时域和频域仿真和设计，观察时域的眼图、 误码率，调整预加重和均衡电路的频域参数，使得信号通道 的物理特性与集成电路和收/发预加重、均衡等相配合，达到 系统性能的最优 ● 有效的PCB的辐射控制与仿真手段，确保系统EMI性能 达标。 现在EDA 市场上已经有一些SI/PI 和EMI/EMC 仿真设计工具，但存在多方面的局限性。我们的PCB 布线工具虽然能解决一定的问

ansys有限元分析报告大作业

有限元大作业 设计题目: 单车的设计及ansys有限元分析 专业班级: 姓名: 学号: 指导老师: 完成日期: 2016.11.23

单车的设计及ansys模拟分析 一、单车实体设计与建模 1、总体设计 单车的总体设计三维图如下，采用pro-e进行实体建模。 在建模时修改proe默认单位为国际主单位(米千克秒 mks) Proe》文件》属性》修改

2、车架 车架是构成单车的基体，联接着单车的其余各个部件并承受骑者的体重及单车在行驶时经受各种震动和冲击力量，因此除了强度以外还应有足够的刚度，这是为了在各种行驶条件下，使固定在车架上的各机构的相对位置应保持不变，充分发挥各部位的功能。车架分为前部和后部，前部为转向部分，后部为驱动部分，由于受力较大，所有要对后半部分进行加固。 二、单车有限元模型 1、材料的选择 单车的车身选用铝合金（6061-T6）T6标志表示经过热处理、时效。 其属性如下： 弹性模量：）（2 N/m 1090E .6

泊松比：0.33 质量密度：）（2 N/m 32.70E + 抗剪模量：）（2N/m 1060E .2+ 屈服强度：） （2N/m 875E .2+ 2、单车模型的简化 为了方便单车的模拟分析，提高电脑的运算效率，可对单车进行初步的简化；单车受到的力的主要由车架承受，因此必须保证车架能够有足够的强度、刚度，抗振的能力，故分析的时候主要对车架进行分析。简化后的车架如下图所示。 3、单元体的选择 单车车架为实体故定义车架的单元类型为实体单元（solid ）。查资料可以知道3D 实体常用结构实体单元有下表。

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疲劳计算，进一部验证产品的设计寿命和可靠性。

课程设计ANSYS有限元分析(最完整)

有限元法分析与建模课程设计报告 学院：机电学院 专业：机械制造及其自动化指导教师：**** 学生：* *** 学号：2012011**** 2015-12-31

摘要 本文通过ANSYS10.0建立了标准光盘的离心力分析模型，采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析，同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。力求较为真实地反映光盘在光驱中实际应力和应变分布情况，为人们进行合理的标准光盘结构设计和制造工艺提供理论依据。 关键词：ANSYS10.0；光盘；应力；应变。

目录 第一章引言 (3) 1.1 引言 (3) 第二章问题描述 (4) 2.1有限元法及其基本思想 (4) 2.2 问题描述 (4) 第三章力学模型的建立和求解 (5) 3.1设定分析作业名和标题 (5) 3.2定义单元类型 (6) 3.3定义实常数 (9) 3.4定义材料属性 (12) 3.5建立盘面模型 (14) 3.6对盘面划分网格 (22) 3.7施加位移边界 (27) 3.8施加转速惯性载荷并求解 (30) 第四章结果分析 (32) 4.1 旋转结果坐标系 (32) 4.2查看变形 (33) 4.3查看应力 (35) 总结 (38) 参考文献 (39)

第一章引言 1.1 引言 光盘业是我国信息化建设中发展迅速的产业之一，认真研究光盘产业的规律和发展趋势，是一件非常迫切的工作。光盘产业发展的整体性强，宏观调控要求高，因此，对于光盘产业的总体部署、合理布局和有序发展等问题，包括节目制作、软件开发、硬件制造、节目生产、技术标准等。 在高速光盘驱动器中，光盘片会产生应力和应变，在用ANSYS分析时，要施加盘片高速旋转引起的惯性载荷，即可以施加角速度。需要注意的是，利用ANSYS施加边界条件时，要将内孔边缘节点的周向位移固定，为施加周向位移，而且还需要将节点坐标系旋转到柱坐标系下。 本文通过ANSYS10.0建立了标准光盘的离心力分析模型，采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析，同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。

ANSYS仿真分析作业解析

桥梁结构仿真分析 学院： 姓名： 学号： 指导教师：

一、下承式钢管混凝土拱桥，跨径90m 。主梁截面布置如下图所示，C50混凝土，吊杆为54根直径7mm 钢丝；钢管混凝土拱圈截面如图所示，直径1m 、壁厚14mm 、内填C50混凝土。吊杆下设置壁厚40cm 的混凝土横梁。设计荷载公路I 级，4车道。按照杆系结构计算结构在移动荷载作用下的弯矩和剪力包络图。 立面布置（单位：m ） 拱圈截面布置 主梁截面布置（单位：cm ） 1、计算模型 1）采用单梁模型； 2）拱圈采用共节点分离模型，单拱轴，钢管、拱轴混凝土和吊杆弹性模量取2倍的各自模量，面积不变； 3）吊杆不考虑初应变。 MA TLAB 求解弯矩和剪力包络图（纵桥向最不利加载）。 2、计算结果 102030405060708090 -2 -1.5 -1 -0.5 00.5 1x 10 4 x (m) M (k N m ) 图1-1 主梁弯矩包络图

图1-2 主梁剪力包络图 3、ANSYS 命令流 FINISH /CLEAR /TITLE,Concrete-Filled Tube Arch Bridge /PREP7 ET,1,82 !辅助单元类型 CYL4,0.5,0.5,0.5,,0.486 !建立一个环面 SMRTSIZE,5 !网格划分 AMESH,ALL SECWRITE,STEEL_TUBE,SECT,,1 !截面存盘（文件名：STEEL_TUBE ；文件拓展名：SECT ；单元类型属性：1） SECTYPE,1,BEAM,MESH, !截面类型和ID 定义 SECOFFSET,CENT,,, !截面偏移，CENT 偏移到质心 SECREAD,'STEEL_TUBE','SECT',,MESH !读入截面 ASEL,ALL !选择所有面 ACLEAR,ALL ADELE,ALL,,,1 !将图形及面、线、点全部删除（只保留截面文件） CYL4,0.5,0.5,0.486 !建立一个圆面 SMRTSIZE,5 AMESH,ALL SECWRITE,CONCRETE,SECT,,1 SECTYPE,2,BEAM,MESH, SECOFFSET,CENT,,, SECREAD,'CONCRETE','SECT',,MESH ASEL,ALL ACLEAR,ALL ADELE,ALL,,,1 K,1,-9 K,2,-9,-0.2 K,3,-7,-0.5 K,4,-6.5,-1.8 K,5,-5,-1.8 K,6,-4.5,-0.5 K,7,-4,-0.25 K,8,4,-0.25 K,9,4.5,-0.25 K,10,5,-1.8 K,11,6.5,-1.8 K,12,7,-0.5 K,13,9,-0.2 K,14,9 A,1,2,3,4,5,6,7,8,9,10,11,12,13,14 SMRTSIZE,5 AMESH,ALL SECWRITE,ZL,SECT,,1 SECTYPE,3,BEAM,MESH, SECOFFSET,CENT,,,

ANSYS有限元分析作业

1.工程背景 房屋刚性独立基础 当建筑物上部结构采用框架结构或单层排架结构承重时，基础常采用方形、圆柱形和多边形等形式的独立式基础，是整个或局部结构物下的无筋或配筋基础。本例以独立坡形基础为例。 2.几何参数及材料 底部：3m*3m，全高：1.8m，上部平台面积：0.6m*0.6m，斜坡坡高：1.2m，坡脚：45°，截面为正方形，选取1/2的单向简化模型。 桩体材料：线弹性材料，弹性模量GPa，泊松比0.2，密度2700kg/m3 土体材料：DP材料，弹性模量25MPa，泊松比0.45，密度2000kg/m3，粘聚力10，摩擦角30°，膨胀角30° 3.建模过程 （1）前处理 1——定义单元类型及材料属性 2——建立平面模型 3——进行网格划分 4——拉伸成体 （2）加载及求解 1——施加约束(整体模型的对称面X=0处施加对称约束，模型底面Y=-2施加全自由度约束，顶面为自由面，其余三个侧面约束其平面外的平动自由度) 2——施加重力荷载 3——施加上部约束 （3）后处理 1——自重荷载下的受力及变形 2——施加约束后的结果 4.命令流 /CLEAR /prep7 et,1,plane182 et,2,solid65 mp,ex,1,2.5e10 !桩的弹性模量 mp,nuxy,1,0.2 !桩的泊松比 mp,dens,1,2700 !桩的密度 mp,ex,2,2.5e7 !土的弹性模量 mp,nuxy,2,0.45 !土的泊松比

mp,dens,2,2000 !土的密度 tb,dp,2 tbdata,1,10,30,30 !粘聚力c为10,摩擦角为30度,膨胀角为30 !keypoints k,1 ！建立模型关键点 k,2,1.5 k,3,1.5,0.3 k,4,0.3,1.5 k,5,0.3,1.8 k,6,0,1.8 k,7,1.5,1.8 k,8,4,1.8 k,9,4,0 k,10,4,-2 k,11,0,-2 *do,i,1,5 ！连接关键点成线 l,i,i+1 *enddo l,1,6 l,3,7 l,5,7 *do,i,7,10 l,i,i+1 *enddo l,2,9 l,11,1 /pnum,line,1 lplot al,1,2,3,4,5,6 ！显示直线编号 al,3,4,8,7 ！绘制直线 al,2,7,9,10,13 ！围成基础面 al,1,13,11,12,14 ！生成土体面 /pnum,area,1 aplot aglue,all ！显示面 nummrg,all ！粘贴各部分 numcmp,all lsel,s,,,2,8,1 ！选择直线 lsel,u,,,3 ！去除L3 lesize,all,0.15 ！设定划分尺寸 lsel,s,,,3 ！选择L3 lesize,all,,,10 ！分10份 lsel,s,,,all ！反选L2-L8 lsel,u,,,2,8,1