CN114528739A - Simulation method for automobile hub fracture failure - Google Patents

Simulation method for automobile hub fracture failure Download PDF

Info

Publication number
CN114528739A
CN114528739A CN202210200465.6A CN202210200465A CN114528739A CN 114528739 A CN114528739 A CN 114528739A CN 202210200465 A CN202210200465 A CN 202210200465A CN 114528739 A CN114528739 A CN 114528739A
Authority
CN
China
Prior art keywords
test
mat
automobile hub
simulation
material card
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.)
Pending
Application number
CN202210200465.6A
Other languages
Chinese (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.)
Chongqing Changan Automobile Co Ltd
Original Assignee
Chongqing Changan Automobile Co Ltd
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 Chongqing Changan Automobile Co Ltd filed Critical Chongqing Changan Automobile Co Ltd
Priority to CN202210200465.6A priority Critical patent/CN114528739A/en
Publication of CN114528739A publication Critical patent/CN114528739A/en
Pending legal-status Critical Current

Links

Images

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
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a true simulation method for automobile hub fracture failure simulation, which comprises the following steps: s1, sampling the rim and the spoke of the automobile hub real object and carrying out a material test to obtain a force-displacement curve under different test working conditions; s2, combining the obtained force-displacement curves under different test working conditions, applying LS-DYNA finite element analysis software, and respectively establishing an MAT _24 material card and an MAT _ ADD _ EROSION material card of a rim and a spoke, wherein the MAT _24 material card of the wheel center is obtained by zooming the MAT _24 material card of the rim or the spoke, and the MAT _ ADD _ EROSION material card of the wheel center is the same as the MAT _ ADD _ EROSION material card of the spoke; s3, establishing a finite element model of the automobile hub; and S4, performing a static crushing test and a dynamic drop test of the hub, endowing the material card established in S2 to the automobile hub finite element model established in S3, and determining an automobile hub simulation model through test calibration. The accurate simulation of deformation and failure of the automobile hub under the collision working condition can be realized.

Description

一种汽车轮毂断裂失效模仿真模拟方法A simulation method of automobile wheel hub fracture failure model

技术领域technical field

本发明涉及有限元仿真分析,具体涉及汽车轮毂断裂失效模仿真模拟方法。The invention relates to finite element simulation analysis, in particular to a simulation method for a fracture failure model of an automobile wheel hub.

背景技术Background technique

中国保险汽车安全指数即C-IASI中包含的正面25%偏置碰撞是车辆以64KM/h的时速、25%的重叠率(驾驶员侧)正面撞击固定刚性壁障,由于碰撞面积小,车身纵梁及防撞横梁等车身关重件并未完全参与碰撞过程,因此车辆在发生碰撞之后,汽车轮胎总成作为重要的传力路径,将碰撞力传递到驾驶舱中,最直接的影响就是造成驾驶舱变形,挤压驾驶员生存空间,最终会造成驾乘人员的伤害。而在碰撞过程中汽车轮毂都会不同程度的发生断裂失效行为。如何精确模拟轮毂在整车碰撞过程中的断裂失效行为,对整车结构设计及碰撞仿真分析显得尤为重要。The frontal 25% offset collision included in the China Insurance Automobile Safety Index (C-IASI) is the frontal collision of the vehicle with a fixed rigid barrier at a speed of 64KM/h and a 25% overlap rate (driver side). Due to the small collision area, the body The important parts of the body such as longitudinal beams and anti-collision beams are not fully involved in the collision process. Therefore, after the vehicle collides, the tire assembly acts as an important force transmission path to transmit the collision force to the cockpit. The most direct impact is It will cause deformation of the cockpit, squeeze the driver's living space, and eventually cause injury to the driver and passenger. In the process of collision, the wheel hub of the car will break and fail to varying degrees. How to accurately simulate the fracture failure behavior of the wheel hub during the collision of the vehicle is particularly important for the structural design of the vehicle and the analysis of the collision simulation.

目前整车碰撞仿真分析中汽车轮毂未设计失效或者采用简单的力失效或时间失效,与实际碰撞过程中轮毂断裂失效存在较大差异,影响整车碰撞仿真分析精度,无法对整车结构及驾乘人员伤害进行评估和预测。At present, the vehicle wheel hub is not designed to fail in the vehicle crash simulation analysis, or simple force failure or time failure is used, which is quite different from the wheel hub fracture failure during the actual collision process, which affects the accuracy of the vehicle crash simulation analysis, and cannot be used for vehicle structure and driving. Occupant injury assessment and prediction.

发明内容SUMMARY OF THE INVENTION

本发明的目的是提供一种汽车轮毂断裂失效模仿真模拟方法,其能够实现对碰撞工况的汽车轮毂变形及失效的准确仿真模拟。The purpose of the present invention is to provide a simulation method for a fracture failure model of an automobile wheel hub, which can realize accurate simulation of the deformation and failure of the automobile wheel hub in a collision condition.

本发明所述的汽车轮毂断裂失效模拟仿真模拟方法,其包括如下步骤:The automobile wheel hub fracture failure simulation simulation method of the present invention comprises the following steps:

S1,对汽车轮毂实物的轮辋和轮辐取样并进行材料测试试验,获取不同试验工况下的力-位移曲线;S1, sample the rims and spokes of the actual automobile wheel hubs and conduct material tests to obtain the force-displacement curves under different test conditions;

S2,结合步骤S1试验获得的不同试验工况下的力-位移曲线,应用LS-DYNA有限元分析软件,分别建立轮辋和轮辐的MAT_24号材料卡和MAT_ADD_EROSION材料卡,轮心的MAT_24号材料卡通过轮辋或轮辐的MAT_24号材料卡缩放得到,轮心的MAT_ADD_EROSION材料卡与轮辐的MAT_ADD_EROSION材料卡相同;S2, combined with the force-displacement curves under different test conditions obtained from the test in step S1, the LS-DYNA finite element analysis software is used to establish the MAT_24 material card and MAT_ADD_EROSION material card for the rim and the spoke respectively, and the MAT_24 material card for the wheel center. Obtained by scaling the MAT_24 material card of the rim or spoke, the MAT_ADD_EROSION material card of the wheel center is the same as the MAT_ADD_EROSION material card of the spoke;

S3,建立汽车轮毂有限元模型,所述汽车轮毂有限元模型包括轮辋、轮辐和轮心三个部分,三个部分的连接处共用节点;S3, establishing a finite element model of an automobile hub, the finite element model of the automobile hub includes three parts, a rim, a spoke and a wheel center, and the joints of the three parts share a node;

S4,进行轮毂静态压溃试验和动态跌落试验,将S2建立的材料卡赋予S3建立的汽车轮毂有限元模型中,通过试验标定,确定汽车轮毂仿真模型。S4, carry out the static crush test and dynamic drop test of the wheel hub, assign the material card established by S2 to the finite element model of the automobile wheel hub established by S3, and determine the simulation model of the automobile wheel hub through test calibration.

进一步,所述材料测试试验包括单向拉伸试验、剪切试验、拉剪试验、R5缺口拉伸试验、R10缺口拉伸试验、中心孔拉伸试验、压缩试验、穿孔试验。Further, the material testing test includes uniaxial tensile test, shear test, tensile shear test, R5 notched tensile test, R10 notched tensile test, central hole tensile test, compression test, and perforation test.

进一步,所述轮辋和轮辐的MAT_24号材料卡的建立具体为:基于获取的单向拉伸工况下的力-位移曲线得到轮辋和轮辐的真实应力应变曲线,将真实应力应变曲线通过硬化分析模型进行外推得到MAT_24号材料卡的硬化曲线,MAT_24号材料卡的其他信息根据汽车轮毂材料进行确定;Further, the establishment of the MAT_24 material card of the rim and the spoke is specifically: based on the obtained force-displacement curve under the uniaxial tensile condition, the true stress-strain curve of the rim and the spoke is obtained, and the true stress-strain curve is analyzed by hardening. The model is extrapolated to obtain the hardening curve of the MAT_24 material card, and other information of the MAT_24 material card is determined according to the material of the automobile wheel hub;

所述轮辋和轮辐的MAT_ADD_EROSION材料卡的建立具体为:采用GISSMO失效模型表征轮辋和轮辐材料断裂失效行为,基于材料测试试验获得的不同试验工况下的力-位移曲线,建立试验测试样品的有限元对标模型,并使得有限元对标模型仿真与试验力-位移曲线吻合度在90%及以上;然后通过有限元对标模型仿真结果得到不同试验工况下的实际应力三轴度及等效塑性失效应变,将不同应力状态下得到的仿真数据进行数据拟合,得到MAT_ADD_EROSION材料卡的断裂失效曲线,MAT_ADD_EROSION材料卡的其他信息根据汽车轮毂材料进行确定。The establishment of the MAT_ADD_EROSION material card of the rim and the spoke is as follows: the GISSMO failure model is used to characterize the fracture failure behavior of the rim and the spoke material, and based on the force-displacement curves under different test conditions obtained from the material test test, the finite limit of the test sample is established. Then, through the simulation results of the finite element benchmarking model, the actual stress triaxiality and so on under different test conditions can be obtained. The plastic failure strain is calculated, and the simulation data obtained under different stress states are fitted to the data to obtain the fracture failure curve of the MAT_ADD_EROSION material card. Other information of the MAT_ADD_EROSION material card is determined according to the material of the automobile wheel hub.

进一步,所述MAT_24号材料卡其他信息包括:材料密度、弹性模量E、泊松比参数,通过查表获得;Further, other information of the MAT_24 material card includes: material density, elastic modulus E, and Poisson's ratio parameters, which are obtained by looking up a table;

所述MAT_ADD_EROSION材料卡的其他信息包括材料失稳曲线、尺寸效应曲线、应力衰减指数;Other information of the MAT_ADD_EROSION material card includes material instability curve, size effect curve, stress decay index;

所述材料失稳曲线通过单向拉伸试验中材料发生颈缩时对应的等效塑性应变,在-2/3压缩应力态到2/3双向拉伸应力态区间下为一定值;The material instability curve is a certain value in the range of -2/3 compression stress state to 2/3 biaxial tensile stress state through the equivalent plastic strain corresponding to the necking of the material in the uniaxial tensile test;

所述尺寸效应曲线通过单向拉伸有限元对标模型,分别建立0.5mm、1mm、2mm、4mm和8mm有限元模型进行仿真与试验的对标,最终分别获得不同尺寸下能够与单向拉伸试验工况下的力-位移曲线吻合度在90%以上,以0.5mm下的等效塑性失效应变为准进行归一化处理,可以得到不同尺寸下对断裂失效曲线的缩放系数,即为尺寸效应曲线;The size effect curve is compared with the unidirectional tensile finite element model, and the 0.5mm, 1mm, 2mm, 4mm and 8mm finite element models are established for simulation and test benchmarking, and finally obtained under different sizes. The fit degree of the force-displacement curve under the tensile test condition is more than 90%, and the equivalent plastic failure strain under 0.5mm is used for normalization, and the scaling factor of the fracture failure curve under different sizes can be obtained, which is size effect curve;

所述应力衰减指数通过手动输入不同值来对标单向拉伸试验工况下的力-位移曲线中发生颈缩之后曲线下掉的幅度,使仿真与试验力-位移曲线吻合度在90%及以上作为应力衰减指数值。The stress decay exponent is used to manually input different values to benchmark the force-displacement curve under the uniaxial tensile test condition with the magnitude of the curve falling off after necking occurs, so that the simulation and test force-displacement curve are in 90% agreement. and above as the stress decay index value.

