CN114781232A - 有限元网格质量自动调整方法、设备和存储介质 - Google Patents

有限元网格质量自动调整方法、设备和存储介质 Download PDF

Info

Publication number
CN114781232A
CN114781232A CN202210682906.0A CN202210682906A CN114781232A CN 114781232 A CN114781232 A CN 114781232A CN 202210682906 A CN202210682906 A CN 202210682906A CN 114781232 A CN114781232 A CN 114781232A
Authority
CN
China
Prior art keywords
target
vertex
grid
mesh
requirement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210682906.0A
Other languages
English (en)
Other versions
CN114781232B (zh
Inventor
武振江
杨建森
王鹏
邹楠
李洪亮
张芳芳
张志国
王建海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CATARC Tianjin Automotive Engineering Research Institute Co Ltd
Original Assignee
CATARC Tianjin Automotive Engineering Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CATARC Tianjin Automotive Engineering Research Institute Co Ltd filed Critical CATARC Tianjin Automotive Engineering Research Institute Co Ltd
Priority to CN202210682906.0A priority Critical patent/CN114781232B/zh
Publication of CN114781232A publication Critical patent/CN114781232A/zh
Application granted granted Critical
Publication of CN114781232B publication Critical patent/CN114781232B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Image Generation (AREA)

Abstract

本发明实施例公开了一种有限元网格质量自动调整方法、设备和存储介质。其中，方法包括：S110、从CAE软件获取车辆零部件有限元网格特征，将质量不合格的三角形网格作为目标网格；S120、如果目标网格的至少一个顶点不满足内角要求，沿至少一个顶点中任一目标顶点到对边的垂线调整目标顶点的位置，得到第一网格，作为后续操作的目标网格；S130、如果目标网格的至少一个顶点不满足高度要求，沿至少一个顶点中任一目标顶点到对边的垂线调整目标顶点的位置；根据调整后的目标顶点对目标网格进行等比例放大，得到第二网格，作为后续操作的目标网格。本实施例实现了网格质量的自动化调整。

Description

有限元网格质量自动调整方法、设备和存储介质
技术领域
本发明实施例涉及有限元网格处理技术，尤其涉及一种有限元网格质量自动调整方法、设备和存储介质。
背景技术
汽车车身结构性能开发过程中，CAE性能仿真是必不可少的环节，而整个CAE仿真过程中汽车结构有限元网格划分占据了60%-80%的时间。为了使划分好的网格满足后续求解器的计算要求，需要对网格进行调整，以满足内角角度、翘曲、横纵比等质量要求。
现有技术中，完成初步的网格划分后，不满足质量要求的网格可以通过CAE软件自动筛选出来，随后通过人工操作拖动网格单元的节点，依靠经验拖动到大致合适的位置来改善质量网格，得到网格划分的最终结果。整个网格质量调整效率低，需投入人员与费用较多，降低了有限元网格划分的效率。
发明内容
本发明实施例提供一种有限元网格质量自动调整方法、设备和存储介质，实现了自动化的有限元网格质量调整流程，提高了网格划分效率。
第一方面，本发明实施例提供了一种有限元网格质量自动调整方法，包括：
S110、从CAE软件获取车辆零部件有限元网格特征，将所述有限元网格特征中质量不合格的三角形网格作为目标网格；
S120、如果所述目标网格的每个顶点均满足内角要求，进入S130；如果所述目标网格的至少一个顶点不满足内角要求，沿所述至少一个顶点中任一目标顶点到对边的垂线调整所述目标顶点的位置，得到第一网格，使调整后的目标顶点在所述第一网格中满足内角要求，并将所述第一网格作为后续操作的目标网格；重新执行S120，直到所述目标网格的每个顶点均满足内角要求；
S130、如果所述目标网格的至少一个顶点不满足高度要求，沿所述至少一个顶点中任一目标顶点到对边的垂线调整所述目标顶点的位置，使调整后的目标顶点满足高度要求；根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第二网格，并将所述第二网格作为后续操作的目标网格；重新执行S130，直到所述目标网格的每个顶点均满足内角要求和高度要求。
第二方面，本发明实施例提供了另一种有限元网格质量自动调整方法，包括：
S210、从CAE软件获取车辆零部件的有限元网格特征，将所述有限元网格特征中质量不合格的四边形网格作为目标网格；
S220、如果所述目标网格满足翘曲要求，进入S230；如果所述目标网格不满足翘曲要求，调整所述目标网格对角线两侧顶点的位置，得到第三网格，使所述第三网格满足翘曲要求，并将所述第三网格作为后续操作的目标网格；
S230、如果所述目标网格的每个顶点均满足内角要求，进入S240；如果所述目标网格的至少一个顶点不满足内角要求，将所述至少一个顶点中的任一顶点作为目标顶点，并确定所述目标网格的两个相邻点构成的对角线；沿所述目标顶点到所述对角线的垂线调整所述目标顶点的位置，得到第四网格，使调整后的目标顶点在所述第四网格中满足内角要求，并将所述第四网格作为后续操作的目标网格；重新执行S230，直到所述目标网格的每个顶点均满足内角要求；
S240、如果所述目标网格的每个顶点均满足最小高度要求，进入S250；如果所述目标网格的至少一个顶点不满足最小高度要求，确定所述至少一个顶点中任一目标顶点的最小高度对应的对边；沿所述目标顶点到所述对边的垂线调整所述目标顶点的位置，得到第五网格，使调整后的目标顶点满足最小高度要求；如果所述第五网格的每个顶点均满足内角要求，将所述第五网格作为后续操作的目标网格，进入S250；如果所述第五网格的至少一个顶点不满足内角要求，根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第六网格，将所述第六网格作为后续操作的目标网格；重新执行S240，直到所述目标网格的每个顶点均满足内角要求和高度要求；
S250、如果所述目标网格不满足横纵比要求，调整所述目标网格最短边和最长边的交点位置，使所述目标网格满足横纵比要求。
第三方面，本发明实施例提供了另一种有限元网格质量自动调整方法，包括：
从CAE软件获取车辆零部件的有限元网格特征；
将所述有限元网格特征中的每个质量不合格的三角形网格作为目标网格，对所述目标网格执行上述实施例提供的S110-S130的操作；同时，
将所述有限元网格特征中的每个质量不合格的四边形网格作为目标网格，对所述目标网格执行上述实施例提供的S210-S250的操作。
第四方面，本发明实施例还提供了一种电子设备，包括：
一个或多个处理器；
存储器，用于存储一个或多个程序，
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如任一实施例所述的有限元网格质量自动调整方法。
第五方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现任一实施例所述的有限元网格质量自动调整方法。
本发明实施例逐个调整三角形网格中不满足质量要求的目标顶点，每个顶点均完成内角调整后再进行逐个进行高度调整；在调整过程中将目标顶点的对边保持不变，通过目标顶点到对边的垂线限定目标顶点的调整轨迹，充分利用垂线与三个顶点的几何关系修改目标顶点的位置，保证调整后的目标顶点符合相应的质量要求，避免了人为调整中随机拖动顶点位置而无法保证调整后的顶点一定满足质量要求的问题。在进行高度调整时，考虑到三角形内角之间的强耦合性，在增大目标顶点的网格高度后对目标网格进行等比例放大，以保持各内角不变，避免了内角调整和高度调整之间的反复跳跃，提高了调整效率。整个调整策略易于实现，且能够保证调整结果满足质量要求，为网格质量调整提供了一种自动化的实现方式。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种板壳结构有限元网格的示意图；
图2是本发明实施例提供的一种三角形网格内角角度的示意图；
图3是本发明实施例提供的一种三角形网格网格高度的示意图；
图4是本发明实施例提供的一种四边形网格翘曲角度的示意图，其中图4(a)为对角线
Figure 222867DEST_PATH_IMAGE001
划分成的两个三角平面的法向量夹角，图4(b)为对角线
Figure 969106DEST_PATH_IMAGE002
划分成的两个三角平面的法向量夹角；
图5是本发明实施例提供的一种四边形网格内角角度的示意图；
图6是本发明实施例提供的一种四边形网格最小高度的示意图；
图7是本发明实施例提供的一种四边形网格横纵比的示意图；
图8是本发明实施例提供的另一种有限元网格质量自动调整方法的流程图；
图9是本发明实施例提供的一种三角形网格的示意图；
图10是本发明实施例提供的一种三角形网格内角调整的示意图，其中，图10(a)为目标顶点的内角小于第一阈值时内角调整的示意图，图10(b)为目标顶点的内角大于第二阈值时内角调整的示意图；
图11是本发明实施例提供的三角形网格三个顶点的网格高度的示意图，其中，图11(a)为顶点b的网格高度h b 的示意图，图11(b)为顶点a的网格高度h a 的示意图，图11(c)为顶点c的网格高度h c 的示意图；
图12是本发明实施例提供的一种三角形网格高度调整的示意图；
图13是本发明实施例提供的一种三角形网格等比例放大的示意图；
图14是本发明实施例提供另一种有限元网格质量自动调整方法的流程图；
图15为本发明实施例提供的一种四边形网格的示意图；
图16是本发明实施例提供的一种四边形网格内角调整的示意图，其中图16(a)为目标顶点的内角小于第五阈值时内角调整的示意图，图16(b)为目标顶点的内角大于第六阈值内角调整的示意图；
图17是本发明实施例提供的一种四边形网格最小高度调整的示意图；
图18是本发明实施例提供的一种四边形网格等比例放大的示意图；
图19是本发明实施例提供的一种四边形网格横纵比调整的示意图；
图20是本发明实施例提供的另一种有限元网格质量自动调整方法的流程图；
图21是本发明实施例提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将对本发明的技术方案进行清楚、完整的描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施例，都属于本发明所保护的范围。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明的描述中，还需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
本发明实施例提供一种有限元网格质量自动调整方法。为了详细说明该方法，优先说明汽车板壳结构零部件有限元网格的质量参数。图1是本发明实施例提供的一种板壳结构有限元网格的示意图。如图1所示，复杂的板壳结构经初步的网格划分后由大量的三角形网格和四边形网格组成，两种网格的质量参数不完全相同。
三角形网格单元的质量参数包括内角角度和网格高度。内角角度指三角形网格各顶点处的内角角度，如图2，三角形网格的内角要求通常为每个顶点的内角大于或等于第一阈值且小于或等于第二阈值。网格高度是指三角形网格中各顶点到对边的垂直距离，如图3所示，三角形网格的高度要求通常为每个顶点到对边的垂直距离大于或等于第三阈值。
四边形网格的质量参数包括；翘曲角度、内角角度、最小网格高度和横纵比。翘曲角度是指四边形网格中每条对角线划分成的两个三角平面的法向量夹角中较小的值；具体的，如图4所示，计算对角线
Figure 730389DEST_PATH_IMAGE002
划分成的两个三角平面的法向量
Figure 790749DEST_PATH_IMAGE003
和法向量
Figure 750614DEST_PATH_IMAGE004
的第一夹角，以及对角线
Figure 871017DEST_PATH_IMAGE001
划分成的两个三角平面的法向量
Figure 119596DEST_PATH_IMAGE003
和法向量
Figure 983647DEST_PATH_IMAGE004
的第二夹角，选择两个夹角中较小的一个作为四边形网格的翘曲角度；四边形网格的翘曲要求通常为为翘曲角度小于或等于第四阈值。内角角度指四边形网格各顶点处的内角角度，如图5所示；四边形网格的内角要求通常为每个顶点的内角大于或等于第五阈值且小于或等于第六阈值。最小高度是指四边形网格中各顶点到两条对边的垂直距离中较小的值，如图6中顶点b到两条对边的高度分别为h2和h3，其中h2的值更小，为最小高度；四边形网格的高度要求通常为每个顶点的最小高度大于或等于第七阈值。横纵比是指四边形网格中的最长边长度L max 与最短边长度L min 的比值，如图7所示；四边形网格的横纵比要求通常为横纵比小于或等于第八阈值。
表1给出一种汽车耐久领域的具体质量要求。
表1
Figure 798019DEST_PATH_IMAGE005
基于以上质量参数的介绍，图8是本发明实施例提供一种有限元网格质量自动调整方法的流程图。该方法适用于在有限元仿真中对汽车板壳结构零部件的三角形网格自动进行质量调整的情况。该方法可以由电子设备执行，如图8所示，该方法具体包括：
S110、从CAE软件获取车辆零部件有限元网格特征，将所述有限元网格特征中质量不合格的三角形网格作为目标网格。
利用CAE软件对车辆零部件的CAD模型进行初步网格划分后，得到该零部件的有限元网格特征，包括组成该零部件的各网格的编号，以及每个网格中各顶点的编号和在整车坐标系中的坐标等。得到有限元网格特征后，CAE软件能够自动筛选出其中不符合质量要求的有限元网格，包括三角形网格和四边形网格。本实施例将不符合质量要求的三角形网格作为目标网格，作为后续质量调整的对象。
此外，经过网格质量筛选完成后，三角形网格会存在多种网格质量问题，例如，仅不满足内角要求、仅不满足高度要求，以及同时不满足内角要求和高度要求；还可能存在多个顶点不满足内角要求、多个顶点不满足高度要求，甚至两种情况的组合状态等。而在质量调整过程中，调整某一顶点的网格高度，会导致该顶点的内角产生变化，反之亦然。因此本实施例设定三角形网格调整策略为优先调整内角，每个顶点均满足内角要求后，再进行高度调整。下面将按照该顺序对目标网格进行质量调整。
S120、如果所述目标网格的每个顶点均满足内角要求，进入S130；如果所述目标网格的至少一个顶点不满足内角要求，沿所述至少一个顶点中任一目标顶点到对边的垂线调整所述目标顶点的位置，得到第一网格，使调整后的目标顶点在所述第一网格中满足内角要求，并将所述第一网格作为后续操作的目标网格；重新执行S120，直到所述目标网格的每个顶点均满足内角要求。
S120完成三角形网格的内角调整。由于CAE软件仅能给出目标网格是否满足质量要求，不能确定具体的质量问题类型或哪个顶点存在质量问题。因此首先计算目标网格每个顶点的内角，判断目标网格是否满足内角要求。如图9所示，将三角形网格的三个顶点a、b、c按照环路顺序进行编号，并且按照编号顺序组成三个边向量V 1、V 2和V 3通过三个变向量的点乘即可计算出
Figure 56700DEST_PATH_IMAGE006
：
Figure 589312DEST_PATH_IMAGE007
（1）
由于计算机在进行空间几何操作中尽量规避除法与反三角函数的计算，因此直接通过上述余弦值判断是否满足内角要求。三角形网格的内角要求为每个顶点的内角大于或等于第一阈值
Figure 991475DEST_PATH_IMAGE008
且小于或等于第二阈值
Figure 863616DEST_PATH_IMAGE009
，因此只需判断每个顶点是否满足
Figure 591400DEST_PATH_IMAGE010
，其中，
Figure 345730DEST_PATH_IMAGE011
。基于不同的判断结果，S120包括以下几种可选实施方式：
第一种可选实施方式，每个顶点均满足内角要求，这时无需进行内角调整，直接进入S130进行高度调整。
第二种可选实施方式，存在一个顶点不满足内角要求，这时以该顶点为目标顶点，沿目标顶点到对边的垂线调整目标顶点的位置。具体的，根据目标顶点内角的大小，包括以下两种情况：
情况一、目标顶点的内角小于所述第一阈值，在所述目标顶点到对边的垂线上确定第一点，使所述第一点与所述目标顶点的两个相邻顶点的连线夹角等于所述第一阈值，并将所述目标顶点调整到所述第一点的位置。具体的，以图10(a)为例，目标顶点a的内角
Figure 551583DEST_PATH_IMAGE012
满足
Figure 543810DEST_PATH_IMAGE013
，说明
Figure 973654DEST_PATH_IMAGE014
，这时在目标顶点a到对边
Figure 684121DEST_PATH_IMAGE015
的垂线（即过目标顶点a的三角形的高所在的直线）上确定一点a’，使得
Figure 920762DEST_PATH_IMAGE016
和
Figure 767495DEST_PATH_IMAGE017
之间的夹角等于
Figure 633820DEST_PATH_IMAGE008
，将目标顶点a调整到a’的位置，由调整后的目标顶点a’和目标顶点的两个相邻顶点b、c构成第一网格
Figure 566004DEST_PATH_IMAGE018
，并将第一网格作为后续操作的目标网格。进一步的，a’的位置可以通过该垂线的方程和顶点b、c的坐标联合求解，或根据从垂线上的多个预设位置中选取最接近“
Figure 113660DEST_PATH_IMAGE016
和
Figure 877217DEST_PATH_IMAGE017
之间的夹角等于
Figure 852126DEST_PATH_IMAGE008
”的条件的位置作为a’的位置，本实施例不作限制。
情况二、目标顶点的内角大于所述第二阈值，在所述目标顶点到对边的垂线上确定第二点，使所述第二点与所述目标顶点的两个相邻顶点的连线夹角等于所述第二阈值，并将所述目标顶点调整到所述第二点的位置。具体的，以图10(b)为例，顶点a的内角
Figure 271606DEST_PATH_IMAGE012
满足
Figure 622953DEST_PATH_IMAGE019
，说明
Figure 241016DEST_PATH_IMAGE020
，这时在所述垂线上确定一点a’，使得
Figure 121247DEST_PATH_IMAGE016
和
Figure 526559DEST_PATH_IMAGE017
间的夹角等于
Figure 478334DEST_PATH_IMAGE009
，将目标顶点a调整到a’的位置, 由调整后的目标顶点a’和目标顶点之外的其他两个顶点b、c构成第一网格
Figure 154166DEST_PATH_IMAGE018
。需要说明的是，通常
Figure 470878DEST_PATH_IMAGE021
，图10(b)中的
Figure 864950DEST_PATH_IMAGE022
，仅用于说明本可选实施方式的方法，在实际应用中通常不会出现。
第三种可选实施方式，存在多个顶点不满足内角要求，这时从所述多个顶点中任选一顶点作为目标顶点，按照第二种可选实施方式提供的方法对目标顶点进行处理，得到一个第一网格，并将第一网格作为后续操作的目标网格；重新执行S120，直到最终的目标网格中每个顶点均满足内角要求。实际应用汇总，在内角要求设置合理的情况下，例如
Figure 354837DEST_PATH_IMAGE023
，每个顶点重复调整的次数不会很多。
S130、如果所述目标网格的至少一个顶点不满足高度要求，沿所述至少一个顶点中任一目标顶点到对边的垂线调整所述目标顶点的位置，使调整后的目标顶点满足高度要求；根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第二网格，并将所述第二网格作为后续操作的目标网格；重新执行S130，直到所述目标网格的每个顶点均满足内角要求和高度要求。
进入S130后，目标网格中每个顶点均满足内角要求，S130完成目标网格的高度调整。首先计算目标网格每个顶点的网格高度，即三个顶点的高，如图11中的h a 、h b 和h c 所示，可以通过海伦公式计算得出，本实施例对具体计算方法不作限制。
得到三个顶点的高度后，判断每个顶点是否满足高度要求。三角形网格的高度要求为每个顶点的网格高度大于或等于第三阈值H1，即判断每个顶点到对边的垂直距离是否满足
Figure 885176DEST_PATH_IMAGE024
。基于不同的判断结果，S120包括以下几种可选实施方式：
第一种可选实施方式，每个顶点均满足高度要求，这时无需进行高度调整，整个三角形单元的质量调整过程结束。
第二种可选实施方式，存在一个顶点不满足高度要求，这时以该顶点为目标顶点，沿目标顶点到对边的垂线调整目标顶点的位置，得到第一网格。具体的，确定所述目标顶点到对边的垂线，在所述垂线上确定上确定第五点，使所述第五点到所述对边的垂直距离等于所述第三阈值；将所述目标顶点调整到所述第五点的位置，由所述调整后的目标顶点，以及目标网格中除目标顶点之外的其它两个顶点构成第一网格。以图12为例，目标顶点a的高度为h 1，不满足高度要求，则在垂线上确定第五点a’，a’到对边的垂直距离=H1，将a调整到a’的位置，由调整后的目标顶点a’和目标顶点之外的其他两个顶点b、c构成第一网格
Figure 638368DEST_PATH_IMAGE018
。所述调整后的目标顶点在所述第一网格中一定是满足高度要求的。进一步的，a’的位置可以通过该垂线的方程和顶点b、c的坐标联合求解，或根据从垂线上的多个预设位置中选取垂直距离最接近H1的位置作为a’的位置，本实施例不作具体限制。得到第一网格后，目标顶点的内角一定会发生变化。由于三角形网格中三个内角之间的耦合性较高，当一个内角发生变化时，其余两个内角均会发生变化，因此第一网格大概率不满足内角要求。为了节省计算资源，本可选实施方式中不再重新对三个内角进行判断，而是直接根据所述调整后的目标顶点（即第五点）对所述目标网格进行等比例放大，得到第二网格。以图13为例，将顶点a调整到第五点a’的位置后，根据调整后的目标顶点a’对目标网格
Figure 50895DEST_PATH_IMAGE025
进行等比例放大，得到第二网格
Figure 282156DEST_PATH_IMAGE026
。等比例放大后，第二网格的三条边与目标网格的三条边平行，因此三个内角的角度相等，一定满足内角要求。
可选的，所述根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第二网格，包括：确定过所述调整后的目标顶点的两条直线，使所述两条直线分别平行于所述目标网格中过所述目标顶点的两条边；分别确定所述对边的延长线与所述两条直线的两个交点，由所述调整后的目标顶点和所述两个交点构成第二网格。仍以图13为例，确定过a’的两条直线，分别平行于目标网格
Figure 932580DEST_PATH_IMAGE025
中过顶点a的两条边
Figure 653412DEST_PATH_IMAGE027
和
Figure 520611DEST_PATH_IMAGE028
；分别确定目标顶点a的对边
Figure 555563DEST_PATH_IMAGE029
的延长线分别与两条直线的交点b’和c’，由调整后的目标顶点a’和所述两个交点b’和c’构成第二网格
Figure 794915DEST_PATH_IMAGE026
。可以看出，
Figure 686648DEST_PATH_IMAGE030
，因此，第二网格
Figure 277029DEST_PATH_IMAGE026
是目标网格
Figure 850093DEST_PATH_IMAGE025
的等比例放大，第二网格
Figure 271847DEST_PATH_IMAGE026
与目标网格
Figure 537743DEST_PATH_IMAGE031
中三个内角的角度相等，一定满足内角要求。
第三种可选实施方式，存在多个顶点不满足高度要求，这时从所述多个顶点中的任选一顶点作为目标顶点，按照第二种可选实施方式提供的方法对目标顶点进行处理，得到一个第二网格，并将所述第二网格作为后续操作的目标网格；重新执行S130，直到最终的目标网格中每个顶点均满足内角要求和高度要求。
S130完成后，最终的目标网格中每个顶点均满足内角要求和高度要求，三角形网格的质量调整结束。
本实施例逐个调整三角形网格中不满足质量要求的目标顶点，每个顶点均完成内角调整后再进行逐个进行高度调整；在调整过程中将目标顶点的对边保持不变，通过目标顶点到对边的垂线限定目标顶点的调整轨迹，充分利用垂线与三个顶点的几何关系修改目标顶点的位置，保证调整后的目标顶点符合相应的质量要求，避免了人为调整中随机拖动顶点位置而无法保证调整后的顶点一定满足质量要求的问题。在进行高度调整时，考虑到三角形内角之间的强耦合性，在增大目标顶点的网格高度后对目标网格进行等比例放大，以保持各内角不变，避免了内角调整和高度调整之间的反复跳跃，提高了调整效率。整个调整策略易于实现，且能够保证调整结果满足质量要求，为网格质量调整提供了一种自动化的实现方式。
图14是本发明实施例提供另一种有限元网格质量自动调整方法的流程图，适用于在有限元仿真中对汽车板壳结构零部件的四边形网格自动进行质量调整的情况。该方法可以由电子设备执行，如图14所示，该方法具体包括：
S210、从CAE软件获取车辆零部件的有限元网格特征，将所述有限元网格特征中质量不合格的四边形网格作为目标网格。
利用CAE软件对车辆零部件的CAD模型进行初步网格划分后，得到该零部件的有限元网格特征，包括组成该零部件的各网格的编号，以及每个网格中各顶点的编号和在整车坐标系中的坐标等。得到有限元网格特征后，CAE软件能够自动识别出其中不符合质量要求的有限元网格，包括三角形网格和四边形网格。本实施例将不符合质量要求的四边形网格作为目标网格，作为后续质量调整的对象。
四边形网格的网格参数包括翘曲、内角角度、最小高度和横纵比。经过网格质量筛选后，四边形网格会存在多种网格质量的问题。由于各种网格参数存在相互耦合的情况，比如调整四边形网格的内角角度会影响最小高度等，因此在调整四边形单元网格质量调整时需要兼顾执行效率，尽量少的反复调整某一网格质量指标。在各种质量参数中，翘曲与其他参数的耦合程度不高，而内角角度对最小高度和横纵影响较大。因此，本实施例设定四边形网格的调整策略为按照“翘曲、内角角度、最小高度和横纵比”的顺序逐依次调整。下面将按照该顺序对目标网格进行质量调整。
S220、如果所述目标网格满足翘曲要求，进入S230；如果所述目标网格不满足翘曲要求，调整所述目标网格对角线两侧顶点的位置，得到第三网格，使所述第三网格满足翘曲要求，并将所述第三网格作为后续操作的目标网格。
S220完成四边形网格的翘曲调整。由于CAE软件仅能给出目标网格是否满足质量要求，不能确定具体的质量问题类型或哪个顶点存在质量问题。因此首先计算目标网格的翘曲角度，判断目标网格是否满足翘曲要求。如图15所示，将四边形网格的四个顶点按照环路顺序编号为a、b、c、d，并按照编号顺序组成四个边向量V 1、V 2、V 3和V 4，内角角度分别为
Figure 615420DEST_PATH_IMAGE032
。图4显示了图15中的四边形网格的翘曲角度的计算过程，如图4所示，计算对角线
Figure 631656DEST_PATH_IMAGE002
划分成的两个三角平面的法向量
Figure 111179DEST_PATH_IMAGE003
和法向量
Figure 344714DEST_PATH_IMAGE004
的第一夹角，以及对角线
Figure 909687DEST_PATH_IMAGE001
划分成的两个三角平面的法向量
Figure 824554DEST_PATH_IMAGE003
和法向量
Figure 424162DEST_PATH_IMAGE004
的第二夹角，选择两个夹角中较小的一个作为四边形网格的翘曲角度。图4中，对角线
Figure 31861DEST_PATH_IMAGE002
对应的夹角最小，则将该夹角作为四边形网格的翘曲角度。四边形网格的翘曲要求为翘曲角度小于或等于第四阈值，基于不同的判断结果，S220包括以下几种可选实施方式：
第一种可选实施方式，目标网格满足翘曲要求，这时无需进行翘曲调整，直接进入S230进行内角调整。
第二种可选实施方式，目标网格不满足翘曲要求，这时调整所述目标网格对角线两侧的顶点位置，得到第三网格，使所述第三网格满足翘曲要求，并将所述第三网格作为后续操作的目标网格。具体的，首先计算所述翘曲角度与第四阈值的差值，并确定用于计算翘曲的对角线，即对应的法向量夹角最小的对角线；然后调整该对角线两侧的顶点位置，使所述对角线切分成的两个三角平面向各自相反的方向旋转所述差值的一半。结合图4，假设当前翘曲角度与第四阈值的差值为
Figure 588525DEST_PATH_IMAGE033
，用于计算翘曲的对角线为
Figure 979187DEST_PATH_IMAGE034
。调整
Figure 902143DEST_PATH_IMAGE002
两侧的顶点b和d的位置，使
Figure 680743DEST_PATH_IMAGE002
切分成的两个三角平面
Figure 187686DEST_PATH_IMAGE031
和
Figure 709934DEST_PATH_IMAGE035
分别绕对角线
Figure 815293DEST_PATH_IMAGE034
，各自向相反的方向旋转
Figure 764795DEST_PATH_IMAGE036
。
S230、如果所述目标网格的每个顶点均满足内角要求，进入S240；如果所述目标网格的至少一个顶点不满足内角要求，将所述至少一个顶点中的任一顶点作为目标顶点，并确定所述目标网格的两个相邻顶点构成的对角线；沿所述目标顶点到所述对角线的垂线调整所述目标顶点的位置，得到第四网格，使调整后的目标顶点在所述第四网格中满足内角要求，并将所述第四网格作为后续操作的目标网格；重新执行S230，直到所述目标网格的每个顶点均满足内角要求。
进入S230后，所述目标网格已满足翘曲要求，S230完成目标网格的内角调整。首先判断目标网格的四个顶点是否均满足内角要求。四边形网格的内角要求为顶点的内角小于或等于第五阈值
Figure 526077DEST_PATH_IMAGE037
且大于或等于第六阈值
Figure 586437DEST_PATH_IMAGE038
，因此判断每个顶点是否满足
Figure 811882DEST_PATH_IMAGE039
，其中，
Figure 932285DEST_PATH_IMAGE040
。基于不同的判断结果，S230包括以下几种可选实施方式：
第一种可选实施方式，每个顶点均满足内角要求，这时无需进行内角调整，直接进入S240进行高度调整。
第二种可选实施方式，存在一个顶点不满足内角要求，这时以该顶点为目标顶点，确定所述目标顶点的两个响铃顶点构成的对角线；沿目标顶点到所述对角线的垂线调整目标顶点的位置。以图16为例，目标顶点b不满足内角要求，确定所述目标网格abcd中目标顶点b的相邻顶点构成的对角线
Figure 915284DEST_PATH_IMAGE034
,沿目标顶点b到对角线
Figure 576073DEST_PATH_IMAGE002
的垂线
Figure 859287DEST_PATH_IMAGE041
调整目标顶点b的位置。具体的，根据目标顶点内角的大小，包括以下两种情况：
情况一、目标顶点的内角小于第五阈值，在所述垂线上确定第六点，使所述第六点分别与所述目标顶点的两个相邻顶点的连线夹角等于第五阈值，并将目标顶点调整到所述第六点的位置。具体的，如图16(a)所示，目标顶点b的内角
Figure 383547DEST_PATH_IMAGE042
满足
Figure 119422DEST_PATH_IMAGE043
，说明
Figure 583901DEST_PATH_IMAGE044
，这时在目标顶点b到对边
Figure 721621DEST_PATH_IMAGE002
的垂线
Figure 918247DEST_PATH_IMAGE041
上确定一点b’，使得
Figure 141418DEST_PATH_IMAGE045
和
Figure 409589DEST_PATH_IMAGE046
之间的夹角等于
Figure 401815DEST_PATH_IMAGE037
，将目标顶点b调整到b’的位置，由调整后的目标顶点b’和目标顶点之外的其它三个顶点a、c和d构成第四网格ab’cd。进一步的，b’的位置可以通过该垂线的方程和顶点a、c的坐标联合求解，或根据从垂线上的多个预设位置中选取最接近“
Figure 769343DEST_PATH_IMAGE045
和
Figure 479810DEST_PATH_IMAGE046
之间的夹角等于
Figure 286092DEST_PATH_IMAGE037
”的条件的位置作为b’的位置，本实施例不作限制。需要说明的是，通常
Figure 365781DEST_PATH_IMAGE047
，图16(a)中的
Figure 435368DEST_PATH_IMAGE048
，仅用于说明本可选实施方式的方法，在实际应用中通常不会出现。
情况二、目标顶点的内角大于所述第六阈值，在所述垂线上确定第七点，使所述第七点到所述目标顶点的两个相邻顶点的连线夹角等于所述第六阈值，并将所述目标顶点调整到所述第七点的位置。具体的，以图16(b)为例，顶点内角满足
Figure 977525DEST_PATH_IMAGE051
，说明
Figure 413186DEST_PATH_IMAGE052
，这时在所述垂线上确定一点b’，使得
Figure 653674DEST_PATH_IMAGE053
和
Figure 869892DEST_PATH_IMAGE054
之间的夹角等于
Figure 955659DEST_PATH_IMAGE055
，将目标顶点b调整到b’的位置，由调整后的目标顶点b’和目标顶点之外的其它三个顶点a、c和d构成第四网格ab’cd，并将第四网格作为后续操作的目标网格。
第三种可选实施方式，存在多个顶点不满足内角要求，这时从所述多个顶点中任选一顶点作为目标顶点，按照第二种可选实施方式提供的方法对目标顶点进行处理，得到一个第四网格，并将所述第四网格作为后续操作的目标网格；重新执行S230，直到最终的目标网格中每个顶点均满足内角要求。实际应用中，在内角要求设置合理的情况下，例如，
Figure 776985DEST_PATH_IMAGE056
，每个顶点调整的次数不会很多。
S240、如果所述目标网格的每个顶点均满足最小高度要求，进入S250；如果所述目标网格的至少一个顶点不满足最小高度要求，确定所述至少一个顶点中任一目标顶点的最小高度对应的对边；沿所述目标顶点到所述对边的垂线调整所述目标顶点的位置，得到第五网格，使调整后的目标顶点满足最小高度要求；如果所述第五网格的每个顶点均满足内角要求，将所述第五网格作为后续操作的目标网格，进入S250；如果所述第五网格的至少一个顶点不满足内角要求，根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第六网格，将所述第六网格作为后续操作的目标网格；重新执行S240，直到所述目标网格的每个顶点均满足内角要求和高度要求。
进入S240后，目标网格满足翘曲要求且每个顶点均满足内角要求，S240完成目标网格的最小高度调整。首先计算目标网格每个顶点的最小高度。如图17所示，分别计算目标顶点b到两条对边
Figure 985112DEST_PATH_IMAGE057
的垂直距离h2和h3，取所述两个垂直距离中的较小值作为目标顶点b的最小高度。图17中h2，作为目标顶点b的最小高度。
得到四个顶点的最小高度后，判断每个顶点的最小高度是否满足最小高度要求。四边形网格的最小高度要求为顶点的最小高度大于第七阈值H2，及判断每个顶点的最小高度是否满足
Figure 891888DEST_PATH_IMAGE058
。基于不同的判断结果，S240包括以下几种可选实施方式：
第一种可选实施方式，每个顶点均满足高度要求，这时无需进行高度调整，直接进入S250进行横纵比调整。
第二种可选实施方式，存在一个顶点不满足高度要求，这时以该顶点为目标顶点，确定所述目标顶点的最小高度对应的对边。以图17为例，目标顶点b不满足最小高度要求，确定b的最小高度h2对应的对边
Figure 297460DEST_PATH_IMAGE059
。确定该对边后，沿所述目标顶点到所述对边的垂线调整所述目标顶点的位置，得到第五网格。具体的，在所述垂线上确定上确定第八点，使所述第八点到所述对边的垂直距离等于所述第七阈值；将所述目标顶点调整到所述第八点的位置，由所述调整后的目标顶点，以及目标网格中除目标顶点之外的其它三个顶点构成第五网格。继续以图17为例，在目标顶点b到对边
Figure 973292DEST_PATH_IMAGE060
的垂线上确定第八点b’，b’到
Figure 352321DEST_PATH_IMAGE060
的垂直距离=H2，将b调整到b’的位置，由调整后的目标顶点b’和目标顶点之外的其它三个顶点a、c和c构成第五网格ab’cd所述调整后的目标顶点在所述第五网格中一定是满足高度要求的。得到第五网格后，目标顶点的内角、以及目标顶点两个相邻顶点的内角会发生变化，除此之外的另外一个内角则保持不变。由于四边形网格中四个内角之间的耦合性相比于三角形较弱，当一个内角发生变化时，仍有很大概率四个顶点均满足内角要求。因此本实施例首先判断所述第五网格的每个顶点是否均内角要求，包括以下两种情况：
情况一、第五网格的每个顶点均满足内角要求，则将所述第五网格作为后续操作的目标网格，进入S250。
情况二、第五网格的至少一个顶点不满足内角要求，这时根据调整后的目标顶点对所述目标网格进行等比例放大，得到第六网格，将所述第六网格作为后续操作的目标网格。等比例放大后，第六网格的四条边与目标网格的四条边是平行的，因此四个内角保持不变，一定满足内角要求。
可选的，所述根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第六网格，包括：确定过所述调整后的目标顶点的两条直线，使所述两条直线分别平行于所述目标网格中过所述目标顶点的两条边；分别在所述两条直线上确定两点，使得所述两点和所述调整后的目标定点构成的两条新边与所述两条边的长度比值相同；根据所述两条新边，生成第六网格。以图18为例，确定过b’的两条直线，分别平行于目标网格abcd中过目标顶点b的两条边
Figure 746393DEST_PATH_IMAGE061
和
Figure 173964DEST_PATH_IMAGE062
；分别在这两条直线上确定两点a’和c’，使得所述两点a’和c’和所述调整后的目标定点b’构成的两条新边
Figure 766619DEST_PATH_IMAGE063
和
Figure 254232DEST_PATH_IMAGE064
与所述两条边
Figure 870021DEST_PATH_IMAGE065
和
Figure 366862DEST_PATH_IMAGE062
的长度比例相同，即
Figure 814023DEST_PATH_IMAGE066
。确定
Figure 738117DEST_PATH_IMAGE067
的过c’点的平行线，确定
Figure 74158DEST_PATH_IMAGE068
的过a’点的平行线，两平行线的交点和a’、c’构成另外两条新边。由所述两条新边和所述另外两条新边公共构成第六网格。第六网格
Figure 905848DEST_PATH_IMAGE069
是目标网格abcd的等比例放大，与目标网格abcd中四个内角的角度相等，一定满足内角要求。
第三种可选实施方式、存在多个顶点不满足高度要求，这时从所述多个顶点中的任选一顶点作为目标顶点，按照第二种可选实施方式提供的方法对目标顶点进行处理，得到一个第六网格，并将所述第六网格作为后续操作的目标网格；重新执行S240，直到最终的目标网格中每个顶点均满足内角要求和最小高度要求。
S250、如果所述目标网格不满足横纵比要求，调整所述目标网格最短边和最长边的交点位置，使所述目标网格满足横纵比要求。
进入S250后，目标网格满足翘曲要求且每个顶点均满足内角要求和高度要求，S250完成目标网格的横纵比要求。首先计算目标网格的横纵比，判断四边形网格是否满足横纵比要求。四边形网格的横纵比要求为横纵比小于或等于第八阈值。基于不同的判断结果，包括以下两种情况：
情况一、四边形网格满足横纵比要求，这时无需进行横纵比调整，整个四边形网格的质量调整结束。
情况二、四边形网格不满足横纵比要求。由于经过翘曲调整、内角调整和最小高度调整后，四边形单元中的最短边和最长边一定是相邻的（如图19所示），则通过调整最短边和最长边的交点位置，使目标网格满足横纵比要求。具体的，沿所述目标网格最短边所在的直线，调整所述目标网格的最长边和所述最短边的交点位置，以增加所述最短边的长度，得到第七网格。由于调整所述交点位置可能会影响内角角度，因此判断第七网格中的每个顶点是否均满足内角要求。基于不同的判断结果，S250包括以下几种实施方式：
第一种实施方式，第七网格的每个顶点均满足内角要求，将第七网格为横纵比调整的结果，整个四边形网格的质量调整结束。
第二种实施方式，第七网格的存在至少一个顶点不满足内角要求，则是通过调整最长边的长度来改变横纵比。可选的，根据所述横纵比要求在所述最长边上确定第三点；确定所述最短边的过所述第三点的平行线；根据所述平行线，调整所述最短边的位置。具体的，根据横纵比要求，计算在最短边不变的情况下最长边需要减小的长度。如图19所示，如果最长边
Figure 676358DEST_PATH_IMAGE068
减小l就可以满足横纵比要求，则在最长边
Figure 505774DEST_PATH_IMAGE070
上确定第三点d’，并确定最短边
Figure 96155DEST_PATH_IMAGE071
的过d’的平行线，根据所述平行线调整最短边的位置。
可选的，所述根据所述平行线，调整所述最短边的位置，包括以下两种情况：
情况一、所述平行线与所述最长边的对边存在交点。当四边形网格的翘曲角度接近于0时，即四边形网格中只存在一个平面内时，所述平行线和最长边的对边一定存在交点。这时由所述交点和所述第三点构成调整后的最短边。以图19为例，平行线与最长边
Figure 731536DEST_PATH_IMAGE068
的对边
Figure 356552DEST_PATH_IMAGE062
存在交点c’，由交点c’和第三点d’构成调整后的最短边
Figure 356869DEST_PATH_IMAGE072
。从图19中可以看出，情况一的处理方法相当于将最短边
Figure 434547DEST_PATH_IMAGE072
沿最长边
Figure 608039DEST_PATH_IMAGE070
及其对边
Figure 821983DEST_PATH_IMAGE062
平行移动，因此四边形的每个内角保持不变，最短边
Figure 757315DEST_PATH_IMAGE073
和最长边
Figure 322289DEST_PATH_IMAGE070
的交点c的最小高度保持不变，则整个目标网格的最小高度也可保持不变。
情况二、所述平行线与所述最长边的对边不存在交点。当四边形网格的翘曲角度不接近于0时，即四边形网格内存在两个平面时，所述平行线和最长边的对边不一定存在交点。这时确定所述平行线上距离所述最长边的对边最近的第四点，由所述第三点和所述第四点构成调整后的最短边。仍以图19为例，平行线与最长边
Figure 237155DEST_PATH_IMAGE068
的对边
Figure 571185DEST_PATH_IMAGE062
不存在交点，则确定所述平行线上距离
Figure 382146DEST_PATH_IMAGE062
最近的第四点c’，由第三点d’和第四点c’构成调整后的最短边
Figure 434415DEST_PATH_IMAGE072
。可以看出，该平行线相当于将最短边
Figure 917087DEST_PATH_IMAGE072
沿最长边
Figure 902361DEST_PATH_IMAGE070
及其对边
Figure 680961DEST_PATH_IMAGE062
平行移动的结果，由于平行线与
Figure 954947DEST_PATH_IMAGE062
无交点，无法直接作为调整后的四边形的一条边，因此将平线性向
Figure 477196DEST_PATH_IMAGE062
投影，投影点c’即为距离
Figure 582555DEST_PATH_IMAGE062
最近的第四点。
S250完成后，最终的目标网格满足翘曲要求和横纵比要求，且每个顶点均满足内角要求和最小高度要求，四边形网格的质量调整结束。
本实施例按照翘曲角度、内角角度、最小高度和横纵比的顺序对四边形网格进行质量调整。在翘曲调整中，通过调整对角线两侧顶点的位置实现两个平面朝相反方向的旋转，以减小翘曲角度。在内角调整和最小高度调整过程中将目标顶点的两个相邻顶点构成的对角线保持不变，通过目标顶点到对角线的垂线限定目标顶点的调整轨迹，充分利用垂线与四个顶点的几何关系修改目标顶点的位置，保证调整后的目标顶点符合相应的质量要求，避免了人为调整中随机拖动顶点位置而无法保证调整后的顶点一定满足质量要求的问题。在进行高度调整时，考虑到四边形内角之间的耦合性，增大目标顶点的最小高度后判断是否满足内角要求，当不满足内角要求时对目标网格进行等比例放大以保持各内角不变，避免了在内角调整和高度调整之间的反复跳跃，提高了调整效率。在进行横纵比调整时，通过调整最短边和最长边的交点增加最短边的长度来调整横纵比；同时考虑到对内角的影响，在内角不满足要求时放弃增加最短边的长度，通过最短边的平行线实现最短边在最长边及其对边之间的平行移动，从而减小最长边的长度来调整横纵比，并能保持内角不变，最小高度也满足高度要求。整个调整策略易于实现，且能够保证调整结果满足质量要求，为网格质量调整提供了一种自动化的实现方式。
图20是本发明实施例提供另一种有限元网格质量自动调整方法，适用于对汽车的板壳结构零部件的有限元网格自动进行质量调整的情况。该方法可以由电子设备执行，参见图20，该方法具体包括：
S310、从CAE软件获取车辆零部件的有限元网格特征；
S320、将所述有限元网格特征中的每个三角形网格作为目标网格，对所述目标网格执行上述任一实施例提供的S110-S130的操作；同时，
S330、将所述有限元网格特征中的每个四边形网格作为目标网格，对所述目标网格执行上述任一实施例提供的S210-S250的操作。
本实施例基于上述任一实施例提供的方法实现，具备上述任一实施例可实现的技术效果。
图21为本发明实施例提供的一种电子设备的结构示意图，如图21所示，该设备包括处理器40、存储器41、输入装置42和输出装置43；设备中处理器40的数量可以是一个或多个，图21中以一个处理器40为例；设备中的处理器40、存储器41、输入装置42和输出装置43可以通过总线或其他方式连接，图21中以通过总线连接为例。
存储器41作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的有限元网格质量自动调整方法对应的程序指令/模块。处理器40通过运行存储在存储器41中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述的有限元网格质量自动调整方法。
存储器41可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储器41可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器41可进一步包括相对于处理器40远程设置的存储器，这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置42可用于接收输入的数字或字符信息，以及产生与设备的用户设置以及功能控制有关的键信号输入。输出装置43可包括显示屏等显示设备。
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现任一实施例的有限元网格质量自动调整方法。
本发明实施例的计算机存储介质，可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子（非穷举的列表）包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器（RAM）、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言—诸如”C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络，包括局域网(LAN)或广域网(WAN)，连接到用户计算机，或者，可以连接到外部计算机（例如利用因特网服务提供商来通过因特网连接）。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (10)

1.一种有限元网格质量自动调整方法，其特征在于，包括：
S110、从CAE软件获取车辆零部件有限元网格特征，将所述有限元网格特征中质量不合格的三角形网格作为目标网格；
S120、如果所述目标网格的每个顶点均满足内角要求，进入S130；如果所述目标网格的至少一个顶点不满足内角要求，沿所述至少一个顶点中任一目标顶点到对边的垂线调整所述目标顶点的位置，得到第一网格，使调整后的目标顶点在所述第一网格中满足内角要求，并将所述第一网格作为后续操作的目标网格；重新执行S120，直到所述目标网格的每个顶点均满足内角要求；
S130、如果所述目标网格的至少一个顶点不满足高度要求，沿所述至少一个顶点中任一目标顶点到对边的垂线调整所述目标顶点的位置，使调整后的目标顶点满足高度要求；根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第二网格，并将所述第二网格作为后续操作的目标网格；重新执行S130，直到所述目标网格的每个顶点均满足内角要求和高度要求。
2.根据权利要求1所述的方法，其特征在于，所述内角要求为内角角度大于或等于第一阈值且小于或等于第二阈值；
所述沿所述至少一个顶点中任一目标顶点到对边的垂线调整所述目标顶点的位置，得到第一网格，使调整后的目标顶点在所述第一网格中满足内角要求，包括：
如果所述至少一个顶点中任一目标顶点的内角小于所述第一阈值，在所述目标顶点到对边的垂线上确定第一点，使所述第一点与所述目标顶点的两个相邻顶点的连线夹角等于所述第一阈值，并将所述目标顶点调整到所述第一点的位置；
如果所述至少一个顶点中任一目标顶点的内角大于所述第二阈值，在所述目标顶点到对边的垂线上确定第二点，使所述第二点与所述目标顶点的两个相邻顶点的连线夹角等于所述第二阈值，并将所述目标顶点调整到所述第二点的位置。
3.根据权利要求1所述的方法，其特征在于，所述根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第二网格，包括：
确定过所述调整后的目标顶点的两条直线，使所述两条直线分别平行于所述目标网格中过所述目标顶点的两条边；
分别确定所述对边的延长线与所述两条直线的两个交点，由所述调整后的目标顶点和所述两个交点构成第二网格。
4.一种有限元网格质量自动调整方法，其特征在于，包括：
S210、从CAE软件获取车辆零部件的有限元网格特征，将所述有限元网格特征中质量不合格的四边形网格作为目标网格；
S220、如果所述目标网格满足翘曲要求，进入S230；如果所述目标网格不满足翘曲要求，调整所述目标网格对角线两侧顶点的位置，得到第三网格，使所述第三网格满足翘曲要求，并将所述第三网格作为后续操作的目标网格；
S230、如果所述目标网格的每个顶点均满足内角要求，进入S240；如果所述目标网格的至少一个顶点不满足内角要求，将所述至少一个顶点中的任一顶点作为目标顶点，并确定所述目标网格的两个相邻点构成的对角线；沿所述目标顶点到所述对角线的垂线调整所述目标顶点的位置，得到第四网格，使调整后的目标顶点在所述第四网格中满足内角要求，并将所述第四网格作为后续操作的目标网格；重新执行S230，直到所述目标网格的每个顶点均满足内角要求；
S240、如果所述目标网格的每个顶点均满足最小高度要求，进入S250；如果所述目标网格的至少一个顶点不满足最小高度要求，确定所述至少一个顶点中任一目标顶点的最小高度对应的对边；沿所述目标顶点到所述对边的垂线调整所述目标顶点的位置，得到第五网格，使调整后的目标顶点满足最小高度要求；如果所述第五网格的每个顶点均满足内角要求，将所述第五网格作为后续操作的目标网格，进入S250；如果所述第五网格的至少一个顶点不满足内角要求，根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第六网格，将所述第六网格作为后续操作的目标网格；重新执行S240，直到所述目标网格的每个顶点均满足内角要求和高度要求；
S250、如果所述目标网格不满足横纵比要求，调整所述目标网格最短边和最长边的交点位置，使所述目标网格满足横纵比要求。
5.根据权利要求4所述的方法，其特征在于，所述根据所述调整后的目标顶点对所述目标网格进行等比例放大，得到第六网格，包括：
确定过所述调整后的目标顶点的两条直线，使所述两条直线分别平行于所述目标网格中过所述目标顶点的两条边；
分别在所述两条直线上确定两点，使得所述两点和所述调整后的目标定点构成的两条新边与所述两条边的长度比值相同；
根据所述两条新边，生成第六网格。
6.根据权利要求4所述的方法，其特征在于，所述调整所述目标网格最短边和最长边的交点位置，使所述目标网格满足横纵比要求，包括：
沿所述目标网格最短边所在的直线，调整所述目标网格的最长边和所述最短边的交点位置，以增加所述最短边的长度，得到第七网格；
如果所述第七网格中存在至少一个顶点不满足内角要求，根据所述横纵比要求在所述最长边上确定第三点；
确定所述最短边的过所述第三点的平行线；
根据所述平行线，调整所述最短边的位置。
7.根据权利要求6所述的方法，其特征在于，所述根据所述平行线，调整所述最短边的位置，包括：
如果所述平行线与所述最长边的对边存在交点，由所述交点和所述第三点构成调整后的最短边；
如果所述平行线与所述最长边的对边不存在交点，确定所述平行线上距离所述最长边的对边最近的第四点，由所述第三点和所述第四点构成调整后的最短边。
8.一种有限元网格质量自动调整方法，其特征在于：
从CAE软件获取车辆零部件的有限元网格特征；
将所述有限元网格特征中的每个质量不合格的三角形网格作为目标网格，对所述目标网格执行如权利要求1-3中任一所述的S110-S130的操作；同时，
将所述有限元网格特征中的每个质量不合格的四边形网格作为目标网格，对所述目标网格执行如权利要求4-7中任一所述的S210-S250的操作。
9.一种电子设备，其特征在于，包括：
一个或多个处理器；
存储器，用于存储一个或多个程序，
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-8中任一所述的有限元网格质量自动调整方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-8中任一所述的有限元网格质量自动调整方法。
CN202210682906.0A 2022-06-17 2022-06-17 有限元网格质量自动调整方法、设备和存储介质 Active CN114781232B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210682906.0A CN114781232B (zh) 2022-06-17 2022-06-17 有限元网格质量自动调整方法、设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210682906.0A CN114781232B (zh) 2022-06-17 2022-06-17 有限元网格质量自动调整方法、设备和存储介质

Publications (2)

Publication Number Publication Date
CN114781232A true CN114781232A (zh) 2022-07-22
CN114781232B CN114781232B (zh) 2022-09-16

Family

ID=82421423

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210682906.0A Active CN114781232B (zh) 2022-06-17 2022-06-17 有限元网格质量自动调整方法、设备和存储介质

Country Status (1)

Country Link
CN (1) CN114781232B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060044309A1 (en) * 2004-08-31 2006-03-02 Satoshi Kanai Tetrahedral mesh generating method for finite-element analysis and finite-element analyzing system using its method
CN108537797A (zh) * 2018-04-09 2018-09-14 哈尔滨理工大学 一种虚拟手术中可变形对象切割仿真的网格变形优化方法
CN113920275A (zh) * 2021-09-30 2022-01-11 广州极飞科技股份有限公司 三角网格构建方法、装置、电子设备及可读存储介质

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060044309A1 (en) * 2004-08-31 2006-03-02 Satoshi Kanai Tetrahedral mesh generating method for finite-element analysis and finite-element analyzing system using its method
CN108537797A (zh) * 2018-04-09 2018-09-14 哈尔滨理工大学 一种虚拟手术中可变形对象切割仿真的网格变形优化方法
CN113920275A (zh) * 2021-09-30 2022-01-11 广州极飞科技股份有限公司 三角网格构建方法、装置、电子设备及可读存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李宗领等: "平面区域有限元三角网格迭代优化方法", 《新型工业化》 *

Also Published As

Publication number Publication date
CN114781232B (zh) 2022-09-16

Similar Documents

Publication Publication Date Title
US11676335B2 (en) Tessellation method using vertex tessellation factors
CN115294258B (zh) 三维模型的展开方法、装置、设备及计算机可读存储介质
US20150002510A1 (en) Sketch-based generation and editing of quad meshes
CN110689492B (zh) 图像边缘平滑方法和装置
US20120029882A1 (en) System and method for generating quadrangulations
US11465361B2 (en) Systems, methods, and media for controlling support structures and build orientation in manufacturing
CN112632818B (zh) 高阶梯度光滑无网格方法、系统、介质、计算机设备
CN114494648B (zh) 基于有限元网格划分的网格调整方法、设备和存储介质
CN114972687B (zh) 基于消除三角形网格对的网格调整方法
US20240221278A1 (en) Ordering of patch selection in tessellation operations
US20240062469A1 (en) Data processing method, apparatus, device, and medium
CN115329628A (zh) 一种基于局部网格自适应剖分的电磁辐射仿真分析方法
CN113797531B (zh) 遮挡剔除实现方法、装置、计算机设备及存储介质
CN114781232B (zh) 有限元网格质量自动调整方法、设备和存储介质
CN111914364A (zh) 基于高阶矩量法与投影的频选天线罩建模方法
CN115081122A (zh) 板壳结构几何中面的自动抽取方法、设备及介质
CN111784796A (zh) 一种地形网格生成方法、装置、设备和介质
CN114842448B (zh) 三维车道线生成方法、装置、电子设备和计算机可读介质
CN117115383B (zh) 照明灯的布置方法、计算机设备和可读存储介质
US8847971B2 (en) Graphic processing device, graphic processing method, and non-transitory medium for storing graphic processing program for generating corrected graphic data
WO2023185917A1 (zh) 一种芯片的模块排布方法及相关设备
CN114037713A (zh) 跨屏绊线勾画方法、装置、电子设备及存储介质
CN118052957A (zh) 一种针对qem算法在三维模型中保持细节特征的优化方法
CN118037992A (zh) 一种三角网局部重构方法、系统、设备及存储介质
CN117011432A (zh) 一种蒙皮动画轻量化方法、装置和存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant