Mechanics 版 (精华区)

发信人: rhc (土藤), 信区: Mechanics
标  题: 有限元分析软件的比较及展望
发信站: BBS 哈工大紫丁香站 (Sat May 21 17:32:28 2005)

有限元分析软件的比较及展望 

大型有限元分析软件ＡＮＳＹＳ的特点


    随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑
物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师
在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温
度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥
梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确
定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。
这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的
。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Elem
ent Analysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实
践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下
几个方面： 
   增加设计功能，减少设计成本； 
   缩短设计和分析的循环周期； 
   增加产品和工程的可靠性； 
   采用优化设计，降低材料的消耗或成本； 
   在产品制造或工程施工前预先发现潜在的问题； 
   模拟各种试验方案，减少试验时间和经费； 
   进行机械事故分析，查找事故原因。 
   在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分
析计算，FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、
60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的
是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的
NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能
最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较
小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PAFE
C、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC
和STARDYNE等公司的产品。 

   以下对一些常用的软件进行一些比较分析： 
1. LSTC公司的LS-DYNA系列软件 
LS-DYNA是一个通用显式非线性动力分析有限元程序，最初是1976年在美国劳伦斯利弗莫尔
国家实验室(Lawrence Livermore National Lab.)由J.O.Hallquist 主持开发完成的，主
要目的是为核武器的弹头设计提供分析工具，后经多次扩充和改进，计算功能更为强大。
此软件受到美国能源部的大力资助以及世界十余家著名数值模拟软件公司（如ANSYS、MSC
.software、ETA等）的加盟，极大地加强了其的前后处理能力和通用性，在全世界范围内
得到了广泛的使用。在软件的广告中声称可以求解各种三维非线性结构的高速碰撞、爆炸
和金属成型等接触非线性、冲击载荷非线性和材料非线性问题。即使是这样一个被人们所
称道的数值模拟软件，实际上仍在诸多不足，特别是在爆炸冲击方面，功能相对较弱，其
欧拉混合单元中目前最多只能容许三种物质，边界处理很粗糙，在拉格朗日——欧拉结合
方面不如DYTRAN灵活。虽然提供了十余种岩土介质模型，但每种模型都有不足，缺少基本
材料数据和依据，让用户难于选择和使用。 

2. MSC.software公司的 DYTRAN软件 
当前另一个可以计算侵彻与爆炸的商业通用软件是MSC.Software Corporation ( MSC公司
) 的MSC.DYTRAN程序。该程序在是在LS-DYNA3D的框架下，在程序中增加荷兰PISCES INTE
RNATIONAL公司开发的PICSES的高级流体动力学和流体——结构相互作用功能，还在PISCE
S的欧拉模式算法基础上，开发了物质流动算法和流固耦合算法。在同类软件中，其高度非
线性、流—固耦合方面有独特之处。 
MSC.DYTRAN的算法基本上可以概况为：MSC.DYTRAN采用基于Lagrange格式的有限单元方法
（FEM）模拟结构的变形和应力，用基于纯Euler格式的有限体积方法（FVM）描述材料（包
括气体和液体）流动，对通过流体与固体界面传递相互作用的流体—结构耦合分析，采用
基于混合的Lagrange格式和纯Euler格式的有限单元与有限体积技术，完成全耦合的流体-
结构相互作用模拟。MSC.DYTRAN用有限体积法跟踪物质的流动的流体功能，有效解决了大
变形和极度大变形问题，如：爆炸分析、高速侵彻。 
但MSC.DYTRAN本身是一个混合物，在继承了LS-DYNA3D与PISCES的优点同时，也继承了其不
足。首先，材料模型不丰富，对于岩土类处理尤其差，虽然提供了用户材料模型接口，但
由于程序本身的缺陷，难于将反映材料特性的模型加上去；其次，没有二维计算功能，轴
对称问题也只能按三维问题处理，使计算量大幅度增加；在处理冲击问题的接触算法上远
不如当前版的LS-DYNA3D全面。 

