本科毕业设计论文_外文翻译对由ansys开发的大型工程模型的降阶.docVIP

对由ansys开发的大型工程模型的降阶 Evgenii B. Rudnyi 和 Jan G. Korvink IMTEK微控技术研究所 弗赖堡大学 Georges-K ohler-Allee,103 D - 79110,德国弗赖堡 { rudnyi,korvink } @imtek.de http://www.imtek.uni-freiburg.de/simulation/ 摘要 工程师能够在ANSYS开发有限元模型运用现代模型降阶技术。我们于如何从ANSYS和C ++上实现执行模型中提取所需的信息与结构力学和热有限元模型的实例计算成本。 商业有限元软件包含两个几乎独立的模块（见图1）[6]。第一个模块用于读取一阶动态系统或二阶系统二进制文件和装载有限元分析软件ANSYS， Mx¨+Ex˙+Kx=Bu. (2.3) y=Cx. 其中M，E和K是三个系统矩阵。第二个模块适用于模型降阶算法式（1.1）或式（2.3），也就是说，它找到一个低维的式V， X = V Z + （2.4）以便让我们在误差范围内保证重现逼近原始状态向量的瞬态行为。 从原方程的子空间投影可以发现，例如（1.2），我们有Er=EV, Ar=AV, Br=B, Cr=CV ，我们支持三种方法来看待二阶系统。当仿照瑞利阻尼阻尼矩阵E=αM+βK，我们可以保存为以α系数和β作为参数的简化模型[7]。在一般情况下，我们可以把二阶系统转换为一阶系统，或者用二阶阿诺尔迪算法[8]。 ANSYS还可以读取、写入原系统矩阵的矩阵市场格式[4]。大量的模型降阶基准已由ANSYS分析制定 [9]。  图1 ANSYS解决方案 2.1与ANSYS的接口技术 生成第一个模块是相当困难的，因为商业有限元程序包中的大多数用户并不是有能力去提取动力系统式（1.1）或式（2.3），因此，这不是一个简单的操作。 ANSYS是一个巨大的包，其行为并不是完全一致的。例如，下面描述的信息是不适用流体动力学模块的FLOTRAN。 ANSYS软件读取二进制EMATfile矩阵与元素，以组装全局系统矩阵，记录文件格式，为ANSYS提供了一个Fortran子程序库[10]。通过例子可以发现更多的ANSYS代码[6]。 ANSYS有一个特殊的命令，称为部分解决带动完全求解 ， 可以不通过真正的解决方案阶段，评价矩阵的元素。这让我们对一个给定的模型拥有有效的EMAT文件。然而，这需要克服以下问题： - EMAT文件中不包含有狄利克雷边界条件或等式约束的信息。他们应分别提取。 - EMAT文件只有对元素矩阵有贡献的载荷向量。如果应用的负载是节点力或加速度，这个信息也应该被单独提取。 - 这矩阵要有必要的组装完整矩阵的元素。 在解决方案阶段，ANSYS可以写一个完整的二进制文件与系统矩阵匹配。当我们开始用ANSYS 5.7时，该文件不包含负载向量（输入矩阵）。而从那时起，已经有许多变化。ANSYS6.0保留所有的原始矩阵，载荷向量，文件中的Dirichlet和等式约束等完整的文件。 ANSYS8.0允许我们只做出装配（相当于一个EMAT的部分解决方案）就书写完整的文件。现在也可以把信息从完整文件转化为Harwell文件信息。因此，自ANSYS8.0以来，它可以有效地使用完整的文件。然而，根据分析类型的完整文件可能包含原来的非刚度矩阵，而不是一个系统矩阵的线性组合。 当前版本的ANSYS有更多的EMAT文件，是构造式（1.1）或式（2.3）的主要来源。我们已经开发的ANSYS宏包含Dirichlet和等式约束和节点力的附加信息。如果不使用完整的文件，这将是困难的，例如，在有加速负载的情况下提取载荷向量， ANSYS的EMAT文件不能充分表示负荷向量，当有多个输入时，如输入式（1.1）或式（2.3）时用户应该： - 删除以前应用的负载，- 申请一个新的负载，- 生成矩阵。 为了改进这一过程，第二个策略也被允许用户不删除以前的负载，在这种情况下，ANSYS会纠正在第一阶段结束的每一个新的负载向量（包含所有先前的载体）。2.2运行模型降阶算法 在克雷洛夫空间的基础上，我们可以通过一个非常有效的计算[11,8]，获得具有优良逼近性质的低维子空间（2.4）。当前版本的ANSYS实现了更多的Arnoldi 算法[11]，以支持多个输入，块大小等于输入数量。 每一步迭代Krylov子空间，都需要我们计算矩阵向量积为一阶系统，例如， （2.5） 其中h是向量，该系统矩阵高维并且可降阶，是一个不明确的计算结果。唯一可行的解决方案是解决如下的线性方程组的