进一步,所述S3中采用四面体单元建立汽车轮毂的有限元模型。Further, in S3, a tetrahedral element is used to establish a finite element model of the automobile wheel hub.

进一步,所述S4中的静态压溃试验具体为:将汽车轮毂固定于底座上,冲头以准静态的速度进行静压轮毂,同时约束冲头除垂向外的其他方向自由度,使得汽车轮毂在静压的作用下发生断裂失效,分别得到静态压溃时基于汽车轮毂有限元模型的仿真最大力值F1和基于汽车轮毂实物的实际试验最大力值F2,若F1和F2的比值为85~115%,且汽车轮毂仿真断裂失效形式与实际试验断裂失效形式相似,则判定汽车轮毂有限元模型对标精度满足要求,反之则判定对标精度不满足要求,返回S2优化MAT_24号材料卡和MAT_ADD_EROSION材料卡后重新进行静态压溃试验。Further, the static crushing test in the S4 is specifically as follows: fixing the automobile wheel hub on the base, the punch presses the wheel hub statically at a quasi-static speed, and constrains the punch in other directions except the vertical and outward degrees of freedom, so that the automobile The wheel hub fractures and fails under the action of static pressure, and the simulated maximum force value F 1 based on the finite element model of the automobile wheel hub and the actual test maximum force value F 2 based on the actual automobile wheel hub during static crushing are obtained respectively. If F 1 and F 2 The ratio is 85 to 115%, and the simulation fracture failure mode of the automobile hub is similar to the actual test fracture failure mode, then it is determined that the benchmarking accuracy of the finite element model of the automobile hub meets the requirements, otherwise it is determined that the benchmarking accuracy does not meet the requirements, and return to S2 to optimize MAT_24 No. material card and MAT_ADD_EROSION material card and re-run static crush test.

进一步,所述S4中的动态跌落试验具体为:将汽车轮毂固定于底座上,约束冲头除垂向外的其他方向自由度,冲头冲击汽车轮毂的速度通过调整冲头质量来控制,使得汽车轮毂在冲头冲击力的作用下发生断裂失效,分别得到动态跌落冲击时基于汽车轮毂有限元模型的仿真最大力值F3和基于汽车轮毂实物的实际试验最大力值F4,若F3和F4的比值为85~115%,且汽车轮毂仿真断裂失效形式与实际试验断裂失效形式相似,则判定汽车轮毂有限元模型对标精度满足要求,反之则判定对标精度不满足要求,返回步骤二优化MAT_24号材料卡和MAT_ADD_EROSION材料卡后重新进行静态压溃试验。Further, the dynamic drop test in the S4 is specifically as follows: fixing the car wheel hub on the base, constraining the punch in other directions of freedom except vertical and outward, and the speed of the punch impacting the car wheel hub is controlled by adjusting the punch quality, so that The car wheel breaks and fails under the impact force of the punch, and the simulated maximum force value F 3 based on the finite element model of the car wheel hub and the actual test maximum force value F 4 based on the real car wheel hub during the dynamic drop impact are obtained respectively. If F 3 The ratio of F 4 to 85% to 115%, and the simulation fracture failure mode of the automobile wheel hub is similar to the actual test fracture failure mode, then it is judged that the benchmarking accuracy of the finite element model of the automobile wheel hub meets the requirements, otherwise, it is determined that the benchmarking accuracy does not meet the requirements, return Step 2: After optimizing the MAT_24 material card and the MAT_ADD_EROSION material card, perform the static crush test again.

本发明与现有技术相比具有如下有益效果。Compared with the prior art, the present invention has the following beneficial effects.

1、本发明所述汽车轮毂断裂失效模拟仿真模拟方法考虑了汽车轮毂不同区域力学性能上的差异,利用前处理软件进行精细化建模。同时利用LS-DYNA分析软件建立MAT_24号材料卡和MAT_ADD_EROSION材料卡,结合设计验证试验对标优化,获取含断裂失效行为的汽车轮毂有限元仿真分析模型,真实的反映了汽车轮毂的力学特性即断裂失效行为,以便于在设计阶段对整车结构进行设计变更。1. The automobile wheel hub fracture failure simulation simulation method of the present invention takes into account the differences in mechanical properties of different regions of the automobile wheel hub, and uses preprocessing software to carry out refined modeling. At the same time, the LS-DYNA analysis software was used to establish the MAT_24 material card and the MAT_ADD_EROSION material card. Combined with the design verification test benchmarking optimization, the finite element simulation analysis model of the automobile wheel hub with fracture failure behavior was obtained, which truly reflected the mechanical characteristics of the automobile wheel hub, namely fracture. Failure behavior to facilitate design changes to the vehicle structure during the design phase.

2、本发明采用了GISSMO模型模拟材料损伤失效,考虑了不同应力状态下的失效行为,真实的反映出汽车轮毂在碰撞工况下断裂失效行为,提升了整车碰撞仿真分析的精度,并且能够有效指导整车结构设计及车体开发。2. The present invention uses the GISSMO model to simulate material damage and failure, considers the failure behavior under different stress states, truly reflects the fracture failure behavior of the automobile wheel hub under the collision condition, improves the accuracy of the vehicle collision simulation analysis, and can Effectively guide the vehicle structure design and body development.

3、本发明所述所述汽车轮毂断裂失效模拟仿真模拟方法流程简单,实施方便,为其他零部件的断裂失效行为的模拟研究提供了思路。3. The method for simulating the fracture failure of the automobile wheel hub according to the present invention has a simple process and is convenient to implement, and provides an idea for the simulation study of the fracture failure behavior of other parts.