3. HKS公司的ABAQUS软件 
ABAQUS是一套先进的通用有限元系统，也是功能最强的有限元软件之一，可以分析复杂的
固体力学和结构力学系统。ABAQUS有两个主要分析模块：ABAQUS/Standard提供了通用的分
析能力，如应力和变形、热交换、质量传递等；ABAQUS/Explicit应用对时间进行显示积分
求解，为处理复杂接触问题提供了有力的工具，适合于分析短暂、瞬时的动态事件，但对
爆炸与冲击过程的模拟相对不如DYTRAN和LS-DYNA3D 

4 ADINA 
ADINA是一个古老的有限元软件, 有一些很老的版本，它们只有基本的计算功能，没有前后
处理。用它算题，必须自己手工建模，现在看来这些实在是太落后了，但是，重要的一点
是它有源代码。有了源码，就可以对程序进行改造，满足特殊的需求。其实国内对ADINA的
改造还是很多的，比如将等带宽存储改为变带宽存储，将元素库从整个程序中分离出来，
可以有选择的将将元素编译连接到程序中。还有的在程序中加入了自己的材料本构关系，
也有在元素库中加进了新的单元等等。经过这些改进，程序的功能得到了扩展，效率得到
了提高，更重要得是在一定程度上具有了自己的知识产权。 

5 ANSYS和NASTRAN 
因为和NASA的特殊关系，msc nastran在航空航天领域有着崇高的地位。而ANSYS则在铁道
，建筑和压力容器方面应用较多。尽管目前, ANSYS已发展了很多版本, 其实它们核心的计
算部分变化不大，只是模块越来越多。比如5.1没有lsdyna，和cad软件的接口，到了5.6还
有疲劳模块等等。其实这些模块并不是ANSYS公司自己搞的，就是把别人的东西买来集成到
自己的环境里。NASTRAN最早是用的for windows 2.0。是nsatran v68集成在femap5里。n
astran的求解器效率比ansys高一些。有一个算例可以说明，20000多个节点，D版的ansys
56建模，用femap7.0转成nastran的dat文件，静力计算及前5阶的线性频率，结果ansys56
在PIII450上所用的时间和D版的nastran707在赛杨400上用的时间相当，内存都是128M，全
部选项都是缺省的，nastran用子空间迭代法求频率，ansys没仔细看，计算的结果倒是没
什么大的差别。 

其他还有一些软件例如sap，algor,cosmos等，只是影响比较小。 

还有一点值得说明, 目前的有限元软件，求出的位移结果都很准，可应力就不太一样了，
这是一个有趣的现象, 大家可以讨论。 

  另外，从发展上来说，国际上数值模拟软件发展呈现出以下一些趋势特征 
a. 由二维扩展为三维 
   早期计算机的能力十分有限，受计算费用和计算机储存能力的限制，数值模拟程序大多
是一维或二维的，只能计算垂直碰撞或球形爆炸等特定问题。随着第三代、第四代计算机
的出现，才开始研制和发展更多的三维计算程序。现在，计算程序一般都由二维扩展到了
三维，如LSDYNA2D和LSDYNA3D，AUTODYN2D和AUTODYNA3D，但也有完全在三维基础上开发的
，如MSC.DYTRAN，就没有二维功能。 
b. 从单纯的结构力学计算发展到求解许多物理场问题 
数值模拟分析方法最早是从结构化矩阵分析发展而来，逐步推广到板、壳和实体等连续体
固体力学分析，实践证明这是一种非常有效的数值分析方法。近年来数值模拟方法已发展
到流体力学、温度场、电传导、磁场、渗流和声场等问题的求解计算，最近又发展到求解
几个交叉学科的问题。例如内爆炸时，空气冲击波使墙、板、柱产生变形，而墙、板、柱
的变形又反过来影响到空气冲击波的传播……这就需要用固体力学和流体动力学的数值分
析结果交叉迭代求解，即所谓“流—固耦合”的问题。 
c. 从单一坐标体系发展多种坐标体系 
   数值模拟软件在开始阶段一般采用单一坐标，或采用拉格朗日坐标或采用欧拉坐标，由
于这两种坐标自身的缺陷，计算分析问题的范围都有很大的限制。为克服这种缺陷，采用
了三种方法，一是两个程序简单组合，如CTH—EPIC，爆炸与侵彻由不同的程序分开计算；
二是在同一程序中采用多种坐标体系，如DYNA3D中早期采用的是拉格朗日坐标，而LSDYNA
3D的最新版除原有类型外，新加了欧拉方法以及拉格朗日与欧拉耦合方法，而最近几年才
发展的DYTRAN则是拉格朗日型的LSDYNA3D(1988版)与欧拉型的PISCES的整合体；三是采用
新的计算方法，如SPH等，SPH法不用网格，没有网格畸变问题，所以能在拉格朗日格式下
处理大变形问题，同时，SPH法允许存在材料界面，可以简单而精确地实现复杂的本构行为
，也适用于材料在高加载速率下的断裂等问题的研究。 
d. 由求解线性工程问题进展到分析非线性问题 
随着科学技术的发展，线性理论已经远远不能满足设计的要求。诸如岩石、土壤、混凝土
等，仅靠线性计算理论就不足以解决遇到的问题，只有采用非线性数值算法才能解决。众
所周知，非线性的数值计算是很复杂的，它涉及到很多专门的数学问题和运算技巧，很难
为一般工程技术人员所掌握。为此，近年来国外一些公司花费了大量的人力和投资，开发
了诸如LSDYNA3D、ABAQUS和AUTODYN等专长于求解非线性问题的有限元分析软件，并广泛应
用于工程实践。这些软件的共同特点是具有高效的非线性求解器以及丰富和实用的非线性
材料库。 
e. 增强可视化的前置建模和后置数据处理功能 
   早期数值模拟计算软件的研究重点在于推导新的高效率求解方法和高精度的单元。随着
数值分析方法的逐步完善，尤其是计算机运算速度的飞速发展，整个计算系统用于求解运
算的时间越来越少，而数据准备和运算结果的表现问题却日益突出。在现在的工程工作站
上，求解一个包含10万个方程的有限元模型只需要用几十分钟。但如果用手工方式来建立
这个模型，然后再处理大量的计算结果则需用几周的时间。可以毫不夸张地说，工程师在
分析计算一个工程问题时有80%以上的精力都花在数据准备和结果分析上[14]。因此目前几
乎所有的商业化数值模拟程序系统都有功能很强的前置建模和后置数据处理模块。在强调
“可视化”的今天，很多程序都建立了对用户非常友好的GUI（图形用户界面—Graphics 
User Interface），使用户能以可视图形方式直观快速地进行网格自动划分，生成有限元
分析所需数据，并按要求将大量的计算结果整理成变形图、等值分布图，便于极值搜索和
所需数据的列表输出。 
f. 与CAD软件的无缝集成 
   与通用CAD软件的集成使用，即在用CAD软件完成结构设计后，自动生成有限元网格并进
行计算，如果分析的结果不符合设计要求则重新进行构造和计算，直到满意为止，从而极
大地提高了设计水平和效率。今天，工程师可以在集成的CAD和数值模拟软件环境中快捷地
解决一个在以前无法应付的复杂工程分析问题。所以当今所有的商业化有限元系统商都开
发了和著名的CAD软件（例如AutoCAD、Pro/ENGINEER、Unigraphics、SolidEdge、SolidW
orks、IDEAS等）的接口。 
g. 工作平台多样化 
   早期的数值分析软件基本上都是在大中型计算机上开发和运行的，后来又发展到以工程
