非线性有限元.docVIP

非线性有限元 　　在20世纪60年代，由于Ed Wilson发布了他的第一个程序，这种激情终于被点燃了。这些程序的第一代没有名字。在遍布世界的许多实验室里，通过改进和扩展这些早期在Berkeley开发的软件，工程师们扩展了新的用途，这些带来了对工程分析的巨大冲击和有限元软件的随之发展。在Berkeley开发的第二代线性程序称之为SAP（structural analysis program）。由Berkeley的工作发展起来的第一个非线性程序使NONSAP，它具有隐式积分进行平衡求解和瞬时问题求解的功能。 　　第一批非线性有限元方法文章的主要贡献者由Argyris（1965），Marcal和King（1967）。不久，大批文章激增，而且软件随之诞生。当时在Brown大学任教的Pedro Marcal，作为第一个非线性商业有限元程序进入市场，与1969年建立了一个公司，程序命名为MARC，目前它仍然是主要软件。大约在同期，John Swanson 为了核能应用在Westinghouse发展了一个非线性有限元程序。为了使ANSYS程序进入市场，他于1969年离开Westinghouse。尽管ANSYS主要是关注非线性材料而非求解完全的非线性问题，但他多年来仍垄断了商业非线性有限元的舞台 　　在早期的商用软件舞台上，另外两个主要人物使David Hibbitt 和Klaus－Jurgen Bathe。Hibbitt与Pedro Marcal 合作到了1972年，后来与其它人合作建立了HKS公司，使ABAQUS商业软件进入市场。因为该程序是能够引导研究人员增加用户单元和材料模型的早期有限元程序之一，所以它对软件行业带来了实质性的冲击。Jurgen Bathe 是在Ed Wilson的指导下在Berkeley获得博士学位的，不久之后开始在MIT任教，这期间他发表了他的程序。这是NONSAP软件的派生品，称之ADINA。 直到大约1990年，商用有限元程序集中在静态解答和隐式方法的动态解答。在20世纪70年代，这些方法取得了非常大的进步，主要贡献来自于Berkeley，起源于Berkeley的研究人员：Thomas J.R. Hughes，Robert Tayor,Juan Simo,Jurgen Bathe,Carlos Felippa,Pal Bergan,Kaspar Willam,Ekerhard Ramm和Michael Ortiz。他们是Berkeley的杰出研究者中的一部分。不容置疑，他们是早期有限元的主要孵化人员。 当代非线性软件的另一个血脉是显式有限元程序。Wilkins（1964）在DOE实验室的工作强烈的影响了早期的显式有限元方法，特别式命名为hydro－codes的软件。 1964年，Costantino在芝加哥的IIT研究院发展了可能是第一个的显式有限元程序（Costantino，1967）。它局限于线性材料和小变形由带状刚度矩阵乘以节点位移计算内部的节点力。它首先在一台IBM7040系列计算机上运行，花费了数百万美元，其数度远远低于一个megaflop（每秒一百万浮点运算）和3200字节RAM。刚度矩阵存储在磁带上，通过观察磁带驱动能够检测计算的过程。当每一步骤完成时，磁带驱动将逆转以便允许阅读刚度矩阵。这些和以后的Control Data机器有类似的性能，如CDC6400和6600，他们是20世纪60年代运行有限元程序的机器。一台CDC6400价值近1000万美元，有32kB内存（存储全部的操作系统和编译器）和大约一个megaflop的真实速度。 1969年，为了实现对空军销售的计划，高级研究人员开发了著名的单元乘单元的技术，节点力的计算不必运用刚度矩阵。因此，发展了名为SAMSON的二维有限元程序，它被美国的武器实验室用了10年，1972年，该程序的功能扩展至结构的完全非线性三维瞬态分析，称之WERCKER。这一工作得到美国运输部敢于幻想的计划经理Lee Ovenshire的基金资助，它在20世纪70年代初期就预言汽车的碰撞实验可能被仿真所代替。 然而，比他预言的时间稍微提前了一点。在当时运行一个300个单元模型的仿真，对于2000万次模拟需要大约30小时的计算机机时，花费大约3万美元，这个相当于助理教授3年的工资。Lee Ovenshire 的计划资助了若干个开拓性的工作：Hughes的接触－冲击工作，Ivor Mclvor 的碰撞工作，以及由Ted Shugar和Carly Ward在Port Hueneme所从事的关于人头颅的模拟研究。但是，大约在1975年，运输部认为仿真太昂贵，决定所有的基金转向实验方面，使这些研究努力令人心痛地停止下来。在Ford公司WRECKER勉强生存了下一个10年。而在Ar