CN109977459A - 一种应用cae技术提高强化载荷谱精度的方法 - Google Patents

一种应用cae技术提高强化载荷谱精度的方法 Download PDF

Info

Publication number
CN109977459A
CN109977459A CN201910109552.9A CN201910109552A CN109977459A CN 109977459 A CN109977459 A CN 109977459A CN 201910109552 A CN201910109552 A CN 201910109552A CN 109977459 A CN109977459 A CN 109977459A
Authority
CN
China
Prior art keywords
components
spectrum
loading spectrum
grade
load
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.)
Granted
Application number
CN201910109552.9A
Other languages
English (en)
Other versions
CN109977459B (zh
Inventor
康一坡
曹正林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun Automotive Test Center Co ltd
FAW Group Corp
Original Assignee
FAW Group Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by FAW Group Corp filed Critical FAW Group Corp
Priority to CN201910109552.9A priority Critical patent/CN109977459B/zh
Publication of CN109977459A publication Critical patent/CN109977459A/zh
Application granted granted Critical
Publication of CN109977459B publication Critical patent/CN109977459B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

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)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

一种应用CAE技术提高强化载荷谱精度的方法,包括下述步骤:一、编辑道路载荷谱,获得多级载荷谱L0;二、建立零部件有限元模型;三、有限元分析获得零部件应力场;四、将应力场、材料力学参数、多级载荷谱L0导入到疲劳软件中,计算零部件损伤D0;五、将多级载荷谱L0替换为1级载荷谱L1,再次计算零部件损伤D1;六、若D1=D0,则L1与L0等价,L1即为1级强化载荷谱;否则调整循环次数为获得1级强化载荷谱L2,1级强化载荷谱应用于台架试验和CAE仿真计算。该方法采用CAE技术,实现了材料S‑N曲线向零部件S‑N曲线转化,解决了零部件S‑N曲线的斜率K值取经验值造成1级强化载荷谱精度较差的问题,同时避免试验测试零部件S‑N曲线,节约了产品开发成本。

Description

一种应用CAE技术提高强化载荷谱精度的方法
技术领域
本发明属于载荷谱测试与编辑技术领域,具体涉及一种应用CAE技术提高强化载荷谱精度的方法。
背景技术
强化载荷谱主要应用于产品开发过程中的台架试验和CAE仿真计算。强化载荷谱的级数越多,试验和计算的周期越长,反之则越短。因此,从试验和计算周期上考虑,强化载荷谱级数越少越好,最好是1级,其对于快速完善结构方案具有重要意义。
载荷谱强化主要用到零部件应力-寿命曲线的斜率K。目前K值获取有两种方法,一是根据经验直接取值,通常取K=5;二是通过试验做出零部件应力-寿命曲线,即S-N曲线,进而获得K值。第一种方法多用于载荷级数较多情况,例如1000级,当载荷谱级数较少时,经验值K则会导致载荷谱精度较差,造成产品过重,不能实现产品的轻量化设计。第二种方法虽然可以用于载荷谱强化为1级的情况,但零部件应力-寿命曲线获取耗费大量人力、物力,试验周期长,且获得的K仅适用于所试验的零部件,对于其他零部件不具有通用性,所以通过试验获取真实零部件的K值意义不大。
目前强化载荷谱遇到的主要问题是,零部件S-N曲线的斜率K值制约着多级载荷谱向1级载荷谱强化时的精度。
发明内容
本发明的目的是提供一种应用CAE技术提高强化载荷谱精度的方法,为获得高精度1级强化载荷谱提供技术支撑。
为实现上述目的,本发明提供一种应用CAE技术提高强化载荷谱精度的方法,包括下述步骤:
步骤一、编辑道路载荷谱,包括去毛刺,删掉小幅值载荷等,然后统计分析载荷谱,获得多级载荷谱L0
步骤二、建立载荷谱所应用零部件的装配有限元模型,施加边界条件及载荷,给定材料弹性模量E和泊松比μ;
步骤三、完成零部件的有限元分析,获得结构应力场;
步骤四、将结构应力场、材料力学参数(包括材料S-N曲线)、多级载荷谱L0导入到疲劳计算软件中,计算零部件的损伤D0
步骤五、多级载荷谱L0替换为1级载荷谱L1,L1包括三个参数:循环次数n、幅值A、均值M,再次计算零部件损伤为D1
步骤六、若D1=D0,则L1与L0等价,L1即为1级强化载荷谱;否则调整循环次数为获得1级强化载荷谱L2,L2包括三个参数:循环次数幅值A、均值M,1级强化载荷谱应用于台架试验和CAE仿真计算。
在所述步骤一中,多级载荷谱L0级数可达1000级以上。
在所述步骤二中,零部件有限元模型应力较大位置应细化网格,以提高仿真精度。搭建的零部件装配有限元模型至少包括施加载荷部件、施加边界条件部件及研究对象。
在所述步骤四中,得到零部件损伤D0同时,也可以获得零部件的S-N曲线,及其斜率值。
步骤五中,所述1级载荷谱L1包括三个参数:循环次数n、幅值A、均值M。循环次数初始值n设为1000000,根据道路载荷谱确定幅值A、均值M,选取道路载荷谱历程中的最大值F1和最小值F1-1,当F=F1-F1-1作用到零部件有限元模型上时,如果零部件应力小于材料强度极限,则A=(F1-F1-1)/2,M=(F1+F1-1)/2;否则,选取道路载荷谱历程中的第二最大值F2和第二最小值F2-1,当F=F2-F2-1作用到零部件有限元模型上,如果零部件应力小于材料强度极限,则A=(F2-F2-1)/2,M=(F2+F2-1)/2;否则,再选取道路载荷谱历程中的第三最大值和第三最小值,并计算其差值F作用到零部件有限元模型时,判断零部件应力是否小于材料强度极限;以上过程需要重复计算至零部件应力小于材料强度极限为止;最终根据道路载荷谱确定1级载荷谱的幅值A、均值M。
本发明利用CAE技术,以及材料S-N曲线的优点,将二者充分结合,从而实现了材料S-N曲线向零部件S-N曲线转化,解决了零部件S-N曲线的斜率K值取经验值造成的1级强化载荷谱精度较差的问题,同时避免因试验测试零部件S-N曲线耗费大量人力、物力、财力,节约了产品开发成本。
附图说明
图1为道路载荷谱;
图2为采用雨流计数获得的多级载荷谱;
图3为零部件的有限元模型;
图4为零部件应力场;
图5为材料应力-寿命曲线;
图6为应用多级载荷谱计算的零部件损伤;
图7为零部件的应力-寿命曲线;
图8为1级强化载荷谱计算的零部件损伤。
具体实施方案
为了使发明要解决的技术问题、技术方案和优点更加清楚,进一步说明本发明的技术方案及优点,下面结合附图及实施例进行详细描述。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本发明提供了一种应用CAE技术提高强化载荷谱精度的方法,包括下述步骤:
一、编辑道路载荷谱,如图1所示,包括去毛刺,删掉小幅值载荷等,然后进行载荷雨流计数统计分析,获得多级载荷谱L0,如图2所示;
二、建立零部件有限元模型,如图3所示,约束零部件的轴向运动以及零部件的轴向旋转运动,根据道路载荷谱确定施加零部件的径向载荷F=10000N,给定材料弹性模量E=2.1e5MPa、泊松比μ=0.3;
三、应用ABAQUS软件进行零部件的有限元分析,获得结构应力场,如图4所示,应力单位为MPa;
四、将结构应力场(如图4)、材料力学参数(表1、如图5)、多级载荷谱L0(如图2)导入到疲劳计算软件FEMFAT中,计算零部件最大损伤为D0=0.59,如图6所示,同时获得零部件的S-N曲线,如图7所示;
表1材料力学参数 单位:MPa
五、依据道路载荷谱经过分析计算,L1为1级载荷谱,包括三个特征参数:循环次数n=1000000、幅值A=5000N、均值M=500N,F=A+M=10000N,检查结构应力最大为62.7MPa,其远小于表1中的材料强度极限500MPa,符合要求;多级载荷谱L0替换为1级载荷谱L1,再次计算零部件损伤,如图8所示,此时计算的损伤为D1=0.17;
六、比较D0和D1,D1≠D0,为此调整循环次数为即得1级强化载荷谱L2,L2包括三个参数:循环次数n=3470588、幅值A=5000N、均值M=5000N,L2应用于台架试验和CAE仿真计算。
所述方法应用CAE技术实现了多级载荷谱向1级载荷谱强化,实现快,精度高,该方法对于完成零部件CAE计算、零部件试验及其结构改进具有重要意义。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (5)

1.一种应用CAE技术提高强化载荷谱精度的方法,其特征在于,包括下述步骤:
步骤一,编辑道路载荷谱,然后统计分析载荷谱,获得多级载荷谱L0
步骤二,建立零部件的装配有限元模型,施加边界条件及载荷;
步骤三,完成零部件的有限元分析,获得结构应力场;
步骤四,将结构应力场、材料力学参数、多级载荷谱L0导入到疲劳计算软件中,计算零部件的损伤D0
步骤五,多级载荷谱L0替换为1级载荷谱L1,L1包括三个参数:循环次数n、幅值A、均值M,再次计算零部件损伤为D1
步骤六,若D1=D0,则L1与L0等价,L1即为1级强化载荷谱;否则调整循环次数为获得1级强化载荷谱L2,L2包括三个参数:循环次数幅值A、均值M,1级强化载荷谱应用于台架试验和CAE仿真计算。
2.根据权利要求1所述的一种应用CAE技术提高强化载荷谱精度的方法,其特征在于,步骤五中,所述1级载荷谱L1包括三个参数:循环次数n、幅值A和均值M,所述循环次数n初始值设为1000000,根据道路载荷谱确定幅值A、均值M,选取道路载荷谱历程中的最大值F1和最小值F1-1
3.根据权利要求2所述的一种应用CAE技术提高强化载荷谱精度的方法,其特征在于,当F=F1-F1-1作用到零部件有限元模型上时,如果零部件应力小于材料强度极限,则A=(F1-F1-1)/2,M=(F1+F1-1)/2;否则,选取道路载荷谱历程中的第二最大值F2和第二最小值F2-1
4.根据权利要求3所述的一种应用CAE技术提高强化载荷谱精度的方法,其特征在于,当F=F2-F2-1作用到零部件有限元模型上,如果零部件应力小于材料强度极限,则A=(F2-F2-1)/2,M=(F2+F2-1)/2;否则,再选取道路载荷谱历程中的第三最大值和第三最小值,并计算其差值F作用到零部件有限元模型时,判断零部件应力是否小于材料强度极限。
5.根据权利要求4所述的一种应用CAE技术提高强化载荷谱精度的方法,其特征在于,需要重复计算至零部件应力小于材料强度极限为止;最终根据道路载荷谱确定1级载荷谱的幅值A、均值M。
CN201910109552.9A 2019-02-11 2019-02-11 一种应用cae技术提高强化载荷谱精度的方法 Active CN109977459B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910109552.9A CN109977459B (zh) 2019-02-11 2019-02-11 一种应用cae技术提高强化载荷谱精度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910109552.9A CN109977459B (zh) 2019-02-11 2019-02-11 一种应用cae技术提高强化载荷谱精度的方法

Publications (2)

Publication Number Publication Date
CN109977459A true CN109977459A (zh) 2019-07-05
CN109977459B CN109977459B (zh) 2022-10-28

Family

ID=67076925

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910109552.9A Active CN109977459B (zh) 2019-02-11 2019-02-11 一种应用cae技术提高强化载荷谱精度的方法

Country Status (1)

