有限元分析50年发展

1、有限元分析50年发展之路计算机辅助工程（CAE）作为一门新兴的学科已经逐渐的走下神坛，成为了各大企业中设计新 产品过程中不可缺少的一环。传统的CAE技术是指工程设计中的分析计算与分析仿真，具 体包括工程数值分析、结构与过程优化设计、强度与寿命评估、运动/动力学仿真，验证未 来工程/产品的可用性与可靠性。如今，随着企业信息化技术的不断发展，CAE软件与CAD/CAM/CAPP/PDM/ERP 一起，已经 成为支持工程行业和制造企业信息化的主导技术，在提高工程/产品的设计质量，降低研究 开发成本，缩短开发周期方面都发挥了重要作用。而CAE技术出现则是要归功于有限元分 析的诞生，在有限元法诞生的早期

2、，几乎所有的CAE软件都是使用有限元法来进行计算求 解。因此，可以说有限元法的发展也间接反映了 CAE软件在这半个世纪的发展历史。1有限元法的诞生每一项新技术的推出都是由于时代的迫切需要，而新技术的出现后也需要经历历史的重 重考验。在上个世纪40年代，由于航空事业的快速发展，对飞机内部结构设计提出了越来 越高的要求，即重量轻、强度高、刚度好，人们不得不进行精确的设计和计算。正是在这一 背景下，有限元分析的方法逐渐的发展起来。早期的一些成功的实验求解方法与专题论文，完全或部分的内容对有限元技术的产生做 出的贡献，首先在应用数学界第一篇有限元论文是1943年Courant R发表的Variatio

3、nal methods for the solution of problems of equilibrium and vibration一文，文中描述了他使用三 角形区域的多项式函数来求解扭转问题的近似解，由于当时计算机尚未出现，这篇论文并没 有引起应有的注意。1956年，M.J.Turner （波音公司工程师），R.W.Clough （土木工程教授），H.C.Martin （航空 工程教授）及LJ.Topp （波音公司工程师）等四位共同在航空科技期刊上发表一篇采用有限元 技术计算飞机机翼的强?的论文，名为Stiffness and Deflection Analysis of Comple

4、x Structures，文中把这种解法称为刚性法（Stiffness），一般认为这是工程学界上有限元法的 开端。1960年，Ray W. Clough教授在美国土木工程学会（ASCE）之计算机会议上，发表另一篇 名为The Finite Element in Plane Stress Analysis的论文，将应用范围扩展到飞机以外之土 木工程上，同时有限元法（Finite Element Method）的名称也第一次被正式提出。由此之后，有限元法的理论迅速地发展起来，并广泛地应用于各种力学问题和非线性 问题，成为分析大型、复杂工程结构的强有力手段。并且随着计算机的迅速发展，有限元法 中人工

5、是难以完成的大量计算工作能够由计算机来实现并快速地完成。因此，可以说计算机 的发展很大程度上促进了有限元法的建立和发展。2由理论到程序的转变2.1 FEA技术的探索起源期有限元法”概念的提出，引出了美国加州大学伯克利分校有限元技术研究小组的最为辉 煌的十年历程。1963年在加州大学Berkeley分校，Edward L.Wilson教授和Ray W. Clough教授为了教授 结构静力与动力分析而开发了 SMIS(Symbolic Matrix Interpretive System)，其目的是为了 弥补在传统手工计算方法和结构分析矩阵法之间的隔阂。1969年，Wilson教授在第一代程 序的

6、基础上开发的第二代线性有限元分析程序就是著名的SAP(Structural analysis program)， 而非线性程序则为NONSAPoWilson 教授的学生 Ashraf Habibullah 于 1978 年创建了 Computer and Structures Inc.(CSI)， CSI的大部分技术开发人员都是Wilson教授的学生，并且Wilson教授也是CSI的高级技术发 展顾问。而SAP2000则是由CSI在SAP5、SAP80、SAP90的基础上开发研制的通用结构分析 与设计软件。同样是1963年，Richard MacNeal博士和Robert Schwendler

