一种网格优化方法、网格生成器、存储介质技术

本发明专利技术公开了一种网格优化方法、网格生成器、存储介质，包括以下步骤：步骤S1、获取网格生成器生成的初始网格，对所述初始网格进行静力分析，得到所述初始网格的应力场，判断应力场是否满足计算需求，若不满足，执行步骤S2；步骤S2、计算所述初始网格的估计误差，得到所述初始网格的网格节点上的度量张量场；步骤S3、计算每个所述网格节点的度量张量场的特征值，获取所述网格节点的约束尺寸，根据约束尺寸生成优化网格；步骤S4、令优化网格为初始网格，重新执行步骤S1进行下一次迭代，直至初始网格满足预设要求。本发明专利技术自动调整网格尺寸，生成满足预设要求的网格，迭代计算，节约了计算资源，提高计算精度的同时提高计算效率。

【技术实现步骤摘要】
一种网格优化方法、网格生成器、存储介质
本专利技术涉及数据处理
，尤其涉及一种网格优化方法、网格生成器、计算机可读存储介质。
技术介绍
随着CAE仿真技术的不断发展，如何同时提升计算效率和提高计算精度，是急需解决的问题。目前可行的方法包括借助并行计算、高性能计算机、提升求解算法的效率等。网格划分也可同时提高计算效率及计算精度。在需要精细计算的区域，如应力集中、大形变区域，使用小尺度网格单元，在不需要精细计算区域，可以使用大尺度网格单元，起到节约计算资源，提升计算效率的目的。但是如何合理的划分的网格、预估需要和不需要精细计算的区域非常依赖工程师的处理经验。如图1所示，为一个几何模型，在网格尺寸为10的约束尺寸下网格生成器生成的四面体网格，其中包含2392个网格单元。图2是基于图1的网格应力值示意图，设置模型材料属性、固定边界条件、施加载荷后获得的应力仿真结果，可以看到应力最大值仅为209Mpa。根据经验及公式分析，该应力数值跟真实值存在较大差距，其网格划分的合理程度有限，仍需优化。专利技术内容为了克服本文档来自技高网...

【技术保护点】
1.一种网格优化方法，其特征在于，包括以下步骤：/n步骤S1、获取网格生成器生成的初始网格，对所述初始网格进行静力分析，得到所述初始网格的应力场，判断所述初始网格的应力场是否满足预设需求，若不满足，则执行步骤S2；/n步骤S2、计算所述初始网格的估计误差，得到所述初始网格的网格节点上的度量张量场；/n步骤S3、计算每个所述网格节点的度量张量场的特征值，获取所述网格节点的约束尺寸；根据所述约束尺寸生成优化网格；/n步骤S4、令所述优化网格为初始网格，重新执行步骤S1进行下一次迭代，直至所述初始网格满足预设要求。/n

【技术特征摘要】
1.一种网格优化方法，其特征在于，包括以下步骤：
步骤S1、获取网格生成器生成的初始网格，对所述初始网格进行静力分析，得到所述初始网格的应力场，判断所述初始网格的应力场是否满足预设需求，若不满足，则执行步骤S2；
步骤S2、计算所述初始网格的估计误差，得到所述初始网格的网格节点上的度量张量场；
步骤S3、计算每个所述网格节点的度量张量场的特征值，获取所述网格节点的约束尺寸；根据所述约束尺寸生成优化网格；
步骤S4、令所述优化网格为初始网格，重新执行步骤S1进行下一次迭代，直至所述初始网格满足预设要求。


2.如权利要求1所述的网格优化方法，其特征在于，所述步骤S2中计算所述初始网格的估计误差，得到所述初始网格的网格节点上的度量张量场包括以下步骤：
步骤S21、计算每个所述网格节点与所有连接节点的应力梯度值，根据所述应力梯度值及连接节点的数量计算平均应力值梯度；所述连接节点为连接所述网格节点的节点；
步骤S22、令表示连接网格节点和连接节点的网格连接线，根据每个网格节点的平均应力值梯度和连接节点平均应力值梯度的，计算网格连接线的估计误差；
步骤S23、计算所述网格连接线的长度变化值，其中为预设最大误差值；
步骤S24、根据每条所述网格连接线的长度变化值，计算每个网格节点的度量张量场。


3.如权利要求2所述的网格优化方法，其特征在于，所述步骤S21包括以下步骤：
步骤S211、令表示连接网格节点和连接节点的网格连接线，网格节点和网格节点之间的应力梯度值为，计算所述网格节点与所有连接该网格节点的连接节点的应力梯度值，其中为网格节点的应力值，为连接节点的应力值；
步骤S212、计算网格节点与所有连接节点的平均应力值梯度，所述平均应力值梯度为所述网格节点与所有连接节点的应力值梯度的平均值。


4.如权利要求3...

【专利技术属性】
技术研发人员：赵佳欣，李会江，张军飞，
申请(专利权)人：广州中望龙腾软件股份有限公司，
类型：发明
国别省市：广东;44