【推荐】偶尔看到的一篇文章，觉得好推荐给大家

                          有限元分析系统的发展现状与展望
   
TYPE：理论指南 TIME：2003-4-21  
AUTHOR：尹柏生 URL：N/A  
HIT：4534  DAYHIT：28  WEEKHIT：4532 　  
   
随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析 新只 叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式往往是不可能的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element Analysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。我国在"九五"计划期间大力推广CAD技术，机械行业大中型企业CAD的普及率从"八五"末的20%提高到目前的70%。随着企业CAD应用的普及，工程技术人员已逐步甩掉图板，而将主要精力投身如何优化设计，提高工程和产品质量，计算机辅助工程分析（CAE，Computer Aided Engineering）方法和软件将成为关键的技术要素。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面：  
  
增加设计功能，减少设计成本；  
缩短设计和分析的循环周期；  
增加产品和工程的可靠性；  
采用优化设计，降低材料的消耗或成本；  
在产品制造或工程施工前预先发现潜在的问题；  
模拟各种试验方案，减少试验时间和经费；  
进行机械事故分析，查找事故原因。  
  
在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA在工程设计和分析中将得到越来越广泛的重视。下图是美国旧金山海湾大桥地震响应计算的有限元分析模型。发展方向及重大进展国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。当今国际上FEA方法和软件发展呈现出以下一些趋势特征：  
  
1. 从单纯的结构力学计算发展到求解许多物理场问题有限元分析方法最早是从结构化矩阵分析发展而来，逐步推广到板、壳和实体等连续体固体力学分析，实践证明这是一种非常有效的数值分析方法。而且从理论上也已经证明，只要用于离散求解对象的单元足够小，所得的解就可足够逼近于精确值。所以近年来有限元方法已发展到流体力学、温度场、电传导、磁场、渗流和声场等问题的求解计算，最近又发展到求解几个交叉学科的问题。例如当气流流过一个很高的铁塔产生变形，而塔的变形又反过来影响到气流的流动……这就需要用固体力学和流体动力学的有限元分析结果交叉迭代求解，即所谓"流固耦合"的问题。  
  
2. 由求解线性工程问题进展到分析非线性问题随着科学技术的发展，线性理论已经远远不能满足设计的要求。例如建筑行业中的高层建筑和大跨度悬索桥的出现，就要求考虑结构的大位移和大应变等几何非线性问题；航天和动力工程的高温部件存在热变形和热应力，也要考虑材料的非线性问题；诸如塑料、橡胶和复合材料等各种新材料的出现，仅靠线性计算理论就不足以解决遇到的问题，只有采用非线性有限元算法才能解决。众所周知，非线性的数值计算是很复杂的，它涉及到很多专门的数学问题和运算技巧,很难为一般工程技术人员所掌握。为此近年来国外一些公司花费了大量的人力和投资开发诸如MARC、ABQUS和ADINA等专长于求解非线性问题的有限元分析软件，并广泛应用于工程实践。这些软件的共同特点是具有高效的非线性求解器以及丰富和实用的非线性材料库。  
  
3. 增强可视化的前置建模和后置数据处理功能早期有限元分析软件的研究重点在于推导新的高效率求解方法和高精度的单元。随着数值分析方法的逐步完善，尤其是计算机运算速度的飞速发展，整个计算系统用于求解运算的时间越来越少，而数据准备和运算结果的表现问题却日益突出。在现在的工程工作站上，求解一个包含10万个方程的有限元模型只需要用几十分钟。但是如果用手工方式来建立这个模型，然后再处理大量的计算结果则需用几周的时间。可以毫不夸张地说，工程师在分析计算一个工程问题时有80%以上的精力都花在数据准备和结果分析上。因此目前几乎所有的商业化有限元程序系统都有功能很强的前置建模和后置数据处理模块。在强调"可视化"的今天，很多程序都建立了对用户非常友好的GUI（Graphics User Interface）,使用户能以可视图形方式直观快速地进行网格自动划分，生成有限元分析所需数据，并按要求将大量的计算结果整理成变形图、等值分布云图，便于极值搜索和所需数据的列表输出。  
  