附图说明Description of drawings

图1是本发明所述汽车轮毂断裂失效模拟仿真模拟方法的流程示意图;Fig. 1 is the schematic flow chart of the simulation simulation method of automobile wheel hub fracture failure according to the present invention;

图2是本发明材料测试样本的取样示意图;Fig. 2 is the sampling schematic diagram of the material test sample of the present invention;

图3是本发明所述汽车轮毂的结构示意图;Fig. 3 is the structural representation of the automobile wheel hub of the present invention;

图4是本发明所述断裂失效曲线示意图;Fig. 4 is the schematic diagram of fracture failure curve of the present invention;

图5是本发明所述汽车轮毂的有限元模型示意图;Fig. 5 is the finite element model schematic diagram of the automobile wheel hub of the present invention;

图6是本发明所述汽车轮毂静态压溃验证试验示意图之一;6 is one of the schematic diagrams of the static crush verification test of the automobile wheel hub according to the present invention;

图7是本发明所述汽车轮毂静态压溃验证试验示意图之二;7 is the second schematic diagram of the static crush verification test of the automobile wheel hub according to the present invention;

图8是本发明所述汽车轮毂动态跌落验证试验示意图。FIG. 8 is a schematic diagram of the dynamic drop verification test of the automobile wheel hub according to the present invention.

具体实施方式Detailed ways

下面结合附图对本发明作详细说明。The present invention will be described in detail below with reference to the accompanying drawings.

参见图1,所示的汽车轮毂断裂失效模拟仿真模拟方法,其包括如下步骤:Referring to Fig. 1, the shown simulation method for simulating failure of automobile wheel hub includes the following steps:

S1,对汽车轮毂实物的轮辋和轮辐取样并进行材料测试试验,具体为:采用CAD软件读取汽车轮毂三维模型数据,根据材料测试试验的测试样本大小,在汽车轮毂三维模型中选取能够作为材料测试样本的区域,参见图2,一般选取汽车轮毂1上较为平整的区域,即以轮辋11和轮辐12区域的部位位置作为取样区域。S1, sample the rims and spokes of the actual automobile wheel hub and conduct material testing tests, specifically: using CAD software to read the three-dimensional model data of the automobile wheel hub, and according to the test sample size of the material test test, select the materials that can be used as the material in the three-dimensional model of the automobile wheel hub. For the area of the test sample, see FIG. 2 , generally a relatively flat area on the automobile wheel hub 1 is selected, that is, the position of the rim 11 and the spoke 12 area is taken as the sampling area.

依据在CAD中确定的取样区域,在汽车轮毂实物中取出包括单向拉伸试验、剪切试验、拉剪试验、R5缺口拉伸试验、R10缺口拉伸试验、中心孔拉伸试验、压缩试验、穿孔试验用的测试样本,然后将取出的测试样本开展材料测试试验,即开展材料本构模型测试试验和材料断裂失效测试试验,获得不同试验工况下的力-位移曲线。According to the sampling area determined in CAD, take out from the real vehicle hub including uniaxial tensile test, shear test, tensile shear test, R5 notched tensile test, R10 notched tensile test, central hole tensile test, compression test , test samples for perforation test, and then carry out the material test test on the test sample taken out, that is, carry out the material constitutive model test test and the material fracture failure test test, and obtain the force-displacement curve under different test conditions.

S2,结合步骤S1试验获得的不同试验工况下的力-位移曲线,采用有限元仿真手段建立测试样本等效的有限元对标模型,进行轮毂材料断裂失效卡片的开发,所述轮毂断裂失效卡片是通过LS-DYNA有限元分析软件中定义关键字MAT_24号材料卡和MAT_ADD_EROSION材料卡弹塑性模型耦合使用,所述分别MAT_24号材料卡和MAT_ADD_EROSION材料卡描述材料的弹塑性行为以及断裂行为。S2, combined with the force-displacement curves under different test conditions obtained in the test in step S1, the finite element simulation method is used to establish an equivalent finite element benchmark model of the test sample, and the development of the hub material fracture failure card is carried out, and the hub fracture failure The card is used by coupling the elastoplastic model with the keywords MAT_24 material card and MAT_ADD_EROSION material card defined in the LS-DYNA finite element analysis software. The MAT_24 material card and MAT_ADD_EROSION material card describe the elastic-plastic behavior and fracture behavior of the material respectively.

受结构特征影响,参见图3,汽车轮毂1中轮辋11、轮辐12以及轮心13的屈服强度存在差异,因此应用LS-DYNA有限元分析软件,分别建立轮辋11和轮辐12的MAT_24号材料卡和MAT_ADD_EROSION材料卡。由于难以在轮心13位置取样,因此轮心13的MAT_24号材料卡通过轮辋11或轮辐12的MAT_24号材料卡缩放得到,缩放系数根据经验幅值,轮心13的断裂失效形式与轮辐12的断裂失效形式相近,因此轮心13的MAT_ADD_EROSION材料卡设置为与轮辐的MAT_ADD_EROSION材料卡相同。Affected by structural characteristics, see Figure 3, there are differences in the yield strength of rim 11, spoke 12 and wheel center 13 in automobile hub 1. Therefore, the LS-DYNA finite element analysis software is used to establish MAT_24 material cards for rim 11 and spoke 12 respectively. and MAT_ADD_EROSION material card. Since it is difficult to take samples at the wheel center 13, the MAT_24 material card of the wheel center 13 is obtained by scaling the MAT_24 material card of the rim 11 or the spoke 12. The scaling factor is based on the empirical amplitude. The fracture failure modes are similar, so the MAT_ADD_EROSION material card of the wheel center 13 is set to be the same as the MAT_ADD_EROSION material card of the spokes.

所述轮辋11和轮辐12的MAT_24号材料卡的建立具体为:基于获取的单向拉伸工况下的力-位移曲线得到轮辋和轮辐的真实应力应变曲线,将真实应力应变曲线通过硬化分析模型进行外推得到MAT_24号材料卡的硬化曲线,所述硬化分析模型包括Swift、HS和VOCE中的至少一种。建立测试样本的有限元对比模型,进行仿真与试验的对标,使得仿真与试验的力-位移曲线吻合度在90%及以上,若吻合度低于90%,则重新通过硬化分析模型将真实应力应变曲线进行外推得到MAT_24号材料卡的硬化曲线。MAT_24号材料卡的其他信息根据汽车轮毂材料进行确定,具体地所述MAT_24号材料卡其他信息包括:材料密度、弹性模量E、泊松比参数,通过查表获得;The establishment of the MAT_24 material card of the rim 11 and the spoke 12 is as follows: based on the obtained force-displacement curve under the uniaxial tensile condition, the true stress-strain curve of the rim and the spoke is obtained, and the true stress-strain curve is analyzed by hardening. The model is extrapolated to obtain the hardening curve of the material card No. MAT_24, and the hardening analysis model includes at least one of Swift, HS and VOCE. The finite element comparison model of the test sample is established, and the simulation and test are compared, so that the force-displacement curve of the simulation and the test agrees with 90% or more. The stress-strain curve was extrapolated to obtain the hardening curve of the MAT_24 material card. The other information of the MAT_24 material card is determined according to the material of the automobile wheel hub. Specifically, the other information of the MAT_24 material card includes: material density, elastic modulus E, and Poisson's ratio parameters, which are obtained by looking up the table;

所述轮辋11和轮辐12的MAT_ADD_EROSION材料卡的建立具体为:采用GISSMO失效模型表征轮辋和轮辐材料断裂失效行为,基于材料测试试验获得的不同试验工况下的力-位移曲线,建立试验测试样品的有限元对标模型,并使得有限元对标模型仿真与试验力-位移曲线吻合度在90%及以上,若吻合度低于90%则重新建立有限元对标模型,直至吻合度满足要求。然后通过有限元对标模型仿真结果得到不同试验工况下的实际应力三轴度及等效塑性失效应变,将不同应力状态下得到的仿真数据进行数据拟合,参见图4,得到MAT_ADD_EROSION材料卡的断裂失效曲线。所述MAT_ADD_EROSION材料卡的其他信息包括材料失稳曲线、尺寸效应曲线、应力衰减指数;所述材料失稳曲线通过单向拉伸试验中材料发生颈缩时对应的等效塑性应变,在-2/3压缩应力态到2/3双向拉伸应力态区间下为一定值;所述尺寸效应曲线通过单向拉伸有限元对标模型,分别建立0.5mm、1mm、2mm、4mm和8mm有限元模型进行仿真与试验的对标,最终分别获得不同尺寸下能够与单向拉伸试验工况下的力-位移曲线吻合度在90%以上,以0.5mm下的等效塑性失效应变为准进行归一化处理,可以得到不同尺寸下对断裂失效曲线的缩放系数,即为尺寸效应曲线;所述应力衰减指数通过手动输入不同值来对标单向拉伸试验工况下的力-位移曲线中发生颈缩之后曲线下掉的幅度,使仿真与试验力-位移曲线吻合度在90%及以上作为应力衰减指数值。The establishment of the MAT_ADD_EROSION material card of the rim 11 and the spoke 12 is as follows: the GISSMO failure model is used to characterize the fracture failure behavior of the rim and the spoke material, and the test sample is established based on the force-displacement curve under different test conditions obtained from the material test test. The finite element benchmarking model is established, and the fit of the finite element benchmarking model and the test force-displacement curve is 90% or above. If the fit is less than 90%, the finite element benchmarking model is re-established until the fit meets the requirements. . Then, the actual stress triaxiality and equivalent plastic failure strain under different test conditions are obtained through the simulation results of the finite element benchmarking model, and the simulation data obtained under different stress states are fitted to the data, as shown in Figure 4, to obtain the MAT_ADD_EROSION material card The fracture failure curve. Other information of the MAT_ADD_EROSION material card includes material buckling curve, size effect curve, and stress decay index; the material buckling curve is the equivalent plastic strain corresponding to the necking of the material in the uniaxial tensile test, at -2 /3 compressive stress state to 2/3 biaxial tensile stress state is a certain value; the size effect curve is based on the uniaxial tensile finite element model, and the 0.5mm, 1mm, 2mm, 4mm and 8mm finite element models are established respectively. The model is compared with the simulation and test, and finally the force-displacement curves of different sizes and the uniaxial tensile test conditions are obtained to be more than 90% consistent, and the equivalent plastic failure strain under 0.5mm shall prevail. By normalizing, the scaling factor of the fracture failure curve under different sizes can be obtained, that is, the size effect curve; the stress decay index can be compared with the force-displacement curve under the uniaxial tensile test condition by manually entering different values. The magnitude of the curve drop after necking occurs in the middle, so that the fit between the simulation and the test force-displacement curve is 90% and above as the stress decay index value.