7、先生联手创办了 MSC公司， 并开发第一个软件程序，名为 SADSAM (Structural Analysis by Digital Simulation of Analog Methods )即数字仿真模拟法结构分析。提到MSC公司，就想到与其有着不解渊源的美国国家太空总署(NASA)，当年美国为了能 够在与前苏联之间的太空竞赛中取得优胜而成立了 NASAo而为了满足宇航工业对结构分析 的迫切需求，NASA于1966年提出了发展世界上第一套泛用型的有限元分析软件 Nastran(NASA STRuctural ANalysis Program)的计划，MSC.Software 则参与了整个

8、 Nastran 程序 的开发过程。1969年NASA推出了其第一个NASTRAN版本，称为COSMIC Nastran。之后 MSC继续的改良Nastran程序并在1971年推出MSC.Nastran。另一个与NASA结缘的是SDRC公司，1967年在NASA的支持下SDRC公司成立，并于1968 年发布了世界上第一个动力学测试及模态分析软件包，1971年推出商业用有限元分析软件 Supertab(后并入I-DEAS软件中，这也就是为什么I-DEAS作为一款设计软件其有限元分析还 如此强大的原因)。1969 年，John Swanson 博士建立了自己的公司 Swanson Analysis

9、 Systems Inc(SASI)。其实 早在1963年John Swanson博士任职于美国宾州匹兹堡西屋公司的太空核子实验室时，就已 经为核子反应火箭作应力分析编写了一些计算加载温度和压力的结构应力和变位的程序，此 程序当时命名为STASYS (Structural Analysis SYStem)。在Swanson博士公司成立的次年，结合 者早期的 STASYS 程序发布了商用软件 ANSYSo 1994 年 Swanson Analysis Systems, Inc.被 TA Associates并购，并宣布了新的公司名称改为ANSYS。2.2 FEA技术的蓬勃发展期进入70年代后，

10、随着有限元理论的趋于成熟，CAE技术也逐渐进入了蓬勃发展的时期， 一方面MSC，ANSYS，SDRC三大CAE公司先后组建，并且致力于大型商用CAE软件的研究 与开发，另一方面，更多的新的CAE软件迅速出现，为CAE市场的繁荣注入了新鲜血液。70年代初当时任教于Brown大学的Pedro Marcal创建了 MARC公司，并推出了第一个商 业非线性有限元程序MARC。虽然在MARC在1999年被MSC公司收购，但其对有限元软件 的发展起到了决定性的推动作用，至今在MSC的分析体系中依然有着MARC程序的身影， 更值得一提的是Pedro Marcal早年也是毕业于Berkeley大学。在早期的商

11、用软件舞台上，还有两位主要人物，他们是David Hibbitt和Klaus J Bathe。David Hibbitt 是 Pedro Marcal 在 Brown 的博士生，David Hibbitt 与 Pedro Marcal 合作到 1972 年，随后Hibbitt与Bengt Karlsson和Paul Sorenson于1978年共同建立HKS公司，推出了 Abaqus软件，使ABAQUS商业软件进入市场。因为该程序是能够引导研究人员增加用户单 元和材料模型的早期有限元程序之一，所以它对软件行业带来了实质性的冲击。2002年HKS 公司改名为ABAQUS，并于2005年被达索公司

12、收购。另外一位对有限元方法做出重大贡献的是Klaus J. Bathe博士。Klaus J. Bathe六十 年代末在Berkeley大学Clough和Ed Wilson博士的指导下攻读博士学位，从事结构动力学 求解算法和计算系统的研究。由于Bathe博士在对结构计算以及SAP软件所做的贡献，Bathe 博士毕业后被MIT聘请到机械与力学学院任教至今。1975年在MIT任教的Bathe博士在NONSAP的基础上发表了著名的非线性求解器ADINA (Automatic Dynamic Incremental Nonlinear Analysis)，而在 1986 年 ADINA R&D Inc.

13、 成 立以前，ADINA软件的源代码是公开的，即著名的ADINA81版和ADINA84版本的fortran源 程序，后期很多有限元软件都是根据这个源程序所编写的。1977年Mechanical Dynamics Inc.(MDI)公司成立，致力于发展机械系统仿真软件，其 软件ADAMS应用于机械系统运动学、动力学仿真分析。后被MSC公司收购，成为MSC分析体 系中一个重要的组成部分。在CAE的历史中另一个神奇的程序是显式有限元程序DYNA，DYNA程序由当时在美国 Lawrence Livermore国家实验室的John Hallquist编写。之所以说DYNA神奇，是因为在 现在我们熟知的众

