65高性能计算应用.pdfVIP

高性能计算应用 65 基于高性能计算的曲轴系统动力学疲劳仿真分析  丁峻宏 邓雪滨 李根国 上海超级计算中心 上海201203  张晓云 上海交通大学机械与动力工程学院 上海200240 摘要： 针对曲轴的动力学仿真和疲劳寿命计算，建立了一种曲轴系统大规模直接计算模型，借助 高性能计算(HPC)技术并利用显式有限元算法实现了动力学模型的直接求解，结果详细描述了包 括曲轴强度和变形在内的曲轴系统动力学特性。采用全寿命分析方法直接对动力学仿真结果进 行曲轴寿命及安全系数的计算，基于应力－时间历程等动态结果进行的疲劳分析结果显得更为 真实。计算过程和结果证明了曲轴系统直接动力学及疲劳仿真分析方法的可行性和有效性。 关键词：曲轴，动力学，有限元，疲劳寿命，高性能计算 引言 模型；然后采用并行性能较好的显式有限元算法进 目前针对曲轴的疲劳分析方法主要包括：传统 行计算，得出曲轴的应力应变历程；最后利用疲劳 准静态疲劳分析方法简单易行，但其仅考虑最大应 分析软件对动力学计算结果进行处理，并对曲轴的 力而且不能反映振动特性和影响；目前应用较多的 强度和寿命分别进行评估。 是将有限元方法和多体系统仿真相结合的动力学分 析方法，其关键技术是模态分析和线性叠加，该方 1.1显式算法基本方程 计算模型中节点在t时刻的运动方程表示为： 法计算开销小，在解决不太复杂的机构（系统）疲 n 劳问题时，能体现出较好的计算效率；第三种是直 接动力学分析方法，可在建立整体曲轴系统的动力 其中 学模型基础上，直接求解得到曲轴的周期动态应力 式中M—总体质量矩阵 历程，继而对其寿命进行评估。由于该方法计算开 P—施加的外力及体力矢量 销相对较大，国内在动力学和疲劳分析方面的应用 F—内力矢量 int 尚不太多。 F—沙漏阻力 hg 近几年来，计算机硬件水平飞速发展，多核处 F —接触力矢量 contact 理器和并行计算的使用日趋广泛，直接动力学疲劳 采用显式中心差分法进行该方程的求解 计算方法也逐渐为人们所接受，在某些工程领域如 (3) 车载等复杂系统动力学疲劳寿命分析等问题的解决 上已有不少应用。本文将借助高性能计算平台，以 某泵机中的曲轴系统为例，详细阐述曲轴直接动力 式中 — 时刻的节点速度矢量 学疲劳分析的理论方法和计算过程。 — 时刻的节点坐标矢量 若已经求得t时刻的节点位置、加速度以及时步