CN103530463B - 一种有限元三维模型网格快速过渡的方法 - Google Patents
一种有限元三维模型网格快速过渡的方法 Download PDFInfo
- Publication number
- CN103530463B CN103530463B CN201310480241.6A CN201310480241A CN103530463B CN 103530463 B CN103530463 B CN 103530463B CN 201310480241 A CN201310480241 A CN 201310480241A CN 103530463 B CN103530463 B CN 103530463B
- Authority
- CN
- China
- Prior art keywords
- several
- section
- divided
- plane
- grid
- 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.)
- Active
Links
Abstract
一种有限元三维模型网格快速过渡的方法,它涉及一种有限元模型网格过渡的方法,以解决目前网格过渡的方法仅限于单方向的过渡手段,即在单方向上采用远端与近端网格尺寸比为2或者3的过渡比例进行多次网格过渡,该手段不能实现大型构件的有限元仿真的问题。方法:一、根据构件不同部位计算数据变化梯度的大小划分为重点区域、过渡区域和非重点区域;二、在重点区域内划分数个立方体或长方体的网格,重点区域A的网格应密集;三、对过渡区域进行过渡网格划分;四、在非重点区域内划分数个立方体或长方体的网格,非重点区域的网格应稀疏,至此整个构件的网格划分完毕。本发明用于大型或超大型构件有限元模拟仿真建模的网格划分。
Description
技术领域:
本发明涉及一种有限元模型网格过渡的方法,具体涉及一种实现有限元三维模型网格快速过渡的方法。
背景技术:
有限元方法是随着电子计算机的发展而迅速发展起来的一种求解数学、物理问题的现代数值计算方法,它可以分析求解静态或动态的力学性能、热传导、电磁场、流体力学等线性或非线性问题。有限元方法用较简单的问题代替复杂问题后再求解,将求解域看成是由许多称为有限元的小的互连子域组成,对小单元进行分析。一般用有限元分析软件来计算实际工程问题,在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器、国防军工、船舶、铁道、石化、能源、科学研究等各个领域的得到广泛使用。用有限元分析方法来代替前期的实验,将这些过程由计算机处理,可以节约大量的人力和物力,缩短研究周期,还可以减少实际工作反复性和盲目性,尤其解决大型或超大型复杂构件无法实现的试验过程,有限元方法更显得具有经济意义和现实意义。
要实现大型构件的有限元模拟仿真,首先解决的就是有限元网格的建立问题,构件的尺寸越大,需要的有限元网格数量就越多,则计算量越大。此外在有限元分析中,网格划分合适与否与计算结果的精度和计算效率息息相关。网格划分越细,计算精度越高,所花费的计算时间越长;反之,计算精度越低,所花费的计算时间越短。而且,网格划分细到一定程度,计算精度变化较小,甚至不发生变化。在计算数据变化梯度较大或关键区域(如应力集中部位),为了较好地反映数据变化规律,需要采用比较密集的小网格。而在计算数据变化梯度较小或对整体影响不大的区域(如应力变化平缓部位),为减小模型规模,则应划分相对稀疏的大网格。这样可以保证计算解精确可靠,又可以提高计算效率。两个区域之间选择映射网格划分方法,即网格过渡。网格过渡应相对平稳,否则,将影响计算结果的准确性甚至计算结果不收敛致使使有限元无法计算。由于建模时既要保证网格的尺寸足够小,又要尽量控制网格数量,不同区域的网格大小差异数倍或几个数量级,这样将造成对过渡区域网格的质量要求较高,给建模带来极大的困难。
目前网格过渡的方法仅限于单方向的过渡,即在单方向上采用远端与近端网格尺寸比为2或3的过渡比进行多次网格过渡,如图6所示。单方向的过渡方法不仅过渡比小,而且操作复杂,无法实现大型构件的有限元仿真过程。
发明内容:
本发明的目的是为解决目前网格过渡的方法仅限于单方向的过渡方法,即在单方向上采用远端与近端网格尺寸比为2或者3的过渡比进行多次网格过渡,单方向的过渡方法不仅过渡比小,而且操作复杂,很难实现大型构件的有限元仿真的问题,而提供的一种有限元三维模型网格快速过渡方法。
本发明的一种有限元三维模型网格快速过渡方法是通过以下步骤实现的:
步骤一、首先对大型复杂构件的三维立体模型进行不同区域的划分:
根据构件不同部位计算数据变化梯度的大小划分为重点区域、过渡区域和非重点区域;
步骤二、对重点区域进行网格划分:
在重点区域内划分数个立方体或长方体的网格,重点区域的网格应密集;
步骤三、对过渡区域进行过渡网格划分:
过渡区域环绕在重点区域的下面、左侧和右侧,过渡区域的截面形状呈凹形,凹形上端面的中心点与凹形下端面左端点之间的斜面为左平分斜面,凹形上端面的中心点至凹形下端面右端点之间的斜面为右平分斜面,左平分斜面与右平分斜面之间的区域内设有由中心点至凹形下端面的数个下斜面,数个下斜面将左平分斜面与右平分斜面之间的夹角等分为数个锐角,凹形上端面与左平分斜面之间的区域内设有数个左斜面,数个左斜面将凹形上端面与左平分斜面之间的夹角等分为数个锐角,凹形上端面与右平分斜面之间的区域内设有数个右斜面,数个右斜面将凹形上端面与右平分斜面之间的夹角等分为数个锐角;
在过渡区域上由内向外划分为数个单元,数个单元从截面看,由内向外依次为第一单元凹形截面、第二单元凹形截面、第三单元凹形截面、第四单元凹形截面、第五单元凹形截面、第六单元凹形截面、第七单元凹形截面和第八单元凹形截面,第一单元凹形截面、第二单元凹形截面、第三单元凹形截面、第四单元凹形截面、第五单元凹形截面、第六单元凹形截面、第七单元凹形截面和第八单元凹形截面、数个下斜面、数个左斜面和数个右斜面将过渡区域划分为数个立体梯形网格,数个立方体或长方体网格由内向外、由小到大逐渐过渡;
步骤四、对非重点区域进行网格划分:
非重点区域环绕在过渡区域的外围,非重点区域的截面形状呈凹形,在非重点区域内划分数个立方体或长方体的网格,非重点区域的网格应稀疏,至此整个构件的网格划分完毕;
各个维度方向上的网格密度不相同,且网格密度由内向外、由小到大逐渐过渡。
本发明的优点是:
一、利用本发明所提供的有限元三维模型网格快速过渡方法,大大提高了计算效率,可完成大型或超大型复杂构件的仿真过程,甚至可以解决目前常规过渡方法因网格数量过多而计算结果不收敛致使有限元无法计算的难题。本发明所提供网格过渡方法较常规的方法可大大节约网格总数,划分效率与过渡比均大大提高。而且当构件尺寸越大,加工总长度越长时,利用本发明所提供的网格划分方法的效率越显著。
二、本发明提供的一种实现有限元网格模型快速三维过渡的方法,针对大型或超大型构件有限元模拟仿真建模的网格划分,实现了超大过渡比的网格过渡。
附图说明:
图1为步骤一中对大型复杂构件划分为重点区域A、过渡区域B和非重点区域C划分示意图;
图2为步骤二中对重点区域A进行网格划分示意图;
图3为步骤三中对过渡区域B进行网格划分示意图;
图4为图3的I局部放大图;
图5为利用本发明的有限元三维模型网格快速过渡的方法对整个构件的重点区域A、过渡区域B和非重点区域C进行网格划分完成的整个构件示意图;
图6为采用常规单方向网格过渡方法进行网格划分示意图。
具体实施方式:
具体实施方式一:结合图1~图5说明本实施方式,本实施方式是通过以下步骤实现的:
步骤一、首先对大型复杂构件的三维立体模型进行不同区域的划分:
根据构件不同部位计算数据变化梯度的大小划分为重点区域A、过渡区域B和非重点区域C;
步骤二、对重点区域A进行网格划分:
在重点区域A内划分数个立方体或长方体的网格,重点区域A的网格应密集;重点区域A为计算数据变化梯度较大或关键区域(如应力集中部位),为了较好地反映数据变化规律,有限元计算对该区域的网格质量要求较高,采用比较密集的网格,所以要着重单独划分重点区域A;
步骤三、对过渡区域B进行过渡网格划分:
过渡区域B环绕在重点区域A的下面、左侧和右侧,过渡区域B的截面形状呈凹形,凹形上端面的中心点O与凹形下端面左端点M之间的斜面为左平分斜面O-M,凹形上端面的中心点O至凹形下端面右端点N之间的斜面为右平分斜面O-N,左平分斜面O-M与右平分斜面O-N之间的区域内设有由中心点O至凹形下端面的数个下斜面O-S,数个下斜面O-S将左平分斜面O-M与右平分斜面O-N之间的夹角等分为数个锐角,凹形上端面与左平分斜面O-M之间的区域内设有数个左斜面O-U,数个左斜面O-U将凹形上端面与左平分斜面O-M之间的夹角等分为数个锐角,凹形上端面与右平分斜面O-N之间的区域内设有数个右斜面O-V,数个右斜面O-V将凹形上端面与右平分斜面O-N之间的夹角等分为数个锐角;
在过渡区域B上由内向外划分为数个单元,数个单元从截面看,由内向外依次为第一单元凹形截面1、第二单元凹形截面2、第三单元凹形截面3、第四单元凹形截面4、第五单元凹形截面5、第六单元凹形截面6、第七单元凹形截面7和第八单元凹形截面8,第一单元凹形截面1、第二单元凹形截面2、第三单元凹形截面3、第四单元凹形截面4、第五单元凹形截面5、第六单元凹形截面6、第七单元凹形截面7和第八单元凹形截面8、数个下斜面O-S、数个左斜面O-U和数个右斜面O-V将过渡区域B划分为数个立体梯形网格,数个立方体或长方体网格由内向外、由小到大逐渐过渡;
步骤四、对非重点区域C进行网格划分:
非重点区域C环绕在过渡区域B的外围,非重点区域C的截面形状呈凹形,在非重点区域C内划分数个立方体或长方体的网格,非重点区域C的网格应稀疏,至此整个构件的网格划分完毕。非重点区域C为计算数据变化梯度较小或对整体影响不大的区域(如应力变化平缓部位),为了最大程度的提高计算效率,减小模型规模,有限元仿真对此区域的网格质量要求较低,则非重点区域C应划分相对稀疏的网格;
各个维度方向上的网格密度不相同,且网格密度由内向外、由小到大逐渐过渡。
位于比较密集的小网格和比较稀疏的大网格之间的过渡网格(即过渡区域B),用于保证重点区域A和非重点区域C的网格保持匹配。
三种区域内的各自均有实体单元,区域之间的单元完全匹配。最终构件的有限元网格模型如图5所示。
具体实施方式二:结合图3说明本实施方式,本实施方式是步骤三中左平分斜面O-M与凹形下端面之间的锐角为45°,右平分斜面O-N与凹形下端面之间的锐角为45°。其它步骤与具体实施方式一相同。
具体实施方式三:结合图4说明本实施方式,本实施方式是步骤三中过渡区域B内侧的网格尺寸与重点区域A的正方形网格尺寸相匹配。其它步骤与具体实施方式一或二相同。
具体实施方式四:结合图4说明本实施方式,本实施方式是步骤四中非重点区域C内侧的网格尺寸与过渡区域B外侧的网格尺寸相匹配。其它步骤与具体实施方式三相同。
应用实例:以一个长宽高尺寸为600mm×1000mm×600mm的大构件的网格划分为例,图1为构件三维造型图和不同区域的划分结果。重点区域A的立方体小网格尺寸为2.5mm×2.5mm×2.5mm,B为过渡区域,非重点区域C的立方体大网格尺寸为40mm×40mm×40mm。结果证明,整个网格划分过程中使用本发明的方法得到的网格数量为29248个,而通过常规单方向网格过渡方法得到的网格数量为928000个,所以本发明所提供网格过渡方法较常规方法可大大节约网格总数,效率提高到96%以上,过渡比达到3375以上。
Claims (4)
1.一种有限元三维模型网格快速过渡方法,其特征在于:所述方法是通过以下步骤实现的:
步骤一、首先对大型复杂构件的三维立体模型进行不同区域的划分:
根据构件不同部位计算数据变化梯度的大小划分为重点区域(A)、过渡区域(B)和非重点区域(C);
步骤二、对重点区域(A)进行网格划分:
在重点区域(A)内划分数个立方体或长方体的网格,重点区域(A)的网格应密集;
步骤三、对过渡区域(B)进行过渡网格划分:
过渡区域(B)环绕在重点区域(A)的下面、左侧和右侧,过渡区域(B)的截面形状呈凹形,凹形上端面的中心点(O)与凹形下端面左端点(M)之间的斜面为左平分斜面(O-M),凹形上端面的中心点(O)至凹形下端面右端点(N)之间的斜面为右平分斜面(O-N),左平分斜面(O-M)与右平分斜面(O-N)之间的区域内设有由中心点(O)至凹形下端面的数个下斜面(O-S),数个下斜面(O-S)将左平分斜面(O-M)与右平分斜面(O-N)之间的夹角等分为数个锐角,凹形上端面与左平分斜面(O-M)之间的区域内设有数个左斜面(O-U),数个左斜面(O-U)将凹形上端面与左平分斜面(O-M)之间的夹角等分为数个锐角,凹形上端面与右平分斜面(O-N)之间的区域内设有数个右斜面(O-V),数个右斜面(O-V)将凹形上端面与右平分斜面(O-N)之间的夹角等分为数个锐角;
在过渡区域(B)上由内向外划分为数个单元,数个单元从截面看,由内向外依次为第一单元凹形截面(1)、第二单元凹形截面(2)、第三单元凹形截面(3)、第四单元凹形截面(4)、第五单元凹形截面(5)、第六单元凹形截面(6)、第七单元凹形截面(7)和第八单元凹形截面(8),第一单元凹形截面(1)、第二单元凹形截面(2)、第三单元凹形截面(3)、第四单元凹形截面(4)、第五单元凹形截面(5)、第六单元凹形截面(6)、第七单元凹形截面(7)和第八单元凹形截面(8)、数个下斜面(O-S)、数个左斜面(O-U)和数个右斜面(O-V)将过渡区域(B)划分为数个立体梯形网格,数个立方体或长方体网格由内向外、由小到大逐渐过渡;
步骤四、对非重点区域(C)进行网格划分:
非重点区域(C)环绕在过渡区域(B)的外围,非重点区域(C)的截面形状呈凹形,在非重点区域(C)内划分数个立方体或长方体的网格,非重点区域(C)的网格应稀疏,至此整个构件的网格划分完毕;
各个维度方向上的网格密度不相同,且网格密度由内向外、由小到大逐渐过渡。
2.根据权利要求1所述一种有限元三维模型网格快速过渡方法,其特征在于:步骤三中左平分斜面(O-M)与凹形下端面之间的锐角为45°,右平分斜面(O-N)与凹形下端面之间的锐角为45°。
3.根据权利要求1或2所述一种有限元三维模型网格快速过渡方法,其特征在于:步骤三中过渡区域(B)内侧的网格尺寸与重点区域(A)的正方形网格尺寸相匹配。
4.根据权利要求3所述一种有限元三维模型网格快速过渡方法,其特征在于:步骤四中非重点区域(C)内侧的网格尺寸与过渡区域(B)外侧的网格尺寸相匹配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310480241.6A CN103530463B (zh) | 2013-10-05 | 2013-10-05 | 一种有限元三维模型网格快速过渡的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310480241.6A CN103530463B (zh) | 2013-10-05 | 2013-10-05 | 一种有限元三维模型网格快速过渡的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103530463A CN103530463A (zh) | 2014-01-22 |
CN103530463B true CN103530463B (zh) | 2017-01-11 |
Family
ID=49932470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310480241.6A Active CN103530463B (zh) | 2013-10-05 | 2013-10-05 | 一种有限元三维模型网格快速过渡的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103530463B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106649903A (zh) * | 2015-11-02 | 2017-05-10 | 株洲时代新材料科技股份有限公司 | 一种货车扭杆式翻转机构有限元模型的混合网格划分方法 |
CN108351906B (zh) * | 2015-11-25 | 2022-05-17 | 西门子工业软件有限公司 | 用于具有晶格结构的部件的建模的系统和方法 |
CN106557638B (zh) * | 2016-12-07 | 2019-11-05 | 重庆科技学院 | 焊接混合双向过渡单元网格模型的建立方法 |
CN111931457B (zh) * | 2020-09-27 | 2021-01-15 | 北京智芯仿真科技有限公司 | 基于混合阶有限元的多层集成电路电磁场计算方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034745A (zh) * | 2011-09-30 | 2013-04-10 | 鞍钢股份有限公司 | 一种模拟型钢切深孔轧制过程的网格生成方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6794766B2 (en) * | 2001-06-29 | 2004-09-21 | General Electric Company | Method and operational strategy for controlling variable stator vanes of a gas turbine power generator compressor component during under-frequency events |
-
2013
- 2013-10-05 CN CN201310480241.6A patent/CN103530463B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034745A (zh) * | 2011-09-30 | 2013-04-10 | 鞍钢股份有限公司 | 一种模拟型钢切深孔轧制过程的网格生成方法 |
Non-Patent Citations (1)
Title |
---|
基于几何分解的六面体网格生成算法的研究;王东风;《中国优秀硕士学位论文全文数据库 信息科技辑》;20111215(第S2期);第40-53页、图0.10-0.12 * |
Also Published As
Publication number | Publication date |
---|---|
CN103530463A (zh) | 2014-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103530463B (zh) | 一种有限元三维模型网格快速过渡的方法 | |
CN104036095B (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
Abedi et al. | An h-adaptive spacetime-discontinuous Galerkin method for linear elastodynamics | |
CN103559418A (zh) | 一种基于故障物理的电子产品任务可靠度计算方法 | |
CN103366085A (zh) | 编织复合材料力学性能的多尺度预测方法 | |
CN105260581A (zh) | 舰船机电控制设备虚拟振动和冲击试验方法 | |
CN106442271A (zh) | 岩心渗透率模拟方法及装置 | |
Peng et al. | A node split method for crack growth problem | |
Xiang et al. | GPU acceleration of CFD algorithm: HSMAC and SIMPLE | |
Jaworski et al. | Composite material microlevel cellular model data transfer and verification by fem analysis software systems | |
CN110362912A (zh) | 介观结构优化方法 | |
CN109388833A (zh) | 一种基于疲劳寿命的弹性元件结构优化设计方法 | |
Akkurt et al. | An efficient edge based data structure for the compressible Reynolds‐averaged Navier–Stokes equations on hybrid unstructured meshes | |
Wang et al. | An Approach to Predicting Fatigue Crack Growth Under Mixed-Mode Loading Based on Improved Gaussian Process | |
CN112685936B (zh) | 一种用于贝壳珍珠母微结构有限元分析的建模方法 | |
Wada et al. | Effective adaptation technique for hexahedral mesh | |
Hsiao et al. | Development of compressible-incompressible link to efficiently model bubble dynamics near floating body | |
Sobotka et al. | Automatic 3D crack placement using the Python API in ABAQUS CAE | |
CN108197398A (zh) | 一种预测基于空间群p4的三维编织复合材料失效的有限元方法 | |
Jamroz et al. | Asynchronous communication in spectral-element and discontinuous Galerkin methods for atmospheric dynamics–a case study using the High-Order Methods Modeling Environment (HOMME-homme_dg_branch) | |
CN117709129B (zh) | 多尺度仿真方法、装置、设备及存储介质 | |
Azevedo | Fracture mechanics using the natural neighbour radial point interpolation method | |
CN101777089A (zh) | 地铁隔断门门扇的有限元分析及优化 | |
Dawes et al. | Acceleration of an unstructured hybrid mesh RANS solver by porting to GPU architectures | |
Deliparaschos et al. | Optimised sensor selection for control: A hardware-in-the-loop realization on FPGA for an EMS system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |