Radioss 稳定性（重复精度）的对比分析.doc

Radioss 稳定性（重复精度）的对比分析 摘要:有限元计算的精度基本上可以概括为3个概念:绝对精度,相对精度,重复精度.在目前硬件软件日趋繁多复杂的背景下如何控制重复精度慢慢受到广大用户的重视.越来越多的工程师面临这样的问题,由于硬件的升级或者软件的升级,甚至计算模型的微小修改导致不稳定的计算结果,有时候甚至导致计算模型难以计算.无锡未知元汽车科技有限公司与无锡超级计算机中心以及Altair共同合作,把Altair公司的求解器RADIOSS安装在我国自制研发的神威4000A上进行一系列56km/h整车正面刚性墙碰撞计算,以此来验证RADIOSS作为求解器稳定的高重复精度的特性. 关键字:有限元计算精度,汽车整车碰撞模拟,并行运算, RADIOSS 简介 当今全世界几乎所有科技密集型产品包括飞机,汽车,电脑,手机等都以CAE(Computer Aided Engineering)为主要技术研发手段,正是因为CAE具有降低试验成本,减少研发周期等突出的特点。以CAE为主要手段的技术研发与其他方法一样有精度的要求,CAE模拟计算主要有绝对精度,相对精度,以及重复精度。 3个精度概念主要区别就是对比的目标值不同,绝对精度对比的目标值是产品试验结果,相对精度是研发工程师挑选出的CAE计算基础结果,重复精度是CAE重复计算结果发生最多的值。本文分析就是通过使用RADIOSS显示非线性对整车56km/h速度碰撞刚性墙模拟计算在不同硬件下的重复精度,最终可以通过重复精度体现RADIOSS作为求解器在有限元领域最复杂的显式非线性计算的稳定性。 背景 高性能计算机(HPC)的发展 高性能计算机发展到今天，现在已经达到了PFLOPS的级别，而且目前世界上还有若干国家或公司正在研制万万亿次浮点计算能力的超级计算机系统。接下来就来回顾一下曾经的超级计算机，峰值性能从百万次到万万次，经历了大概半个世纪的时间。1964年，CDC 6600系统的计算能力首次达到了MetaFLOP/s（每秒百万次浮点计算），1985年，Cray-2超级计算机系统则是突破GigaFLOPS（每秒十亿次浮点计算）。这些超级计算机系统由于年代久远，查阅其最初价格十分的困难，国外专家根据之后的CDC、Cray-2超级计算机系列的价格，推算出其当时的价格分别是在7百万美元和1千6百万美元。考虑通货膨胀因素换算到2010，则是4千9百万美元和3千2百万美元。超级计算机性能再次提升千倍则是到了1996年，ASCI Red超级计算机系统的计算能力首次达到TeraFLOPS（每秒万亿次浮点计算），价格在5千5百万美元左右。折合为2010年的7千6百万美元。2008年，IBM Roadrunner（走鹃）再次将计算性能提升千倍，首次超过PFLOPS（每秒千万亿次浮点计算），其费用为1亿美元，折合到2010年为1亿1百万美元。2011年底，日本的K Computer，首次超过10PFLOPS的浮点计算能力。将超级计算机带到了一个10PFLOPS的级别，而根据现在透露的消息，中国天河以及美国IBM的超级计算机，都计划在2015年进入100PFLOPS的级别。 近些年来，中国在超级计算机研究的也有着飞速的发展，尤其是在2010年登上TOP500首位的天河一号A型，更是打破了美国长年占据首位的传统。 RADIOSS的发展与主要特点 有限元方法思想的萌芽可以追溯到18世纪末，欧拉在创立变分法的同时就曾用与现代有限元相似的方法求解轴力杆的平衡问题，但那个时代缺乏强大的运算工具解决其计算量大的困难。Courant(1943)用最小势能原理和现代有限元法中的线性三角元求解st Venant弹性扭转问题，但未能引起足够重视。波音飞机工程师Turner，Clough等人在1956年首次将有限元法用于飞机机翼的结构分析，吹响了有限元的号角，有限元这一名称在1960年正式提出。Wilkins(1964)在DOE实验室的工作强烈地影响了早期的显式有限元方法，特别式命名为hydro－codes的软件, RADIOSS就是属于hydro－codes家族的求解器。他有如下主要的特征： 3维的拉格朗日网格计算方法 显式的时间积分方案有效的减少了计算时间 简化了有限元积分模型 相比隐性分析对输入/输出以及存储方面要求低 非交互式方法 惩罚式修复方法的接触 高矢量化的计算 RADIOSS是的线性和非线性仿真有限元求解器。它可以用来模拟结构、流体、流固耦合、金属薄板冲压和机械系统。这一功能强大的，多学科的解决方案使制造商能够最大限度地提高设计的耐久性、噪音和振动性能、耐撞性、安全性和可制造性，以使新产品更快地投向市场。 有限元计算的稳定性主要由三种精度概念决定:绝对精度,相对精度以及重复精度。 绝对精度顾名思义就是有限元计算结