14、多软件中，都可以发现DYNA的踪迹，因此LS-DYNA系列也被公认为显式 有限元程序的鼻祖。下面我们来细数一下由DYNA所演变出来的有限元程序：在20世纪80年代，DYNA程序首先被法国ESI公司商业化，命名为PAM-CRASH，现已成 为了 ESI的明星产品。除此之外，ESI公司还有多个被人熟知的软件，如铸造软件ProCAST， 钣金软件PAM-STAMP，焊接软件SYSWELD，振动噪声软件VA One，空气动力学软件 CFD-FASTRAN，多物理场软件CFD-ACE+等等。1988 年, John Hallquist 自己创建 LSTC(Livermore Software Techn

15、ology Corporation)公司，发行和扩展DYNA程序商业化版本LS-DYNA0同样是1988年，MSC在DYNA3D的框架下开发了 MSC.Dyna并于1990年发布第一个版本，随 后于1993年发布了著名的MSC.Dytran。另外，ANSYS收购了 Century Dynamics公司，把该公司以DYNA程序开发的高速瞬态动力分 析软件AUTODYN纳入到ANSYS的分析体系中。并且在1996年，ANSYS与LSCT公司合作推出 了 ANSYS/LS-DYNAo1984年，ALGOR公司成立于，总部位于宾州的匹兹堡，ALGOR公司在购买SAP5源程序 和 vizicad 图象处

16、理软件后，同年推出 ALGOR FEAS (Finite Element Analysis System)o 随着有限元技术的日趋成熟，市场上不断有新的公司成立并推出CAE软件，1983年AAC公 司成立，推出COMET程序，主要用于噪声及结构噪声优化分析等领域。随后Computer Aided Design Software Inc推出提供线性静态、动态及热分析的 PolyFEM软件包。1988年 Flomerics公司成立，提供用于空气流及热传递的分析程序。同时期还有多家专业性软件公 司投入专业CAE程序的开发。由此，CAE的分析已经逐渐的扩展到了声学、热传导以及流体 等更多的领域。在早期

17、有限元技术刚刚提出时，其应用范围仅在航空航天领域，且研究的对象也只局限 在线性问题与静力分析。而经过近十年的发展研究，有限元技术的应用范围已经囊括了力学、 热、流体、电磁的自然界四大基本物理场，并且已经发展到多场耦合技术(如表1所示)。 可以说有限元技术经过十年的研究发展，其应用范围与研究对象发生了翻天覆地的变化。2.3 FEA技术的成熟壮大期上世纪90年代至今是CAE技术的成熟壮大时期，这一时期的CAE领域呈现出了大鱼吃 小鱼的市场局面，大的软件公司为了提升自己的分析技术、拓宽自己的应用范围寻找机会收 购、并购小的、专业的软件商，因此CAE软件本身的功能得到了极大的提升。MSC公司作为最早成

18、立的CAE公司，先后通过开发、并购，已经把数个CAE程序集成到 其分析体系中。目前MSC公司旗下拥有10几个产品，如Nastran、patran、Marc、Adams、 Dytran和Easy 5等，覆盖了线性分析、非线性分析、显式非线性分析以及流体动力学问题 和流场耦合问题。另外，MSC公司还推出了多学科方案(MD)来把以上的诸多产品集成为了 一个单一的框架解决多学科仿真问题。ANSYS公司通过一连串的并购与自身壮大后，把其产品扩展为ANSYS Mechanical系列， ANSYS CFD (FLUENT/CFX)系列，ANSYS ANSOFT 系列以及 ANSYS Workbench和

19、EKM 等。由此 ANSYS塑造了一个体系规模庞大、产品线极为丰富的仿真平台，在结构分析、电磁场分析、 流体动力学分析、多物理场、协同技术等方面都提供完善的解决方案。SDRC把其有限元程序Supertab并入到I -DEAS中，并加入耐用性、NVH、优化与灵敏度、 电子系统冷却、热分析等技术，且将有限元技术与实验技术有机地结合起来，开发了实验信 号处理、实验与分析相关等分析能力。而在2001年SDRC公司被EDS所收购，并将其与UGS 合并重组,SDRC的有限元分析程序也演变成了 NX中的I-deas NX Simulation，与NX Nastran 一起成为了 NX产品生命周期中的仿真分析

20、中的重要组成部分。说到NX Nastran，大家都会想到另一个以 Nastran为名的有限元软件 MSC.Nastran。 MSC.Nastran与NX Nastran可谓是同根同源，皆是由NASA推出的Nastran程序的源代码发 展出来的。下面我们可以简单的介绍下MSC.Nastran与NX Nastran的由来。在当时开发Nastran程序的不止MSC 一家公司，还有另外两家公司也推出了 Nastran程序的 商业版，1972年UAI公司发布基于COSMIC NASTRAN的UAI Nastran软件，1985年CSAR公 司发布了基于COSMIC NASTRAN的CSAR Nastra

21、n软件。当时是市场上有这3家公司共同经营 Nastran 软件。而在1999年，MSC收购了 UAI和CSAR，成为市场上惟家提供Nastran商业代码的供应 商。而在此后的几年，独自享有源代码的MSC Nastran软件价格不断上涨，但是其功能和服 务却没有得到相应的提升，从而引发大量客户的抱怨，为此NASA则向美国联邦贸易委员会 （FTC）提出了申诉。美国FTC判MSC Nastran垄断”，MSC Nastran源代码须公开，而这一决定也引来了 UGS 公司加入到Nastran的市场中来。而后，UGS根据MSC所提供的源代码、测试案例、开发工 具和其他技术资源开发出了 NX Nastra

22、n。至此，源于NASA的Nastran 一分为二，齐头并进， 为用户带来了更多的新技术与服务。进入21世纪后，早期的三大软件商MSC、ANSYS、SDRC的命运各不相同，SDRC被EDS收购 后与UGS进行了重组，其产品I-DEAS已经逐渐的淡出了人们的视线；MSC自从Nastran被 反垄断拆分后一蹶不振，2009年7月被风投公司STG收购，前途至今还不明朗；而ANSYS 则是最早出现的三大巨头中最为强劲的一支，收购了 Fluent、CFX、Ansoft等众多知名厂商 后，逐渐的塑造了一个体系规模庞大、产品线极为丰富的仿真平台。而在CAE市场的其他厂商也发生了不少的并购和重组，一些新近的厂商

23、也逐渐在崭露头角。 如并入达索SIMULIA的ABAQUS，能否如SolidWorks 一样借助达索的强劲在CAE市场中打出 一片天地；以前后处理而进入CAE领域的Altair公司，其Hypermesh软件自诞生之日起就 备受业界的关注，而围绕前后处理建立起来的HyperWorks软件，也已经成为了现在市场上 很有竞争力的软件，近几年来收入也持续上涨；LMS也是一个比较有特点的CAE软件公司， 其软件的分析集1D、3D、“试验”于一身，不仅可以加速虚拟仿真，还能使仿真结果更准 确可靠；COMSOL则是以多物理场耦合仿真开辟出了一片新天地，为其发展、更为CAE技术的发展拨开迷雾。另外，在市场中占

24、有一定份额的还有如前后处理软件ANSA、Truegrid，流体仿真软件Fluent （被 ANSYS 收购）、CFX （被 ANSYS 收购）、Phoenics、NUMECA、Star-CD，铸造仿真软件ProCAST、FLOW-3D、MAGMA SOFT 等一批专业 CAE 分析软件。3国内有限元法的发展之路我国的力学工作者为有限元方法的初期发展做出了许多贡献，其中比较著名的有:陈伯屏（结 构矩阵方法），钱令希（余能原理），钱伟长（广义变分原理），胡海昌（广义变分原理）， 冯康（有限单元法理论）。遗憾的是由于当时环境所致，我国有限元方法的研究工作受到阻 碍，有限元理论的发展也逐渐与国外拉开了

25、距离。20世纪60年代初期，我国的老一辈计算科学家较早地将计算机应用于土木、建筑和机 械工程领域。当时黄玉珊教授就提出了 “小展弦比机翼薄壁结构的直接设计法”和“力法一 应力设计法”；而在70年代初期，钱令希教授提出了 “结构力学中的最优化设计理论与方 法的近代发展”。这些理论和方法都为国内的有限元技术指明了方向。1964年初崔俊芝院士研制出国内第一个平面问题通用有限元程序，解决了刘家峡大坝 的复杂应力分析问题。20世纪60年代到70年代，国内的有限元方法及有限元软件诞生之 后，曾计算过数十个大型工程，应用于水利、电力、机械、航空、建筑等多个领域。20世纪70年代中期，大连理工大学研制出了 J

