三维有限元六面体网格几何自适应再生成方法

基于几何自适应的六面体网再生成方法,提出了一种以栅格法为基本方法,对有限元网格再划分过程中的网格再划分标准、几何形状的获取以及新旧网格之间物理量的传递等关键问题进行了研究.重点介绍了基于几何自适应的六面体网格再生成算法,首先对旧网格实体模型进行识别,根据实体模型几何特征建立加密源点信息场;然后采用栅格法生成核心网格并对...     

基于有限元弹性变形运算的全六面体网格生成方法研究

以映射法为基础并结合网格划分经验,提出了原理简单的六面体网格生成新办法.该方法根据物体轮廓选择初始网格,设定表面结点强制变形到目的曲面,经由有限元弹性计算确定内部节点的位置.在检查全体单元质量以后,调整畸形单元从而生成目的网格.通过为一个复杂的马头门模型构建全六面体网格,最后证明了本文所述方法的可行性.     

基于二重网格的动网格松弛法改进

相比传统的弹簧法等方法,基于球松弛算法的动网格松弛法在复杂边界大变形条件下可以得到质量更高的边界网格以及更大的极限变形量,但该方法在时间效率上还有提升的空间.引入二重网格,采用动网格松弛法进行稀疏网格的网格变形,将边界位移传递到整个网格计算域;再利用二重网格映射,将稀疏网格位移映射到原有计算网格的节点上.算例表明,改进...     

三维六面体有限元网格自动划分中的一种单元转换优化算法

本文采用十节点曲边四面体转换为六面体网格,并采用非线性约束优化算法取Ｌａｐｌａｃｉａｎ光滑处理算法有效地提高六面体单元的质量,实现了对任意实体的六面体网格自动划分。     

基于双重自适应的六面体网格再生成方法

提出了一种以栅格法为基本方法,基于几何特征和物理场量双重自适应的六面体网格再生成方法。首先,依据旧网格的表面曲率和几何特征,采用基于栅格法的几何自适应网格再生成方法,生成密度受控的基础网格;然后,将旧网格的物理场量传递到基础网格中;最后,采用有限元误差估计方法对新网格单元的计算误差进行估计,对误差较大的单元进行加密,减...     

一种新的准结构网格生成方法

在有限元分析中，高质量的结构网格可以有效地提高有限元分析的精度，但结构网格的几何适应性差，针对复杂边界的二维计算模型，现有的方法很难自动生成高质量的结构网格；而非结构网格几何适应性很好，但存在计算效率低和精度差等问题。提出了一种新的准结构网格生成方法，能够实现复杂区域的网格自动生成并且具有高网格质量。该方法首先对计算区域运用Delaunay三角剖分技术生成粗背景网格；然后利用背景网格，使用优化的Voronoi图生成过渡的蜂巢网格；最后，通过中心圆方法对蜂巢网格单元进行结构网格剖分。分析NACA0012翼型数值模拟结果表明，提出的新准结构网格生成方法能够对边界复杂的模型自动生成高质量的网格，并且通过三种不同拓扑类型网格计算结果相互对比及与实验结果对比，证明准结构网格具有高计算精度。     

复杂三维组合曲面的有限元网格生成方法

提出一种基于映射法的复杂三维组合曲面的有限元网格全自动生成方法。通过引入虚边界解决了闭合曲面在参数域中边界不完整的问题；通过调节虚边界提高了复杂组合曲面网格生成的质量。改进二维多边形区域的裁减算法，解决了闭合曲面在参数空间中的边界环形成问题。对曲面片公共边界进行统一离散化处理，以满足有限元网格的相容性要求。以边界表示(B—Rep)数据结构为基础，实现了组合曲面全自动网格剖分的总成算法．改进了曲面网格剖分布点算法，并结合局部连接、诊断交换等技术，优化了网格的整体质量。     

基于转换模板的三维实体全六面体网格生成方法

以四面体一六面体基本转换模板为基础,提出一系列具有伸缩性的扩展转换模板．可根据需要选择不同模板及其组合,将四面体分解为不同数量、不同密度过渡形式的六面体单元。这样,初始四面体网格不需要划分得很细,生成的六面体单元数量也可以通过采用不同规格的扩展转换模板而得到控制。提出了基于CAD几何造型的边界节点坐标修正方法．从而使边界网格更好地拟合几何模型边界。     

岩石、混凝土类材料断裂破坏有限元数值模拟中的有效方法

岩石、混凝土类材料断裂破坏有限元数值模拟中的网格重划,依据单元畸变和裂缝介质间的单元干涉作为网格重划判据,采用几何体重构技术把几何实体分解成能在ANSYS上实现六面体网格划分的几个部分,利用体积判断法确定新结点在旧单元的单元编号,在场量传递上采用基于解析性质的等参有限元逆变换,把旧网格场量信息传递到新网格中。本文对ANSYS进行二次开发,实现了三维网格重划,网格重划采用单元畸变和界面干涉两个判据,在网格再划分前进行几何体重构,提取变形后的点线面信息重新生成实体,充分利用AN-SYS的函数和体积判断法找到新结点在旧网格中的位置,在新旧网格间的场量传递中采用基于解析逆等参单元法。在平台上实现了三维有限元网格重划技术,最后利用方料的单轴压缩断裂模拟计算检验了传递前后等效塑性应变分布用载荷信息的变化,证明了所开发系统的正确性。     