工作站（EWS，Engineering Work Station）上，它们的共同特点都是采用UNIX操作系统。
PC机的出现使计算机的应用发生了根本性的变化，工程师渴望在办公桌上完成复杂工程分
析的梦想成为现实。但是早期的PC机采用16位CPU和DOS操作系统，内存中的公共数据块受
到限制，因此当时计算模型的规模不能超过1万阶方程。Microsoft Windows操作系统和32
位的Intel Pentium处理器的推出，为PC机用于有限元分析提供了必需的软件和硬件支撑平
台。因此当前国际上著名的有限元程序研究和发展机构都纷纷将他们的软件移值到Window
s平台上。最新高档PC机的求解能力已和中低档的EWS不相上下。 
为了将在大中型计算机和EWS上开发的有限元程序移值到PC机上，常常需要采用Hummingbi
rd公司的一个仿真软件Exceed。这样做的结果比较麻烦，而且不能充分利用PC机的软硬件
资源。所以最近有些公司，例如ANSYS、MSC.software等开始在Windows平台上开发有限元
程序，大多采用了OpenGL图形编程软件，同时还有在PC机上的Linux操作系统环境中开发的
有限元程序包。 
h. 软件开发强强联合 
   由于数值软件的开发是一项长期而艰巨的任务，开发一个通用软件是十分困难的，各家
开发的软件由于应用背景的不同而各有千秋，随着数值模拟软件商业化的进展，各数值模
拟软件公司为扩大市场，追求共同的利润，出现了强强联合的局面。典型的如ANSYS与LSD
YNA3D联合，MSC.software软件公司对ABAQUS、LS DYNA3D及PISCES等的购买。 
　　 

   再谈一下国内的发展情况和前景 
   1979年美国的SAP5线性结构静、动力分析程序向国内引进移植成功，掀起了应用通用有
限元程序来分析计算工程问题的高潮。这个高潮一直持续到1981年ADINA非线性结构分析程
序引进，一时间许多一直无法解决的工程难题都迎刃而解了。大家也都开始认识到有限元
分析程序的确是工程师应用计算机进行分析计算的重要工具。但是当时限于国内大中型计
算机很少，大约只有杭州汽轮机厂的Siemens7738和沈阳鼓风机厂的IBM4310安装有上述程
序，所以用户算题非常不方便，而且费用昂贵。PC机的出现及其性能奇迹般的提高，为移
植和发展PC版本的有限元程序提供了必要的运行平台。可以说国内FEA软件的发展一直是围
绕着PC平台做文章。在国内开发比较成功并拥有较多用户（100家以上）的有限元分析系统
有大连理工大学工程力学系的FIFEX95、北京大学力学与科学工程系的SAP84、中国农机科
学研究院的MAS5.0和杭州自动化技术研究院的MFEP4. 等。但正如上面所述，国外很多著名
的有限元分析公司已经从前些年对PC平台不屑一顾转变为热衷发展，对国内FEA程序开发者
来说发展PC版本不再具有优势，而以后应该从下面几方面加以努力： 
1. 研究开发求解非固体力学和交叉学科的FEA程序经过几十年的研究和发展，用于求解固
体力学的有限元方法和软件已经比较成熟，现在研究的前沿问题是流体动力学、可压缩和
不可压缩流体的流动等非固体力学和交叉学科的问题。由于国内没有类似功能的商品化软
件，所以国外的软件就卖得非常贵。为了破这种垄断局面，我们必须发展有自主版权、用
于分析流体等非固体力学和交叉学科的软件。因为流体力学问题远比固体复杂得多，而且
很少有现成的软件可以借鉴，所以需要投入大量的人力和经费。这就必须有国家和大型企
业集团来支持。 
2. 开发具有中国特色的自动建模技术和GUI开发建模技术和GUI的投入比前述课题要少得多
，但却可以大大提高FEA软件的性能和用户接受程度，从而起到事半功倍的效果。国内不少
人在这方面做了很多工作，但是由于当时PC机上的图形支撑环境有限，所以开发的效果都
不甚理想。Windows中提供了OpenGL图形标准，为在PC机上应用可视化图形技术开发GUI提
供了强有力的工具。OpenGL是当今国际上公认的高性能图形和交互式视景处理标准，应用
它开发出来的三维图形软件深受专业技术人员的钟爱，目前世界上占主导地位的计算机公
司都采用了这一标准。正如前面所述，近年来国外有的FEA程序已抛开仿真软件，直接在W
indows平台上开发有限元程序。杭州自动化技术研究院1997-1999年采用OpenGL图形标准和
相应的Visual C++等编程工具，在PC机上成功地开发了一套可视化有限元程序包。它能直
观地通过对"菜单"、"窗口"、"对话框"和"图标"等可视图形画面和符号的操作，自动建立
有限元分析模型，并以交互方法式实现计算结果的可视化处理，因而可大大提高有限昂分
析的效率和精确性，也便于用户学习和掌握。 
3. 与具有我国自主版权的CAD软件集成前面已经讲过，当今有限元方法的一个重要特点是
和CAD软件的无缝集成。作为我国自行开发的FEA程序，首先要考虑和我国自主版权的CAD软
件集成。因为有限元分析主要用于形状比较复杂的零部件，所以要和具有三维造型功能和
CAD软件集成，使设计和分析紧密结合、融为一体。 