Country Link
CN (1) CN109977459B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111914358A (zh) * 2020-07-08 2020-11-10 中国第一汽车股份有限公司 发动机冲击作用下的变速器壳体极限承载能力预报方法
CN111950094A (zh) * 2020-07-01 2020-11-17 燕山大学 一种多载荷疲劳寿命约束的拓扑优化方法
CN111950170A (zh) * 2020-07-03 2020-11-17 中国第一汽车股份有限公司 获得高精度麦弗逊式前悬架转向节台架试验载荷的方法
CN112948975A (zh) * 2021-03-11 2021-06-11 北京理工大学 一种仿真与试验载荷融合的载荷谱分析方法
CN113255081A (zh) * 2021-06-18 2021-08-13 中国第一汽车股份有限公司 一种电动汽车减速器的载荷谱的构建方法
CN113281069A (zh) * 2021-06-16 2021-08-20 中国第一汽车股份有限公司 提高汽车橡胶衬套耐久试验载荷精度的方法、装置及介质

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101393079A (zh) * 2008-11-06 2009-03-25 上海交通大学 汽车车身结构疲劳寿命预测系统
CN101718651A (zh) * 2009-11-18 2010-06-02 上海理工大学 强化和损伤共同作用下剩余强度和寿命的评价方法
CN102680270A (zh) * 2012-06-07 2012-09-19 哈尔滨工程大学 不平衡激励力及不平衡弯矩联合作用时设备对船体激励载荷的间接测量方法
CN103853879A (zh) * 2013-12-20 2014-06-11 淮阴工学院 一种表征组合路况作用下车辆结构疲劳损伤的网状图法
CN104833536A (zh) * 2014-02-12 2015-08-12 大连理工大学 一种基于非线性累积损伤理论的结构疲劳寿命计算方法
CN105718633A (zh) * 2016-01-15 2016-06-29 重庆长安汽车股份有限公司 一种底盘件的载荷分析方法
CN106066921A (zh) * 2016-06-14 2016-11-02 奇瑞汽车股份有限公司 一种整车路面振动噪声的cae仿真预测方法
CN106644464A (zh) * 2016-11-18 2017-05-10 南京工业大学 一种基于载荷谱分析的轧机传动系统关键零部件的疲劳寿命预警方法
CN107505130A (zh) * 2017-09-19 2017-12-22 中国第汽车股份有限公司 车架总成扭转疲劳加载试验台架及试验方法
CN107644149A (zh) * 2017-09-27 2018-01-30 交通运输部公路科学研究所 考虑载荷分布及损伤一致性的多通道均衡关联新方法
CN108427844A (zh) * 2018-03-16 2018-08-21 北京工业大学 考虑温度和随机振动载荷的加筋板结构疲劳寿命计算方法
CN108491627A (zh) * 2018-03-22 2018-09-04 东北大学 一种机械零部件结构的可靠性分析方法
CN109033709A (zh) * 2018-08-30 2018-12-18 电子科技大学 基于非线性疲劳损伤累积理论的构件疲劳寿命评估方法
CN109030016A (zh) * 2018-06-11 2018-12-18 中国北方车辆研究所 基于载荷谱的综合传动装置耐久性评估方法
CN109060497A (zh) * 2018-07-18 2018-12-21 天津大学 一种考虑低于疲劳极限加载的疲劳损伤分析方法

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101393079A (zh) * 2008-11-06 2009-03-25 上海交通大学 汽车车身结构疲劳寿命预测系统
CN101718651A (zh) * 2009-11-18 2010-06-02 上海理工大学 强化和损伤共同作用下剩余强度和寿命的评价方法
CN102680270A (zh) * 2012-06-07 2012-09-19 哈尔滨工程大学 不平衡激励力及不平衡弯矩联合作用时设备对船体激励载荷的间接测量方法
CN103853879A (zh) * 2013-12-20 2014-06-11 淮阴工学院 一种表征组合路况作用下车辆结构疲劳损伤的网状图法
CN104833536A (zh) * 2014-02-12 2015-08-12 大连理工大学 一种基于非线性累积损伤理论的结构疲劳寿命计算方法
CN105718633A (zh) * 2016-01-15 2016-06-29 重庆长安汽车股份有限公司 一种底盘件的载荷分析方法
CN106066921A (zh) * 2016-06-14 2016-11-02 奇瑞汽车股份有限公司 一种整车路面振动噪声的cae仿真预测方法
CN106644464A (zh) * 2016-11-18 2017-05-10 南京工业大学 一种基于载荷谱分析的轧机传动系统关键零部件的疲劳寿命预警方法
CN107505130A (zh) * 2017-09-19 2017-12-22 中国第汽车股份有限公司 车架总成扭转疲劳加载试验台架及试验方法
CN107644149A (zh) * 2017-09-27 2018-01-30 交通运输部公路科学研究所 考虑载荷分布及损伤一致性的多通道均衡关联新方法
CN108427844A (zh) * 2018-03-16 2018-08-21 北京工业大学 考虑温度和随机振动载荷的加筋板结构疲劳寿命计算方法
CN108491627A (zh) * 2018-03-22 2018-09-04 东北大学 一种机械零部件结构的可靠性分析方法
CN109030016A (zh) * 2018-06-11 2018-12-18 中国北方车辆研究所 基于载荷谱的综合传动装置耐久性评估方法
CN109060497A (zh) * 2018-07-18 2018-12-21 天津大学 一种考虑低于疲劳极限加载的疲劳损伤分析方法
CN109033709A (zh) * 2018-08-30 2018-12-18 电子科技大学 基于非线性疲劳损伤累积理论的构件疲劳寿命评估方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
樊宝锋 等: ""一种新的循环荷载作用下考虑损伤累积应力-应变关系"", 《建筑结构》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111950094A (zh) * 2020-07-01 2020-11-17 燕山大学 一种多载荷疲劳寿命约束的拓扑优化方法
CN111950170A (zh) * 2020-07-03 2020-11-17 中国第一汽车股份有限公司 获得高精度麦弗逊式前悬架转向节台架试验载荷的方法
CN111914358A (zh) * 2020-07-08 2020-11-10 中国第一汽车股份有限公司 发动机冲击作用下的变速器壳体极限承载能力预报方法
CN111914358B (zh) * 2020-07-08 2022-07-19 中国第一汽车股份有限公司 发动机冲击作用下的变速器壳体极限承载能力预报方法
CN112948975A (zh) * 2021-03-11 2021-06-11 北京理工大学 一种仿真与试验载荷融合的载荷谱分析方法
CN112948975B (zh) * 2021-03-11 2022-02-22 北京理工大学 一种仿真与试验载荷融合的载荷谱分析方法
CN113281069A (zh) * 2021-06-16 2021-08-20 中国第一汽车股份有限公司 提高汽车橡胶衬套耐久试验载荷精度的方法、装置及介质
CN113255081A (zh) * 2021-06-18 2021-08-13 中国第一汽车股份有限公司 一种电动汽车减速器的载荷谱的构建方法
CN113255081B (zh) * 2021-06-18 2022-09-02 中国第一汽车股份有限公司 一种电动汽车减速器的载荷谱的构建方法