4. 与CAD软件的无缝集成当今有限元分析系统的另一个特点是与通用CAD软件的集成使用 即,在用CAD软件完成部件和零件的造型设计后，自动生成有限元网格并进行计算，如果分析的结果不符合设计要求则重新进行造型和计算，直到满意为止，从而极大地提高了设计水平和效率。今天，工程师可以在集成的CAD和FEA软件环境中快捷地解决一个在以前无法应付的复杂工程分析问题。所以当今所有的商业化有限元系统商都开发了和著名的CAD软件（例如Pro/ENGINEER、Unigraphics、SolidEdge、SOLIDWORKS、IDEAS、Bentley和AutoCAD等）的接口。  
  
5. 在Wintel平台上的发展早期的有限元分析软件基本上都是在大中型计算机（主要是Mainframe）上开发和运行的，后来又发展到以工程工作站（EWS，Engineering Work Station）为平台，它们的共同特点都是采用UNIX操作系统。PC机的出现使计算机的应用发生了根本性的变化，工程师渴望在办公桌上完成复杂工程分析的梦想成为现实。但是早期的PC机采用16位CPU和DOS操作系统，内存中的公共数据块受到限制，因此当时计算模型的规模不能超过1万阶方程。Microsoft Windows操作系统和32位的Intel Pentium处理器的推出为将PC机用于有限元分析提供了必需的软件和硬件支撑平台。因此当前国际上著名的有限元程序研究和发展机构都纷纷将他们的软件移值到Wintel平台上。下表列出了用ADINA V7.3版在PC机的Windows NT环境和SGI工作站上同时计算4个工程实例所需要的求解时间。从中可以看出最新高档PC机的求解能力已和中低档的EWS不相上下。　  
　　　　　　　  
ADINA V7.3版几种机器上的运行时间比较  
计算问题的特征--方程数NEQ--运行时间（秒） PentiumII450MHZ PentiumIII450MHZ--工程工作站  
非线性三维接触  
（计算一步）--51159--219--188--219  
线性静力分析--58707--262--234--249  
动力分析频率计算--31205--311--258--287  
三维流动分析--36162--218--183--208  
  
为了将在大中型计算机和EWS上开发的有限元程序移值到PE机上，常常需要采用Hummingbird公司的一个仿真软件Exceed。这样做的结果比较麻烦，而且不能充分利用PC机的软硬件资源。所以最近有些公司，例如IDEAS、ADINA和R&D开始在Windows平台上开发有限元程序，称作"Native Windows"版本，同时还有在PC机上的Linux操作系统环境中开发的有限元程序包。  
  
国内发展情况和前景1979年美国的SAP5线性结构静、动力分析程序向国内引进移植成功，掀起了应用通用有限元程序来分析计算工程问题的高潮。这个高潮一直持续到1981年ADINA非线性结构分析程序引进，一时间许多一直无法解决的工程难题都迎刃而解了。大家也都开始认识到有限元分析程序的确是工程师应用计算机进行分析计算的重要工具。但是当时限于国内大中型计算机很少，大约只有杭州汽轮机厂的Siemens7738和沈阳鼓风机厂的IBM4310安装有上述程序，所以用户算题非常不方便，而且费用昂贵。PC机的出现及其性能奇迹般的提高，为移植和发展PC版本的有限元程序提供了必要的运行平台。可以说国内FEA软件的发展一直是围绕着PC平台做文章。在国内开发比较成功并拥有较多用户（100家以上）的有限元分析系统有大连理工大学工程力学系的FIFEX95、北京大学力学与科学工程系的SAP84、中国农机科学研究院的MAS5.0和杭州自动化技术研究院的MFEP4. 0等。但正如上面所述，国外很多著名的有限元分析公司已经从前些年对PC平台不屑一顾转变为热衷发展，对国内FEA程序开发者来说发展PC版本不再具有优势，而应该从下面几方面加以努力：  
  