说明: 来自于网上,但我作了一定概括总结, 可以算原创吧 
wgao 编辑于 2004-03-31 18:30



大型有限元分析软件ＡＮＳＹＳ的特点
________________________________________
王友海　颜慧军　胡长胜

　　ＡＮＳＹＳ程序是美国ＡＮＳＹＳ公司研制的大型有限元分析 (ＦＥＡ)软件 ,自 1 
970年ＪｏｈｎＳｗａｎｓｏｎ博士洞察到计算机模拟工程应该商品化 ,创建了ＡＮ ＳＹ
Ｓ公司以来 ,ＡＮＳＹＳ程序已发展成为全球范围一个多用途的设计分析软件。
ＡＮＳＹＳ程序是一个功能强大的设计分析及优化软件包。与其它有限元分析软件如ＳＡ
Ｐ或ＮＡＳ ＴＲＡＮ等相比 ,它有以下特点 :(1 )ＡＮＳＹＳ是完全的ＷＩＮＤＯＷＳ程
序 ,从而使应用更加方便 ;(2 )产品系列由一整套可扩展的、灵活集成的各模块组成 ,因
而能满足各行各业的工程需要 ;(3)它不仅可以进行线性分析 ,还可以进行各类非线性分析
 ;(4)它是一个综合的多物理场耦合分析软件 ,用户不但可用其进行诸如结构、热、流体流
动、电磁等的单独研究 ,还可以进行这些分析的相互影响研究 ,例如 :热—结构耦合 ,磁
—结构耦合以及电—磁—流体—热耦合等。本文将以ＡＮＳＹＳ/Ｓｔｒｕｃｔｕｒａｌ(
结构 )模块为例 ,详细研究该软件的功能及特点。
1　结构静、动力分析
1 1　结构静力分析
ＡＮＳＹＳ程序中结构静力分析 ,用来求解外载荷引起的位移、应力、和力。静力分析适
合于求解
惯性及阻尼的时间相关作用对结构响应的影响不显著的问题。ＡＮＳＹＳ程序中静力分析
同样能包括非线性 ,如塑性、蠕变、膨胀、大变形、大应变及接触面等。有关非线性内容
后面将详细叙述。
1 2　结构动力分析
结构动力分析用来求解随时间变化的载荷对结构或部件的影响。ＡＮＳＹＳ程序可以求解
下列类型
的动力分析问题 :瞬态动力、模态、谐波响应及随机振动。
1 2 1　瞬态分析
瞬态分析 (也称时间—历程分析 )用于确定结构承受随时间变化载荷时的动力响应。ＡＮ
ＳＹＳ求
解瞬态动力问题有三种方法 :全瞬态动力分析方法 ,凝聚法和模态叠加法。
1 2 2　模态分析

 
图 1　皮带轮模态分析

