考虑参数不确定性的转向架构架结构强度分析
智鹏鹏1 李永华2 陈秉智2 1.大连交通大学机械工程学院,大连,116028
2.大连交通大学机车车辆工程学院,大连,116028 : , D
摘要 考虑工程实际中设计参数的不确定性 基于 最优试验设计和可靠性理论对转向架构架进
. UIC510G3G1994 ,
行结构强度分析 根据 标准对构架进行结构强度分析 确定受力较大部位的设计参
. APDL , D ,
数 运用 语言建立构架的参数化模型 并对设计参数进行 最优试验设计 进而计算构架的结. , , ,构强度 利用最小二乘法对试验数据进行拟合 建立多项式响应面函数 并对其进行方差分析 检验响
. , MonteGCarlo
应面的精度 基于高精度响应面函数 采用 拉丁超立方抽样方法对构架进行结构强度的. : , .可靠性分析 研究结果表明 设计参数的变化对结构强度的影响较大 且各设计参数的交互作用较小
0.9843, .
转向架构架的结构可靠度为 表明了不确定参数对构架结构强度和安全性的影响程度 基于D , .最优试验设计的响应面法工作效率高 可为其他复杂结构的不确定性分析提供参考:D ; ; ;关键词 最优试验设计 转向架构架 响应面函数 不确定性分析:U272.1中图分类号
DOI:10.3969/j.issn.1004 132X.2019.01.004 ( ) ( OSID):
开放科学 资源服务 标识码
段对可能出现的各种问题进行安全评判和设计参
.
数的修改
,
目前 已有学者对转向架构架的结构强度进
. [3] ANSYS Workbench行了研究 李万莉等 利用
UIC615G4G2003
软件根据 标准对新型轨道运料
. [4]车转向架构架进行了静强度分析 史艳民等 对
4跨座式单轨交通作业车转向架构架进行了基于
Goodman种超常载荷工况的静强度和基于 疲劳
. [5]极限图的疲劳强度分析 李金城等 根据欧洲EN13749 70%
标准对 低地板有轨电车动力转向
、 .架构架进行了静强度 疲劳强度和模态分析 周
[6] [2] ANSYS PDS ,
元 和卢耀辉等 利用 的 模块 对
,构架静强度可靠性进行分析 得到了转向架构架
.LI [7]参数对可靠性的灵敏度 等 利用改进的正
, ANSYS交设计去拟合响应面 并利用 软件对高速动车组转向架构架焊接接头进行了疲劳可靠性
, .
分析 证明了改进正交设计方法的可行性 上述
,研究多以确定性的参数为属性建立模型 分析转
,向架构结构强度是否满足标准和要求 未能考虑
. PDS工程实际中参数的不确定性 基于 模块的
, ,不确定性分析虽然计算精度高 但是计算效率低
.对复杂程度较高或模型较大的结构不适用
,本文以某货车转向架构架为研究对象 综合
,考虑工程实际中设计参数的不确定性 对转向架
,构架进行结构强度的可靠性分析 进而评估设计
.参数的不确定性对转向架构架结构强度的影响
转向架构架设计参数的确定
转向架构架的设计参数对其整体的结构性能
,
有着重要的影响 是设计阶段对结构进行安全评
. ,
定的基础 为确定关键的设计参数 提高转向架
,构架结构强度分析的效率 本文根据国际铁路联
UIC510G3G1994
盟标准 对转向架构架进行结构
.
强度分析
Shell181 ,转向架构架整体采用 壳单元 转臂
Solid185 ,
座采用 实体单元进行网格划分 节点总
128819 , 121978 .
数为 个 单元总数为 个 在轴箱
Rbe3 ,吊挂的两个圆孔采用 单元模拟螺栓连接
Beam188 ,同时在两圆孔之间采用 单元模拟螺栓
Spring COMBIN14转臂座之间采用 单元 模拟弹
,
簧连接 转臂座实体网格与侧梁壳单元采用刚性
, CP_STRUC.元模拟连接 单元类型选择 构架有
1.
限元模型见图
UIC510G3G1994 ,
根据 标准规定 转向架构架
5
结构强度的评定为 种超常载荷工况下结构的静
.
强度 本次分析仅施加第一超常工况下的载荷的
, 366kN .
情况 即在中心盘处施加 的载荷 在一
4系弹簧单元和转臂座弹簧单元刚性连接处的 个
:、 ,位置分别施加如下约束YZ方向的位移约束
2 3; , 3; 、 、
即和 Z方向的位移约束 即 XYZ方向
, 1、2 3; 、 ,的位移约束 即 和 XZ方向的位移约束
1 3. 2.
即 和 具体约束情况见图
Q345E,转向架构架钢板材质为 屈服强度为345 MPa,转向架构架的最大应力取有限元分析
,结果的最大节点等效应力 构架的应力云图见图3. 3a ,
由图 可知 转向架构架的最大应力出现在
, 181.9 MPa,
中心销座孔 为 小于材料的屈服强
. ,
度 通过分析转向架构架的结构应力可知 构架横梁和侧架的上下盖板与侧板对其结构应力的影
,
响较大 可将其厚度作为结构强度不确定性分析
.
的重要设计参数
不确定性的分析方法
2.1
响应面法
(response surface methodology,响应面法
RSM)
是一种利用试验设计数据求解多元方程组并[ 8G9].进行多元回归分析的经验统计建模技术 采用RSM
法进行结构强度分析能够较大地提高分析效, .
率 有利于大型复杂模型的不确定性分析 基于响
4.应曲面的转向架构架结构强度分析流程见图
2.2 D
最优试验设计D
最优试验设计在选择试验运行方面具有
, .较强的预测能力 适用于每一个特定的研究 它
2~24 ,
可以设计 个因素 能够最大限度地减小与
23
(a)
正面轴侧图 (b)
反面轴侧图 (c)
主视图 (d)
侧视图
,该模型系数估计相关的方差 对于高度受限的设[ 10G11].计是有用的 本文选择该方法对转向架构架
.
进行试验设计
不论因变量与自变量之间存在何种回归关,
系 可设其回归模型为
y β1f1( x )+ β2f2( x )+ ( x )+ ε (1)
= + βmf
α α α m α α , ,
式中x 为给定的因子区域X中一点 若因子空间为p维
α
, , ( ,, , ); x p x =欧氏空间 则 为 维向量 xα1 xα2 x
α a αp
f1( ), f2( ), , ( ) ; ,
x x f x 为连续函数β1,β2, βm
α α m α
24
; .为m个待定系数ε 为服从正态分布的随机变量
α
,, ,若试验方案由N个试验点xx x 组
1 2 N
, (1) :
成 则可以得到式 的参数估计é f1( x ) f2( x ) ( x ) f ù
1 1 m 1
f1( ) f2( ) ( ) x x f x
( )= 2 2 m 2 2) Fx
⋮ ⋮ ⋮ f1( ) f2( ) ()
ë x x f x û
N N m N
信息矩阵
Σ N
A F′xFx ( ) ( )= f′( xx ) f( ) (3)
=
α α
α= 1 D
最优设计就是使得信息矩阵的值极大化
.D
的一种设计 最优试验设计可以使回归系数
,, ,
的估计值bb b 的方差所构成的密集椭球12 p
,
体的体积最小化 同时可以使模型回归预测值方[ 12].差最大值达到最小
D ,基于 最优试验设计提供的初始数据点 采用二次多项式拟合设计参数与响应值的函数关
,
系 并通过最小二乘法估计获得各项的待定系
.
数 为判断拟合的响应面函数能否准确表示转向
,
架构架的结构特征需要对其进行方差分析(analysisofvariance,ANOVA) .
和精度的检验
1 ,,
表 给出了五因素的方差分析表 其中n为试验
.
总次数通过方差分析能够得到对响应值有显著影响
、
的因素 各因素之间的交互作用以及显著影响因
, .素的水平 进而确定设计参数选择的合理性
.响应面函数精度的测定常采用r检验法 设
,计参数x的各个观测点与回归方程越靠近 即r2
1, ,
越接近 表明回归直线和观测点越接近 响应面
, 13]. [拟合程度越好 精度越高
2.3
结构可靠性评估根据机械结构可靠性评估的应力 强度干涉, ,
理论 假定设计参数均服从正态分布 其表达式为
1 x - μX
f ( x )= 2πexp(- ( ) 2) (4) X σ σ
X
X ,、 ,
式中μX σ 为分布参量 分别表示随机变量X的均值和
X
; ( ∞,∞).
标准差x的取值范围为 -
(4),
根据式 转向架构架结构可靠性评估的极
:
限状态方程定义如下
( (5) R = f x )- S X ,( f x ) ; .式中 为构架的最大应力S为材料屈服强度
X
(5),
根据式 可靠度可定义为
n′ (6) Z =
N′
, ( ≤0) ;式中n′为应力不大于强度 R 的抽样次数 N′为抽
.样的总次数
,对转向架构架进行可靠性评估 实际是在考
, ≤0虑各设计参数不确定性的条件下 求出R 的
,
概率分布特征 评估设计参数的不确定性对结构
.
强度的影响程度
, MCLHS为提高计算效率 本文采用 对不确
, 14]: [定性设计参数进行随机采样基本步骤为① 将 每个随机变量的分布函数按照等概率原则
;②
分成互不重叠的区间 根据每个随机变量的概
, ;③率密度函数 从每个区间任意选择一个值 重
(1) (2),
复步骤 和步骤 直到选择了每个随机变
;④
量 将x 的n′个值和x 的n′个值进行随机
i j≠ i
.
组合
3
转向架构架结构强度不确定性分析
3.1 选择初始样本点
1
根据第 节分析得到的转向架构架的重要设
,
计参数 将板壳厚度和材料的弹性模量作为不确
, APDL
定的设计参数 运用 语言对有限元分析模
, 2.型进行参数化 各设计参数的统计特征见表
表2 设计参数统计特征
Tab.2 Statisticalcharacteristicsofdesignparameters
D ,采用 最优试验设计选取初始样本点 并
,将其代入到参数化的有限元分析模型中 得到试
.D验设计后转向架构架的最大应力值 最优设
3.
计试验及响应值见表
表3D 最优试验设计及响应值
Tab.3 DGoptimalexperimentaldesignandresponsevalue
,取任意两个设计参数为坐标轴XY和 根据3 MATLAB2015b
表 通过 编程绘制设计参数的
, 5. 5
样本空间 见图 图 能够直观地反映设计点在
. 5空间的位置以及重复试验的样本点 由图 可
25
,31 ,
知 组设计点均匀地分布在样本空间 且某一
,
设计点存在重复设计 其目的在于提升试验结果
, .的可靠度与精确度 减小误差
(a)设计参数t1、t2 对应的样本空间
(b) 、
设计参数t3 t4 对应的样本空间
图5 设计参数的样本空间
Fig.5 Samplespaceofdesignparameters
3.2 响应面的建立
3 31
根据表 得到的 组包含仿真模型信息的
, ,数据点 利用最小二乘法求解多项式的系数 建立
:转向架构架最大应力的响应曲面函数
Smax =- 3793.98784-78.248t1 + 16.56 t2 - 176.23653t3 + 41.6372t4 + 0.05329 - 0.75589t1t2 +
E
3.48750t1t3 - 1.19447t1t4 + 1.04286×10 - 4E - 2.07642t2t3 + 0.79469t2t4 - 6.224×10 - 5t 2E -
2.76799t3t4 - 8.8535×10 - 5t 3E - 7.0162×10 - 5t 4E + 1.41053t1 + 0.85715t2 + 6.92059t3 + 0.46283t4 -
2 2 2 2
1.23773×10 - 7E2 (7)为了更好地观测不确定设计参数和构架结构
,6 7
强度之间的关系 图 和图 分别示出了设计参数
.对构架结构强度的三维响应曲面图和等高线图
6 ,
由图 可知 设计参数的变化对转向架构架
, 6a ,的结构强度具有一定的影响 由图 可以看出侧架侧板和横梁侧板厚度的增加均会导致构架所
, .
受应力的减小 与实际情况相符 在相同的厚度
, ,
变化区间范围内 相比横梁侧板 由侧架侧板厚度
,变化产生的应力减小程度更大 表明侧架侧板对
. 6b ,构架结构强度的敏感性较强 由图 可以看出 26
(a)设计参数t1、t2 的响应曲面 .横梁上下盖板对构架结构强度的敏感性较强 设计参数与构架结构强度的等高线形状为椭圆说明
,两个设计参数的交互作用显著 圆形说明交互作[ 15]. 7 ,
用不显著 由图 可知 设计参数之间的交互
, ,作用显著性较低 表明在转向架构架设计时 可以
.忽略交互作用不显著的设计参数的变化
3.3
响应面的精度检验响应面的精度是保证试验设计和响应面函数
,
有效性的基础 是进一步利用该模型进行分析的
.
前提 通过对响应面精度的检验可以获得响应面
,函数的准确性及选择设计参数的显著性 有助于
.
判定选取的设计参数是否合理 本文利用ANOVA , 4.
对该模型进行分析 结果见表
4,
表 中FP值和 值均代表相关系数的显著, , , .性F值越大 P值越小 说明相关系数越显著
4 , ,
通过表 可知 该模型的拟合精度较高 且设计参
、、、 ,
数tttt 均非常显著 说明设计参数的选择1234
, .
是合理的 能够反映结构强度的变化 参数E表
,
现不显著说明该设计参数对结构强度影响不
. ,
大 各设计参数之间的交互作用显著性较低 说
,明各设计参数的相关性较小 与等高线图的分析
. 0.9948,结果吻合 响应面函数的测定系数r2为
.
说明拟合程度较高
,8为了更好地观测响应面的拟合精度 图 和9
图 分别给出了试验值和预测值的相关图和残差
. 8 9 ,
分布图 由图 和图 可知 所有的设计点均在45° , 2.0 ,
对角线附近 且残差基本控制在 范围内说明转向架构架响应面函数的预测值和试验值相
.
当接近
图9 试验值和预测值的残差图
Fig.9 Theresidualvaluesoftestandpredictedvalues
3.4
转向架构架结构可靠性评估
(4), MATLAB2015b
基于式 利用 软件使设
10000 ,
计参数均产生 个随机数 用以进行可靠性
, 10.评估 不确定设计参数的概率分布特征见图
10 ,根据图 所产生的不确定设计参数 运用MCLHS 10000 ,
对其进行 次抽样 并将结果代入(7),
式 得到转向架构架最大应力的频率分布直方
, 11. 11 ,
图 见图 由图 可知 设计参数的不确定性
,使转向架构架的结构强度出现了一定的波动 主
150~300 MPa,
要分布在 同时该波动导致转向架构架的最大应力出现超出材料屈服强度的情
.
况 相比传统意义上仅按照标准对构架进行确定
,
性的静强度分析而言 不确定性分析能够更好地
.评估构架设计的合理性和安全性
(7) (5),
将式 得到的计算结果代入式 并结合(6), ,式 得到构架结构强度的可靠性评估结果 见
12 13.
图 和图
12 10000
图 展示了进行 次计算后转向架构
27
(a)
产生的随机分布 (b)
产生的随机分布 (c)
产生的随机分布 (d)
产生的随机分布
, 200 MPa架的强度抽样历史 其中构架应力在 左
.
右居多 空心圆圈代表应力超过构架材料的屈服
,
强度 说明在这些设计参数范围内构架是不可靠
, 1.
的 安全系数小于
13 , ,
由图 可知 当抽样次数较少时 设计参数
,
的随机组合较少 构架所受应力均低于材料的屈
, 1. ,服强度 结构可靠度为 当抽样次数增加 构架
, 500 ,的可靠度下降 在 次左右时趋于收敛 得到转
0.9843,向架构架的结构可靠度为 相比确定性的
,
结构强度分析 考虑设计参数的不确定性使得计
, .算结果更加贴近工程实际 计算结果更加精确
结论
(1) UIC510G3G1994
根据 对转向架构架进行
,
结构强度分析 得到了构架受力较大处的位置为
、 ,
构架横梁 侧架的上下盖板和侧板 并将其厚度确
.定为结构强度分析的不确定设计参数
(2) D ,采用 最优试验设计并结合响应面法
,对转向架构架结构强度进行不确定性分析 获得
,
了高精度的响应面 直观地展示了设计参数的变化对构架结构强度的影响和各设计参数之间的交
.
互作用 与传统的直接对设计参数进行抽样的强
, ,度计算比较 本文方法提高了分析效率 同时保证
,
了分析结果的准确性 适合其他大型复杂结构的
.
不确定性分析
(3) ,
基于高精度响应面的结构可靠性评估 进
———考虑参数不确定性的转向架构架结构强度分析 智鹏鹏 李永华 陈秉智
.一步验证了考虑不确定设计参数的重要性 转向
0.9843,架构架的结构可靠度为 表明了在考虑设
,计参数不确定的条件下 构架的结构强度有超过
, ,,材料屈服强度的可能性 即结构失效 因此 为确
,保结构设计的安全性 在设计过程中考虑参数的
.不确定性更加符合工程实际的需要
:
[1] , .
王洪娇 杨翊仁 地铁车辆转向架构架静强度分析
[J]. ,2016,19的一般方法研究 城市轨道交通研究(11):43G45.
WANG Hongjiao,YANG Yiren.ConventionalApG proachtotheStaticStrengthAnalysisofMetroVeG hicleBogieFrame[J].Urban MassTransit,2016, 19(11):43G45.
[2] , , ,.
卢耀辉 曾京 邬平波 等 铁道车辆转向架构架可
[J]. ,2010,31靠性参数灵敏度分析 中国铁道科学(1):111G115.
LU Yaohui,ZENGJing,WU Pingbo,etal.ParaG metricSensitivityAnalysisfortheBogieFrameReG liabilityofRailway Vehicle[J].ChinaRailwaySciG ence,2010,31(1):111G115.
[3] , , ,.
李万莉 严俊 朱福民 等 新型轨道运料车转向架
[J]. ,2014,12构架结构分析 中国工程机械学报(4):336G341. LIWanli,YANJun,ZHUFumin,etal.Structural AnalysisonBogieFrameforNew RailVehicle[J]. ChineseJournalofConstruction Machinery,2014, 12(4):336G341.
[4] , , ,.
史艳民 缪炳荣 陈建政 等 跨座式单轨交通作业
[J]. ,车转向架构架静强度及疲劳分析 机车电传动2015(5):65G68.
SHIYanmin,MIAO Bingrong,CHENJianzheng, etal.StaticStrengthandFatigue AnalysisofBogie Frame for StraddleGtype Monorail Transportation Working Vehicle[J].Electric Drivefor LocomoG tives,2015(5):65G68.
[5] , , .70%
李金城 李芾 杨阳 低地板有轨电车动力转向
[J]. ,2016(3):87G91.架构架强度分析 机车电传动LIJincheng,LIFu,YANG Yang.Strength Analysis for70% LowGfloorTram Car MotorBogieFrame[J]. ElectricDriveforLocomotives,2016(3):87G91.
[6] . ANSYS
周元 基于 的高速客车转向架构架可靠性
[D]. : ,2012.
研究 成都 西南交通大学
ZHOU Yuan.Reliability ResearchofBogieFramefor High Speed Passenger Car Based on ANSYS[D]. Chendu:SouthwestJiaotongUniversity,2012.
[7] LIY,QIN Q,WANG Y,etal.WeldedJointFatigue Reliability AnalysisonBogieFrameofHighGspeedEG lectric MultipleUnit[J].JournalofShanghaiJiaotong
University(Science),2017,22(3):365G370.
[8] ANGIN D,TIRYAKI A E.Application of ReG sponseSurface Methodologyand Artificial Neural NetworkonPyrolysisofSafflowerSeedPressCake [J].Energy Sources,Part A:Recovery,UtilizaG tion,and Environmental Effects,2016,38(8): 1055G1061.
[9] , , ,.
张伟 毛建国 魏特特 等 基于响应曲面法的煤油
[J]. ,2016,27发动机气门正时优化 中国机械工程(7):933G938.
ZHANG Wei,MAOJianguo,WEITete,etal.OpG timizationof Value Timing ofa Kerosene Engine Basedon RSM [J].China MechanicalEngineering, 2016,27(7):933G938.
[10] DJURISJ,VASILJEVICD,JOKICS,etal.ApG plicationofDGOptimalExperimentalDesign MethG odtoOptimizetheFormulationofO/W Cosmetic Emulsions[J].InternationalJournalof Cosmetic Science,2014,36(1):79G87.
[11] SALVO P,SMAJDA R,DINIV,etal.ADGOptiG malDesignto ModelthePerformancesofDressG ingsand Devices for Negative Pressure Wound Therapy[J].JournalofTissueViability,2016,25 (2):83G90.
[12] , , .DG
柯作楷 唐绍裘 杜海清 最优实验设计及其在
[J]. ,1987(2):陶瓷材料研究中的应用 陶瓷研究
10G16.
KEZuokai,TANGShaoqiu,DU Haiqing.ExperiG mentalDesignofDGOptimizationandItsApplicaG tionin Ceramic Materials[J].Ceramic Studies Journal,1987(2):10G16.
[13] .MATLAB [M].张德丰 概率与数理统计分析 北
: ,2010:245G247.
京 机械工业出版社
ZHANG Defeng.MATLABProbabilityand MathG ematicalStatistics Analysis[M].Beijing:China MachinePress,2010:245G247.
[14] MIJ,LIYF,YANGYJ,PENG W,etal.ReliabiliG tyAssessmentofComplexElectromechanicalSystems underEpistemicUncertainty[J].ReliabilityEngineerG ing & SystemSafety,2016,152:1G15.
[15] MONTGOMERY DC.DesignandAnalysisofExG periments[M].New York:John Wiley & Sons, 1976:30G35.
( )
编辑 袁兴玲
: , ,1989 , . CAE
作者简介 智鹏鹏 男 年生 博士研究生 研究方向为
、 、 .EGmail:关键技术车辆结构可靠性不确定性分析与优化
( ), ,1971 ,.
zhipeng17@yeah.net.李永华 通信作者 女 年生 教授
、研究方向为轨道车辆现代化设计方法 车辆结构的疲劳可靠性分
、 、 RAMS .析 机械产品数字仿真与优化设计 质量与 工程 发表论50 .EGmail:yonghuali@163.com.文 余篇