一种基于双前沿推进的边界层网格生成算法

高雷诺数粘性流动模拟对边界层内的网格正交性有特殊要求.对于复杂外形,这类问题的网格自动化生成十分困难.面向该问题,提出一种双前沿推进思想,并形成一种面向复杂几何外形的边界层网格全自动生成算法.结合多种网格技术处理局部几何特征以保证边界层网格的质量.双前沿推进思想同时适用于多块结构网格和混合网格的边界层网格生成.多个模型...     

A high‐order element based adaptive mesh refinement strategy for three‐dimensional unstructured grid

Adaptive mesh refinement (AMR) shows attractive properties in automatically refining the flow region of interest, and with AMR, better prediction can be obtained with much less labor work and cost compared to manually remeshing or the global mesh refinement. Cartesian AMR is well established; however, AMR on hybrid unstructured mesh, which is heavily used in the high‐Reynolds number flow simulation, is less matured and existing methods may result in degraded mesh quality, which mostly happens in the boundary layer or near the sharp geometric features. User intervention or additional constraints, such as freezing all boundary layer elements or refining the whole boundary layer, are required to assist the refinement process. In this work, a novel AMR strategy is developed to handle existing difficulties. In the new method, high‐order unstructured elements are first generated based on the baseline mesh; then the refinement is conducted in the parametric space; at last, the mesh suitable for the solver is output. Generating refined elements in the parametric space with high‐order elements is the key of this method and this helps to guarantee both the accuracy and robustness. With the current method, 3‐dimensional hybrid unstructured mesh of huge size and complex geometry can be automatically refined, without user intervention nor additional constraints. With test cases including the 2‐dimensional airfoil and 3‐dimensional full aircraft, the current AMR method proves to be accurate, simple, and robust.     

Adaptive panel methods and their integration with CAD

The work outlined below presents simple but effective adaptive meshing algorithms for boundary integral methods modelling inviscid flows (panel method) using the IGES standard for describing geometry. By using certain IGES entities in describing the boundary, CAD-derived geometry may be used such that the geometric integrity of the boundary is maintained after an adaptive redistribution of the mesh. Three types of error estimators are tested and all are shown to produce a more accurate representation of the flow phenomena for the same number of panels as compared with a uniform mesh distribution.     

适应复杂外形的三维粘性混合网格生成算法

提出一类适应复杂外形的粘性混合网格生成算法。表面网格由前沿推进三角形曲面网格程序获得,边界层布置各向异性的三棱柱体网格,远物面区域采用Delaunay方法生成四面体网格。针对模型的复杂几何特征,综合采用了各种网格处理技术,以保证边界层网格的质量,并避免算法失效问题。网格实例及计算结果表明了本文算法的实用性及和效性。     

基于背景网格变形的动态网格移动方法

基于Delaunay背景网格插值技术的动态网格生成方法无需迭代计算,效率较高。但对复杂构形大幅运动的动边界问题,尤其当边界大幅转动时,背景网格极易交叉重叠。重新生成背景网格和重新定位网格节点信息不仅费时而且会导致网格质量的严重下降。本文提出改进的基于背景网格的动态网格变形方法,通过在初始Delaunay背景网格中添加辅助点,生成一层新的背景网格和新的映射关系;采用ball-vertex弹簧法驱动新背景网格的变形,进而牵动目标网格的变形。算例表明,本文提出的动态网格变形方法对所关心区域的网格具有良好保形性,边界可转动更大角度而不会出现网格交叉重叠问题,总体上提高了动态网格更新的效率和质量。     

Shape morphing and topology optimization of fluid channels by explicit boundary tracking

An integrated shape morphing and topology optimization approach based on the deformable simplicial complex methodology is developed to address Stokes and Navier‐Stokes flow problems. The optimized geometry is interpreted by a set of piecewise linear curves embedded in a well‐formed triangular mesh, resulting in a physically well‐defined interface between fluid and impermeable regions. The shape evolution is realized by deforming the curves while maintaining a high‐quality mesh through adaption of the mesh near the structural boundary, rather than performing global remeshing. Topological changes are allowed through hole merging or splitting of islands. The finite element discretization used provides smooth and stable optimized boundaries for simple energy dissipation objectives. However, for more advanced problems, boundary oscillations are observed due to conflicts between the objective function and the minimum length scale imposed by the meshing algorithm. A surface regularization scheme is introduced to circumvent this issue, which is specifically tailored for the deformable simplicial complex approach. In contrast to other filter‐based regularization techniques, the scheme does not introduce additional control variables, and at the same time, it is based on a rigorous sensitivity analysis. Several numerical examples are presented to demonstrate the applicability of the approach.