1. 研究开发求解非固体力学和交叉学科的FEA程序经过几十年的研究和发展，用于求解固体力学的有限元方法和软件已经比较成熟，现在研究的前沿问题是流体动力学、可压缩和不可压缩流体的流动等非固体力学和交叉学科的问题。由于国内没有类似功能的商品化软件，所以国外的软件就卖得非常贵。为了破这种垄断局面，我们必须发展有自主版权、用于分析流体等非固体力学和交叉学科的软件。因为流体力学问题远比固体复杂得多，而且很少有现成的软件可以借鉴，所以需要投入大量的人力和经费。这就必须有国家和大型企业集团来支持。  
  
2. 开发具有中国特色的自动建模技术和GUI开发建模技术和GUI的投入比前述课题要少得多，但却可以大大提高FEA软件的性能和用户接受程度，从而起到事半功倍的效果。国内不少人在这方面做了很多工作，但是由于当时PC机上的图形支撑环境有限，所以开发的效果都不甚理想。Windows中提供了OpenGL图形标准，为在PC机上应用可视化图形技术开发GUI提供了强有力的工具。OpenGL是当今国际上公认的高性能图形和交互式视景处理标准，应用它开发出来的三维图形软件深受专业技术人员的钟爱，目前世界上占主导地位的计算机公司都采用了这一标准。正如前面所述，近年来国外有的FEA程序已抛开仿真软件，直接在Windows平台上开发有限元程序。杭州自动化技术研究院1997-1999年采用OpenGL图形标准和相应的Visual C++等编程工具，在PC机上成功地开发了一套可视化有限元程序包。它能直观地通过对"菜单"、"窗口"、"对话框"和"图标"等可视图形画面和符号的操作，自动建立有限元分析模型，并以交互方法式实现计算结果的可视化处理，因而可大大提高有限昂分析的效率和精确性，也便于用户学习和掌握。  
  
3. 与具有我国自主版权的CAD软件集成前面已经讲过，当今有限元方法的一个重要特点是和CAD软件的无缝集成。作为我国自行开发的FEA程序，首先要考虑和我国自主版权的CAD软件集成。因为有限元分析主要用于形状比较复杂的零部件，所以要和具有三维造型功能和CAD软件集成，使设计和分析紧密结合、融为一体。   
   
  

使用 ANSYS 进行 CFD 流体力学计算的一些技巧
   
TYPE：使用技巧 TIME：2003-5-23  
AUTHOR：金泰木 URL：四方机车车辆厂客车产品开发部  
HIT：2737  DAYHIT：16  WEEKHIT：2734 　  
   
关于计算流体力学主要有以下几个主要问题大家比较关心：（原稿：金泰木，四方机车车辆厂客车产品开发部）  
1． 关于瞬态计算的问题  
2． 关于建模的问题  
3． 关于网格化的问题  
4． 关于动画显示的问题  
5． 关于交变载荷的问题  
  
