CN105631112B - 一种厚板类零件的热变形仿真分析与建模方法 - Google Patents
一种厚板类零件的热变形仿真分析与建模方法 Download PDFInfo
- Publication number
- CN105631112B CN105631112B CN201510993887.3A CN201510993887A CN105631112B CN 105631112 B CN105631112 B CN 105631112B CN 201510993887 A CN201510993887 A CN 201510993887A CN 105631112 B CN105631112 B CN 105631112B
- Authority
- CN
- China
- Prior art keywords
- thermal deformation
- thick plate
- model
- simulation analysis
- value
- 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.)
- Expired - Fee Related
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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种适用于厚板类零件的热变形仿真分析与建模的方法,在设定温度载荷条件下,利用有限元仿真分析方法对厚板类零件进行热变形分析,然后根据仿真分析结果建立厚板类零件的热变形快速计算模型。其特征在于,基于厚板类零件长度方向和厚度方向的热变形协调规律,通过提取长度方向上设定位置点的热变形值和中性面上厚度方向的热变形值,利用多项式拟合与多元线性回归的方法,建立中性面上厚度方向的热变形值与长度方向上设定位置点的热变形值之间的函数映射关系。本发明能够通过对厚板类零件的热变形仿真分析,实现通过长度方向设定位置点的热变形值预测中性面上厚度方向的热变形值的目的。
Description
技术领域
本发明涉及零件的变形仿真分析与建模领域,特别是一种适用于厚板类零件的热变形仿真分析与建模的方法。
背景技术
近年来,我国制造业的迅猛发展对零件的精度提出了越来越高的要求,人们在关注零件几何变形的同时,开始关注零件的热变形。但是目前传统的方法主要集中于对零件几何变形的分析,提出的技术方法也多针对于如何控制零件的几何变形,而针对零件热变形的分析与研究方法相对较少,因此如何对零件的热变形进行分析与建模,已成为我国制造业亟待解决的一大关键性技术问题。
发明内容
本发明为解决公知技术中存在的技术问题而提供了一种厚板类零件的热变形仿真分析与建模方法。
本发明为解决公知技术中存在的技术问题所采取的技术方案是:一种适用于厚板类零件的热变形仿真分析与建模方法,包括厚板类零件的热变形仿真分析与热变形建模方法。
所述厚板类零件的热变形仿真分析方法是基于通用有限元分析软件Ansys进行的,,包括如下步骤:
步骤一,几何模型的简化;
步骤二,添加材料属性并进行网格划分;
步骤三,载荷的施加;
步骤四,约束的定义;
步骤五,进行有限元求解;
步骤六,提取有限元分析结果,包括在有限元分析完成后,提取厚板类零件需要建立数学模型方向的热变形结果,以及垂直于需要建立热变形数学模型方向的热变形结果。
所述步骤一,几何模型的简化,包括对厚板类零件有限元分析影响较小的结构的简化。
所述步骤三,载荷的施加,包括在厚板类零件的相同部位施加不同大小的温度载荷,以及施加零件表面的传热系数和环境温度。
所述厚板类零件的热变形建模方法,包括利用多项式拟合与多元线性回归的方法,建立关于模型热变形的数学模型。其中,首先利用多项式拟合的方法,根据提取出的厚板类零件需要建立数学模型方向的热变形结果,建立相同阶次的不同加载情况下的热变形数学模型;其次,利用多元线性回归的方法,建立提取出的垂直于需要建立热变形数学模型方向的热变形结果与第一次所建立的数学模型系数之间的函数关系。
附图说明
图1为厚板类零件的热变形仿真分析与建模方法的流程示意图;
图2为一种厚板类零件的CAD模型示意图;
图中:1、模型左侧面预先选定点1;2、模型左侧面预先选定点2;3、模型左侧面预先选定点3;4、模型右侧面固定支撑;5、第一处施加温度载荷位置;6、第二处施加温度载荷位置。
图3为图2所示厚板类零件的有限元网格模型示意图;
图4为图2所示厚板类零件施加23℃以及51℃温度载荷情况下的温度场、Z向以及Y向结果云图示意图;
图中:a)施加23℃温度载荷时的温度场分布;b)施加23℃温度载荷时的Z向热变形;c)施加23℃温度载荷时的Y向热变形;d)施加51℃温度载荷时的温度场分布;e)施加51℃温度载荷时的Z向热变形;f)施加51℃温度载荷时的Y向热变形。
图5为图2所示厚板类零件所建立的Y向热变形数学模型的效果验证示意图。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
请参阅图1,一种厚板类零件的热变形仿真分析与建模方法的具体流程,包括确定厚板类零件所需建立热变形模型的方向;对厚板类零件进行有限元仿真分析;利用所述热变形仿真分析的结果,对于厚板类零件的两个相互垂直方向的热变形,通过其中一个方向的热变形仿真分析结果,利用多项式拟合及多元线性回归的方法,建立与之垂直方向热变形的数学模型。
现根据图1所示流程,对图2所示的一种厚板类零件进行热变形仿真分析与建模,步骤如下:
首先,确定图2所示的厚板类零件所需要建立热变形数学模型的方向为Y方向。
其次,进行有限元仿真分析,具体步骤如下:
步骤一,几何模型简化,去掉对有限元分析影响较小的结构;
步骤二,模型的材料属性按照QT500-7的属性,包括弹性模量1730GPa,泊松比0.3,密度7300kg/m3,热膨胀系数9.1×10-61/K,热导率47W/(m·K),添加完成后进行网格划分;
步骤三,载荷的施加,包括在第一处施加温度载荷位置与第二处施加温度载荷位置依次施加23℃、25℃、29℃、34℃、39℃、46℃、51℃的恒定温度载荷,并对不同的加载条件进行编号依次为“01加载条件”,“02加载条件”以此类推;施加表面传热系数10W/(m2·K),环境温度施加为20℃;
步骤四,约束的定义,对模型右侧面固定支撑处X、Y与Z方向的自由度进行限制,图3所示为该厚板类零件的有限元网格模型;
步骤五,进行有限元求解,图4中列出了有限元分析完成后,施加23℃以及51℃温度载荷情况下的温度场、Y向以及Z向结果云图;
步骤六,提取有限元分析结果,包括对Y、Z两个相互垂直方向的热变形的提取,为了便于描述,提取图2所示的模型左侧面上预先选定的3点的Z向热变形,以及模型上表面中线的Y向热变形。
再者,根据不同加载情况下提取的Y向热变形结果,利用多项式拟合的方法建立相同阶次的零件模型上表面中线的Y向热变形数学模型,
δyi=pi·z (1)
式中pi=(pi1 pi2 pi3),i=1,2…,7,z=(z3 z2 z)T,δyi为第“i个加载条件”下该厚板类零件上表面中线的Y向热变形,pi为表达式的系数,z为坐标值;再利用多元线性回归的方法,建立提取出的预先选定的3点的Z向热变形结果与式(1)系数之间的函数关系,
pi=A·δzi+B (2)其中,δzi=(δzi1 δzi2 δzi3)T,i=1,2…,7,δzi为提取出第“i个加载条件”下预先选定的3点的Z向热变形结果;联立式(1)与式(2),最终可以得到该零件模型上表面中线的Y向热变形数学模型,
δy=(A·δz+B)·z (3)
式中,δy为模型上表面中线的Y向热变形,δz为提取出的预先选定的3点的Z向热变形结果,z=(z3 z2 z)T,z为坐标值,A是3×3的系数矩阵,
B是3×1的常数向量,
B=[1.3806×10-10 -4.6663×10-7 3.8320×10-6]T
最后,对所建立的数学模型进行验证,按照上文有限元仿真分析步骤,在第一处及第二处载荷施加位置施加41℃的温度载荷并进行有限元分析,分析完成后将施加23℃,41℃,51℃温度载荷时提取出的预先选定的3点的Z向热变形,代入式(3)中,得到3种情况下该种厚板类模型的上表面中线Y向变形情况,并与直接从有限元分析软件中提取出的Y向变形情况进行比对,对比结果如图5所示。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护范围之内。
Claims (1)
1.一种厚板类零件的热变形仿真分析与建模方法,包括如下步骤:首先在设定温度载荷条件下,利用有限元仿真分析方法对厚板类零件进行热变形分析,有限元分析时对模型右侧面固定支撑处X、Y与Z方向的自由度进行限制,然后根据仿真分析结果建立厚板类零件的热变形快速计算模型,基于厚板类零件长度方向和厚度方向的热变形协调规律,通过提取长度方向上设定位置点的热变形值和中性面上厚度方向的热变形值;
利用多项式拟合与多元线性回归的方法,建立中性面上厚度方向的热变形值与长度方向上设定位置点的热变形值之间的函数映射关系,具体方式为:根据不同加载情况下提取的Y向热变形结果,利用多项式拟合的方法建立相同阶次的零件模型上表面中线的Y向热变形数学模型δyi=pi·z,式中pi=(pi1 pi2 pi3),i=1,2…,7,z=(z3 z2 z)T,δyi为第“i个加载条件”下该厚板类零件上表面中线的Y向热变形,pi为表达式的系数,z为坐标值;再利用多元线性回归的方法,建立提取出的预先选定的3点的Z向热变形结果与δyi=pi·z系数之间的函数关系,pi=A·δzi+B,其中,δzi=(δzi1 δzi2 δzi3)T,i=1,2…,7,δzi为提取出第“i个加载条件”下预先选定的3点的Z向热变形结果,联立δyi=pi·z和pi=A·δzi+B,得到该零件模型上表面中线的Y向热变形数学模型δy=(A·δz+B)·z,式中,δz为模型上表面中线的Y向热变形,δz为提取出的预先选定的3点的Z向热变形结果,z=(z3 z2 z)T,z为坐标值,A是3×3的系数矩阵,B是3×1的常数向量,B=[1.3806×10-10-4.6663×10-7 3.8320×10-6]T,从而实现通过长度方向设定位置点的热变形值预测中性面上厚度方向的热变形值的目的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510993887.3A CN105631112B (zh) | 2015-12-25 | 2015-12-25 | 一种厚板类零件的热变形仿真分析与建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510993887.3A CN105631112B (zh) | 2015-12-25 | 2015-12-25 | 一种厚板类零件的热变形仿真分析与建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105631112A CN105631112A (zh) | 2016-06-01 |
CN105631112B true CN105631112B (zh) | 2021-08-24 |
Family
ID=56046041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510993887.3A Expired - Fee Related CN105631112B (zh) | 2015-12-25 | 2015-12-25 | 一种厚板类零件的热变形仿真分析与建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631112B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0655597A (ja) * | 1992-08-10 | 1994-03-01 | Sekisui Chem Co Ltd | 射出成形プロセスシミュレーション方法及びその装置 |
CN101629801A (zh) * | 2009-08-18 | 2010-01-20 | 上海理工大学 | 数控磨床导轨热误差确定方法 |
CN101943896A (zh) * | 2010-07-16 | 2011-01-12 | 浙江大学 | 数控机床误差的轨迹再生补偿方法 |
CN102672527A (zh) * | 2012-05-03 | 2012-09-19 | 四川大学 | 数控机床进给系统全工作行程热误差补偿方法及其实施系统 |
CN102831274A (zh) * | 2012-08-30 | 2012-12-19 | 中南林业科技大学 | 铝合金厚板淬火应力场分布的修正方法 |
CN103926874A (zh) * | 2013-12-03 | 2014-07-16 | 江苏齐航数控机床有限责任公司 | 数控机床热误差补偿建模温度测点组合的选择优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016012122A (ja) * | 2014-06-03 | 2016-01-21 | キヤノン株式会社 | フォトマスク、光学素子アレイの製造方法、光学素子アレイ |
-
2015
- 2015-12-25 CN CN201510993887.3A patent/CN105631112B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0655597A (ja) * | 1992-08-10 | 1994-03-01 | Sekisui Chem Co Ltd | 射出成形プロセスシミュレーション方法及びその装置 |
CN101629801A (zh) * | 2009-08-18 | 2010-01-20 | 上海理工大学 | 数控磨床导轨热误差确定方法 |
CN101943896A (zh) * | 2010-07-16 | 2011-01-12 | 浙江大学 | 数控机床误差的轨迹再生补偿方法 |
CN102672527A (zh) * | 2012-05-03 | 2012-09-19 | 四川大学 | 数控机床进给系统全工作行程热误差补偿方法及其实施系统 |
CN102831274A (zh) * | 2012-08-30 | 2012-12-19 | 中南林业科技大学 | 铝合金厚板淬火应力场分布的修正方法 |
CN103926874A (zh) * | 2013-12-03 | 2014-07-16 | 江苏齐航数控机床有限责任公司 | 数控机床热误差补偿建模温度测点组合的选择优化方法 |
Non-Patent Citations (1)
Title |
---|
基于接触分析的立式加工中心热态特性研究;李少龙;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20130930(第09期);第9-63页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105631112A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108062432A (zh) | 一种激光选区熔化过程的数值模拟方法 | |
Miao et al. | Study on the effects of changes in temperature-sensitive points on thermal error compensation model for CNC machine tool | |
TWI309846B (en) | System and method for using first-principles simulation to control a semiconductor manufacturing process | |
CN102073762B (zh) | 一种高精度模具设计方法及模具制造管理系统 | |
Hu et al. | A new cohesive crack tip symplectic analytical singular element involving plastic zone length for fatigue crack growth prediction under variable amplitude cyclic loading | |
CN105631122B (zh) | 一种机床大件的热变形仿真分析与建模方法 | |
CN102169518A (zh) | 精铸涡轮叶片模具型腔精确定型方法 | |
US8370774B2 (en) | Constructing mapping between model parameters and electrical parameters | |
Iftikhar et al. | Turbine blade manufacturing through rapid tooling (RT) process and its quality inspection | |
Lee et al. | Shape optimization of the workpiece in the forging process using equivalent static loads | |
TWI689724B (zh) | 晶片溫度計算方法及晶片溫度計算裝置 | |
CN104881514A (zh) | 一种面向3d打印的壳状部件轻量化建模系统 | |
CN107357946A (zh) | 基于热压膜成型技术的壳状牙科器械制作工艺的验证方法 | |
CN106529005A (zh) | 一种基于有限元仿真的超高强度钢螺纹类工件局部回火感应线圈的设计方法 | |
CN107153724B (zh) | 基于迭代算法的芯片温度分析方法 | |
CN105631112B (zh) | 一种厚板类零件的热变形仿真分析与建模方法 | |
Mayavaram et al. | Optimization of bearing lengths in aluminum extrusion dies | |
CN108595895A (zh) | 一种用于预测铝合金大型构件残余应力的方法及系统 | |
Xiong et al. | Numerical methods to improve the computational efficiency of thermal analysis for the die casting process | |
CN101976291B (zh) | 一种热交换器板片的制作方法 | |
JP2016031628A (ja) | 電池パックに関する温度予測方法及び温度予測装置 | |
CN106326574A (zh) | 一种计及焊接变形和残余应力的极限强度数值计算方法 | |
CN106442407B (zh) | 一种增材制造产品的可靠性验证方法 | |
CN113297758A (zh) | 一种大型复杂筋板件预锻成形初始毛坯优化设计方法 | |
CN110955999B (zh) | 轮胎三维温度场仿真方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210824 Termination date: 20211225 |
|
CF01 | Termination of patent right due to non-payment of annual fee |