26、EFIX有限元软件，航空工业部研制了 HAJIF系列程序。80年代中期，北京大学的袁明武教授通过对国外SAP软件的移植和重大改 造，研制出了 SAP-84；北京农业大学的李明瑞教授研发了 FEM软件；建筑科学研究院在国 家“六五”攻关项目支持下，研制完成了“BDP-建筑工程设计软件包”；中国科学院开发了 FEPS、SEFEM；航空工业总公司飞机结构多约束优化设计系统YIDOYU等一批自主程序。 上世纪90年代以来，大批国外CAE软件涌入国内市场，遍及国内的各个领域，国外的专家 则深入到大学、院所、企业与工厂，展示他们的CAE技术、系统功能及使用技巧，因此使得 国内自主研发CAE软件受到强烈打压

27、。同时，有关管理部门在对直接为先进装备制造业服务 的CAE软件核心技术的认识上产生了偏差：CAE既不属于基础科学，又不属于科技攻关，故 而失去了必要的支持，使其发展举步维艰，以至于在上世纪的最后十几年国内CAE自主创新 的步伐已经非常缓慢，也逐渐的拉开了与国外CAE软件的距离。进入21世纪后，虽然国外CAE软件占据市场主流的现状短时间内已经无法撼动，但国内自 主知识产权CAE软件逐渐市场化，获得了一定的发展：北京飞箭软件有限公司推出的FEPG、 郑州机械研究所推出的紫瑞CAE、大连的大工安道公司的CAE软件Adopt.Smart；湖南大学 与吉林大学开发了针对汽车结构的KMAS分析系统；华中科

28、技大学针对铸造成型开发的华铸 CAE软件；清华大学、上海交大在注塑成型CAE领域也推出了相应的分析软件。虽然国内CAE自主研发之路历经艰辛，但是广大专家学者用锲而不舍的战斗精神顽强地生存 下来。尤其是在近几年，数字化产品设计的概念逐渐深入人心，国内高校技术研究和应用水 平不断提高，有限元技术已经为广大企业所认可。随着国家对发展自主CAE平台已经愈发重 视，国内CAE的研究已经逐渐走出低迷状态，获得了一定的发展，而且值得注意的是，有限 元技术不再仅仅停留在高校中，而是更多的走向了企业。同时，更多使用方便、操作简单的 专用分析软件也得到了广泛应用。4多种技术、同一目的发展到今天，CAE软件不仅仅只

29、有有限元法一种基本算法，目前已经发展出了包括有限差分 法、有限元法、有限体积法等多种数学算法。而究其本质都是采用微积分的方法对离散方程 进行求解，从而得出所求的结果。三种方法的比较如表2所示。多种求解方法使CAE技术得到的长足的发展，而目前有限元法覆盖的领域已经非常广泛，并 已大量的应用于结构力学，结构动力学，热力学，流体力学等仿真分析，并且向着多物理场 耦合分析的方向发展着，进而逐步的推动着CAE技术的发展。5结语CAE技术诞生至今已半个多世纪，随着计算机技术与软件技术的不断发展壮大，CAE的 技术手段与应用范围已经不可同日而语，随之而来也诞生了数以百计的CAE软件商。如今的 CAE市场是一

30、个群雄割据的年代，据不完全统计全球有超过200种仿真分析的软件在被企业 所使用着，昔日的巨头命运也各不相同，更有着无数的新贵崛起进入人们的视野。种类繁多的CAE软件为人类探索未知提供了工具，而CAE软件已经逐渐的被更多的工程人员 所接受，在CAE软件中进行着产品的仿真、分析与优化，同时工程和产品的自主创新也为 CAE发展提供了强大的动力，相信CAE技术将会随着工程科学技术的飞速发展而迎来一个更 为辉煌灿烂的明天分享PATRAN的一些小技巧！ ！、在Patran里如何Move 一组Points的位置，而不改变这组Points的ID编号？ Group/Transform/Translate 的功能，这样不 但编号不会变，连property跟边界条件都会保留2、Patran如