CN106777630A - 用于机械零件的模态有限元仿真精度的确定方法 - Google Patents
用于机械零件的模态有限元仿真精度的确定方法 Download PDFInfo
- Publication number
- CN106777630A CN106777630A CN201611114544.6A CN201611114544A CN106777630A CN 106777630 A CN106777630 A CN 106777630A CN 201611114544 A CN201611114544 A CN 201611114544A CN 106777630 A CN106777630 A CN 106777630A
- Authority
- CN
- China
- Prior art keywords
- finite element
- machine components
- element simulation
- mode
- determination method
- 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.)
- Pending
Links
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/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种用于机械零件的模态有限元仿真精度的确定方法,包括下列步骤:基于机械零件的材料,制备材料样条并进行样条实验,得到材料的多组弹塑性数据;运用试验设计方法对材料的多组弹塑性数据进行分类选择,得到多组材料参数;应用多组材料参数对机械零件进行相应的模态有限元仿真,得到与每组材料参数相应的仿真结果;根据模态有限元仿真中建立的有限元模型,进行相同边界条件的物理实验,得到物理实验的实验结果;对仿真结果和实验结果进行正态分布分析,根据仿真结果和实验结果的正态分布均值进行计算,得到模态有限元仿真精度。与现有技术相比,本发明具有仿真精度定义准确以及接近实际情况等优点。
Description
技术领域
本发明涉及机械有限元仿真领域,尤其是涉及一种用于机械零件的模态有限元仿真精度的确定方法。
背景技术
“工业4.0”作为应对科技变革与工业格局重塑的重要工业国策,其实现的关键一环在于设计与制造的智能化。而基于计算机模型的数字化设计能力和虚拟仿真能力是智能化设计与制造的基础和保证。数字化设计技术作为一种集成化产品设计和开发手段,在高性能计算环境下运用分析模型代替实际物理系统,使得产品设计和制造过程具有智能、高效、安全、受环境约束少等优点,大大提升了新技术研发能力,降低了新产品开发成本和开发周期。在汽车、航天、航空、武器系统以及车辆等战略性现代工程领域的研究中,需要深度解析和构建具有领域特征的复杂系统以支持设计与开发,而有效地运用数字化设计技术对复杂产品的性能稳定性、质量可靠性、以及环境适应性等指标的保障和提升有巨大促进作用。
然而,虚拟仿真技术只是真实物理系统的一个有限的映像,在其建模过程中的各类抽象、前提、假设、有限的输入参数,以及因物理现象理解有限而可能出现的模型选择错误,通常很难完全从计算机仿真模型得到和物理实验结果完全一样的预测结果。另外,连续区间里的离散解也可能导致数值误差及模型求解的收敛问题。所以在使用数学或者计算机仿真模型的时候,仿真模型预测结果的可信度是虚拟样机技术的一个关键问题。
虚拟仿真模型验证的主要目标是评估数字化仿真模型在拟定应用下对实际物理实验结果的预测能力,同时结合物理实验来量化仿真模型的误差。验证过程通常需要通过对比仿真模型输出与对应的物理实验测量结果,并基于对比结果对仿真模型一致性做出决策。然而现有的验证结果往往是基于单次的仿真结果与实验结果进行比对,比对结果不准确,无法对仿真精度有一个精确的定义。
发明内容
本发明的目的是针对上述问题提供一种用于机械零件的模态有限元仿真精度的确定方法。
本发明的目的可以通过以下技术方案来实现:
一种用于机械零件的模态有限元仿真精度的确定方法,所述方法包括下列步骤:
1)基于机械零件的材料,制备材料样条并进行样条实验,得到材料的多组弹塑性数据;
2)运用试验设计方法对步骤1)得到的材料的多组弹塑性数据进行分类选择,得到多组材料参数;
3)应用步骤2)得到的多组材料参数对机械零件进行相应的模态有限元仿真,得到与每组材料参数相应的仿真结果;
4)根据步骤3)中模态有限元仿真中建立的有限元模型,进行相同边界条件的物理实验,得到物理实验的实验结果;
5)对步骤3)得到的仿真结果和步骤4)得到的实验结果进行正态分布分析,根据仿真结果和实验结果的正态分布均值进行计算,得到模态有限元仿真精度。
所述模态有限元仿真精度具体为:
所述步骤2)具体为:
21)剔除材料的多组弹塑性数据中的噪音点,得到有效的多组弹塑性数据;
22)对步骤21)中得到的有效的多组弹塑性数据进行正态分布处理;
23)根据试验设计方法中的数据减缩规则,对正态分布处理后的数据进行减缩处理,得到多组材料参数。
所述步骤3)具体为:
31)建立机械零件的有限元模型;
32)对机械零件的有限元模型进行前模态分析计算,得到该有限元模型的前六阶模态。
所述步骤32)具体为:
321)将机械零件的有限元模型定义为粘性阻尼单自由度状态下的系统,得到系统受到的阻尼力fd:
其中,fd为有系统受到的阻尼力,c为系统阻尼,为系统振动速度;
322)建立该有限元模型的振动方程:
其中,m为系统质量,为系统加速度,x(t)为系统位移,k为系统的刚度,f(t)为时间t下的系统载荷;
323)定义该有限元模型为忽略阻尼的自由振动并进行简谐运动,得到振动方程:
(k-ω2m)·x(t)=0
对振动方程进行求解得到ω1,ω2,…,ω6即为该有限元模型的前六阶模态。
所述有限元模型包括正六面体网格。
所述样条实验具体为:对材料样条进行低速拉伸。
所述弹塑性数据包括弹性模量和应力应变曲线。
所述边界条件包括施加的激励条件、施加的激励方式和梁端固定边界。
所述材料的弹塑性数据的组数在20-30组之间。
与现有技术相比,本发明具有以下有益效果:
(1)利用仿真结果和实验结果的正态分布均值计算仿真精度,与传统的单次仿真和实验进行结果比对相比,可以对仿真精度有准确的定义,同时由于定义的仿真精度是基于正态分布的均值,因而结果准确,有参考价值。
(2)制备材料样条并进行样条实验,得到材料的多组弹塑性数据,考虑了机械零件的材料的不确定性分布,使得最终得到的模态有限元仿真精度更加接近真实情况。
(3)运用试验设计方法对于得到的材料的多组弹塑性数据进行了选择,剔除了其中的噪音点,并对有效的多组弹塑性数据进行正态分布处理和减缩处理,既提高了仿真的精度又避免了进入仿真的数据过多。
(4)在有限元仿真过程中,建立的有限元模型为正六面体网格,既保证了建模的立体性也避免了建模复杂而导致的计算量过大,实用性强。
(5)通过对材料样条进行低速拉伸完成样条实验从而得到样条的弹塑性数据,由于拉伸为低速,保证了得到的弹塑性数据准确可靠。
(6)物理实验和仿真采用相同的边界条件,保证了仿真和物理实验的一致性,从而保证了仿真精度定义的准确性。
附图说明
图1为本发明的方法流程图;
图2为试验设计方法的相关示意图,其中(2a)为原始的多组弹塑性数据,(2b)为对有效的多组弹塑性数据进行正态分布处理后的示意图,(2c)为对正态分布处理后的数据进行减缩处理的示意图;
图3为仿真结果和实验结果的正态分布图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本实施例提供了一种用于简易悬梁臂的模态有限元仿真精度的确定方法,包括下列步骤:
1)确定悬梁臂的材料,制备材料样条并进行样条实验,得到多组材料的弹性模量和应力应变曲线等弹塑性数据,本实施例中,根据拉伸实验国家标准,采用低速拉伸条件,最终实施了20组拉伸实验;
2)如图2所示,运用试验设计方法(DOE)对步骤1)得到的材料的多组弹塑性数据进行分类选择,得到多组材料参数,具体为:
21)剔除材料的多组弹塑性数据中的噪音点,得到有效的多组弹塑性数据;
22)对步骤21)中得到的有效的多组弹塑性数据进行正态分布处理;
23)根据试验设计方法中的数据减缩规则,对正态分布处理后的数据进行减缩处理,得到多组材料参数,本实施例中,减缩后的材料参数为6组;
3)应用步骤2)得到的6组材料参数对机械零件进行相应的模态有限元仿真,得到与每组材料参数相应的仿真结果,即得到与6组材料参数所对应的悬臂梁前六阶模态频率值:
31)建立机械零件的有限元模型;
32)对机械零件的有限元模型进行前模态分析计算,得到该有限元模型的前六阶模态:
321)将机械零件的有限元模型定义为粘性阻尼单自由度状态下的系统,得到系统受到的阻尼力fd:
其中,fd为有系统受到的阻尼力,c为系统阻尼,为系统振动速度;
322)建立该有限元模型的振动方程:
其中,m为系统质量,为系统加速度,x(t)为系统位移,k为系统的刚度,f(t)为时间t下的系统载荷;
323)定义该有限元模型为忽略阻尼的自由振动并进行简谐运动,得到振动方程:
(k-ω2m)·x(t)=0
对振动方程进行求解得到ω1,ω2,…,ω6即为该有限元模型的前六阶模态。
4)根据步骤3)中模态有限元仿真中建立的有限元模型,建立相同材料的悬臂梁试验台,进行相同边界条件(如相同的激励条件、相同的激励方式和相同的梁端固定边界)的物理实验,得到物理实验的实验结果;
5)对步骤3)得到的仿真结果和步骤4)得到的实验结果进行正态分布分析,如图3所示,根据仿真结果和实验结果的正态分布均值进行计算,得到模态有限元仿真精度,具体为:
Claims (10)
1.一种用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述方法包括下列步骤:
1)基于机械零件的材料,制备材料样条并进行样条实验,得到材料的多组弹塑性数据;
2)运用试验设计方法对步骤1)得到的材料的多组弹塑性数据进行分类选择,得到多组材料参数;
3)应用步骤2)得到的多组材料参数对机械零件进行相应的模态有限元仿真,得到与每组材料参数相应的仿真结果;
4)根据步骤3)中模态有限元仿真中建立的有限元模型,进行相同边界条件的物理实验,得到物理实验的实验结果;
5)对步骤3)得到的仿真结果和步骤4)得到的实验结果进行正态分布分析,根据仿真结果和实验结果的正态分布均值进行计算,得到模态有限元仿真精度。
2.根据权利要求1所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述模态有限元仿真精度具体为:
仿真精度=100%-仿真偏差
3.根据权利要求1所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述步骤2)具体为:
21)剔除材料的多组弹塑性数据中的噪音点,得到有效的多组弹塑性数据;
22)对步骤21)中得到的有效的多组弹塑性数据进行正态分布处理;
23)根据试验设计方法中的数据减缩规则,对正态分布处理后的数据进行减缩处理,得到多组材料参数。
4.根据权利要求1所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述步骤3)具体为:
31)建立机械零件的有限元模型;
32)对机械零件的有限元模型进行前模态分析计算,得到该有限元模型的前六阶模态。
5.根据权利要求4所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述步骤32)具体为:
321)将机械零件的有限元模型定义为粘性阻尼单自由度状态下的系统,得到系统受到的阻尼力fd:
其中,fd为有系统受到的阻尼力,c为系统阻尼,为系统振动速度;
322)建立该有限元模型的振动方程:
其中,m为系统质量,为系统加速度,x(t)为系统位移,k为系统的刚度,f(t)为时间t下的系统载荷;
323)定义该有限元模型为忽略阻尼的自由振动并进行简谐运动,得到振动方程:
(k-ω2m)·x(t)=0
对振动方程进行求解得到ω1,ω2,…,ω6即为该有限元模型的前六阶模态。
6.根据权利要求4所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述有限元模型包括正六面体网格。
7.根据权利要求1所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述样条实验具体为:对材料样条进行低速拉伸。
8.根据权利要求1所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述弹塑性数据包括弹性模量和应力应变曲线。
9.根据权利要求1所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述边界条件包括施加的激励条件、施加的激励方式和梁端固定边界。
10.根据权利要求1所述的用于机械零件的模态有限元仿真精度的确定方法,其特征在于,所述材料的弹塑性数据的组数在20-30组之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611114544.6A CN106777630A (zh) | 2016-12-07 | 2016-12-07 | 用于机械零件的模态有限元仿真精度的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611114544.6A CN106777630A (zh) | 2016-12-07 | 2016-12-07 | 用于机械零件的模态有限元仿真精度的确定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106777630A true CN106777630A (zh) | 2017-05-31 |
Family
ID=58874767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611114544.6A Pending CN106777630A (zh) | 2016-12-07 | 2016-12-07 | 用于机械零件的模态有限元仿真精度的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777630A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108089457A (zh) * | 2017-11-29 | 2018-05-29 | 北京航空航天大学 | 一种基于在线有限元仿真的过程质量控制方法 |
CN109902326A (zh) * | 2018-06-14 | 2019-06-18 | 中国煤炭科工集团太原研究院有限公司 | 一种有限元仿真实验效果测评方法 |
CN113361162A (zh) * | 2021-06-03 | 2021-09-07 | 郑州大学 | 一种计算撞振模型节点位移的方法、装置 |
CN116011124A (zh) * | 2023-03-20 | 2023-04-25 | 西安航天动力研究所 | 一种振动试验系统的仿真建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004823A (zh) * | 2010-11-11 | 2011-04-06 | 浙江中科电声研发中心 | 一种扬声器振动和声学特性的数值模拟方法 |
CN103279657A (zh) * | 2013-05-21 | 2013-09-04 | 北京航空航天大学 | 一种基于工程经验的产品加速退化试验方案设计方法 |
CN105447228A (zh) * | 2015-11-11 | 2016-03-30 | 中国人民解放军国防科学技术大学 | 一种地球扰动引力赋值模型适用性能评估方法 |
CN105825025A (zh) * | 2016-03-30 | 2016-08-03 | 清华大学 | 板状零件加工变形仿真预测方法 |
-
2016
- 2016-12-07 CN CN201611114544.6A patent/CN106777630A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102004823A (zh) * | 2010-11-11 | 2011-04-06 | 浙江中科电声研发中心 | 一种扬声器振动和声学特性的数值模拟方法 |
CN103279657A (zh) * | 2013-05-21 | 2013-09-04 | 北京航空航天大学 | 一种基于工程经验的产品加速退化试验方案设计方法 |
CN105447228A (zh) * | 2015-11-11 | 2016-03-30 | 中国人民解放军国防科学技术大学 | 一种地球扰动引力赋值模型适用性能评估方法 |
CN105825025A (zh) * | 2016-03-30 | 2016-08-03 | 清华大学 | 板状零件加工变形仿真预测方法 |
Non-Patent Citations (3)
Title |
---|
朱跃: ""基于分层思想的复杂机械结构模型修正及确认技术研究"", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
赵士亮: ""粗糙集理论中若干问题的研究及应用"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
钟伟斌: ""基于有限元分析的轮毂轴承单元轴铆合工艺优化"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108089457A (zh) * | 2017-11-29 | 2018-05-29 | 北京航空航天大学 | 一种基于在线有限元仿真的过程质量控制方法 |
CN109902326A (zh) * | 2018-06-14 | 2019-06-18 | 中国煤炭科工集团太原研究院有限公司 | 一种有限元仿真实验效果测评方法 |
CN109902326B (zh) * | 2018-06-14 | 2023-04-07 | 中国煤炭科工集团太原研究院有限公司 | 一种有限元仿真实验效果测评方法 |
CN113361162A (zh) * | 2021-06-03 | 2021-09-07 | 郑州大学 | 一种计算撞振模型节点位移的方法、装置 |
CN113361162B (zh) * | 2021-06-03 | 2022-11-22 | 郑州大学 | 一种计算撞振模型节点位移的方法、装置 |
CN116011124A (zh) * | 2023-03-20 | 2023-04-25 | 西安航天动力研究所 | 一种振动试验系统的仿真建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | Optimizing fluid–structure interaction systems with immersogeometric analysis and surrogate modeling: Application to a hydraulic arresting gear | |
Vahdati et al. | Mechanisms for wide-chord fan blade flutter | |
CN103150446B (zh) | 临近空间飞行器模型验证与仿真实现方法 | |
CN106777630A (zh) | 用于机械零件的模态有限元仿真精度的确定方法 | |
Biancolini et al. | Static aeroelastic analysis of an aircraft wind-tunnel model by means of modal RBF mesh updating | |
CN105260581A (zh) | 舰船机电控制设备虚拟振动和冲击试验方法 | |
CN105678015B (zh) | 一种高超声速三维机翼的非概率可靠性气动结构耦合优化设计方法 | |
CN106934136A (zh) | 一种细长型飞行器的振动试验仿真系统及方法 | |
CN115438598A (zh) | 基于一般时间根方尺度的雷诺应力湍流模型的数值方法 | |
Gao et al. | On the dispersion mechanism of the flutter boundary of the AGARD 445.6 wing | |
Chen et al. | Towards high-accuracy deep learning inference of compressible turbulent flows over aerofoils | |
Zangeneh | Development of a new algorithm for modeling viscous transonic flow on unstructured grids at high Reynolds numbers | |
CN102629283A (zh) | 一种转动部件对挠性动力学影响的仿真分析方法 | |
Yabili et al. | Unsteady navier-stokes simulation of low-reynolds stall flutter | |
US11080446B2 (en) | Method to regulate clock frequencies of hybrid electronic systems | |
CN106529086A (zh) | 用于机械结构件的静载荷强度有限元仿真精度的确定方法 | |
CN115618501A (zh) | 基于数据融合修正的亚跨气动特性获取方法、系统及装置 | |
Bohbot et al. | Computation of the flutter boundary of an airfoil with a parallel Navier-Stokes solver | |
CN112214843B (zh) | 一种用于风洞试验机翼模型有限元刚度修正方法及装置 | |
Imamura et al. | Immersed boundary method for high Reynolds number compressible flows around an aircraft configuration | |
Duan et al. | A high order overset FR/CPR method for dynamic moving grids | |
Sugaya et al. | Unsteady turbulent flow simulations on moving Cartesian grids using immersed boundary method and high-order scheme | |
Miyaji et al. | On Accuracy of Prediction of Flutter Boundaries on Unstructured Grids | |
Silva | A unified approach for computational aeroelasticity | |
JP7331877B2 (ja) | システム計算装置、及びシステム計算プログラム |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |
|
RJ01 | Rejection of invention patent application after publication |