S3,参见图5,采用前处理软件Hypermesh建立汽车轮毂有限元模型,所述汽车轮毂有限元模型包括轮辋、轮辐和轮心三个部分,三个部分的连接处共用节点。所述汽车轮毂有限元模型所述模型采用四面体单元建模,由于汽车轮毂造型及结构的复杂,单元网格大小控制在2-4mm。由于轮毂结构复杂并不规整,六面体体单元的网格划分难度较大,不利于工程化应用,采用四面体单元进行划分,通过采用不同单元积分类型对比六面体单元计算结果,选择仿真计算结果对比差异最小的四面体单元积分类型。同时为确保仿真分析结果的准确性,厚度方向上至少分布两层体网格单元;S3, referring to FIG. 5, the preprocessing software Hypermesh is used to establish a finite element model of an automobile wheel hub. The finite element model of the automobile wheel hub includes three parts, a rim, a spoke and a wheel center, and the joints of the three parts share a node. The finite element model of the automobile wheel hub The model is modeled with tetrahedral elements. Due to the complexity of the shape and structure of the automobile wheel hub, the mesh size of the element is controlled at 2-4mm. Due to the complex and irregular structure of the hub, the meshing of hexahedral elements is difficult, which is not conducive to engineering applications. Tetrahedral elements are used for division, and the calculation results of hexahedral elements are compared by using different element integration types, and the simulation results are selected to compare the differences. The smallest tetrahedral element integration type. At the same time, in order to ensure the accuracy of the simulation analysis results, at least two layers of volume mesh elements are distributed in the thickness direction;

S4,进行轮毂静态压溃试验和动态跌落试验,将S2建立的材料卡赋予S3建立的汽车轮毂有限元模型中,通过试验标定,确定汽车轮毂仿真模型。S4, carry out the static crush test and dynamic drop test of the wheel hub, assign the material card established by S2 to the finite element model of the automobile wheel hub established by S3, and determine the simulation model of the automobile wheel hub through test calibration.

所述静态压溃试验具体为:建立汽车轮毂冲头静压仿真模型,包括赋有MAT_24号材料卡和MAT_ADD_EROSION材料卡的汽车轮毂有限元模型、底座有限元模型和冲头有限元模型。参见图6和图7,将汽车轮毂1固定于底座2上,冲头3以准静态的速度进行静压汽车轮毂1,同时约束冲头3除垂向外的其他方向自由度,使得汽车轮毂1在静压的作用下发生断裂失效,分别得到静态压溃时基于汽车轮毂有限元模型的仿真最大力值F1=286kN和基于汽车轮毂实物的实际试验最大力值F2=273kN,F1和F2的比值为104.8%,判定汽车轮毂有限元模型对标精度满足要求。The static crush test is specifically: establishing a static pressure simulation model of an automobile wheel hub punch, including an automobile wheel hub finite element model, a base finite element model and a punch finite element model with MAT_24 material cards and MAT_ADD_EROSION material cards. Referring to Figure 6 and Figure 7, the automobile wheel hub 1 is fixed on the base 2, and the punch 3 statically presses the automobile wheel hub 1 at a quasi-static speed, while constraining the punch 3 in other directions except the vertical and outward degrees of freedom, so that the automobile wheel hub 1 Fracture failure occurs under the action of static pressure, and the simulation maximum force value F 1 = 286kN based on the finite element model of the automobile wheel hub and the actual test maximum force value F 2 = 273kN based on the actual automobile wheel hub during static crushing are obtained. The ratio of F 2 and F 2 is 104.8%, and it is judged that the benchmarking accuracy of the finite element model of the automobile wheel hub meets the requirements.

参见图5,冲头3静压位置位于汽车轮毂1靠近轮辐12和轮心13一侧,进行轮辐12和轮心13强度对标及断裂失效模式的对标。参见图6,头3静压位置位于汽车轮毂1远离轮辐12和轮心13一侧,进行轮辋11强度对标及断裂失效模式的对标。Referring to FIG. 5 , the static pressure position of the punch 3 is located on the side of the automobile wheel hub 1 close to the wheel spokes 12 and the wheel center 13 . Referring to FIG. 6 , the static pressure position of the head 3 is located on the side of the automobile wheel hub 1 away from the spoke 12 and the wheel center 13 , and the rim 11 strength benchmarking and fracture failure mode benchmarking are performed.

所述动态跌落试验具体为:参见图8,将汽车轮毂1固定于底座2上,约束冲头3除垂向外的其他方向自由度,冲头3冲击汽车轮毂1的速度通过调整冲头3质量来控制,使得汽车轮毂1在冲头3冲击力的作用下发生断裂失效,分别得到动态跌落冲击时基于汽车轮毂有限元模型的仿真最大力值F3=365kN和基于汽车轮毂实物的实际试验最大力值F4=355kN,若F3和F4的比值为102.8%,判定汽车轮毂有限元模型对标精度满足要求,反之则判定对标精度不满足要求。The dynamic drop test is specifically as follows: referring to FIG. 8 , the automobile wheel hub 1 is fixed on the base 2, the freedom of the punch 3 is restricted in other directions except the vertical direction, and the speed of the punch 3 impacting the automobile wheel hub 1 is adjusted by adjusting the punch 3. To control the quality, the automobile hub 1 breaks and fails under the impact force of the punch 3, and the simulated maximum force value F 3 = 365kN based on the finite element model of the automobile hub during dynamic drop impact and the actual test based on the actual automobile hub are obtained respectively. The maximum force value F 4 =355kN, if the ratio of F 3 and F 4 is 102.8%, it is determined that the benchmarking accuracy of the finite element model of the automobile wheel hub meets the requirements, otherwise it is determined that the benchmarking accuracy does not meet the requirements.

以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention shall be included in the protection of the present invention. within the range.

Claims (7)

