CN107563071A - 一种综合环境应力的仿真方法 - Google Patents
一种综合环境应力的仿真方法 Download PDFInfo
- Publication number
- CN107563071A CN107563071A CN201710799881.1A CN201710799881A CN107563071A CN 107563071 A CN107563071 A CN 107563071A CN 201710799881 A CN201710799881 A CN 201710799881A CN 107563071 A CN107563071 A CN 107563071A
- Authority
- CN
- China
- Prior art keywords
- result
- stress
- analysis
- random vibration
- model
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及有限元仿真技术,其公开了一种综合环境应力的仿真方法,能够分析在综合环境下的仿真模型的应力分布,使得仿真结果更符合实际情况。该方法包括:a.创建几何模型,定义不同温度下的材料参数,切分模型并划分网格,生成有限元模型;b.进行热‑结构耦合的预应力分析:即施加温度及压力载荷,得到高温所产生的热应力和压力载荷产生的静应力结果的综合应力,作为后续随机振动分析的预应力;c.将所述预应力加载到有限元模型中,进行带预应力的模态分析;d.施加振动激励,进行随机振动分析,得到随机振动应力结果;e.将随机振动应力结果和高温所产生的热应力、压力载荷产生的静应力结果相叠加。本发明适用于综合环境应力仿真分析。
Description
技术领域
本发明涉及有限元仿真技术,具体涉及一种综合环境应力的仿真方法。
背景技术
随着科学技术的迅猛发展,仿真已成为各种复杂系统研制工作的一种必不可少的手段。由于很多产品都工作在高温、高压、振动等综合环境中,其可靠性的研究也愈发的重要。通过仿真,我们可以缩短研发周期、改进生产过程、降低成本并且更为真实的了解到产品受外界环境的影响,其中应力仿真分析就是极其重要的一部分。
目前,许多研究在对结构体进行应力仿真分析时,往往在高温或高压等单一环境下进行,然而结构体通常处于高温、高压、振动等综合环境之下,比单一环境更为复杂,同时在综合环境之下的各种因素会相互影响,为了使仿真结果更符合实际,得出更为真实的应力分布结果,必须找到更加准确的仿真方法。
ANSYS是一个大型通用的有限元仿真软件,其拥有强大的实体建模技术以及对各种物理场的处理分析能力。通过ANSYS软件对仿真模型进行分析,能够得到模型上的应力分布。然而ANSYS软件并未提供在综合环境下的应力仿真,这使得对于处在综合环境下的结构体而言,其仿真结果不符合实际,制约了我们对其可靠性的研究。
为了得出仿真模型在综合环境下的应力分布,使仿真分析结果更符合实际,针对ANSYS研究综合应力仿真方法是必要的。
发明内容
本发明所要解决的技术问题是:提出一种综合环境应力的仿真方法,能够分析在综合环境下的仿真模型的应力分布,使得仿真结果更符合实际情况。
本发明解决上述技术问题所采用的技术方案是:
一种综合环境应力的仿真方法,包括以下步骤:
a.创建几何模型,定义不同温度下的材料参数,切分模型并划分网格,生成有限元模型;
b.进行热-结构耦合的预应力分析:即施加温度及压力载荷,得到高温所产生的热应力和压力载荷产生的静应力结果的综合应力,作为后续随机振动分析的预应力;
c.将所述预应力加载到有限元模型中,进行带预应力的模态分析;
d.施加振动激励,进行随机振动分析,得到随机振动应力结果;
e.将随机振动应力结果和高温所产生的热应力、压力载荷产生的静应力结果相叠加,从而得到综合环境下的应力仿真结果。
作为进一步优化,步骤a中,所述材料参数包括:材料密度、热传导系数和比热容。
作为进一步优化,步骤a中,在模型划分网格时采用六面体网格,以此减少单元数量,加快收敛速度,提高分析的精度,减少数值误差。
作为进一步优化,步骤e中,采用载荷工况的组合技术对随机振动应力结果与高温所产生的热应力以及压力载荷产生的静应力结果进行叠加,具体包括:
e1.将得到的随机振动分析结果文件读入数据库,将载荷工况1指向随机振动分析结果的应力载荷步上;
e2.创建载荷工况1的工况文件,并将后处理设置为缺省状态,则工况文件1的内容即为随机振动应力结果;
e3.将得到的高温和压力载荷形成的静应力结果文件读入数据库,将载荷工况2指向静应力结果文件中的第一个载荷步上;
e4.将载荷工况文件1的工况文件读入数据库,将结果写入静应力结果文件的第二个载荷步;
e5.将当前数据库的结果与载荷工况2的结果进行叠加,即得到综合环境应力的结果,并将叠加的结果写入第三个载荷步中进行保存。
本发明的有益效果是:
分析了高温状况下的热应力、随机振动应力、压力载荷产生的静应力,并采用载荷工况的组合技术对此三个应力进行叠加,从而获得综合环境下的应力分析结果,使得仿真结果更符合实际情况。
附图说明
图1为本发明实施例中的综合环境应力的仿真方法流程图;
图2为ANSYS仿真模型的结构示意图;
图3为施加压力载荷得到的应力分布图;
图4为施加温度载荷得到的应力分布图;
图5为综合环境下的应力分布图。
具体实施方式
本发明旨在提出一种综合环境应力的仿真方法,能够分析在综合环境下的仿真模型的应力分布,使得仿真结果更符合实际情况。
下面结合附图及实施例对本发明的方案作进一步的描述:
如图1所示,本实施例中的综合环境应力的仿真方法包括以下实现步骤:
S01:生成有限元模型:
本步骤中,通过创建几何模型,定义不同温度下的材料参数,切分模型并划分网格,生成有限元模型,其中,所述材料参数包括:材料密度、热传导系数和比热容。在模型划分网格时采用六面体网格,以此减少单元数量,加快收敛速度,提高分析的精度,减少数值误差。
S02:进行热-结构耦合的预应力分析:
本步骤中,施加温度及压力载荷,得到高温所产生的热应力和压力载荷产生的静应力结果的综合应力,作为后续随机振动分析的预应力。
S03:将所述预应力加载到有限元模型中,进行带预应力的模态分析:
本步骤中,将步骤S02获得的预应力(即高温所产生的热应力和压力载荷产生的静应力结果的综合应力)加载到有限元模型中,进行带预应力的模态分析。
S04:施加振动激励,进行随机振动分析,得到随机振动应力结果;
S05:将随机振动应力结果和高温所产生的热应力、压力载荷产生的静应力结果相叠加,从而得到综合环境下的应力仿真结果:
本步骤中,可以采用载荷工况的组合技术对随机振动应力结果与高温所产生的热应力以及压力载荷产生的静应力结果进行叠加,具体包括:
(1)将得到的随机振动分析结果文件读入数据库,将载荷工况1指向随机振动分析结果的应力载荷步上;
(2)创建载荷工况1的工况文件,并将后处理设置为缺省状态,则工况文件1的内容即为随机振动应力结果;
(3)将得到的高温和压力载荷形成的静应力结果文件读入数据库,将载荷工况2指向静应力结果文件中的第一个载荷步上;
(4)将载荷工况文件1的工况文件读入数据库,将结果写入静应力结果文件的第二个载荷步;
(5)将当前数据库的结果与载荷工况2的结果进行叠加,即得到综合环境应力的结果,并将叠加的结果写入第三个载荷步中进行保存。
实施例:
在本实施例中,采用ANSYS软件分别进行单一环境应力仿真分析和综合环境应力仿真分析,仿真模型结构如图2所示,在该模型基础上,施加2MPa的压力载荷,通过仿真分析获得最大应力为27MPa,其应力分布图如图3所示,为单一的压力环境下的应力分布;
在仿真模型结构图2的基础上,施加600度的温度载荷时,通过仿真分析获得最大应力为57MPa,其应力分布图如图4所示,为单一的温度环境下的应力分布。
在仿真模型结构图2的基础上,同时施加2MPa的压力载荷和、600度的温度载荷以及振动载荷时,通过仿真分析获得最大应力为58MPa,其应力分布如图5所示,为压力、温度以及振动综合环境下的应力分布。
可以看出,对于处在综合环境下的结构体而言,利用单一环境去模拟实际的工作状态,以此得出的应力分布结果显然不符合实际情况,并且,在综合环境下的应力分布也并非是各个单一环境下应力分布结果的简单叠加,这是由于综合环境下的各种因素相互影响,只有综合考虑影响应力仿真的各个方面,仿真结果才更接近实际情况。
Claims (4)
1.一种综合环境应力的仿真方法,其特征在于,包括以下步骤:
a.创建几何模型,定义不同温度下的材料参数,切分模型并划分网格,生成有限元模型;
b.进行热-结构耦合的预应力分析:即施加温度及压力载荷,得到高温所产生的热应力和压力载荷产生的静应力结果的综合应力,作为后续随机振动分析的预应力;
c.将所述预应力加载到有限元模型中,进行带预应力的模态分析;
d.施加振动激励,进行随机振动分析,得到随机振动应力结果;
e.将随机振动应力结果和高温所产生的热应力、压力载荷产生的静应力结果相叠加,从而得到综合环境下的应力仿真结果。
2.如权利要求1所述的一种综合环境应力的仿真方法,其特征在于,步骤a中,所述材料参数包括:材料密度、热传导系数和比热容。
3.如权利要求1所述的一种综合环境应力的仿真方法,其特征在于,步骤a中,在模型划分网格时采用六面体网格,以此减少单元数量,加快收敛速度,提高分析的精度,减少数值误差。
4.如权利要求1所述的一种综合环境应力的仿真方法,其特征在于,步骤e中,采用载荷工况的组合技术对随机振动应力结果与高温所产生的热应力以及压力载荷产生的静应力结果进行叠加,具体包括:
e1.将得到的随机振动分析结果文件读入数据库,将载荷工况1指向随机振动分析结果的应力载荷步上;
e2.创建载荷工况1的工况文件,并将后处理设置为缺省状态,则工况文件1的内容即为随机振动应力结果;
e3.将得到的高温和压力载荷形成的静应力结果文件读入数据库,将载荷工况2指向静应力结果文件中的第一个载荷步上;
e4.将载荷工况文件1的工况文件读入数据库,将结果写入静应力结果文件的第二个载荷步;
e5.将当前数据库的结果与载荷工况2的结果进行叠加,即得到综合环境应力的结果,并将叠加的结果写入第三个载荷步中进行保存。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710799881.1A CN107563071A (zh) | 2017-09-07 | 2017-09-07 | 一种综合环境应力的仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710799881.1A CN107563071A (zh) | 2017-09-07 | 2017-09-07 | 一种综合环境应力的仿真方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107563071A true CN107563071A (zh) | 2018-01-09 |
Family
ID=60979504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710799881.1A Pending CN107563071A (zh) | 2017-09-07 | 2017-09-07 | 一种综合环境应力的仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107563071A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109101700A (zh) * | 2018-07-20 | 2018-12-28 | 电子科技大学 | 一种变间隙式结构在仿真中设置接触对的方法 |
CN109190279A (zh) * | 2018-09-18 | 2019-01-11 | 中国人民解放军海军航空大学 | 一种温度振动加速耦合效应模型的构建方法 |
CN110309617A (zh) * | 2019-07-10 | 2019-10-08 | 电子科技大学 | 传感器上温度、压力、振动三种载荷直接耦合的仿真方法 |
CN112861280A (zh) * | 2021-01-29 | 2021-05-28 | 三一石油智能装备有限公司 | 一种泵头体设计方法、加工方法及压裂泵 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105912825A (zh) * | 2016-05-25 | 2016-08-31 | 济南大学 | 一种基于ansys的叶轮结构热应力影响下的模态分析方法 |
-
2017
- 2017-09-07 CN CN201710799881.1A patent/CN107563071A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105912825A (zh) * | 2016-05-25 | 2016-08-31 | 济南大学 | 一种基于ansys的叶轮结构热应力影响下的模态分析方法 |
Non-Patent Citations (3)
Title |
---|
吕晓东: "火焰筒浮动壁的随机振动特性研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
奚浩晨: "伺服机构热-力耦合仿真与电子舱热仿真分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
王艳荣 等: "基于有限元仿真的机电组件的可靠性分析", 《2015年全国机械行业可靠性技术学术交流会暨第五届可靠性工程分会第二次全体委员大会论文集》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109101700A (zh) * | 2018-07-20 | 2018-12-28 | 电子科技大学 | 一种变间隙式结构在仿真中设置接触对的方法 |
CN109190279A (zh) * | 2018-09-18 | 2019-01-11 | 中国人民解放军海军航空大学 | 一种温度振动加速耦合效应模型的构建方法 |
CN110309617A (zh) * | 2019-07-10 | 2019-10-08 | 电子科技大学 | 传感器上温度、压力、振动三种载荷直接耦合的仿真方法 |
CN110309617B (zh) * | 2019-07-10 | 2022-05-13 | 电子科技大学 | 传感器上温度、压力、振动三种载荷直接耦合的仿真方法 |
CN112861280A (zh) * | 2021-01-29 | 2021-05-28 | 三一石油智能装备有限公司 | 一种泵头体设计方法、加工方法及压裂泵 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107563071A (zh) | 一种综合环境应力的仿真方法 | |
CN103559418B (zh) | 一种基于故障物理的电子产品任务可靠度计算方法 | |
CN106485032A (zh) | 一种计算叶盘振动可靠性的双重极值响应面法 | |
CN106873987A (zh) | 一种基于Modelica的多人协同建模系统 | |
CN108663948A (zh) | 一种航空发动机控制系统数值仿真分析平台的设计方法 | |
CN103942091A (zh) | Matlab自定义模型和psasp联合仿真的励磁系统仿真方法及系统 | |
CN104809273A (zh) | 一种描述蠕变变形的方法 | |
CN109684743A (zh) | 基于ansys脚本语言的航空电子产品振动云仿真自动接口方法 | |
CN104376157A (zh) | 大型汽轮发电机组轴系动静态特性分析计算系统 | |
Lu et al. | Modeling and simulation of aircraft assembly line based on quest | |
Blahoš et al. | Parallel harmonic balance method for analysis of nonlinear dynamical systems | |
CN107291992B (zh) | 一种适用沙漠地区电子装备综合环境试验仿真系统及方法 | |
CN102279900B (zh) | 小型涡轮发动机涡轮虚拟试验系统 | |
Nely et al. | Analysis of the capabilities of software products to simulate the behavior of dynamic fluid flows | |
CN106777630A (zh) | 用于机械零件的模态有限元仿真精度的确定方法 | |
CN105677995A (zh) | 一种基于全网格配点理论的模糊稳态热传导问题数值求解方法 | |
CN112861374A (zh) | 基于预控制器的多物理耦合仿真处理方法、装置和设备 | |
CN110334459B (zh) | 一种输电塔塔线体系快速精细化建模系统及方法 | |
CN106372359A (zh) | 一种设备虚拟样机系统的通用架构 | |
CN108846211A (zh) | 一种基于Matlab的直流控制保护功能模块建模的方法和装置 | |
CN104484509B (zh) | 穿孔式高增益光电系统级封装结构的热疲劳优化设计方法 | |
CN107727350A (zh) | 微纳卫星矢量振动试验方法 | |
Hai Linh et al. | Digital Hardware Realization of Forward and Inverse Kinematics for a Five‐Axis Articulated Robot Arm | |
CN104268320A (zh) | 适用于卫星敏感载荷的新型隔振器刚度矩阵快速预估方法 | |
Lagloire et al. | Single platform integration environment for turbine rotor design and analysis |
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: 20180109 |
|
RJ01 | Rejection of invention patent application after publication |