基于道路谱的多体仿真迭代疲劳载荷技术_沙大亮
基于刚弹耦合多体动力学仿真的沙滩车下横臂疲劳寿命预测
; s ( N Iy nt at e1e o , i st d ren rts cr sn ys s 2 s) dh ldi . w crn i fe eot rtad i i — a lia eo c,c e hhe hp eercf e uleg a m t r u i u
s o ees ni at f up ti l A V id n ip p r ae nsf aeN S R NT e i / t se t l t o s se,o o l T oei t s a e b  ̄do t r A T A .h na s rh ap ’ s n fa s n h ow
1 0i… m i、 t t . r - ‘ ~ t 坩 ‘, +' ~ t i j、 i、 、 ¨ f m “
中图分类 号 :H1 ,P 9 . 文献标 识码 : T 6T 3 1 9 A
l l 商 了I 下横臂是沙滩车的关键 部件之 . 其工作条件 比较恶劣 , 要 承 各种交变载荷 , 使其产 生损坏 , 甚至断裂。 这就需要在保持合 的几何结构基础上提高下横臂 的抗疲劳 能 ,要在设训‘ 初期 , 利用计算机模拟仿真 , 对结构 的疲 劳寿命进 行分 忻预测 , 有纠对 性 的进行设计 的改进 , 具有较夫 的现实意 义 构件在规定 的使用期 限内能安全使用 。 基于有限元分析预测疲劳寿命首先根据载倚和几何结构计
32 下横 臂有 限元分 析 .
冈结构的阻尼很小 , 固有频率和振 型的影响很小 , 对 可忽略
不计 取载荷 厂 为 0则求解下横臂结构的动力学特性方 为: . () ,
M Kx 0 + : ( 、 1
( 只要模态截断合理 ) 子结构( 单冗 ) 超 的动力学特征是 由苴
横拉杆载荷谱采集与疲劳仿真分析
横拉杆载荷谱采集与疲劳仿真分析作者:朱雄丁玲王朝斌吴旭来源:《汽车科技》2018年第05期摘要:本文以某型汽车多连杆悬架横拉杆为模型进行载荷谱时域历程采集、静强度校核、疲劳仿真分析。
在横拉杆上选择合适位置粘贴应变片,组全桥测量轴向应变。
在拉压力试验机上标定出横拉杆轴向受力与测点应变之间的线性关系。
根据可靠性试验规范采集一个完整循环的横拉杆载荷谱,为疲劳分析提供力信号输入。
建立横拉杆有限元模型,对比实测应变与仿真输出对应点应变,修改验证模型,保证有限元模型的准确性。
以实测载荷谱为输入对横拉杆进行疲劳仿真分析,验证横拉杆是否满足可靠性要求。
关键词:横拉杆;载荷谱;有限元;疲劳仿真中图分类号:U467.3 文献标识码:A 文章编号:1005-2550(2018)05-0064-03Abstract: In this paper, the load spectrum time history acquisition, static strength checking and fatigue simulation analysis of a multi-link suspension rack rod are carried out. Strain gauges were pasted on the track rod to measure the axial strain. The linear relationship between the axial force of the track rod and the strain at the measuring point is determined on the tensile testing machine. According to the reliability test specification, a full cycle load spectrum of track rod is collected to provide force signal input for fatigue analysis. The finite element model of thetrack rod is established, the corresponding strain of the measured strain and the simulation are comparedto ensure the accuracy of the finite element model. Based on the measured load spectrum, the fatigue analysis is carried out to verify whether the track rod meets the durability requirements.Key Words: track rod; load spectrum; finite element; fatigue simulation1 前言多连杆悬架横拉杆在汽车行驶过程中将路面激励传递到车身,长时间承受着交变载荷作用,可靠性需要得到充分验证。
汽车道路模拟试验路谱迭代
汽车道路模拟试验路谱迭代“实车道路采谱试验就是为了得到汽车在实际道路行驶中的载荷(应变、加速度、力等信息),在该车的实际运用地区的公路以及试验场进行的实车道路试验。
实车道路试验在汽车开发过程中占有十分重要的地位,通过道路试验可以分别评价汽车的耐久性、舒适性和安全性等个方面,同时考察各个系统和总成的性能。
”实车道路采谱试验就是为了得到汽车在实际道路行驶中的载荷(应变、加速度、力等信息),在该车的实际运用地区的公路以及试验场进行的实车道路试验。
实车道路试验在汽车开发过程中占有十分重要的地位,通过道路试验可以分别评价汽车的耐久性、舒适性和安全性等个方面,同时考察各个系统和总成的性能。
道路试验是汽车开发过程中不可或缺的重要阶段,它包括在高速公路、普通路面、恶劣道路以及各种特殊路面上的测试,是一种检验汽车性能的有效手段。
由于西方国家的路面条件与我国实际情况存在较大差异,因而难以参考国外引进的试验规范和试验路面谱,例如福特公司的JerryZ. Wang和Mark W. Muddiman等人曾于1996年至1997年对中国用户道路载荷谱与福特公司在美国和比利时的试车场道路载荷谱进行了比较研究,发现在国外某种道路路面上不会发生故障的零部件却在国内出现刚度强度问题。
另外我国幅员辽阔,各地道路情况差异较大,因而也有必要对典型地区道路载荷谱进行分析,找出其与试车场道路载荷谱对应关系,可为制定适合我国的试验谱系及规范提供理论依据和有效参数。
将地区道路等效成试车场道路不同路段混合而成的组合路段,即得到地区道路与试验场道路载荷谱的当量关系,就可在试车场按一定比例混合各种路面来再现目标用户地区道路载荷输入,进一步扩展外推后,便可了解较长里程后的损伤情况,达到加速试验的目的。
JerryZ. Wang和Mark W. Muddiman等人曾于1996年至1997年对中国用户道路载荷谱与福特公司在美国和比利时的试车场道路载荷谱进行了比较研究,发现在国外某种道路路面上不会发生故障的零部件却在国内出现刚度强度问题。
道路模拟试验中道路载荷谱的选择方法
图 2 测量轮 ( WF T) 和测试系统 Fig. 2 WF T and testing system
图 3 EV P 道路载荷谱时间历程 Fig. 3 Time history of loading spectra acquired in EV P road
通道
019- [ 35 ] 020- [ 36 ] 021- [ 37 ] 022- [ 38 ] 023- [ 39 ] 024- [ 40 ] 025- [ 41 ] 026- [ 42 ] 027- [ 43 ] 028- [ 44 ] 029- [ 45 ] 030- [ 46 ] 031- [ 47 ] 032- [ 48 ] 033- [ 49 ]
第 1 期
彭 为 ,等 :道路模拟试验中道路载荷谱的选择方法
·9 ·
通道的 ,包括力 (力矩) 、位移 、应变等 。在数据筛选 时 ,为保证合理性并方便分 ,可只考虑与实际疲 劳损伤关系最为密切的应变信号 。先对每个循环 试验数据的 30 个应变通道进行损伤计算 ,每个循
环得到一个 30 ×1 的损伤向量 。接下来再对这 15 个损伤向量进行分析筛选 。各循环损伤向量如表 1 所示 (限于篇幅 ,表中仅列出 5 个循环 ,每个循环 列出 15 个通道) 。
收稿日期 : 2003 - 11 - 10 作者简介 : 彭 为 (1971 - ) ,男 ,安徽霍邱人 ,博士生 ,研究
方向为汽车系统动力学与控制. 靳晓雄 ,教授 ,博 士生导师.
荷的情况 。 近年来 ,上海大众等汽车生产厂家开始利用试
验场和试验台来对整车及零部件在实际使用中的 载荷进行当量模拟 ,在相对短的时间内完成疲劳性 能试验 。由于条件限制 ,目前在进行道路模拟试验 时 ,一般以道路试验实测的加速度或应变信号作为 迭代目标信号 ,这种方法存在精度不够高 、迭代收 敛时间长等缺点 。如以实测车轮轴头的力信号作 为迭代目标信号则可克服以上缺点 ,大大提高精度 和效率 。
基于道路载荷谱的车身疲劳寿命改进研究
整 车室 内道 路 模 拟试 验 分 为轮 耦 合 型道路 模 拟 试验 和轴 耦 合型 道 路模 拟 试 验 两种 类 型 。失 效
位 置 的损 伤 贡献 分 析 结 果 表 明 , 身 失 效 位 置 处 车 占比 9 % 以上 的损 伤 发生在 以垂 向激励 为 主 的各 8 种 粗糙 路 面上 , 因此 本 课 题 中采 用 轮 耦 合 方 式 的
目
—
架试 验加 载信 号 驱 动 台架 , 各 监 测 点 处 采 集 应 在 变响 应信 号 , 后 与 实 际道 路 上 采 集 的 信 号进 行 然 比较 , 图 3所 示 , 算 出误 差 信 号 , 此 修 正 驱 如 计 据
动信 号并 再 次 完 成 加 载 , 此 循 环 下 去 。当迭 代 如 点 响应 与 目标 之 间 的 误 差 达 到 允 许 范 围之 内 时 ,
9 % 以 上 的粗 糙 路 面 采 集 的 信 号 进 行 迭 代 , 流 8 其
程 图见 图 2 。
1 代1 r嘉 模 _ 嚣 拟I . 迭- 1 - ' 膝 等
图 1 车 身 开 裂 研 究 技 术 路 线
2 1 典型 路面 道路 谱 采集 .
研 究 过 程 中采 用 虚拟 试 验 台 的形 式来 建 立 目 标 车辆 的耐 久 性 虚 拟 分 析 模 型 , 虚 拟 模 型 的边 该
界输 入 条件 为轮 轴 六 分 力 信 号 和轮 轴 位 移 信 号 , 而 此信 号必 须 通 过 实 车 的道 路 载 荷 谱 采 集 得 到 。
此外 , 虑 到后期 室 内道 路 模 拟 试 验 和 C E耐久 考 A
分 析 的需要 , 车 场 道 路 试 验 采 集 的 信 号 包 括 车 试
基于道路谱的整车疲劳分析
MAX
设置分析 参数
保存Group信息文件
保存Channel MAX信息文件 保存分析任 文件
提交分析
疲劳分析操作流程图
疲劳计算 实车路试对比
三、FEMFAT软件在长安开发 架疲劳仿真
扭力梁 架疲劳仿真
稳定杆 架疲劳仿真
悬架系统 架疲劳仿真
2 模型检查 计算TB状态模型的模态,对模型进行检查
3 静力分析,采用惯性释放方法,计算各通道6方向的静力结果
某款汽车接附点编号示意图
二、基于道路谱的车身疲劳分析
5.在FEMFAT中计算基于道路谱的车身疲劳损伤
新建任
导入分析 模型
生成通道临时文件
后处理
按照材料 分Group
导入材料并指 定Group材料
速度及在Z向
的位移 根据此原则路普采集试验一般安装以 传感器用于相关信号采集 1 变片 标定出车 X Y方向力 标定出减震器 或弹簧 Z方向力 2 位移传感器 测量 心Z向位移 3 速度传感器 测量 心处X Y Z方向 速度 测量车身 毂包处Z方向
速度 4 GPS 测量试验车纬度 经度 海拔 车速等
某款汽车数据采集路线示意图
某款汽车选择的9段特征路面
3. 多体动力学仿真
根据选择的特征路面,进行载荷谱迭代,载荷分解及工况仿真计算
二、基于道路谱的车身疲劳分析
4.计算单位载荷 的 力分布
1 计算单位 力的模型准备
1 一般用TB状态的模型
2 焊点焊缝的处理
3 接附点编号
4 配重保持 多体模型一
5 按照材料对模型进行分组
踏板 架耐久仿真
车 弯曲疲劳 架耐久仿真
三、FEMFAT软件在长安开发体系中的应用
基于道路谱的整车疲劳分析
2 模型检查 计算TB状态模型的模态,对模型进行检查
3 静力分析,采用惯性释放方法,计算各通道6方向的静力结果
某款汽车接附点编号示意图
二、基于道路谱的车身疲劳分析
5.在FEMFAT中计算基于道路谱的车身疲劳损伤
新建任
导入分析 模型
生成通道临时文件
后处理
按照材料 分Group
导入材料并指 定Group材料
某款汽车传感器布置示意图
二、基于道路谱的车身疲劳分析
2.载荷信号采集及数据后处理
1 采集 根据路试道路规划进行多次 少2次以 典型路面载荷采集 2 数据减缩 根据相对损伤值进行数据减缩,确保减缩后的数据 原数据的相对 损伤值相当 3 路普选择 1 该路普多数通道最大值均大于另外的路谱 2 该路谱多数通 道相对损伤值均大于另外的路谱 3 该路谱 存在 合理数据
速度及在Z向
的位移 根据此原则路普采集试验一般安装以 传感器用于相关信号采集 1 变片 标定出车 X Y方向力 标定出减震器 或弹簧 Z方向力 2 位移传感器 测量 心Z向位移 3 速度传感器 测量 心处X Y Z方向 速度 测量车身 毂包处Z方向
速度 4 GPS 测量试验车纬度 经度 海拔 车速等
某款汽车数据采集路线示意图
某款汽车选择的9段特征路面
3. 多体动力学仿真
根据选择的特征路面,进行载荷谱迭代,载荷分解及工况仿真计算
二、基于道路谱的车身疲劳分析
4.计算单位载荷 的 力分布
1 计算单位 力的模型准备
1 一般用TB状态的模型
2 焊点焊缝的处理
3 接附点编号
4 配重保持 多体模型一
5 按照材料对模型进行分组
基于道路谱的整车疲劳分析 ——长安汽车
车辆道路载荷谱的Hilbert边际谱分析方法及应用
车辆道路载荷谱的Hilbert边际谱分析方法及应用车辆道路载荷谱是指车辆在道路上行驶时产生的载荷信号,它对道路的损伤和寿命具有重要影响。
因此,准确地估计道路上车辆的载荷谱显得十分重要。
而Hilbert边际谱分析方法是一种有效的处理信号的方法,能够从时域和频域两个方面对载荷信号进行分析,因此能够应用于车辆道路载荷谱的分析。
Hilbert边际谱分析方法是指通过利用信号在时频域的特性进行边际谱分析,并通过计算峭度来识别信号的特征。
它的基本思想是将信号分解为振幅调制和频率调制两个部分,然后可以对这两个部分进行进一步分析。
对于车辆道路载荷谱,可以将其视为由多个频率分量组成的信号,采用Hilbert边际谱分析方法可以分析出各个频率分量对道路的影响,从而准确估计道路的损伤和寿命。
在实际应用中,Hilbert边际谱分析方法还可以结合相关分析和聚类分析等统计方法,对载荷谱进行进一步分析。
例如,可以将同一路段内不同时间的载荷数据进行关联分析,找出不同车辆的共性和差异性,为道路维修和改造提供科学依据。
此外,还可以采用聚类分析方法对不同道路类型、不同车型和不同行驶状态下的载荷谱进行聚类分析,得到相应的聚类中心,为道路结构设计和材料选择提供决策支持。
总之,Hilbert边际谱分析方法是一种高效准确的分析车辆道路载荷谱的方法,可为道路维护和改造提供科学依据,为道路结构设计和材料选择提供决策支持。
这种方法的应用还需要进一步深入研究和发展,以适应不断增长的道路交通和车辆技术需求。
Hilbert边际谱分析方法的应用不仅局限于车辆道路载荷谱分析,还可以拓展至其他领域。
例如,在医学领域,可以将人体脉搏信号视为一种载荷信号,采用Hilbert边际谱分析方法对其进行分析,从而得出有关人体健康状况的信息。
在工业领域,可以对机器设备的振动信号进行分析,找出设备故障的原因和位置,提高生产效率和设备利用率。
然而,Hilbert边际谱分析方法也存在一些问题和限制。
汽车零部件可靠性测评技术的发展与应用
10.16638/ki.1671-7988.2021.012.035汽车零部件可靠性测评技术的发展与应用宁世儒,张冠勇,庞方超(中汽研汽车检验中心(天津)有限公司,天津300300)摘要:随着汽车零部件可靠性验证已成为产品开发过程中越来越重要的一个环节,并占用了整车研发过程中的大量时间,因此迫切需要研发出更为高效的汽车零部件可靠性测评技术。
文章从编辑道路载荷谱、搭建疲劳耐久性台架试验以及虚拟仿真试验技术三个方面着手进行汽车零部件可靠性测评技术研究及应用现状阐述,并对未来车辆可靠性测评技术的发展方向进行展望,为实现汽车零部件产品开发流程的最优化,降低研发费用以及增强产品竞争力提供重要的理论支撑。
关键词:可靠性强化试验技术;道路载荷谱;可靠性台架试验;虚拟仿真试验技术中图分类号:U467 文献标识码:A 文章编号:1671-7988(2021)12-113-04Development and Application of Reliability StrengtheningTest Technology for Auto PartsNING Shiru, ZHANG Guanyong, PANG Fangchao( China Automotive Research Institute Automotive Inspection Center (Tianjin) Co., Ltd., Tianjin 300300 )Abstract: As the reliability verification of auto parts has become an increasingly important part of product development and takes up a lot of time in the vehicle parts development process, there is an urgent need to develop more efficient reliability evaluation technologies for auto parts. This paper starts to conduct auto parts reliability evaluation research and application status elaboration from three aspects: editing the road load spectrum, building reliability bench test and virtual simulation test technology, and looks forward to the future vehicle reliability evaluation technology, and provides important theoretical support for realizing the optimization of the development process of auto parts products, reducing R&D costs and enhancing product competitiveness.Keywords: Reliability enhancement test technology; Road load spectrum; Reliability bench test; Virtual simulation test technologyCLC NO.: U467 Document Code: A Article ID: 1671-7988(2021)12-113-04前言近年来,国内汽车市场结束了连续29年的增长奇迹,市场对于汽车产品的可靠性和安全性提出越来越高的标准和要求[1-2]。
疲劳载荷谱
疲劳载荷谱
疲劳载荷谱是指衡量或模拟结构物在其寿命内受到的疲劳载荷的统计数据。
它描述了结构物在不同时间点或时间段内所受到的载荷大小和频率。
疲劳载荷谱的构建通常需要进行大量的实验测量或使用计算机模拟技术。
在实验测量中,可以通过使用应变计、加速计、压力传感器等传感器对结构物所受到的加载进行监测和记录。
通过多次重复加载,并将测量结果进行统计分析,就可以得到相应的疲劳载荷谱。
在计算机模拟中,可以通过有限元分析等方法对结构物在不同工况下的载荷进行模拟。
通过对结构物所受到的载荷进行统计分析,就可以得到相应的疲劳载荷谱。
疲劳载荷谱对于结构物寿命评估和疲劳强度设计具有重要意义。
它可以用于确定结构物在预期的使用寿命内是否会发生疲劳破坏,并为疲劳强度设计提供基础数据。
疲劳载荷谱的准确性和合理性对于结构物的安全性和可靠性具有关键影响。
基于多轴道路模拟激励谱的摩托车车架虚拟试验方法
基于多轴道路模拟激励谱的摩托车车架虚拟试验方法邹喜红;熊锋;袁冬梅;易鹏【期刊名称】《农业工程学报》【年(卷),期】2014(000)015【摘要】摩托车车架的疲劳可靠性是摩托车整车最重要的性能之一,为了对摩托车车架的疲劳可靠性进行准确高效的分析和评价,提出了基于时频域误差加权和力-位移混合控制的摩托车车架多轴道路模拟试验方法,结合采集的实际行驶道路载荷谱和开发的车架多轴道路模拟试验台,采用时域和频域误差加权系数均为0.5的力-位移混合控制多轴道路模拟试验使摩托车车架实际行驶道路载荷谱模拟迭代精度达到93%,并提取了多轴道路模拟激励谱。
基于摩托车车架多轴道路模拟试验系统,利用HYPERWORKS有限元分析软件和MSC.ADAMS多体动力学仿真软件,建立了摩托车车架刚柔耦合多轴道路模拟虚拟试验平台,以多轴道路模拟激励谱为输入进行了仿真分析,并采用多轴道路模拟试验进行了结果验证,从而建立了基于多轴道路模拟激励谱的摩托车车架虚拟试验方法。
结果表明,仿真结果和试验结果在时频域中趋势和幅值都吻合很好,时域曲线几乎重合,只是频域曲线在13 Hz和24 Hz附近幅值略有差异,但频域曲线均方根值误差均在10%之内,基于多轴道路模拟激励谱的摩托车车架虚拟试验方法能够对摩托车车架疲劳可靠性进行高效准确的考核。
%Motorcycle frame is the key load-bearing part of vehicle, its fatigue reliability, stiffness, strength, quality and other factors play a decisive role in vehicle's control stability and occupant comfort. New frame products must pass the reliability test before put into market. However, current frame fatigue reliability tests mainly use road test and testingground test, which are time-energy consuming and the results are easily affected by environmental factors; Indoor bench tests at present mainly use single axial single incentive or single axial multi incentive mode. The motorcycle frame installation and constraints on the vehicle are very complex, the actual driving loads are complicated and changeable, thus multi-axial loads on motorcycle frame due to road excitation during actual driving could not be simulated or simulation precision is low, prone to over or under test. While the computer simulation technology development provides a new solution to dynamic performance analysis and fatigue reliability evaluation of vehicles and their parts, which including complete 3D entity model building online, creating a virtual laboratory, a virtual running environment and realizing virtual test and evaluation. For accurate and efficient analysis and evaluation of motorcycle frame fatigue reliability, a multi-axes road simulation test method for motorcycle frame based on time-frequency domain error weighting and force-displacement hybrid control was put forward. A motorcycle frame multi-axes road simulation test bench was built with American MTS road simulation test device, which was arranged with two actuator input channels in horizontal axial and vertical axial. As the motorcycle frame multi-axes road simulation test system is a multiple input multiple output system, the remote parameter matrix control strategy was applied to iterate and reproduce actual driving load spectrum. As the frame front fork and shock absorber were substituted by rigid frame fixture, in the horizontal direction, the frame dynamic response was more effective to theforce loading, so the horizontal actuator used force control mode and loaded force spectrum. In the vertical direction, the system used inertial loading and frame dynamic response was more effective to displacement loading, so the vertical actuator was selected displacement control mode and loaded displacement spectrum. The motorcycle frame actual road load spectra were collected and multi-axes road simulation test bench was developed. With the simulation test bench, the multi-axes road simulation test was conducted based on force-displacement hybrid control under time and frequency domain with error weighting coefficient of 0.5. The simulation iteration precision of frame actual road load spectrum achieved 93%. Therefore the multi-axes road simulation excitation spectrum were extracted. The HYPERWORKS and ADAMS were applied in building the virtual rigid-flexible coupled multi-axes road simulation test platform with reference of the actual motorcycle frame multi-axes road simulation test system. Simulation analysis was conducted with the multi-axes road simulation excitation spectrum as input and the outcomes were verified by actual multi-axes road simulation test, thus the virtual test method for motorcycle frame based on multi-axes road simulation excitation spectrum were established. The results show that, the curve trend and amplitude of simulation results coincide well with test results in time and frequency domain, only slight amplitude deviation exists near frequency value 13 and 24 Hz, while the frequency curve RMS error is within 10%. Thus the virtual test method based on multi-axes road simulation excitation spectrum could be well applied to motorcycle frame fatigue reliability assessment.【总页数】7页(P39-45)【作者】邹喜红;熊锋;袁冬梅;易鹏【作者单位】重庆理工大学汽车零部件制造及检测技术教育部重点实验室,重庆400054; 重庆建设摩托车股份有限公司,重庆 400054;重庆理工大学汽车零部件制造及检测技术教育部重点实验室,重庆 400054;重庆理工大学机械工程学院,重庆 400054;重庆市科学技术研究院,重庆 400054【正文语种】中文【中图分类】U467.3【相关文献】1.摩托车车架多轴向多激励道路模拟试验方法研究 [J], 邹喜红;熊锋;余勇;王锐利2.基于道路模拟激励的汽车下摆臂多轴疲劳分析 [J], 钱立军;吴道俊;杨年炯;祝安定;王挺3.AMT执行机构多轴向多激励道路模拟试验方法研究 [J], 邹喜红;熊锋;余勇;刘栓起;彭吉刚4.基于远程参数控制的汽车驱动桥虚拟多轴道路模拟试验系统 [J], 邹喜红;罗洋;袁冬梅;柳春林;刘瑜5.基于道路模拟激励谱的AMT执行机构疲劳寿命分析 [J], 杨平;石晓辉因版权原因,仅展示原文概要,查看原文内容请购买。
基于多体模型仿真的载荷谱虚拟迭代技术分析
基于多体模型仿真的载荷谱虚拟迭代技术分析
邵建;董益亮;肖攀;王俊翔
【期刊名称】《重庆理工大学学报(自然科学版)》
【年(卷),期】2010(024)012
【摘要】以在路试场测取的轮心加速度和车辆内力载荷为输入,以Adams所创建的多体模型为载体,通过虚拟迭代的方法反求得到外部驱动载荷,驱动多体模型仿真,得到各部件的受力谱.某乘用车的虚拟迭代结果证明该方法有效、准确.
【总页数】4页(P84-87)
【作者】邵建;董益亮;肖攀;王俊翔
【作者单位】重庆长安汽车股份有限公司,汽车工程研究院,重庆,401120;重庆长安汽车股份有限公司,汽车工程研究院,重庆,401120;重庆长安汽车股份有限公司,汽车工程研究院,重庆,401120;重庆长安汽车股份有限公司,汽车工程研究院,重
庆,401120
【正文语种】中文
【中图分类】U463
【相关文献】
1.基于整车动力学模型虚拟迭代仿真的转向节载荷谱提取 [J], 汪随风;武振江;杨建森
2.基于虚拟迭代方法的后扭力梁载荷谱提取 [J], 龚春辉;葛文韬;段龙杨;王祖建
3.基于多体模型仿真的载荷谱虚拟迭代技术分析 [J], 邵建;董益亮;肖攀;王俊翔
4.一种加速多体动力学模型虚拟迭代的载荷谱编制方法 [J], 葛文韬;黄晖;刘洲;邱
星;贾慧芳
5.基于虚拟迭代的装载机后处理支架载荷谱获取方法 [J], 付春雨;曾超;刘宏杰;刘伦伦;程市
因版权原因,仅展示原文概要,查看原文内容请购买。
一种加速多体动力学模型虚拟迭代的载荷谱编制方法
第31卷第12期中国机械工程V o l .31㊀N o .122020年6月C H I N A M E C HA N I C A LE N G I N E E R I N Gp p.1504G1511一种加速多体动力学模型虚拟迭代的载荷谱编制方法葛文韬㊀黄㊀晖㊀刘㊀洲㊀邱㊀星㊀贾慧芳江铃汽车股份有限公司产品开发中心,南昌,330200摘要:为加速多体动力学模型的虚拟迭代,提出了一种能够完整保留长里程路面载荷谱损伤值㊁幅值特征和频率特征的编制方法,采用路面时域信号分割以及路面片段组合优化等手段,使得比利时路单次循环路面里程由2.61k m 缩减至1.49k m ,生成了用于虚拟迭代的加速谱.对比了下摆臂及转向节分别在比利时路原始谱和加速谱作用下的疲劳损伤,结果显示缩减前后损伤分布一致,摆臂和转向节的最大损伤比值分别为1.081和1.205,表明了加速谱能够替代原始谱.运用该方法对耐久规范中的越野路及山路长里程路面进行了同样的里程缩减处理,并对摆臂和转向节进行了疲劳分析预测,分析结果显示其寿命能够满足疲劳性能目标,并通过了实车道路耐久验证.所提出的载荷谱编制方法为加速多体动力学模型的虚拟迭代提供了一种解决途径.关键词:多体动力学;虚拟迭代;时域信号分割;里程缩减;疲劳损伤中图分类号:U 463D O I :10.3969/j .i s s n .1004 132X.2020.12.015开放科学(资源服务)标识码(O S I D ):L o a dS p e c t r u m D e v e l o p m e n tM e t h o d o l o g y f o rA c c e l e r a t i n gD yn a m i c sM o d e lV i r t u a l I t e r a t i o n G E W e n t a o ㊀HU A N G H u i ㊀L I UZ h o u ㊀Q I U X i n g ㊀J I A H u i f a n gP r o d u c tD e v e l o p m e n t&T e c h n i c a l C e n t e r ,J i a n g l i n g M o t o r sC o .,L t d .,N a n c h a n g,330200A b s t r a c t :I no r d e r t o a c c e l e r a t e v e h i c l e d y n a m i c sm o d e l v i r t u a l i t e r a t i o n ,a l o a d s p e c t r u md e v e l o pGm e n tm e t h o d o l o g y w h i c hc o u l d r e s e r v e l o n g m i l e s t o n e l o a ds p e c t r u md a m a g e ,a m pl i t u d e c h a r a c t e r i s Gt i c s a n d f r e q u e n c y c h a r a c t e r i s t i c sw a s p r o p o s e d .T h e m i l e s t o n eo fB e l g i a nb l o c ks i n g l e l o o p w a sr e Gd u c e d t o 1.49k mf r o mo r i g i n a l 2.61k m w i t h a p p l i c a t i o n s o f r o a d t i m e d o m a i n s e gm e n t a t i o n a n d c o m Gb i n a t o r i a l o p t i m i z a t i o n ,a n a c c e l e r a t e d s p e c t r u mu s e d f o r d yn a m i c sm o d e l v i r t u a l i t e r a t i o nw a s c r e a t e d .T h e f a t i g u e d a m a g e o f l o wc o n t r o l a r ma n dk n u c k l ew a s p r e d i c a t e du n d e r o r i g i n a l s pe c t r u ma n d a c c e l Ge r a t e ds p e c t r u m o fB e l g i a nb l o c k .A n a l y s i sr e s u l t ss h o wt h a tt h ed a m a ged i s t r i b u t i o nr e m a i n st h e s a m eb e f o r e a n d a f t e r t h e r e d u c t i o n ,t h em a x i m u md a m a ge r a t i o sf o r l o wc o n t r o l a r ma n dk n u c k l e a r e 1.081a n d 1.205r e s p e c t i v e l y ,w h i c hd e m o n s t r a t e t h a t t h e a c c e l e r a t e d s p e c t r u m m a y r e p l a c e t h e o r ig i Gn a l s p e c t r u m.Th e o f f Gr o a d a n dm o u n t ai n r o a d i n d u r a b i l i t y t e s t p r o c e d u r ew a s r e d u c e d b yt h e s a m e r e Gd u c t i o nm e t h o d o l o g y .T h e f a t i g u e l i f e o f l o wc o n t r o l a r ma n dk n u c k l ew a s p r e d i c t e d .A n a l ys i s r e s u l t s s h o wt h a t t h e f a t i g u e l i f em a y s a t i s f y t h e t a r g e t ,a n d t h e v e h i c l e c o m p l e t e s t h e d u r a b i l i t y v a l i d a t i o no f p r o v i n gg r o u n ds u c c e s s f u l l y .T h e l o a ds p e c t r u m d e v e l o p m e n tm e t h o d o l o g y h e r e i n p r o v i d e saw a y to a c c e l e r a t e v e h i c l e d yn a m i c sm o d e l v i r t u a l i t e r a t i o n .K e y w o r d s :m u l t i Gb o d y d y n a m i c s ;v i r t u a l i t e r a t i o n ;t i m ed o m a i ns e g m e n t a t i o n ;m i l e s t o n e r e d u c Gt i o n ;f a t i g u e d a m a ge 收稿日期:20181217㊀㊀修回日期:202001080㊀引言应用多体动力学模型虚拟迭代方法来获得零部件疲劳分析载荷边界已经在行业内广泛应用[1G3],虚拟迭代用载荷谱通常是按照整车道路耐久性试验规范在试验场通过载荷谱采集试验获得的,规范中一般有单个循环里程较长的比利时路㊁山路及越野路等路面,如若对上述长里程路面全部进行虚拟迭代,则时间跨度长且人力资源消耗大,严重影响产品开发时间周期和成本,故有必要对道路载荷谱进行编制,以达到加速的目的.目前常用的载荷谱编制方法为时域损伤编辑4051 Copyright©博看网 . All Rights Reserved.法[4G5],主要应用于整车室内道路模拟试验中,通过设定应变信号的阈值来删除时域载荷谱中的小损伤片段,从而实现对整车室内道路模拟的加速;编辑后的各段信号之间采用半个正弦波或直线拼接成连续的时域载荷谱信号,通过压缩载荷谱时间减少台架试验时间.但此方法的难点在于阈值的选取,阈值过大使得台架试验无法复现道路耐久试验的故障,阈值过小不利于压缩室内道路模拟试验时间.本文采用路面时域信号分割及路面片段组合优化的编制手段来对长里程路面进行里程缩减处理,通过增加大损伤值路面片段的循环次数来替换小损伤片段,从而实现对多体动力学模型虚拟迭代的加速.文献中主要是从时域损伤的角度来确保原始谱和加速谱之间的等效性[6],而本文从损伤㊁幅值和频域特征多个角度来确保加速谱可替换原始谱.1㊀伪损伤理论在工程应用中,如分析汽车底盘及车身部件在坏路路面激励下的响应特征,需要采集很多通道的道路载荷谱数据.在这些多通道数据中,除应变通道外,还有大量的力㊁位移和加速度信号,由于这些通道为非损伤变量,因此引入伪损伤计算来描述这些通道[7G8],利用时域载荷谱的雨流计数以及合适的SGN曲线,按照M i n e r线性累积损伤理论计算得到伪损伤值.载荷G寿命之间的曲线(SGN曲线)用幂函数公式表示为S b N f=C(1)式中,S为载荷幅值;N f为对应S的寿命;b为损伤斜率, b的选择反映被分析零部件的类型,根据经验一般焊接部件b=3,典型汽车零部件b=5[9],光滑表面的零部件b=7;C为与材料㊁应力比㊁加载方式相关的参数.累积损伤D按线性累积损伤理论计算:D=ðg i=1n i N f i(2)式中,D为累积伪损伤;n i为雨流计数统计第i级载荷S i 作用下的循环次数;N f i为第i级载荷S i作用下的寿命;g 为雨流计数时单元格的数量.将式(1)代入式(2)得D=1Cðg i=1n i S b i(3)采集的载荷谱数据可以是应变信号,也可以是加速度㊁位移和力信号.此方法是一种简便且有效计算各通道伪损伤值的方法,并不用于计算某一个特定零件的寿命,伪损伤值本身并无明确意义,主要是用于计算载荷谱数据缩减前后各通道的相对损伤系数.2㊀载荷谱采集载荷谱采集是通过借助各类型传感器直接获得坏路路面对车辆的激励信号,信号一般包括应变㊁位移㊁加速度㊁力㊁力矩㊁车速及经纬度等.本文中的皮卡前悬架类型为扭杆弹簧双叉臂独立悬架,后悬架类型为非独立板簧悬架.采集信号的类型及通道数分配如表1所示,轮心六分力测量见图1,采用拉线传感器测量前悬架减振器位移,见图2.表1㊀采集信号类型及通道分配T a b.1㊀Ty p e o f a c q u i s i t i o n s i g n a l a n d c h a n n e l d i s t r i b u t i o n 信号类型通道数量采集对象六分力28轮心六分力及轮速加速度16轴头㊁车桥位移4前后减振器位移应变16扭杆弹簧㊁板簧前卷耳支架等图1㊀轮心六分力测量F i g.1㊀M e a s u r e m e n t o f s i xGc o m p o n e n t f o r c e o fw h e e l c e n t e r图2㊀前悬架减振器位移测量F i g.2㊀M e a s u r e m e n t o f f r o n t s u s p e n s i o nd a m p e r d i s p l a c e m e n t 3㊀载荷谱编制此车型的整车道路耐久规范中共包含长波路㊁短波路㊁搓板路㊁扭曲路㊁凸块路㊁比利时路㊁越野路及山路等8种不同的坏路,其中比利时路㊁越野路及山路单个循环里程较长,分别为2.61k m㊁2.85k m和11.51k m,故为加速多体动力学模型的虚拟迭代,需对上述3种路面进行里程缩减处理,其余路面均不做处理.里程缩减处理主要步骤为载荷谱时域信号分割㊁分割后各片段时域损伤值比较以及各片段组合优化,如图3所示.5051一种加速多体动力学模型虚拟迭代的载荷谱编制方法 葛文韬㊀黄㊀晖㊀刘㊀洲等Copyright©博看网 . All Rights Reserved.图3㊀道路载荷谱缩减技术路线图F i g .3㊀T e c h n i q u e r o u t e o f r o a d l o a d s pe c t r u mr e d u c t i o n 3.1㊀载荷谱时域信号分割及损伤值比较根据已有项目工程经验,所能接受的最长单次虚拟迭代时间约为40s ,因此将整个比利时路面载荷谱平均分割为9个路面片段,每个路面片段时间为42s,如图4所示.图4㊀比利时路面分割处理F i g .4㊀B e l g i a nb l o c k s s p l i t t i n gpr o c e s s 由于实测轮心三向力(纵向㊁侧向和垂向)信号直接代表了路面不平度对车辆的载荷激励,故按图3所示伪损伤方法分别计算了9个路面片段轮心三向力的伪损伤值之和,计算结果如图5所示.图6所示为9个路面片段左后轮轮心垂向力累积伪损伤值对比.由图5和图6可以看出,编号为3㊁7㊁8㊁9的路面片段比其他5个片段损伤值小,因此选取编号为1㊁2㊁4㊁5㊁6的路面片段作为下步组合优化的对象.3.2㊀各片段组合优化雨流计数法是一种典型的双参计数法,计数原理与实际工作载荷对金属零件的循环应力应变原理相似,能表达载荷谱数据的幅值㊁均值及循环次数等信息,并以雨流矩阵的形式显示出来.由于道路载荷谱采集通常只采集单个循环的载荷图5㊀路面片段损伤值对比F i g .5㊀D a m a g e c o m pa r i s o no f e a c h r o a d s e c t i on 图6㊀累积伪损伤值对比F i g .6㊀C u m u l a t i v e p s e u d od a m a g e c o m pa r i s o n 谱数据,而整车道路耐久试验需在此路面上循环若干次,故需将单个循环的数据通过雨流循环外推法[10]得到实际里程数,即用于各片段组合优化的目标雨流矩阵.比利时路是由石块随机分布组成的环形道路,平均车速约为25k m /h ,车辆通过此路面时会受到宽频随机振动的激励,由于轮心三向力(纵向㊁侧向和垂向)直接代表了路面激励传递给车辆的直接载荷激励,因此各片段直接依据这些通道的伪损伤值来进行组合优化.各片段组合优化的核心问题是在加速谱保留原始谱损伤值㊁幅值域及频率域特征的基础上,求得大损伤片段替换小损伤片段后的循环次数,从而缩短单次循环的里程.大损伤片段循环次数的求解按下式进行:E =D t a r -D o pt D o pt =B K D t a r =r 1r 2⋮r méëêêêêêùûúúúúú㊀B =b 11b 21 b n 1b12b 22 b n 2⋮⋮⋮b 1mb 2mb n m éëêêêêêùûúúúúúK =k 1,k 2, ,k n ()T㊀E =(e 1,e 2, ,e m )T 式中,D t a r 为目标雨流矩阵的伪损伤值;D o p t 为各片段组合优化雨流矩阵的伪损伤值;E 为D t a r 和D o p t 之间的误差值矩阵;B 为分割后各路面片段的损伤值矩阵;K 为经过缩减后各路面片段的循环次数矩阵;r j 为经过雨流外推后第j 个通道的目标损伤值,j =1,2, ,m ;b i j 为第i 个片段的第j 个通道的损伤值,i =1,2, ,n ;k i 为第i 个片段的循环次数,i =1,2, ,n ;e j 为第j 个通道的误差值;n 为路面片段数量;m 为载荷谱数据的通道数.6051 中国机械工程第31卷第12期2020年6月下半月Copyright©博看网 . All Rights Reserved.使用T e c W a r e软件对各路面片段进行组合优化,使得误差的平方值S E最小[11],即m i n S ES E=E T E=ðm j=1(r j-ðn i=1k i b i j)2同时满足以下约束条件:k iɪk m i n i,k m a x i[]式中,S E为误差的平方值;k m i n i为最小循环次数;k m a x i为最大循环次数.各片段循环次数的结果如表2所示.由表2可知,4个片段的比利时路循环次数为0,在零部件的疲劳分析中可直接省略这些片段.由表3可知,单次循环时间由378s缩短至210s,单次循环路面里程由2.61k m缩减至1.49k m.对耐久规范中的越野路及山路同样进行了里程缩减处理,缩减前后单个循环的时间和里程对比如表3所示,缩短了整车虚拟迭代的里程,加速效果显著.表2㊀各路面片段循环次数T a b.2㊀R e p e a t c o u n t o f e a c h r o a d s e c t i o n表3㊀缩减前后路面单个循环时间和里程对比T a b.3㊀S i n g l e l o o p r o a d t i m e a n dm i l e s t o n e c o m p a r i s o nb e f o r e a n da f t e r r e d uc t i o n路面名称缩减前缩减后时间(s)里程(k m)时间(s)里程(k m)比利时路3782.612101.49越野路5502.853601.86山路98511.514505.333.3㊀原始谱与加速谱对比分析对路面进行时域信号分割㊁损伤值比较及片段组合优化处理后,得到了用于虚拟迭代的加速谱.为避免加速谱改变原始谱的损伤值㊁幅值和频率域特征,需对原始谱和加速谱进行对比分析,目前在行业内主要从损伤域㊁幅值域和频率域进行对比分析[12G13].由于篇幅限制,本文仅对比利时路面的原始谱和加速谱进行比较.(1)损伤域对比.对载荷谱进行编制的首要原则就是保证缩减前后的损伤等效.损伤域对比主要从损伤值保留比例[14]D R P(d a m a g e r e t e n t i o n p e r c e n t a g e,D R P)和各通道的相对损伤系数[14] R D R(r e l a t i v e d a m a g e r a t i o,R D R)来评价:D R P=ðm j=1ðn i=1k i b i jðm j=1r jR D Rj=ðni=1b i jr j㊀㊀j=1,2, ,m表4所示为原始谱和加速谱轮心三向力各通道的伪损伤值,经计算,损伤值保留比例为105%,满足90%~110%的接受范围[14],缩减前后各通道的相对损伤系数均在1.0左右,满足0.5~2.0的接受范围[14].为保证缩减结果的有效性,对轴头加速度㊁减振器位移等响应信号同样计算了其损伤值保留比例和各通道的相对损伤系数,结果如表5所示,损伤值保留比例为103%,各通道的相对损伤系数均在1.0左右.表4和表5的计算结果表明该方法使得加速谱完整地保留了原始谱的损伤值.表4㊀原始谱与加速谱轮心三向力的相对损伤系数T a b.4㊀R e l a t i v e d a m a g e r a t i o o fw h e e l c e n t e r t h r e eGc o m p o n e n tf o r c e s o f o r ig i n a l s p e c t r u ma n da c c e l e r a t e d s p e c t r u m通道名称原始谱伪损伤值加速谱伪损伤值相对损伤系数左前轮纵向力F X11.50541.50941.0027左前轮侧向力F Y10.07390.07020.9499左前轮垂向力F Z116.946217.09591.0088右前轮纵向力F X21.50411.54281.0257右前轮侧向力F Y20.07670.07200.9387右前轮垂向力F Z215.009814.44250.9622左后轮纵向力F X32.99613.24481.0830左后轮侧向力F Y30.13520.13551.0022左后轮垂向力F Z3224.6851231.39231.0298右后轮纵向力F X43.18893.38661.0619右后轮侧向力F Y40.12560.12350.9832右后轮垂向力F Z4210.0618227.49931.0830表5㊀原始谱与加速谱响应信号的相对损伤系数T a b.5㊀R e l a t i v e d a m a g e r a t i o o f r e s p o n s e s i g n a l s o f o r i g i n a l s p e c t r u ma n da c c e l e r a t e d s p e c t r u m通道名称原始谱伪损伤值加速谱伪损伤值相对损伤系数左前轴头垂向加速度9.67ˑ10-131.05ˑ10-121.0913右前轴头垂向加速度7.73ˑ10-137.08ˑ10-130.9153左后车桥垂向加速度4.51ˑ10-134.84ˑ10-131.0729右后车桥垂向加速度3.54ˑ10-133.38ˑ10-130.9542左前减振器位移6.56ˑ10-106.53ˑ10-100.9948右前减振器位移5.20ˑ10-105.10ˑ10-100.9817左后减振器位移3.77ˑ10-103.80ˑ10-101.0083右后减振器位移6.66ˑ10-107.01ˑ10-101.0531左前扭杆弹簧力矩7.58ˑ10-47.57ˑ10-40.9986右前扭杆弹簧力矩5.11ˑ10-45.05ˑ10-40.9882前稳定杆扭矩6.691ˑ10-36.342ˑ10-30.9478左侧转向横拉杆力8.28008.42581.0176右侧转向横拉杆力8.99309.40391.04567051一种加速多体动力学模型虚拟迭代的载荷谱编制方法 葛文韬㊀黄㊀晖㊀刘㊀洲等Copyright©博看网 . All Rights Reserved.㊀㊀(2)幅值域对比.幅值域从雨流计数和累积循环次数两个方面进行对比分析,左后轮垂向力雨流计数对比如图7a 和图7b 所示,缩减前后雨流矩阵的形状基本相同.左后轮垂向力累积循环次数对比如图7c 所示,可以看出原始谱与加速谱的累积循环次数基本吻合,说明该方法使得加速谱完整地保留了原始谱的幅值域特征.(a)原始谱雨流(b)加速谱雨流(c)累积循环次数对比图7㊀原始谱与加速谱雨流计数及累加循环次数对比F i g .7㊀R a i n f l o wa n d c u m u l a t i v e c y c l e c o m pa r i s o nb e t w e e n o r i g i n a l s p ec t r u ma n da c c e l e r a t ed s pe c t r u m ㊀㊀(3)频率域对比.加速谱不仅要满足缩减前后的损伤等效,而且也要满足频率域上的等效.不同零部件对不同频率的敏感程度不一致,不能由于载荷谱缩减处理而改变其频率,进而改变零部件的失效模式.频率域从功率谱密度(p o w e r s pe c Gt r a l d e n s i t y,P S D )进行对比分析.原始谱中车速通道的平均值为25k m /h ,构成加速谱的路面片段1㊁2㊁4㊁5㊁6的车速依次为26k m /h ㊁24.79k m /h㊁24 43k m /h ㊁24.51k m /h 和24.04k m /h ,均值为24.75k m /h .加速谱各路面片段的时域信号按照表2所示的循环次数进行拼接处理,原始谱的时域信号按照耐久规范中的循环次数进行拼接处理,然后分别求取原始谱和加速谱的P S D .图8所示为原始谱和加速谱的左前轮心垂向力P S D 分析对比,可看出能量主要集中在在1.1H z 和11H z,分别表征了簧上质量跳动的低频成分和簧下质量跳动的高频成分[15].在各个频率上两种载荷谱的P S D 分布趋势几乎一致,且具有相同形状的P S D 分布曲线,表明了加速谱没有改变原始谱的频率域特征,实现了缩减前后的频率域特征的等效.由于加速谱删除了部分路面片段,所以加速谱在各频率上的平均能量较原始谱高,导致加速谱的P S D 曲线较原始谱的P S D 曲线整体上移[16].图8㊀垂向力原始谱与加速谱P S D 对比F i g .8㊀V e r t i c a l f o r c eP S Dc o m p a r i s o nb e t w e e no r i gi n a l s p e c t r u ma n da c c e l e r a t e d s pe c t r u m 4㊀缩减效果验证虚拟迭代方法[17]以易测位置处的信号(如减振器位移信号)作为目标信号,通过计算图9所示的整车多体动力学模型的传递函数及其逆函数,反求得到其激励信号,使得模型仿真结果接近载荷谱采集的目标信号,最终获得下摆臂和转向节分别在原始谱和加速谱作用下接附点的载荷谱,然后分别计算下摆臂和转向节在两种载荷谱作用下的疲劳寿命.图10所示为第5片段比利时路左侧下摆臂前点与车架连接处X ㊁Y ㊁Z 方向的时域力信号.8051 中国机械工程第31卷第12期2020年6月下半月Copyright©博看网 . All Rights Reserved.图9㊀刚柔耦合整车多体动力学模型F i g .9㊀R i g i d Gf l e x i b l e c o u pl e dM B Sm o d e l o f f u l l v e h i c le 图10㊀比利时路左侧下摆臂前点与车架连接处载荷F i g.10㊀L e f t l o w e r a r mi n t e r n a l f o r c ew i t h f r a m e u n d e rB e l gi a nb l o c k ㊀㊀下摆臂和转向节的疲劳分析需要单位应力结果㊁接口点时域载荷以及材料的特性曲线.采用惯性释放法,在下摆臂和转向节各接口点施加单位载荷,利用有限元软件计算下摆臂和转向节在单位载荷下的应力结果.由于车辆道路耐久试验属于高周疲劳工况,故本文应用应力寿命分析方法(即S GN 法),采用雨流计数法对应力幅值进行计数统计,应用M i n e r 线性疲劳损伤法[18]分别计算下摆臂和转向节在原始谱和加速谱作用下的疲劳损伤.图11所示为下摆臂在原始谱和加速谱作用下的疲劳分析结果.图12所示为转向节在原始谱和加速谱作用下的疲劳分析结果.为更好地对比疲劳分布区域和损伤值大小,在下摆臂和转向节中各选取3个区域进行对比分析,下摆臂3个区域的节点号分别为16997㊁11263和10206,其中损伤值最大位置出现在节点16997处,详细结果如图11所示.对比结果如表6所示,所选取的3个节点在原始谱和加速谱作用下的疲劳损伤比值范围为1.005~1.081.转向节用于疲劳损伤对比的节点号分别为176327㊁6378和141395,其中176327节点损伤(a)原始谱(b)加速谱图11㊀下摆臂在原始谱和加速谱作用下的疲劳分析F i g .11㊀L o w e r a r mf a t i g u e a n a l ys e s u n d e r o r i g i n a l a n da c c e l e r a t e d s p e c t r u m 表6㊀下摆臂疲劳分析结果对比T a b .6㊀C o m p a r i s o no f f a t i g u e a n a l ys i s r e s u l t s o f l o w e r a r m 对比节点损伤值原始谱加速谱比值169970.0756710.0817771.081112630.0515320.0523061.015102060.0487730.0490341.005值最大,详细结果如图12所示.对比结果如表7所示,所选取的3个节点在原始谱和加速谱作用下的疲劳损伤比值范围为1.156~1.205.综上所述,下摆臂和转向节在原始谱和加速谱作用下,疲劳损伤分布趋势一致且损伤值也很接近,表明加速谱能够替代原始谱用于加快虚拟迭代的进程.最后,按照道路耐久规范分析了下摆臂和转向节的疲劳寿命,图13a 所示为下摆臂的分析结果,疲劳风险集中在节点10206附近,损伤值为0 214.图13b 所示为转向节的分析结果,疲劳风险集中在节点140416附近,损伤值为0.075.两者的损伤值均小于设计目标值1 0,表明其寿命能够满足疲劳性能目标.下摆臂和转向9051 一种加速多体动力学模型虚拟迭代的载荷谱编制方法葛文韬㊀黄㊀晖㊀刘㊀洲等Copyright©博看网 . All Rights Reserved.(a)原始谱(b)加速谱图12㊀转向节在原始谱和加速谱作用下的疲劳分析F i g .12㊀K n u c k l e f a t i g u e a n a l y s e s u n d e r o r i gi n a l a n da c c e l e r a t e d s p e c t r u m 表7㊀转向节疲劳分析结果对比T a b .7㊀C o m p a r i s o no f f a t i g u e a n a l ys i s r e s u l t s o f k n u c k l e 对比节点损伤值原始谱加速谱比值1763270.0063920.0076331.19463780.0006000.0007231.2051413950.0019320.0022331.156节在后期实车道路耐久试验中均没有出现裂纹及磨损等失效问题,如图14所示,通过了实车道路耐久验证.5㊀结论(1)本文提出了一种路面时域信号分割及路面片段组合优化的载荷谱编制方法,通过增加大损伤值片段的循环次数来替换小损伤片段,实现了长里程路面的里程缩减,比利时路单个循环里程由2.61k m 缩减至1.49k m ,越野路单个循环里程由2.85k m 缩减至1.86k m ,山路单个循环里程由11.51k m 缩减至5.33k m ,加速了多体动力学模型虚拟迭代,且加速谱保留了原始谱损伤值㊁幅值特征和频率特征.(2)下摆臂及转向节在比利时路原始谱和加(a)下摆臂(b)转向节图13㊀疲劳分析结果F i g .13㊀F a t i g u e a n a l ys i s r e s u l t s (a)下摆臂㊀㊀(b)转向节图14㊀实车道路耐久试验结果F i g .14㊀V e h i c l e r o a dd u r a b i l i t yt e s t r e s u l t s 速谱两种载荷作用下的疲劳损伤分布一致且损伤值接近,下摆臂和转向节的最大损伤比值分别为1.081和1.205,表明了加速谱能够替代原始谱且加快了虚拟迭代的进程.运用此方法对耐久规范中的越野路及山路等长里程路面同样进行了里程缩减处理.(3)对摆臂和转向节进行了疲劳分析预测,结果显示其寿命能够满足疲劳性能目标,并通过了实车道路耐久验证.所述载荷谱编制方法为加速多体动力学模型的虚拟迭代提供了一种解决途径,具有较高的工程应用价值.参考文献:[1]㊀刘俊,刘亚军,张少辉,等.基于虚拟迭代及有限元理论的某中型货车驾驶室疲劳寿命研究[J ].中国机械工程,2018,29(3):1588G1595.L I UJ u n ,L I U Y a j u n ,Z H A N GS h a o h u i ,e t a l .F a t i gu e0151 中国机械工程第31卷第12期2020年6月下半月Copyright©博看网 . All Rights Reserved.A n a l y s i s o f aM e d i u mT r u c kC a bB a s e d o nV i r t u a l I t e rGa t i o n a n dF i n i t eE l e m e n tT h e o r y[J].C h i n aM e c h a n i c a lE n g i n e e r i n g,2018,29(3):1588G1595.[2]㊀V E MU R IA,T A L E K A R N,A V U T A P A L L I B.R o a dL o a d sf o rD u r a b i l i t y A n a l y s i s U s i n g V i r t u a lI t e r a t i o n s[J].S A ET e c h n i c a l P a p e r,2018G01G0567.[3]㊀吴泽勋,张林波,孟凡亮,等.基于虚拟迭代的轿车车身耐久性虚拟试验方法[J].计算机辅助工程,2014,23(6):37G40.WUZ e x u n,Z HA N GL i n b o,M E N GF a n l i a n g,e t a l.V i r t u a lT e s t M e t h o do fC a rB o d y D u r a b i l i t y B a s e do nV i r t u a l I t e r a t i o n[J].C o m p u t e rA i d e dE n g i n e e rGi n g,2014,23(6):37G40.[4]㊀于佳伟,郑松林,冯金芝,等.某轿车前副下摆臂服役载荷模拟试验加速方法研究[J].机械工程学报,2016,52(22):112G120.Y UJ i a w e i,Z H E N G S o n g l i n,F E N GJ i n z h i,e ta l.R e s e a r c ho n A c c e l e r a t e d T e s t i n g M e t h o df o rt h eS e r v i c eGs i m u l a t i o n F a t i g u e T e s t f o r A u t o m o t i v eF r o n tS u bGf r a m e[J].J o u r n a lo f M e c h a n i c a lE n g iGn e e r i n g,2016,52(22):112G120.[5]㊀L O N D H E A,K A N G D ES,K A R T H I K E Y A N K.D e r i v i n g t h e C o m p r e s s e d A c c e l e r a t e d T e s t C y c l ef r o m M e a s u r e dR o a dL o a dD a t a[J].S A ET e c h n i c a lP a p e r,2012G01G0063.[6]㊀钱立军,吴道俊,杨年炯,等.基于室内道路模拟技术的整车加速性试验的研究[J].汽车工程,2011,33(2):91G96.Q I A NL i j u n,WU D a o j u n,Y A N G N i a n j i o n g,e t a l.A R e s e a r c ho nA c c e l e r a t e dD u r a b i l i t y T e s tB a s e d o nI n d o o rR o a dS i m u l a t i o nT e c h n o l o g y[J].A u t o m o t i v eE n g i n e e r i n g,2011,33(2):91G96.[7]㊀王国军.M S C.F a t i g u e疲劳分析实例指导教程[M].北京:机械工业出版社,2009:39G42.W A N G G u o j u n.M S C.F a t i g u eE x a m p l eG u i d e f o rF aGt i g u e A n a l y s i s[M].B e i j i n g:C h i n a M a c h i n e P r e s s,2009:39G42.[8]㊀吴道俊.车辆疲劳耐久性分析㊁试验与优化关键技术研究[D].合肥:合肥工业大学,2012.WUD a o j u n.R e s e a r c h o n t h eK e y T e c h n o l o g y o fA n a lGy s i sT e s t a n dO p t i m i z a t i o n f o rV e h i c l e F a t i g u eD u r a b i lGi t y[D].H e f e i:H e f e iU n i v e r s i t y o fT e c h n o l o g y,2012.[9]㊀郑松林,周亚捷,冯金芝,等.一种新的考虑伪损伤保留的道路载荷模拟试验加速方法[J].机械强度,2017,39(2):410G416.Z H E N GS o n g l i n,Z H O UY a j i e,F E N GJ i n z h i,e t a l.AN e w S p e c t r u m D e v e l o p m e n t M e t h o d o l o g y B a s e d o nP s e u d oD a m a g e R e s e r v i n g f o rR o a d L o a dS i m u l a t i o nT e s t[J].J o u r n a lo f M e c h a n i c a lS t r e n g t h,2017,39(2):410G416.[10]㊀L E EYL,P A N J,HA T HAWA Y R B.F a t i g u eT e s t i n g a n d A n a l y s i s T h e o r y a n d P r a c t i c e[M].A m s t e r d a m:E l s e v i e r,2004:250G254.[11]㊀Z H A N G H W,G U I LJ,F A NZJ.A N e w M e t h o d t o A c c e l e r a t e R o a d T e s t S i m u l a t i o n o n M u l t iGa x i a lT e s tR i g[J].S A ET e c h n i c a lP a p e r,2017G01G0200.[12]㊀Y UJW,Z H E N G S L,F E N G JZ,e ta l.N e w M e t h o d o l o g y f o rD e t e r m i n a t i o n o f L o a dS p e c t r a f o rt h eV e h i c l eA c c e l e r a t e dD u r a b i l i t y T e s t i n g A s s o c iGa t e d w i t h t h e T i m e C o r r e l a t e d F a t i g u e D a m a g eA n a l y s i sM e t h o d[J].I n t e r n a t i o n a l J o u r n a l o fA uGt o m o t i v eT e c h n o l o g y,2017,18(3):547G560.[13]㊀郑国锋.典型汽车部件载荷谱及加速耐久性编辑与实验方法研究[D].广州:华南理工大学,2017.Z H E N G G u o f e n g.S t u d i e so nL o a dS p e c t r u m a n dA c c e l e r a t e d D u r a b i l i t y T e s t M e t h o d st o T y p i c a lA u t o m o t i v eP a r t s[D].G u a n g z h o u:S o u t h C h i n aU n i v e r s i t y o fT e c h n o l o g y,2017.[14]㊀K UMA RPP,P R A K A A S H J,P A L A N I S AMY K.O p t i m i z a t i o n o f P r o v i n g G r o u n dD u r a b i l i t y T e s tS e q u e n c eB a s e do nR e l a t i v eD a m a g eS p e c t r u m[J].S A ET e c h n i c a l P a p e r,2018G01G0101.[15]㊀张新宇.车辆结构试验场疲劳耐久试验规范及台架试验研究[D].北京:清华大学,2014.Z HA N G X i n y u.R e s e a r c ho nP r o v i n g G r o u n dT e sGt i n g S t a n d a r da n d R i g T e s t i n g o fD u r a b i l i t y P e rGf o r m a n c e o fV e h i c l eS t r u c t u r e[D].B e i j i n g:T s i n gGh u aU n i v e r s i t y,2014.[16]㊀花菲菲,郑松林,冯金芝,等.基于时域损伤编辑法的整车道路模拟试验谱编制[J].机械强度,2015,37(5):964G971.HU AF e i f e i,Z H E N GS o n g l i n,F E N GJ i n z h i,e t a l.S p e c t r u m C o m p l i c a t i o nf o r V e h i c l e R o a dS i m u l aGt i o nT e s t B a s e d o nT i m eC o r r e l a t e dF a t i g u eA n a l yGs i s w i t h D a m a g e E d i t i n g M e t h o d[J].J o u r n a lo fM e c h a n i c a l S t r e n g t h,2015,37(5):964G971.[17]㊀D A N N B A U E R H,M E I S E M,G A T T R I G N E R O,e ta l.I n t e g r a t i o n V i r t u a l T e s t M e t h o d sa n dP h y s i c a lT e s t i n g t o A s s u r e A c c u r a c y a n dt o R eGd u c eE f f o r ta n d T i m e[J].S A E Te c h n i c a lP a p e r,2018G01G0567.[18]㊀王忠校,闫鑫,郭茂林,等.载货汽车驾驶室疲劳仿真方法研究[J].汽车技术,2014(7):54G58.W A N G Z h o n g x i a o,Y A N X i n,G U O M a o l i n,e t a l.R e s e a r c h o nD u r a b i l i t y S i m u l a t i o nA p p r o a c ho fT r u c kC a b[J].A u t o m o t i v eT e c h n o l o g y,2014(7):54G58.(编辑㊀王旻玥)作者简介:葛文韬,男,1989年生,工程师.研究方向为多体动力学载荷提取㊁疲劳试验及仿真分析.EGm a i l:g e w e n t a o0517@126.c o m.1151一种加速多体动力学模型虚拟迭代的载荷谱编制方法 葛文韬㊀黄㊀晖㊀刘㊀洲等Copyright©博看网 . All Rights Reserved.。
摩托车实际道路载荷谱采集及车架疲劳仿真分析
摩托车实际道路载荷谱采集及车架疲劳仿真分析
闫豆豆;徐刚;李先文;涂奎;孙于胤
【期刊名称】《机械科学与技术》
【年(卷),期】2022(41)4
【摘要】摩托车车架的疲劳性能决定着摩托车整车所能达到的最大服务寿命,影响着摩托车行驶安全性。
为此,较为完整的摩托车虚拟耐久性分析过程被引入。
本文在完成实际道路载荷谱采集工作的基础上,利用搭建的摩托车多体动力学模型,将等效后的道路载荷谱通过虚拟迭代再现;对摩托车车架应用准静态叠加法进行了疲劳仿真分析。
结果表明:摩托车多体动力学模型精度较高,车架疲劳薄弱位置预测合理,摩托车虚拟耐久性分析过程基本满足工程应用需求。
【总页数】7页(P619-625)
【作者】闫豆豆;徐刚;李先文;涂奎;孙于胤
【作者单位】同济大学汽车学院;隆鑫通用动力股份有限公司
【正文语种】中文
【中图分类】U483
【相关文献】
1.基于实际道路激励谱的摩托车车架随机响应分析
2.摩托车车架疲劳载荷谱的编制
3.摩托车车架疲劳载荷谱的编制
4.摩托车车架疲劳试验载荷谱的研究与编制
5.基于道路载荷谱的副车架疲劳分析及优化
因版权原因,仅展示原文概要,查看原文内容请购买。
基于疲劳-蠕变载荷等效转换的涡轮盘载荷谱编制及寿命预测
基于疲劳-蠕变载荷等效转换的涡轮盘载荷谱编制及寿命预测肖阳;徐可君;秦海勤;贾明明
【期刊名称】《推进技术》
【年(卷),期】2020(41)10
【摘要】针对雨流计数法在峰谷值提取时进行等值压缩,忽略保载时间的问题进行了改进,提出了一种基于损伤曲线的疲劳-蠕变载荷等效转换方法。
利用非线性疲劳损伤累积函数和损伤等效原则,建立了不同应力水平、不同保载时间下疲劳-蠕变载荷与疲劳载荷之间的等效换算模型。
利用涡轮盘材料试验数据,计算了不同循环加载条件下的等效换算比,得到了其随保载时间的变化规律。
利用改进的雨流计数法,编制了航空发动机高压涡轮盘载荷谱,并将其与寿命-时间分数预测法相结合,得到了涡轮盘剩余寿命。
结果表明,改进的雨流计数法综合考虑了疲劳-蠕变耦合损伤对涡轮盘寿命造成的影响,相比于传统雨流计数法,寿命预测误差降低了15.02%,验证了该方法的有效性。
【总页数】9页(P2316-2324)
【作者】肖阳;徐可君;秦海勤;贾明明
【作者单位】海军航空大学青岛校区航空机械工程与指挥系
【正文语种】中文
【中图分类】V232.3
【相关文献】
1.基于 nsoft的塔式起重机疲劳载荷谱编制及疲劳剩余寿命估算
2.多轴变幅载荷下基于载荷支配模式的缺口件疲劳寿命预测方法
3.装载机驱动桥试验载荷谱编制及整车疲劳寿命预测
4.某轿车摆臂疲劳试验载荷谱编制及寿命预测
因版权原因,仅展示原文概要,查看原文内容请购买。
基于虚拟迭代技术的汽车喇叭支架疲劳分析
基于虚拟迭代技术的汽车喇叭支架疲劳分析王良模;赵野;徐娟;袁刘凯;张汤赟【摘要】针对某汽车喇叭支架避振片断裂问题,采用虚拟迭代技术与有限元相结合的方法进行支架的疲劳寿命分析.在对半车架和驾驶室的有限元模态分析的基础上,获得了半车架和驾驶室的中性文件,采用刚柔耦合建模法建立了驾驶室、车架和虚拟台架多体动力学模型,利用虚拟迭代方法提取了喇叭支架安装点处的载荷谱并对支架进行了疲劳分析.结果表明:支架避振片处最大损伤值为31.356,将发生疲劳破坏,破坏区域与实际断裂部位一致;采取增加一片避振片的改进措施,改进后的喇叭支架避振片最大损伤值远小于1.000,满足疲劳寿命要求;在模型准确的前提下,由虚拟迭代方法获得的载荷谱可有效地应用于汽车零部件的疲劳分析中.【期刊名称】《江苏大学学报(自然科学版)》【年(卷),期】2015(036)002【总页数】6页(P130-135)【关键词】支架;疲劳分析;虚拟迭代;载荷谱;刚柔耦合【作者】王良模;赵野;徐娟;袁刘凯;张汤赟【作者单位】南京理工大学机械工程学院,江苏南京210094;南京理工大学机械工程学院,江苏南京210094;南京理工大学机械工程学院,江苏南京210094;南京依维柯汽车有限公司,江苏南京210028;南京依维柯汽车有限公司,江苏南京210028【正文语种】中文【中图分类】U461.7汽车零部件的疲劳寿命已成为评价汽车产品品质的重要因素之一,零部件的疲劳问题越来越受到各大主机厂的关注.疲劳寿命评估的传统方法是在试验场进行道路试验,这种方法直接准确,但是会消耗大量人力物力[1].CAE(computer aided engineering)方法的基础是获得车身或零部件连接点处可信的载荷谱或者加速度谱.采用多体动力学仿真提取载荷谱时一般是以编制的路面位移输入[2],由于输入非真实的道路载荷谱,提取的载荷缺乏可信度,而且轮胎特性的准确模拟一直是个难题[3],需投入大量资源.利用六分力传感器可以获得较精确的轴头六分力响应信号,但六分力传感器的造价昂贵.如何以低成本获取合理的零部件载荷谱,已成为制约CAE疲劳分析的重要因素.虚拟迭代技术基于实车道路试验采集容易得到的信号,通过建立多体动力学模型迭代求解,可以得到各连接点的力,再进行疲劳寿命预测和结构优化,满足产品可靠性设计目标[4],具有成本低、周期短等优点.目前Magna公司正在推广虚拟迭代技术,国内厂商正逐步将此技术应用于汽车产品设计开发上[5-6].笔者针对某轻卡喇叭支架的避振片断裂问题,对该避振片进行疲劳分析,重点介绍刚柔耦合多体动力学模型的建立及虚拟迭代载荷提取的过程,在此基础上对喇叭支架进行疲劳分析,分析避振片损伤位置与实际断裂位置吻合情况,验证虚拟迭代技术应用在汽车零部件疲劳分析中的有效性.1 虚拟迭代原理简介虚拟迭代的实质是根据已知的响应,迭代出系统的驱动.可以通过试验场实车测试得到相关部件响应信号(悬置下方的测量信号),通过FEMFATLab进行虚拟迭代反求出所加激励处的外部激励,将反求出的外部激励重新激励系统便可提取出所关心部位的载荷.迭代的整个过程:创建白噪声,获得响应,求取传递函数,加载期望信号求解初始驱动,获得虚拟迭代响应信号,响应信号与期望信号对比,重复迭代.文献[7]对此作了比较详尽的阐述.驱动可以是外力或者绝对位移,响应可以是加速度、应变、相对位移或者系统内力.利用白噪声激励系统,获得相应响应,得到系统的传递函数:式中:u0(s)为系统输入;y0(s)为相应的系统响应.从而可得F-1(s),根据试验得到的期望信号yDesired(s)可以获得初始驱动信号:用初始驱动信号激励系统,可以获得相应的响应y1(s).因Adams中建立的样机模型是个非线性系统,而传递函数是线性的,所以y1(s)与yDesired(s)会有一定的误差,需要反复修正驱动信号(激励),使响应信号逼近实测值,最终得到近似等效的激励输入.用F-1(s)继续作用于响应误差(yDesired(s)-y1(s)),用修正参数a加以修正并与初始驱动叠加,以此类推,第n次的迭代可表示为信号比较包含时域曲线比较、功率谱曲线比较和相对损伤值比较等.时域曲线主要比较相位幅值和趋势吻合度;功率谱曲线比较幅值和趋势吻合度;相对损伤值接近1.0.2 虚拟迭代过程2.1 多体动力学模型建立多体动力学模型的精度是影响迭代结果的重要因素,而刚柔耦合多体动力学已经成为车辆系统动力学发展的一个趋势[8].车架的运动形式是否正确将会直接影响迭代结果,为了提高迭代的准确性,采用柔性体建模方法建立车架以及驾驶室的多体动力学模型;同时建立模型过程中对结构进行简化,这一简化要以符合结构主要的力学特性为前提.由于后半部分车架与驾驶室并非直接相连,驾驶室运动的主要直接输入点在前半部分车架,采用半车架模型进行迭代.喇叭的质量约为1 kg,质量相对较轻,对迭代结果影响较小.模型中对喇叭总成进行柔性化,图1a为喇叭支架前4阶模态振型,分别为避振片1阶弯曲(41.4 Hz)、避振片1阶扭转(57.7 Hz)、避振片2阶弯曲(111.3 Hz)和避振片2阶扭转(152.4 Hz),采用固定副连接模拟喇叭支架的螺栓连接.驾驶室形状复杂,有很多钣金件,各零件之间通过焊接、铆接等方式进行连接[9],利用Hypermesh进行网格划分处理.驾驶室中性文件中包含有很多局部模态,图1b为整个驾驶室参与变形的前4阶模态振型,性质分别为左门参与较多的驾驶室弯曲变形(27.9 Hz)、右门参与较多的弯曲变形(28.7 Hz)、1阶扭转(53.5 Hz)以及驾驶室地板前部弯曲变形(60.1 Hz).取车架的前59阶模态,驾驶室取前46阶模态,模型中乘员采用质量点模拟,以刚性单元连接到相应部位.图1 支架及驾驶室前4阶整体振型悬置刚度是影响模型准确度的重要参数,车架与驾驶室之间有6个悬置,悬置的刚度曲线由拉伸、压缩试验获得,通过Gforce函数来模拟衬套的非线性刚度特性. 当所研究的多体动力学模型系统的固有频率分布贴近期望信号的频率分布时迭代过程会更平稳.车架以下部分可视作台架,当迭代发散或者驱动对应的响应出现毛刺时可适当增加台架的质量,提高系统的稳定性,使迭代更加容易进行.根据车辆7自由度模型的简化思想,对模型的4个激励点的z向(垂直方向)、2个y向(侧向)和1个x向(纵向)施加外部激励约束.模型中下方4个激励点总共限制了4个自由度;车架为柔性体,可视为小变形,侧向及纵向的3个激励约束限制了车架沿x,y向的运动以及绕z向的转动.在6个悬置上、下端以建立加速度响应;在喇叭支架连接点处建立提取载荷所需的力响应.完整的多体动力学模型如图2所示,在6个悬置处分别建立相应的Request,其中左前悬置下端的z向、右前悬置下端的x,y,z向、左后悬置下端的z向和右后悬置下端的y,z向这7个通道为迭代信号通道.多体动力学系统迭代模型总共设有7个外部激励点和13个响应点. 图2 车架与驾驶室多体动力学模型2.2 数据采集与处理试验数据曲线由在定远试验场进行的实车路试采集而来.在车架与驾驶室之间的6个悬置点上、下端各布置1个3向加速度传感器,悬置下端传感器采集的信号将作为迭代目标信号(期望信号);悬置上端的测量信号作为监测信号.分别在砂石路、搓板路、鹅卵石路、扭曲路等路面上进行信号采集.最终采集的信号总时长为640 s,如图3所示,对应的总测试路程约为6 km.对采集的信号进行单位转换、将测量信号按256.0 Hz进行重采样、滤波、截取等处理.图3 7通道期望信号所研究对象喇叭支架是安装在车架右前端靠近右前悬置处,右前悬置处的信号迭代结果精确与否将显著影响所要提取的载荷谱的准确度.将左前悬置下端的z向、右前悬置下端的x,y,z向、左后悬置下端的z向和右后悬置下端的y,z向加速度信号筛选出来作为迭代目标信号.2.3 迭代过程及结果评价基于FEMFAT-Lab相关软件进行迭代.将多体动力学求解文件导入,定义输入、输出通道,输出的7条响应顺序应与期望信号的排列顺序相一致.经过创建白噪声、求取模型传递函数、选择期望信号、产生初始驱动和选择迭代步数等过程,即可进行迭代,初步选择迭代步数为8步.图4仅显示了右前悬置下端z向信号在初步的8次迭代结束之后的迭代情况,通过与相应的期望信号在时域内对比可见迭代响应信号与期望信号已经比较接近.各路段以及各过渡路段的迭代信号都在逼近期望信号,未有出现发散现象.其他6个通道的迭代效果也很好,将迭代步数增加至15步.图4 右前悬置下端z向信号第8次迭代结果迭代结果通常采用的评价方法:通过比较迭代出来的响应信号与试验采集信号的时域曲线图可以直观地看到信号迭代的准确程度,从而直观地评价结果的准确性;相对损伤值分布图从损伤的角度反映了迭代出来的响应曲线是否与期望曲线一致;功率谱密度曲线可以从能量分布的角度比较二者的差异.以相对损伤值为例评判此次迭代结果的准确性:图5a为迭代过程中各通道响应信号与期望信号之间的相对损伤值变化趋势图,迭代进行了15次,由图5a可见各通道的相对损伤值趋于平稳并且接近数值1.000,说明了迭代结果的可信度很高;图5b是15次迭代结束之后各通道的迭代响应信号与期望信号的相对损伤值,左前 z 向、右前 x,y,z向、左后z向、右后y,z向7个通道的相对损伤值分别为1.006,1.064,0.999,1.007,1.007,0.987 和1.006.各通道相对损伤都很接近1.000,迭代结果较准确.时域信号越吻合、相对损伤值越接近1.000,说明响应信号与期望信号越吻合,即相应部位的变形越接近实测情况.本研究的对象喇叭支架是安装在车架右前悬置附近,对应的右前悬置3个方向(x,y,z向)的迭代信号与试验获得的期望信号的相对损伤值分别为1.064,0.999和1.007,迭代结果可以接受.图5 相对损伤值分布3 载荷提取与疲劳分析3.1 载荷提取进行疲劳寿命分析首先要了解构件的载荷情况,载荷随时间变化的历程称为载荷谱.结构所承受的真实载荷往往是不规则的,在计算时必须对这些不规则载荷进行处理,得到计算载荷谱,而对这类事件的分析普遍应用循环计数法[10].迭代完成之后,从生成的结果文件中可以直接读取到支架安装点处的载荷.提取喇叭支架连接点处的3个方向的力作为支架疲劳分析所需的载荷谱,如图6所示.对载荷谱进行适当滤波,因低频部分对疲劳计算的影响很小,将0.1 Hz以下的低频信号滤掉.图6 喇叭支架安装点处提取的载荷3.2 疲劳损伤机理及累积损伤模型疲劳失效是荷载反复作用的结果,材料在应力或应变的反复作用下发生了性能变化[11].工程中常用的疲劳分析方法有3种:名义应力法、局部应力应变法和损伤容限法.通常,疲劳过程可以观察到以下4个阶段:裂纹成核、微观裂纹扩展、宏观裂纹扩展和最终断裂.工程应用中,通常将零件在裂纹成核和微观裂纹扩展期间的寿命长度称为裂纹萌生阶段,而将零件在宏观裂纹扩展期间的寿命长度称为裂纹扩展阶段,对从萌生到扩展的过渡阶段无法做出精确的定义.利用道林(Dowling)1998年提出的光滑试样线弹性断裂力学可以估算出裂纹萌生尺寸:式中:ΔSe为疲劳极限应力范围;ΔKth为R=-1时的临界强度因子的范围,R为应力比.零件的损伤可以用趋向可接受最大裂纹长度的裂纹长度累积量来表示.曼森和哈尔福德(Manson and Halford)在1981年提出了裂纹扩展方程:式中:n为达到裂纹长度a时外加载荷的循环次数;a0为初始裂纹长度;Nf为在达到最终断裂的裂纹长度af时外加载荷的循环次数;αf为根据经验确定的指数,可根据式(6)进行计算:累积损伤(D)是瞬时裂纹长度与最终裂纹长度的比值:这个损伤公式表示:当疲劳损伤D大于或等于1(即a≥af)时,就会发生疲劳失效.3.3 疲劳分析及结果对比运用有限元法进行零部件疲劳分析时主要有模态叠加法和惯性释放法2种.惯性释放法的原理相对简单,没有考虑共振的影响;模态法为任何以振动为主或者动态加载的部件的疲劳寿命分析提供了技术分析手段[12].如果结构中某部位的载荷可能激发共振模态,那么线性叠加将可能无法准确地计算由此产生的应力时间历程.选择模态叠加法对支架进行疲劳分析,在HyperWorks中对模型进行前处理,获得模态应力和模态坐标作为耐久性分析的输入,采用E-N法计算喇叭支架避振片的疲劳累积损伤,避振片的材料为65Mn.图7为循环1 400次、对应的道路试验里程约为8 400 km的计算结果,避振片最大损伤值为31.356,远大于1.000,该区域位于避振片上端螺栓孔下方附近位置,将产生疲劳失效,这与喇叭避振片的实际断裂部位一致,验证了虚拟迭代结果的正确性.图7 计算结果与样品实际断裂位置对比4 结构改进根据仿真分析和实际断裂的情况,考虑到避振片属于标准件,单个避振片的结构不可随意改动,因此,通过增加1片避振片来提高支架的使用寿命.用同样的方法对改进后的支架进行疲劳分析,循环1 400次的计算结果显示避振片处最大的损伤值为0.011,循环5 000次的损伤值为0.038,远小于1.000;将改进后的支架安装在整车上,在定远试验场进行疲劳道路试验,结果表明:支架通过了试验,满足疲劳寿命的要求.5 结论1)通过建立柔性车架和驾驶室,提高了多体动力学模型的精确度,使得迭代的结果更加可信;基于提取到的载荷谱得到的避振片疲劳分析结果与实际完全吻合,验证了虚拟迭代结果的正确性.2)研究对象喇叭为典型的汽车挂件之一,使用的方法具备通用性,对于车架上的其他挂件的疲劳载荷提取具有一定的借鉴意义;对于驾驶室疲劳分析所需的载荷可通过提取驾驶室悬置上的力获得.参考文献(References)【相关文献】[1]沈宏杰,周鋐.汽车零部件道路模拟加载谱研究[J].汽车工程,2010,32(2):159-162,142.Shen Hongjie,Zhou Hong.A study on the loading spectra for the road simulation of vehicle components[J].Automotive Engineering,2010,32(2):159-162,142.(in Chinese) [2]Cui An,Liu Zihui,Zeng Qingyang,et al.The layout and fatigue life analysis of welding spots for the cab body in white of a commercial vehicle[C]∥Electronic&Mechanical Engineering and Information Technology.Piscataway,USA:IEEE,2011:2089-2093.[3]Taheri S,Sandua C,Taheri S,et al.A technical survey on Terramechanics modelsfor tire-terrain interaction used in modeling and simulation of wheeled vehicles[J].Journal of Terramechanics,2015,57:1-22.[4]Fischer P,Witteveen W.Intergrated MBS-FE-Durability analysis of truck frame components by modal stresses[C]∥ADAMS User Meeting,2000:1-15.[5]邵建,董益亮,肖攀,等.基于多体模型仿真的载荷谱虚拟迭代技术分析[J].重庆理工大学学报:自然科学,2010,24(12):84-87.Shao Jian,Dong Yiliang,Xiao Pan,et al.Virtual iteration technology of load spectrum based oil multiple body simulation[J].Journal of Chongqing University of Technology:Natural Science,2010,24(12):84-87.(in Chinese)[6]钱立军,吴道俊,杨年炯,等.基于道路模拟激励的汽车下摆臂多轴疲劳分析[J].汽车工程,2012,34(3):249-254.Qian Lijun,Wu Daojun,Yang Nianjiong,et al.Multiaxial fatigue analysis of vehicle lower control arm based on simulated road excitation[J].Automotive Engineering,2012,34(3):249-254.(in Chinese)[7]Helmut D,Matthias M,Otmar G,et al.Integrating virtual test methods and physical testing to assure accuracy and to reduce effort and time[C]∥SAE Technical Paper A:SAE Publication Group,Paper Number:2006-01-3563.[8]Lu Yaohui,Zeng Jing,Wu Pingbo,et al.Modeling of rigid-flexible coupling systemdynamics for railway vehicles with flexible bogie frame[C]∥Fourth International Conference on Innovative Computing,Information and Control,2009:1355-1360.[9]王良模,陈东益,邹小俊,等.某轻型卡车正面A柱撞击仿真与改进[J].江苏大学学报:自然科学版,2014,35(2):132-136,143.Wang Liangmo,Chen Dongyi,Zou Xiaojun,et al.Simulation and improvement of a light truck in A-pillar impact[J].Journal of Jiangsu University:Natural Science Edition,2014,35(2):132-136,143.(in Chinese)[10]高东明,王德成,王光辉,等.草地切根机切刀的疲劳寿命分析[J].江苏大学学报:自然科学版,2012,33(3):283-287.Gao Dongming,Wang Decheng,Wang Guanghui,et al.Fatigue life analysis of root cutting machine blade[J].Journal of Jiangsu University:Natural Science Edition,2012,33(3):283-287.(in Chinese)[11]任伟平,李小珍,李俊,等.公轨两用钢桁桥轨道横梁与整体节点连接头的疲劳荷载[J].中国公路学报,2007,20(1):79-84.Ren Weiping,Li Xiaozhen,Li Jun,et al.Fatigue load of welded connections between rail across beams and integral joints of combined steel truss bridges[J].China Journal of Highway and Transport,2007,20(1):79-84.(in Chinese)[12]Gaier C,Dannbauer H.A multiaxial fatigue analysis method for ductile,semi-ductile,and brittle materials[J].The Arabian Journal for Science and Engineering,2008,33(1):223-235.。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
y Desired 比较,如果满足条件,则迭代结束, un 1 即
满足实际路谱采集时对样车的激励。并可此次的仿 真结果作为迭代的最终结果,提取相应零部件的疲 劳载荷。如果不满足条件,则通过相应信号 yn 1 与 期望信号 y Desired 的差异( yDesied yn1 )以及逆传函 对驱动信号 un 1 进行修正得到驱动信号 u n , 再次驱 动整车样机进行仿真,之道获得 yn 1 逼近 y Desired , 满足判定要求。为了使迭代过程不发散,能够尽快 收敛,首次的驱动信号 u0 通过逆传递 F 和由试验 采集得到的期望路谱 y Desired 获得,如下式:
图 4 路谱信号
的轮心处的外部激励 (驱动信号, 如轮心垂向位移、 水平向力等) 。将整个整车样机模型看作一个系统, 路面在轮心处的激励作为输入,各个传感器的响应 作为输出。利用拉普拉斯变换可建立系统的输入、 输出关系,即传递函数。然后对传递函数求逆得到 逆传递函数,作为由输出反求得到输入的关系。 Adams 建立的整车样机模型是一个非线性系统,而 传递函数是线性的,需要反复修正驱动信号(激 励) ,使响应信号逼近实测值,最终得到准确的激 励输入,这就是虚拟迭代的原理。
是建立输出相应信号通道。在底盘及车身各个信号 采集点建立相标位置应与各个传感器的安装零件 及坐标位置一致。
图 4 样车 Adams 疲劳多体模型图
将经过修正和调试正确的整车模型导出为一 定格式的文件(*.adm) ,用于虚拟迭代时的调用。 3. 道路载荷谱虚拟迭代原理 道路载荷谱虚拟迭代的思想是通过试验场实 车测试得到相关部件响应信号(如弹簧位移、轮心 加速度等) ,通过迭代反求技术获得难以实测获取
1
图 8 Grave 和 Cobblestone 路面各个响应通道相对损伤
参考文献
1、沙大亮,极限工况多体动力学仿真及载荷分解作业指导书 长安汽车研究院 CAE 工程所 2、陈军,MSC.ADAMS 技术与工程分析实例所
u0 F 1 * y Desied
4、迭代结果评价:
(3)
作者简介:沙大亮,男,出生于 1980 年 10 月,硕士,重庆大学机械设 计及理论专业,长安汽车工程研究总院 CAE 工程所。研究方向:道路谱 疲劳载荷分析,车辆多体动力学仿真及载荷分析,联系电话
迭代信号与路谱实测值的比较包括三个方面,
13638342625,邮箱:sdl29@
基于道路谱的多体仿真迭代疲劳载荷技术
作者:沙大亮 卲建 陈建华 王俊翔
单位:重庆长安汽车工程研究院 CAE 工程所
摘要: 对于疲劳载荷谱的获取一直是疲劳分析中极为重要的一个部分, 直接决定疲劳分析的结果是否准确。 本文介绍一种实用性强、精度高、经济可靠的疲劳载荷获取方法。这种技术是用应标定技术获取真实道路 谱,通过驱动 Adams 建立的车辆虚拟样机进行迭代,使车辆虚拟模型的相应通道逼近实测道路谱信号,从 而实现虚拟样机重现实际道路过程, 进而获取车身及其零部件道路谱疲劳载荷, 为疲劳分析提供载荷输入。 关键字:道路谱 ADAMS 迭代 疲劳载荷
2 整车虚拟样机建立 整车虚拟样机的质量直接影响到能否获取可 信了车身及底盘零部件的疲劳载荷。整车模型与实 际样车差距较大将导致迭代所需的传递函数的不 能反映真实车辆状态,线性性也不能满足要求,甚 至会导致仿真迭代不能正常进行。首先,根据样车 底盘参数(硬点、零部件质量、弹簧刚度、减振器 阻尼以及衬套的刚度数据等) 和配重数据, 用 Adams 建立该车的多体动力学模型。模型中不含有轮胎, 这是由于迭代的位移和力均作用在轮心处,作为整 车的驱动信号,直接驱动整车模型运动。其次在道 路测试之前,应该测量车辆的轮胎载荷、重心位置 和悬架限位器间隙等参数,并根据测量结果对模型 进行修正。第三是建立驱动信号通道。在多体模型 的轮心处建立垂向的位移 (Point Motion)驱动作为 垂向的位移驱动信号(如多通道的还需建立横向和 纵向的力驱动) 。位移数据以 Spline 形式通过函数 输入,例如 CUBSPL(time,0,SPLINE_DZ_FL, 0)。最后
图 1 道路谱多体虚拟迭代流程
获取整车及零部件真实状态的疲劳载荷。以此作为
图 2 标定试验
图 3 信号标定曲线
道路谱采集的数据是在现场动态测试中通过 数据采集设备对测试仪器检拾到的动态模拟信号 进行离散化、数字化、转换和变换后得到的一列具 有等时间间隔特征的数字信号。为了满足载荷分解 的需要,应对路谱采集试验获得的各类数据进行科 学的变换处理,包括数据有效性处理、信号处理、 数据缩减。
其中:
(2)
图 6 GRAVE 路面悬架位移时域、频域对比
F 1 ——逆传函,由(1)式求逆获得。
yDesied ——期望信号,即采集信号。 yn 1
——第 n-1 次的输出信号。
u n、 u n1 ——第 n-1、n 次的驱动信号。
总的过程是,第 n-1 次的驱动信号驱动整车样 机获得第 n-1 次的输出相应信号 yn 1 ,将 yn 1 与
该技术是通过标定样车进行道路试验,采集得 到的路谱信号(如加速度、位移、应变等) ,利用 多体动力学软件(Adams)建立车辆虚拟样机,然 后通过迭代软件 FemFat-Lab 调用 ADAMS 的求解器 对样机模型进行迭代计算,并通过各个信号通道控 制,使车辆虚拟样机状态逼近样车实际状态,进而
输入进行疲劳寿命预测和结构优化,以满足产品可 靠性设计目标值。该项技术涉及道路谱采集处理、 整车虚拟样机的建立、道路谱虚拟迭代三部分技 术。 1 道路谱采集及处理 载荷谱采集使用力标定法,通过标定技术确定 零部件的外来与应变之间的线性关系,建立力分量 与应变之间的标定系数矩阵。由于零部件结构形式 及采集外力分量的不同,存在各外力分量对各向应 变之间相互干扰,这样就需要要进行解耦。应变片 的粘贴位置以及组桥形式就需要相应的发生变化, 这是该测量方法技术能力的主要体现之一。
图 5 虚拟迭代原理图
如图 5 所示,首先生成噪声信号,以噪声信号
作为输入通道信号,根据迭代通道不同,噪声也分 为四通道和多通道两种。将噪声输入多体动力学模 型驱动模型进行仿真求解,获取相应输出通道信 号。然后通过拉式变换获得车辆多体动力学模型的 传递函数,如下面公式。 传递函数:
时间域、功率谱密度和相对损伤值。前两项是在时 域和频域对相应信号的幅值吻合度进行评价,主要 是通过曲线趋势及峰值对比进行主观评价 (如图 6、 7) 。这里的“相对损伤值”指标,是忽略了 S-N 曲 线、表面质量、加工质量等与部件相关的性能的影 响(因为分析所用的是同一个部件) ,而只是对载
F ( s) y ( s) / u ( s)
(1)
荷的疲劳性能进行比较。该方法在这里只是对相应 信号和期望信号进行相对比较。 经验表明时间域和功率谱方面吻合得不够好 则相对损伤值会有较大的偏差,只有在时间域和功 率谱方面都吻合得较好,才可能保证相对损伤值满 足要求,一般损伤值取应接近 1(如图 8) 。
输入的载荷类型基本是固定的(如四通道的轮 心的垂向位移) ;在选取输出项(即用于迭代的信 号)时需要选取与输入相关性较好的响应信号,如 弹簧变形和轮心垂向加速度,与输入信号在垂向都 有较好的相关性,利于虚拟迭代。 传递函数生成后,就可以传递函数的逆建立驱 动信号和输出信号之间的关系,如下式表示:
un un1 F 1 * ( y Desied yn1 )