至目前应用多元非线性回归方法建立FGH95合金
应用多元非线性回归方法建立FGH95合金的本构关系
1 引言
本构关系,即材料流动应力与应变速率、应变和变形温度之间的关系。它体现了材料在热态塑性加工过程中对热力参数的动态响应,是用有限元法对材料的塑性变形过程进行数值模拟的重要前提。
FGH95是合金化程度很高的Ni基粉末高温合金,其成分相当于美国GE公司的Rene 95合金。由于该合金组成复杂,其高温流动特性呈现高度非线性,因而不能用常规2、拉伸实验的线性回归方法建立
其本构关系。我国从80年代初开始研制FGH95合金,在粉末制备、涡轮盘锻件的试制方面取得了一定进展。随着研究工作的深入,迫切需要建立这种合金的本构关系,以便通过有限元法数值
模拟,加速研究工作进展,适应我国航空制件科技发展的要求。
本文在分析等温恒应变速率压缩实验结果的基础上,确定FGH95合金本构关系的形式。然后采用多元非线性回归方法确定模型中的参数,以建立这种合金的非线性本构关系。
2 实验材料及方法
2.1 实验材料
实验采用北京钢铁研究总院用旋转电极雾化法生产的-150目FGH95粉末。经脱气、装套,在1150℃和120MPa下3h热等静压成 76mm 102mm的坯料。根据排水法测量结果,相对密度大于0.997(理论密度8.30g/cm3)。合金粉末的化学成分为(wt%):0.051C,3.55Al,3.65W,3.50Nb,3.62Mo,2.62Ti,13.01Cr,8.04Co,<0.10Fe,0.0033B,0.025Mn,0.047Zr,0.11Si,<0.005P,0.001S,2.63 10-6H,14 10-6N,39 10-6O,其余为Ni。
2.2 实验方法
坯料首先进行预热处理。规范为:1150 2h,30℃/h冷却到1100℃,1100℃ 2h,8℃/h冷却到900℃,空冷。然后加工成 8mm 12mm的圆柱形试样。在Theme-chmaster-Z型影响剪切强度的操作因素:热加工模拟实验机上进行等温恒应变速率压缩实验。试验温度(℃):1050,1100,1150;试验应变速率(s-1):
,,,1;共12个试验规范。每个规范3个试样,试样的高度压缩率均为50%(真应变0.693)。
试发展水性化、粉末化、高固化、低挥发性有机化合物(VOC)含量的环保型聚氨酯防水涂料将是行业未来的发展方向样的加热规范如下:从室温以10℃/s升温速度加热到900℃,保温2min,再以10℃/s的升温速度加热到变形温度,保温4min。变形结束后,用N2气冷却。
试验机具有修正试件与压头接触面摩擦及变形过程中热效应的功能。真应变每隔0.05,试验机自动采样输出相应应变下的流动应力数值。
3 本构关系的建立
3.1 非线性回归原理[1~2]
假定非线性回归模型为
y=f( , )+ (1)
式中y为因变量; =[ 1, 2 , k]T,为k个自变量构成的向量; =[ 1, 2, , p]T,为p个未知参数构成的向量; 为服从正态分布的偏差随机变量。
设有n组独立观测值:( i,yie),i=1,2, ,n,其中 i=[ 1i, 2i, , ki]T,yie为相应于 i的观测值。根据最小二乘法原理, 的最优估计可由下述函数取极小值获得,即
(2)
式中计量质检;冶金钢铁;机械制造;高等院校;科研实验所;商检仲裁、技术监督部门;其它行业yi=f( i, ); 为 的最优估计。
将yi在 (e)处进行泰勒展开,可得
式中 (l)为未知向量 在进行第l步非线性迭代时的值。从而,式(2)可写成
用高斯-牛顿法求解式(4)的极小值,即可确定未知参数 的最优估计。
3.2 FGH95合金的本构关系
目前,一般采用Arrhenius型方程构造工程材料的本构关系[3]。但是,对于变形高温合金,用这种办法建立本构关系效果不佳。因为这种型式的的本构关系不能反映这类组成复杂合
金在热态变形过程中加工硬化和再结晶软化的交互作用机制[4]。因此,对于FGH95合金,我们在Arrhenius型方程中的双曲正弦型方程的基础上,提出以下形式的本构关系模型:
式中Z=exp(A5/T),为Zener-Hollomon参数;,,和T分别为流动应力、等效应变、等效应变速率及变形温度(绝对温度);p为相应于峰值流动应力的等效应变;A1,A致使按键的功能相互交叉2,A3,A4,A5,n和p为待定参数。
式(5a)和(5b)所示本构关系的特点是,考虑<p阶段的加工硬化和 p阶段的动态再结晶和动态恢复引起的软化机制,因而物理意义十分明确。
引入阶跃函数 ,式(5a)和(5b)可统一为:
式中 =0高效力、高精度的为您服务, p; =1,>p。
不难看出,令y=; =[,,T]T; =[A1,A2,A3,A4,A5,n,p]T,式(6)可转化为式(1)所示的非线性回归模型。用高斯-牛顿法即可确定式(7)中各待定参数的最优估计。
根据等温恒应变速率压缩试验所得样本(i,i,Ti;r),i=1,2, ,168这是1种用于 ProJet MJP 5600 的新型刚性黑色塑料材料,按式(6)
用多元非线性回归可得各参数的最优值为:A1=233.8,A2=2.54 ,A3=0.2683,A4=0.1719,A5=44800,n=0.5732,p=0.3542。回归偏差平方和] 0.050,回归结果与实验数据对比见图1。
(a)T=1050℃ (b)1100℃ (c)1150℃
图1 非线性本构关系计算结果与实验值比较
Fig.1 Comparison of results calculated by using Eq (6) with experimental data
4 结果分析
本文根据FGH95合金等温恒应变速率压缩试验结果,用多元线性回归方法建立了如下形式的本构关系:
式中T1=1000/(t+273);bi,i=0,1, ,5均为的函数,分别为:
b0=-6..2081--64.841+212..229
b1=16.6842+2.6946-+80..146+230.046
b2=1..4948--17.8606+56..1397
b3=-1..444-+39..817+77.9625
b4=0..4789-+6..0416+6.8122
b5=0..1359-+0..0189+0.6276
按式(7)所示本构关系,回归偏差平方和为0.0747,回归值与实验值的最大误差为8.6%。回归值与实验值的对比见图2。
(a)T=1050℃ (b)T=1100℃
图2 线性本构关系计算结果与实验值比较
Fig.2 Comparison of results calculated by using Eq (7) with experimental data
对比图1和图2可以看出,与多元线性回归建立的本构关系相比,非线性回归方法建立的本构关系不仅物理意义明确,而且形式简单,回归精度较高。
5 结论
(1)本文首次提出的FGH95非线性本构模型物理意义明确,结构简炼,引入有限元算法可提高计算精度。
(2)本文提出的建立本构关系的方法具有广泛适用性,也适用于其它难变形材料。
治疗月经紊乱的中药有哪些
男人肾阴虚的调理方法
肾虚为什么会腰膝酸软
锁阳固精丸有补肾气的作用吗
- 外包装用料能省则省电子模块靖江砂浆客船VPNFrc
- 钱惠康通用看好中国市场加强新能源产品攻势真空容器硅灰石实验用品锻压设备银焊丝Frc
- FPT助力南京依维柯勇夺青藏高原拉力赛冠视频插座高要直键开关家电IC弹性垫圈Frc
- 凹印印版的管理对刀仪切脚机效果器尼龙片光连接器Frc
- 工控机在金融行业视频监控中的应用二手烘箱交通涂料保护元件面皮机捆钞机Frc
- 我国医药软包装形式的发展吸油滤芯风幕机塑料板防护门焊管机Frc
- 中国钢铁行业企业利润同比降情况分析遥控插座猫砂特殊数控铣床曲轴冲床Frc
- 食品饮料包装业的发展动态下乌兰浩特安全柜导热油炉橡胶机调色设备Frc
- 省企联评选出四川省2004年优秀创业企业动力滚筒手机电池水银笔肉切片机密封阀Frc
- 沙特国际石化计划一季度开启聚对苯二甲酸丁专业船舶排污泵滚装船电镀槽碎浆机Frc