1. An automobile hub fracture failure simulation method is characterized by comprising the following steps:
s1, sampling the rim and the spoke of the automobile hub object and performing a material test to obtain force-displacement curves under different test working conditions;
s2, combining the force-displacement curves under different test conditions obtained by the test of the step S1, applying LS-DYNA finite element analysis software, and respectively establishing an MAT _24 material card and an MAT _ ADD _ EROSION material card of a rim and a spoke, wherein the MAT _24 material card of the wheel center is obtained by scaling the MAT _24 material card of the rim or the spoke, and the MAT _ ADD _ EROSION material card of the wheel center is the same as the MAT _ ADD _ EROSION material card of the spoke;
s3, establishing an automobile hub finite element model, wherein the automobile hub finite element model comprises a rim, a spoke and a wheel center, and joints of the three parts share a node;
and S4, performing a static crushing test and a dynamic drop test on the hub, endowing the material card established in S2 to the automobile hub finite element model established in S3, and determining an automobile hub simulation model through test calibration.
2. The automobile hub fracture failure simulation method according to claim 1, characterized in that: the material test comprises a unidirectional tensile test, a shear test, a tensile-shear test, an R5 notch tensile test, an R10 notch tensile test, a center hole tensile test, a compression test and a perforation test.
3. The automobile hub fracture failure simulation method of claim 2, wherein the establishment of the MAT _24 material cards of the rim and the spoke specifically comprises the following steps: obtaining a real stress-strain curve of a rim and a spoke based on the obtained force-displacement curve under the unidirectional tension working condition, extrapolating the real stress-strain curve through a hardening analysis model to obtain a hardening curve of the MAT _24 material card, and determining other information of the MAT _24 material card according to the automobile hub material;
the establishment of the MAT _ ADD _ EROSION material card of the rim and the spoke is specifically as follows: adopting a GISSMO failure model to represent the fracture failure behavior of rim and spoke materials, establishing a finite element benchmarking model of a test sample based on force-displacement curves under different test working conditions obtained by a material test, and enabling the fit degree of the finite element benchmarking model simulation and the test force-displacement curve to be 90% or more; and then obtaining the three-axis degree of actual stress and equivalent plastic failure strain under different test working conditions through a finite element standard model simulation result, performing data fitting on simulation data obtained under different stress conditions to obtain a fracture failure curve of the MAT _ ADD _ EROSION material card, and determining other information of the MAT _ ADD _ EROSION material card according to the test data.
4. The automobile hub fracture failure simulation method according to claim 3, characterized in that: the MAT _24 material card comprises other information: the material density, the elastic modulus E and the Poisson ratio parameters are obtained by looking up a table;
other information of the MAT _ ADD _ EROSION material card comprises a material instability curve, a size effect curve and a stress decay index;
the material instability curve passes through corresponding equivalent plastic strain when the material is necked in a uniaxial tension test, and is a certain value in a range from a-2/3 compressive stress state to a 2/3 biaxial tension stress state;
the size effect curve is subjected to simulation and test benchmarking by respectively establishing 0.5mm, 1mm, 2mm, 4mm and 8mm finite element benchmarking models through a unidirectional stretching finite element benchmarking model, finally, the matching degree of force-displacement curves which can be matched with the unidirectional stretching test under the working condition under different sizes is respectively obtained to be more than 90%, normalization processing is carried out by taking equivalent plastic failure strain under 0.5mm as a reference, and the scaling coefficient of fracture failure curves under different sizes can be obtained, namely the size effect curve;
the stress attenuation index is characterized in that different values are manually input to compare the falling amplitude of a curve after necking in a force-displacement curve under the working condition of a standard unidirectional tensile test, so that the coincidence degree of simulation and test force-displacement curves is 90% or more and is used as the stress attenuation index value.
5. The automobile hub fracture failure simulation method according to claim 1 or 2, characterized in that: and in the step S3, a finite element model of the automobile hub is established by adopting tetrahedral units.
6. The automobile hub fracture failure simulation method according to claim 1 or 2, wherein the static crushing test in the S4 is specifically as follows: the automobile hub is fixed on the base, the punch is used for carrying out static pressure on the hub at a quasi-static speed, and meanwhile, the freedom degrees of the punch in other directions except the vertical direction are restrained, so that the automobile hub is broken and fails under the action of static pressure, and the simulation maximum force value F based on the automobile hub finite element model is obtained when static crushing is respectively obtained1And practical test maximum force value F based on automobile hub real object2If F is1And F2The ratio of the automobile hub simulation fracture failure mode to the standard precision is 85-115%, if the automobile hub simulation fracture failure mode is similar to the actual test fracture failure mode, the standard precision is judged to meet the requirement, otherwise, the standard precision is judged not to meet the requirement, and the static crushing test is carried out again after the MAT _24 material card and the MAT _ ADD _ EROSION material card are optimized in the second step.
7. The automobile hub fracture failure simulation method according to claim 1 or 2, wherein the dynamic drop test in the S4 specifically comprises: on being fixed in the base with automobile wheel hub, restraint drift other direction degrees of freedom except vertical, the speed that the drift impacted automobile wheel hub is controlled through adjustment drift quality for automobile wheel hub is dashing towardsThe head is broken and failed under the action of impact force, and the simulation maximum force value F based on the finite element model of the automobile hub during dynamic falling impact is respectively obtained3And practical test maximum force value F based on automobile hub real object4If F is3And F4The ratio of the automobile hub simulation fracture failure mode to the actual test fracture failure mode is 85-115%, if the automobile hub simulation fracture failure mode is similar to the actual test fracture failure mode, the standard alignment precision of the automobile hub finite element model is judged to meet the requirement, otherwise, the standard alignment precision is judged not to meet the requirement, and the automobile hub simulation fracture failure mode returns to S2 to optimize the MAT _24 material card and the MAT _ ADD _ EROSION material card and then the static crushing test is carried out again.
CN202210200465.6A 2022-02-28 2022-02-28 Simulation method for automobile hub fracture failure Pending CN114528739A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210200465.6A CN114528739A (en) 2022-02-28 2022-02-28 Simulation method for automobile hub fracture failure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210200465.6A CN114528739A (en) 2022-02-28 2022-02-28 Simulation method for automobile hub fracture failure