一、 关于第一个问题的解答：  
计算瞬态设置参数与稳态不同，主要设置的参数为：  
1． FLDATA1,SOLU,TRAN,1 设置为瞬态模式  
2． FLDATA4,TIME,STEP,0.02, 自定义时间步时间间隔0.02秒  
3． FLDATA4,TIME,TEND,0.1, 设置结束时间0。1秒  
4． FLDATA4,TIME,GLOB,10, 设置每个时间步多少次运算  
5． fldata4a,time,appe,0.02 设置记录时间间隔  
6．SET,LIST,2 查看结果  
7．SET,LAST 设为最后一步  
8．ANDATA,0.5, ,2,1,6,1,0,1 动态显示结果  
以上为瞬态和稳态不同部分的设置和操作，特别是第五步。为了动态显示开始到结束时间内气流组织的情况，还是花了我们很多时间来找到这条命令。如果你是做房间空调送风计算的，这项对你来说非常好，可以观察到从开空调机到稳定状态的过程。  
二．关于建模的问题  
大家主要关心的建模问题是模型的导入和导出，及存在的一些问题。这些问题主要体现在：  
1． AUTOCAD建模导出后的格式与ANSYS兼容的只有SAT格式。PROE可以是IGES格式或SAT格式。当然还有其它格式，本人使用的限于正版软件，只有上述两种格式。SAT格式可由PROE中导出为IGES格式。ANSYS默认的导入模型为IGES格式的图形模型。  
2． 使用AUTOCAD一般绘制界面比较复杂的拉伸体非常方便。如果是不规则体，用PROE和ANSYS都比较方便，当然本人推荐用ANSYS本身的建模功能。对于PROE，因为它的功能强大，本人推荐建立很复杂的模型如变截面不规则曲线弯管（如血管）。  
3． 导入过程中会出现默认选项和自定义选项，一般本人推荐使用自定义选项，以避免一些操作带来的问题。有时出现显示只有线而没有面颜色的问题，可以用命令：  
/FACET,NORML 来解决这个问题。  
三．关于网格化的问题。  
网格化对结果影响很大，如果网格化不合理，出现的结果会不准确，或者计算时不收敛。更甚者，网格数量太大，减慢求解速度。对计算流体力学来说，实际应用中三维问题偏多，计算量一般非常大，由于ANSYS采用的是有限元，所以同有限差分比较来说，收敛慢，内存需要量大。但这并不是说水平不如有限差分的流体计算软件。ANSYS的计算结果直观性较好，特别对涡流的处理很形象很准确（其它软件往往看不到该有的涡流，给人的感觉太粗糙）。当然对于稍大的模型，就有点力不从心的感觉。  
首先运行前合理设置内存需求量，充分估计本人计算机的计算能力，本人不愿意看到黄牛拉火车的现象（可能拉不动）。合理划分网格可以大大节省机时，还可以避免自动网格划分中带来的过度疏密。本人建议首先采用手动设置的网格划分，可以大大减少网格数量。对于长大模型，推荐采用比例长度缩小模型，网格化后再比例放大（因为四面体网格不能使用SWEEP）。对于小尺寸或者风口位置，应该尽可能地加密，可以很好的实现收敛。网格化有很多技巧，需要在实践中自己摸索。  
四．动画显示  
在计算完以后，我们最关心的是是用图形方式查看计算的结果，当然最好是动画显示（如粒子流），可以非常形象地显示分析后的效果。ANSYS提供了这方面的强大功能，我们常用的有以下几种：  
1． 切片  
2． 粒子流 ANFLOW,10,0, , ,1.03840361777,0.375E-02,0.2  
下面是切片显示设置的命令流：  
/TYPE,1,8 8指的是Q-SLICE Z_BUFFER  
/CPLANE,1 当前工作面WORKSPACE  
……  
PLNSOL,V,SUM,0, PLNSOL为查看云图  
如果是粒子流，则命令为：  
TRPOIN,P51X,P51X 在工作面上定义多个跟踪点  
。。。。。。  
ANFLOW,10,0, , ,1.03840361777,0.375E-02,0.2 通过跟踪点的离子流线  
对于三维粒子流，如果设置为切片模式，切片会遮挡流线，所以应该设置为：  
TRPOIN,P51X,P51X 在工作面上定义多个跟踪点  
/GRAPHICS,POWER  
/TYPE,1,0 消隐类型  
/CPLANE,0 正常视图  
/SHADE,1,1  
ANFLOW,10,0, , ,1.03840361777,0.375E-02,0.2 通过跟踪点的离子流线  
以上各项的菜单位置为：  
1. Main菜单下的General PostProc子菜单plot result菜单下flow trace\Defi Trace PT+。  
2. Utility菜单plotctrls子菜单下style\hidden-line options\下，/TYPE设置为nohidden；/CPLANE设置为nomal to view。  
ANSYS未提供切片的动态显示，如果你的结果显示速度可以的话，编写一个程序可以实现切片的动画。其中的一个循环命令流见下：  
WPOF,,, .1 步进量为0.1  
PLNSOL,V,SUM,0 显示矢量切片  
/WAIT,0.1 暂停0.1秒  
只要控制开始到结束的步数，用循环操作来实现。本人用VB编写了一个生成移动的命令流，运行时输入移动方向和移动距离、步进量、暂停时间，确定后生成两个宏文件（向前和返回）。在程序中调用即可。