论文部分内容阅读
CAD与CAE一体化已成为工程及科学计算领域中的一大研究热点。目前商业CAE分析软件普遍采用有限单元法,其分析模型是通过网格生成模块对CAD模型离散而获得的。CAD与CAE对应模型表达方式不同,前者是连续的参数模型,后者为近似的离散网格模型。因此,虽然很多CAD软件中添加了CAE分析功能,但还不能从根本上实现两者的无缝集成。基于边界积分方程的边界面法具有许多优点,例如一般只需要离散问题域的表面且可采用不连续单元。在边界面法中,对边界的数值积分和场变量的插值都是在几何模型边界曲面的二维参数空间里进行的。其几何数据直接由参数曲面计算得到,而不是通过分段多项式插值近似,从而能够避免几何误差,有利于实现CADCAE的一体化。本文基于边界面法,在CADCAE一体化的程序设计实现和理论方面做出一系列的探索。主要完成了以下研究工作:(1)为实现对复杂三维实体的表面网格自动生成,本文利用C++语言构建了一个实体表面网格自动生成框架。该框架分为公共几何接口模块、拓扑修复、网格尺寸场、曲线离散、网格数据管理和网格生成方法六个部分。本文设计了一个统一的网格数据管理器,以便于程序的调整、维护及扩充。在该框架中可方便地添加和实现各种二维网格生成方法,可将所得的网格单元转换为计算单元,且可导出到其他模块(比如体网格生成模块)。(2)为实现对复杂三维实体的表面三角形网格自动生成,本文提出了一种将推进波前法与Delaunay方法相结合的方式来生成网格。该方法利用推进波前法生成面内点,之后采用Delaunay内核插点算法将该点加入到网格中。本文添加了对前沿的识别,将前沿分为活动前沿及非活动前沿,只对满足要求的活动前沿插入新点。创建了基于曲面邻近特征的自适应尺寸场,并将其作为背景网格,使得所生成的网格能较好地反映曲面邻近特征。该方法能够生成较高质量的表面网格,且可以保证网格生成的收敛。同时,实现了栅格法,使得对于几类曲面特征的网格生成规模大大降低。(3)为实现对具有焊缝的钢架结构的曲面网格生成,鉴于边界面法的特点,本文首先直接基于连续曲面模型,完成了短边、窄面和不光滑边界等3类曲面缺陷几何特征的拓扑修复。拓扑修复中涉及的所有操作皆为虚操作,不改变原有模型的几何定义。并且提出了基于栅格法的混合网格生成方式,对不同的子区域采用不同的网格生成方法,最后将子区域上的网格合并得到最终的曲面网格。与传统的网格生成方式对比,大大降低了网格规模。(4)为求解边界积分方程中的弱奇异及近奇异积分,提出了一种通用的自适应单元细分方法。该算法通过一系列以源点为球心且半径不断减小的球与单元相交来细分单元,能够对靠近源点的积分片进行自动加密细分,并在生成积分片的过程中额外添加点的迁移算法来保证每个积分片具有较“好”的形状及大小。该算法适用于任意形状的单元,而且源点可以在任意位置。使用本文的方法所获得的计算精度有显著的提高。而且,在要求获得相同数量级计算精度的情况下,该算法需要少得多的高斯积分点,因而大大提高了计算效率。(5)在多域边界面法实现中,提出了一个域编号序列优化算法来减小总体组装矩阵的带宽。由于矩阵的稀疏结构与未知量在整体方程组系统中的顺序直接相关,因此采用了一种排序策略来获得最佳的块结构。该域序号优化算法的优点是使得矩阵中非零块尽可能地靠近主对角线。在该算法中,首先通过考虑域的自由度及连通性产生一个或者多个层次结构。其次,对于每个连续的层次结构,分别计算这些层次结构的带宽。最后,根据具有最小带宽的层次结构对域进行重新排序。数值算例表明采用该算法对整体矩阵进行LU分解的时间显著减少,同时内存开销也显著降低。