(虚线表示未变形形态 )
当需要结构的自然频率时 ,模态分析是很有用的 (图 1 )。模态分析用于抽取结构的自然
频率和振型 ,结构的基本模态和频率信息能有助于归结动力响应特征 ,该分析结果也有助
于确定用于后继瞬态动力分析的模态数目和积分时间步长。ＡＮＳＹＳ程序允许作预应力
模态分析及在大变形分析后作模态分析。
1 2 3　谐波响应分析
谐波响应分析用于求解非线性结构承受正弦变化载荷的响应。该分析用于研究随时间谐波
变化的
载荷的影响 ,如房屋、支座、旋转机构底座。
1 2 4　响应谱分析
响应谱分析可用于求解冲击载荷条件下的结构响应 ,该分析类型使用模态分析的结果 ,再
加上已
知的响应谱数据 ,计算每一个固有频率点在结构中发生的真实位移和应力 ,典型的响应谱
分析应用是地震分析 ,它用于研究地震对管道系统、塔和桥梁等结构的影响。
1 2 5　随机振动分析
随机振动分析是一种谱分析 ,用于研究结构对随机激励的响应 ,如喷气发动机或火箭引擎
产生的那些激励。

2　ＡＮＳＹＳ/ＬＳ ＤＹＮＡ显示求解
2 1　ＡＮＳＹＳ/ＬＳ ＤＹＮＡ发展情况
ＡＮＳＹＳ/ＬＳ ＤＹＮＡ是一个通用显示非线性动力
分析有限元程序 ,可以求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等接
触非线性、冲击载荷非线性和材料非线性问题。
2 2　应用领域
ＡＮＳＹＳ/ＬＳ ＤＹＮＡ的应用可分为国防和民用两大类 ,主要有 :汽车、飞机、火车
、轮船等运输工具的碰撞分析 ;金属成型、金属切割 ;汽车零部件的机械制造 ;塑料成型
、玻璃成型 ;生物力学 ;地震工程 ;消费品、建筑物、乘员、高速结构等的安全性分析 ;
点焊、铆焊、螺栓连接 ;液体—结构相互作用 ;运输容器设计 ;爆破工程的设计分析 ;战
斗部结构设计分析 :内弹道发射对结构的动力响应分析 ;终点弹道的爆炸驱动和破坏效应
分析 ;侵彻过程与爆炸成坑模态分析 ;军用设备和结构设施受碰撞和爆炸冲击加载的结构
动力分析 ;介质 (包括空气、水和地质材料等 )中爆炸及对舰船和结构作用的全过程模拟
分析 ;军用新材料 (包括炸药、复合材料、特种金属等 )的研制和动力特性分析 ;超高速
碰撞模拟分析 ;战地上有生力量的毁伤效应分析等等。
2 3　分析能力
ＡＮＳＹＳ/ＬＳ ＤＹＮＡ具有广泛的分析功能 ,可以模拟许多二维、三维结构的物理特
性 ,诸如 :非线性动力分析、热分析、失效分析、裂纹扩展分析、接触分析、二维静力分
析、任意拉格朗日—欧拉(ＡＬＥ)分析、流体—结构相互作用分析、实时场分析、多物理
场分析等等。
2 4　材料模型
ＡＮＳＹＳ/ＬＳ ＤＹＮＡ程序目前有 1 0 0余种金属和非金属材料模型可供选择 ,如弹
性、弹塑性、超弹性、泡末、玻璃、地质、土壤、混凝土、流体、复合材料、炸药及起爆
燃烧、刚性及用户自定义材料 ,并可考虑材料失效、损伤、粘性、蠕变、与温度相关、与
应变率相关等性质。
2 5　单元库
ＡＮＳＹＳ/ＬＳ ＤＹＮＡ程序的单元类型众多 ,有二维、三维单元 ,薄壳、厚壳、体、
梁单元 ,ＡＬＥ、Ｅｎｌｅｒｉａｎ、Ｌｎｇｒａｎｇｉａｎ单元等。各类单元又有多种
理论算法可供选择 ,单元积分采用沙漏粘性阻尼以克服零模式 ,单元计算速度快 ,节省存
储量并且精度都达到二阶 ,可以满足各种实体结构、薄壁结构和流体—固体耦合结构有限
元网格划分的需要。