Also Published As

Publication number Publication date
CN109977459B (zh) 2022-10-28

Similar Documents

Publication Publication Date Title
CN109977459A (zh) 一种应用cae技术提高强化载荷谱精度的方法
CN110555229B (zh) 一种无网格固体力学仿真方法、电子设备及存储介质
CN110926973B (zh) 循环荷载下混凝土损伤积累分析方法
CN103366085A (zh) 编织复合材料力学性能的多尺度预测方法
CN109446731B (zh) 一种基于abaqus的岩土工程数值模拟方法
CN111062162B (zh) 一种岩土材料精确本构模型的数值建模与应用方法
CN104915483A (zh) 一种地震作用下建筑物基础稳定性验算方法
CN110516365A (zh) 一种测试螺栓连接刚度的方法
CN114330042A (zh) 基于sn曲线的程序载荷谱编制方法、系统和存储介质
CN106295015B (zh) 一种渐开线直齿圆柱齿轮副的齿廓修形方法及与其配套的专用参数化cad系统
CN110110406B (zh) 一种基于Excel计算平台实现LS-SVM模型的边坡稳定性预测方法
CN110705137A (zh) 一种应力幅值和均值的确定方法以及装置
CN107480355A (zh) 发动机安装接头承载能力验证方法
CN109190208B (zh) 一种张拉整体结构的主动索预应力张拉量确定方法
CN109388833A (zh) 一种基于疲劳寿命的弹性元件结构优化设计方法
CN110245455B (zh) 一种单向碳纤维增强复合材料的设计方法
CN105653772B (zh) 一种pbx炸药颗粒压制成型数值模拟自动化建模方法及系统
CN110688784A (zh) 一种高强双角钢组合截面构件受力性能的有限元分析方法
CN112100873A (zh) 水工建筑物承载力的确定方法
CN105005210B (zh) 机电一体化仿真系统及使用其的方法
CN103942368B (zh) 一种激光切割机床的结构设计方法
CN112307575B (zh) 一种评判装配结果的复合材料壁板零件合格性的检查方法
Sobotka et al. Automatic 3D crack placement using the Python API in ABAQUS CAE
CN109165438A (zh) 碳纤维增强复合材料工件的疲劳寿命预测方法
CN109357957B (zh) 一种基于极值窗口的疲劳监测计数方法

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231220

Address after: No. 8899, Dongfeng Street, Changchun automobile economic and Technological Development Zone, Jilin Province, 130011

Patentee after: China Faw Co.,Ltd.

Patentee after: CHANGCHUN AUTOMOTIVE TEST CENTER Co.,Ltd.

Address before: No. 8899, Dongfeng Street, Changchun automobile economic and Technological Development Zone, Changchun City, Jilin Province, 130011

Patentee before: China Faw Co.,Ltd.