Publications (1)

Publication Number Publication Date
CN114528739A true CN114528739A (en) 2022-05-24

Family

ID=81626034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210200465.6A Pending CN114528739A (en) 2022-02-28 2022-02-28 Simulation method for automobile hub fracture failure

Country Status (1)

Country Link
CN (1) CN114528739A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116227087A (en) * 2023-03-29 2023-06-06 南京理工大学 A Prediction Method of Tensile Properties of Welded Joints of Thin-walled Parts Based on Finite Element Method
CN116384258A (en) * 2023-05-29 2023-07-04 浙江大学 A wheel-tire integrated wheel impact dynamics simulation method
CN117313238A (en) * 2023-09-22 2023-12-29 中国汽车工程研究院股份有限公司 Multi-physical quantity calibration method based on trolley subsystem

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105701301A (en) * 2016-01-14 2016-06-22 北京汽车股份有限公司 Vehicle forward intersection detection method and device
CN111125960A (en) * 2019-12-26 2020-05-08 中国汽车工程研究院股份有限公司 A Parameter Optimization Method for GISSMO Material Failure Model
CN111523183A (en) * 2020-04-20 2020-08-11 湖南大学 A Simulation Modeling Method for Mechanical Properties and Fracture Failure of Welded Joints
CN111680358A (en) * 2020-05-09 2020-09-18 中国第一汽车股份有限公司 Collision simulation method for automobile aluminum alloy section parts
CN111753373A (en) * 2020-06-24 2020-10-09 中国第一汽车股份有限公司 Collision simulation method of automobile foamed plastic material

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105701301A (en) * 2016-01-14 2016-06-22 北京汽车股份有限公司 Vehicle forward intersection detection method and device
CN111125960A (en) * 2019-12-26 2020-05-08 中国汽车工程研究院股份有限公司 A Parameter Optimization Method for GISSMO Material Failure Model
CN111523183A (en) * 2020-04-20 2020-08-11 湖南大学 A Simulation Modeling Method for Mechanical Properties and Fracture Failure of Welded Joints
CN111680358A (en) * 2020-05-09 2020-09-18 中国第一汽车股份有限公司 Collision simulation method for automobile aluminum alloy section parts
CN111753373A (en) * 2020-06-24 2020-10-09 中国第一汽车股份有限公司 Collision simulation method of automobile foamed plastic material

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
吴昊 等: "C-NCAP工况下铝合金动态力学与失效特性仿真分析", 《新技术新工艺》, no. 9, 25 September 2020 (2020-09-25), pages 53 - 57 *
吴昊;孟宪明;田杰斌;方锐;黄亚峰;张赛;郑崇嵩;: "C-NCAP工况下铝合金动态力学与失效特性仿真分析", 新技术新工艺, no. 09, 25 September 2020 (2020-09-25) *
曾龙: "高强钢点焊接头失效行为及仿真模型研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》, no. 7, 15 July 2015 (2015-07-15), pages 022 - 417 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116227087A (en) * 2023-03-29 2023-06-06 南京理工大学 A Prediction Method of Tensile Properties of Welded Joints of Thin-walled Parts Based on Finite Element Method
CN116384258A (en) * 2023-05-29 2023-07-04 浙江大学 A wheel-tire integrated wheel impact dynamics simulation method
CN116384258B (en) * 2023-05-29 2023-09-29 浙江大学 A hub-tire integrated wheel impact dynamics simulation method
CN117313238A (en) * 2023-09-22 2023-12-29 中国汽车工程研究院股份有限公司 Multi-physical quantity calibration method based on trolley subsystem

Similar Documents

Publication Publication Date Title
CN114528739A (en) Simulation method for automobile hub fracture failure
CN111753373B (en) Collision simulation method of automobile foamed plastic material
CN111060396B (en) Material mechanical property calibration method based on Ls-Dyna No. 187 material card
EP2191928B1 (en) Method of, computer readable medium and system for determination of spot weld failure using finite element analysis with stress and location sensitivity scale factors
US7818147B2 (en) Analysis data judging apparatus, simulation system and simulation program
US9043187B2 (en) Method of material modeling for crash test dummy finite element models
CN110377966A (en) Truck door crushing strength simulation method
CN101241521A (en) Modeling Method of Comprehensive Performance Index of Car Body Based on Support Vector Machine
CN113033040B (en) Accurate modeling method for vehicle flexible connection
CN113420371A (en) Method for improving simulation analysis precision of misuse working condition of McPherson front suspension subframe
Bojanowski Verification, validation and optimization of finite element model of bus structure for rollover test
Shcherba et al. Developing of phenomenological damage model for automotive low-carbon structural steel for using in validation of Euroncap frontal impact
Zhu et al. Application of AM60B magnesium alloy material model to structural component crush analysis
Balla et al. Incommensurate vibro-acoustic performance due to in-process blank holder force variation during deep drawing process
Sato impact
Zhang et al. Study on the failure behaviour of Advanced High Strength Steel considering stamping history
Gao et al. Impact simulation of Hydro-formed front end vehicle structure
CN114912317B (en) A high-speed collision material card for automobile steel and its application method
Böttcher et al. Consideration of manufacturing effects to improve crash simulation accuracy
CN114218704A (en) Failure analysis and judgment method for high-strength steel for automobile
Samuel Numerical validation of carbon fibre composite crush absorbers via finite element analysis for dynamic impact energy absorption
Mo et al. A study of multibody modeling and calibration method for FMVSS 208 rollover test
Untaroiu et al. Development and validation of a headform impactor finite element model with application to vehicle hood design for pedestrian protection
Bistedt et al. Explicit dynamic analysis of Drive Over Curb load case
Avalle et al. Stochastic crash analysis of vehicle models for sensitivity analysis and optimization

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