3　结构屈曲分析
在ＡＮＳＹＳ中 ,稳定分析分为两类 :线性 (特征值 )屈曲分析和非线性屈曲分析。
3 1　线性屈曲分析
线性或特征值屈曲分析考虑了应力刚化效应 ,这种效应会导致结构在承受应力后抵抗横向
载荷能
力降低 ,当压应力增加时 ,结构抵抗横向能力减小。在某一载荷水平下 ,这种负的应力刚
度超过线
性结构刚度 ,造成结构屈曲。屈曲发生的点称为分叉点 (图 2 ),由于力—变形曲线达到该
点后可能
沿两条不同途径前进 ,所以当超过分叉点时 ,结构将屈曲或者在不稳定状态下承受载荷。

3 2　非线性屈曲分析
要精确地确定屈曲载荷 ,应该使用非线性屈曲分析。非线性屈曲分析中得出的极限载荷通
常比线性屈曲分析确定的分叉点低 (见图 2 )。这是由于非线性屈曲能考虑真实结构中存
在的初始缺陷 ,以及几何和材料的各种非线性。ＡＮＳＹＳ程序在大变形分析中是把弧长
法和Ｎｅｗｔｏｎ Ｒａｐｈｓｏｎ法结合起来修正结构单元的方位 ,从而求出屈曲临界载
荷。
当单独使用增量Ｎｅｗｔｏｎ Ｒａｐｈｓｏｎ法时 ,刚度矩阵可能会变成奇异矩阵。弧长
法使Ｎｅｗｔｏｎ Ｒａｐｈｓｏｎ平衡迭代沿一条弧收敛到其平衡路径。因此 ,可避免矩
阵在那些奇异点处变为奇异矩阵 ,并控制收敛
性 (图 3)。
 
图 2　分叉点或线性屈　　图 3　弧长法和Ｎｅｗｔｏｎ 

曲与极限载荷屈曲的比较　　　Ｒａｐｈｓｏｎ法结合使用

4　结构非线性结构
非线性导致结构或部件的响应随外载荷不成比例变化 ,实际上 ,所有结构本质上是非线性
的 ,只是在对分析影响很小时常常被忽略。ＡＮＳＹＳ程序可求解静态和瞬态非线性问题
。非线性静态分析将载荷分解成一系列增量的载荷步 ,并且在每一载荷步内进行一系列线
性逼近以达到平衡。每次线
性逼近需要对方程进行依次求解 (称为平衡迭代 )。类似地 ,非线性瞬态问题可被分解为
连续的随时间变化的载荷增量 ,在每步进行平衡迭代 ,然而瞬态情况也可能包括惯性效应
的时间积分。
在非线性瞬态分析中 ,动力平衡方程用Ｎｅｗ ｍａｒｋ时间积分求解。瞬态分析被分为离
散时间点。任意两个连续时间点之差称为积分时间步长。ＡＮ ＳＹＳ程序具有自动定义时
间步长的能力 ,它根据响应频率和非线性程度 ,增加或减小积分时间步长。在保证精度的
前提下 ,使得所需的时间步数量少。在静态和瞬态分析中 ,ＡＮＳＹＳ程序可考虑多种非
线性的影响 ,这些非线性可分为三类 :材料、几何和单元非线性。
4 1　材料非线性
当应力和应变不成比例时 ,存在材料非线性。ＡＮＳＹＳ程序可模拟各种非线性材料性质
 ,塑性、
多线性弹性和超弹性的特点是存在非线性的应力 应变关系。而粘塑性、蠕变和粘弹性的特
点是其应
变与其它因素 (如时间、温度和应力 )有关 ,非线性材料性质用Ｎｅｗｔｏｎ Ｒａｐｈｓ
ｏｎ法解决。
4 2　几何非线性
当结构位移显著地改变其刚度时 ,则被视为几何非线性。ＡＮＳＹＳ程序可解决这几类几
何非线性
效应 :大应变、大变形、应力刚化和旋转软化。大应变几何非线性解决大的局部变形问题
。它可作为结构变形而出现。材料中的应变和转角数量没有假定 ,程序通过调整反映几何
变化的单元形状来解决大应变问题。大变形表示由于结构变形导致的单元空间防卫变化而
引起的总结构刚度变化。这时假定发生小应变和大转动。程序利用更新单元方位作为结构
变形来解决大变形问题。应力刚化解决由应力状态引起的结构刚度的增大或减小。在该软
件中 ,旋转软化是另一种非线性效应 ,它往往是非常重要的。应力刚化解决内应力引起的
刚度变化 (通常增加 ),而旋转软化是针对由物体变形引起的刚度减小 ,诸如气轮机叶片在
平面的旋转。通常应力刚化和旋转软化共同用于旋转体分析。
4 3　单元非线性
非线性单元是其本身具有非线性行为的单元 ,而与其它单元无关。典型表现为由于状态变
化而引起刚度的突变 (诸如接触单元由开放变为关闭 )。
单元非线性提供了总体非线性不可能实现的各种功能。ＡＮＳＹＳ单元程序库包括下列非
线性单元 :一般的点面接触单元 ;刚性对柔性接触单元 ;界面单元 ;加筋实体单元 ;非线
性阻尼 (弹簧 );翘曲壳单元 ;组合单元 ;控制单元。
5　静力和运动学分析
运动学是力学的一个分支 ,它涉及运动而不考虑力或质量 ,可描述两种物体运动 :刚性或
柔性体。作为大变形和有限 (大的 )转动分析能力的一部分 ,ＡＮＳＹＳ程序可以分析大
型三维柔体运动 ,如挖掘机复杂的连杆结构。ＡＮＳＹＳ程序的下列特点使它能处理发生
发运动的结构 :Ｎｅｗｔｏｎ Ｒａｐｈｓｏｎ求解方法 :结构单元承受大转动的能力 ;代
表旋转铰的三维单元 ;代表线性调节器的三维单元。ＡＮ ＳＹＳ单元已设计得能使其在二
维或三维空间作无限制空间运动 ,该单元可用于构造曲柄连杆机构 ,诸如汽车挡风玻璃上
的雨刷 ,其中一部分结构 (曲柄 )的旋转运动导致另一部分结构 (清洗器 )的往复运动。
为表示整个模型的响应 ,必须使曲柄在分析中至少完成一次完整的旋转 ,具有无限旋转功
能的二维弹性梁可用于曲柄部分的建模 ,这就使得ＡＮＳＹＳ程序能追踪曲柄运动 ,并确
定整个模型运动。三维铰链单元用来表示连接模型两部分的铰链或销钉 ,这种单元表达多
种效应 ,诸如接头的柔性(或刚性 )、摩擦、阻尼和一定控制功能。
ＡＮＳＹＳ程序在许多领域都有其独到之处 ,这里不一一介绍。总之 ,该软件是目前国际
上最流
行、最实用、最全面的大型有限元分析程序 ,对它的深入研究必会对我国的经济建设起到
促进作用。

--

※ 来源:·哈工大紫丁香 http://bbs.hit.edu.cn·[FROM: 202.118.226.129]