CN112966407B - Simulation method for rubber riveting composite connection collision working condition - Google Patents

Simulation method for rubber riveting composite connection collision working condition Download PDF

Info

Publication number
CN112966407B
CN112966407B CN202110120476.9A CN202110120476A CN112966407B CN 112966407 B CN112966407 B CN 112966407B CN 202110120476 A CN202110120476 A CN 202110120476A CN 112966407 B CN112966407 B CN 112966407B
Authority
CN
China
Prior art keywords
riveting
viscose
glue
mechanical property
rivet
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.)
Active
Application number
CN202110120476.9A
Other languages
Chinese (zh)
Other versions
CN112966407A (en
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.)
China Automotive Technology and Research Center Co Ltd
CATARC Tianjin Automotive Engineering Research Institute Co Ltd
Original Assignee
China Automotive Technology and Research Center Co Ltd
CATARC Tianjin Automotive Engineering Research Institute 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 China Automotive Technology and Research Center Co Ltd, CATARC Tianjin Automotive Engineering Research Institute Co Ltd filed Critical China Automotive Technology and Research Center Co Ltd
Priority to CN202110120476.9A priority Critical patent/CN112966407B/en
Publication of CN112966407A publication Critical patent/CN112966407A/en
Application granted granted Critical
Publication of CN112966407B publication Critical patent/CN112966407B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/26Composites

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Insertion Pins And Rivets (AREA)
  • Connection Of Plates (AREA)

Abstract

The invention provides a simulation method for a rubber riveting composite connection collision working condition, which comprises the following steps of S1: testing the mechanical property of pure viscose glue and calibrating simulation parameters, testing the mechanical property of riveting and calibrating simulation parameters, testing the mechanical property of rubber-riveting composite connection and correcting the simulation parameters; s2: correcting the combined riveting performance of the same material and different plate thicknesses; s3: a viscose defect simulation method is used for simulating and analyzing the collision working condition of the glue-rivet composite connection of the whole vehicle structure. Compared with the prior art, the simulation method for the collision working condition of the glue-riveting composite connection provided by the invention is applied to the whole vehicle structure, and the simulation method for the viscose defect is invented, so that the precision of a glue-riveting simulation model under the collision working condition can be greatly improved, and the accurate prediction of the deformation of the glue-riveting composite connection structure in the collision process can be realized.

Description

Simulation method for rubber riveting composite connection collision working condition
Technical Field
The invention belongs to the technical field of metal material connection performance analysis, and particularly relates to a glue-rivet composite connection collision condition simulation method.
Background
The steel-aluminum hybrid vehicle body structure is mainly applied to the multi-material hybrid vehicle body, and the steel-aluminum hybrid vehicle body structure faces the connection problem between steel aluminum and aluminum, wherein the SPR (self-piercing riveting) process is widely applied. Because the requirements of steel-aluminum connection on the performances of fatigue resistance, corrosion resistance, sealing property and the like are high, the composite connection technology, particularly the glue-rivet composite connection, is widely applied to a steel-aluminum mixed vehicle body. The glue-riveting composite connection has good connection performance, the rigidity performance of the automobile body can be greatly improved by applying the glue-riveting composite connection to the automobile body, and the lightweight design level of the automobile body can be effectively improved.
The influence of the glue-riveting composite connection on collision performance is uncertain, in a front cabin collapse energy absorption area, the deformation of a key energy absorption structure is greatly influenced by the mechanical performance and failure parameters of glue rivets, for example, an aluminum longitudinal beam is extruded, a lap plate is arranged on the side surface of the aluminum longitudinal beam, the lap plate and the longitudinal beam are in a glue-riveting composite connection mode, the glue-riveting composite connection can increase the deformation resistance of the lap area, and the crushing energy absorption of the longitudinal beam is influenced; in the side collision or side column collision working condition, in some specific joint parts, such as joints of the B column and a threshold, the B column and an upper boundary beam, a seat cross beam and a middle channel and the like, the high strength is needed to bear the impact of collision, the mechanical property and failure parameters of the glue-rivet composite connection play a role in lifting the side of the whole vehicle, the accurate glue-rivet simulation method can help to realize reasonable arrangement of the position and number of the glue-rivet composite connection, accurate prediction of the collision performance is realized, and the light work can be effectively unfolded.
The existing glue-riveting composite connection technology focuses on researching how to improve the static strength and the fatigue resistance of glue-riveting connection; or the rivets for riveting the finite element model are subjected to fine modeling, the grid size is too small, and the calculation efficiency is reduced sharply when the rivets are applied in a whole vehicle collision model; a numerical simulation method of glue-riveting connection is urgently needed in vehicle structure crashworthiness development, structural deformation prediction can be carried out on glue-riveting composite connection, and calculation efficiency is not reduced.
Disclosure of Invention
In view of this, the invention aims to provide a simulation method for a rubber-riveted composite connection collision working condition to solve the problems that the modeling precision of rubber-riveted composite connection is inaccurate, the collision performance prediction cannot be realized and the like in the vehicle collision.
In order to achieve the purpose, the technical scheme of the invention is realized as follows: the method comprises the following steps:
s1: the pure viscose glue mechanical property test and simulation parameter calibration specifically comprises the following steps:
s11: processing a pure viscose test sample;
s12: testing the mechanical property of the pure viscose sample;
s13: building a pure viscose simulation model;
s14: calibrating mechanical property and failure parameters of the pure viscose simulation model;
s2: and (4) testing the riveting mechanical property and calibrating simulation parameters.
S21: processing a riveting sample (a shearing working condition and a cross stretching working condition);
s22: testing the riveting mechanical properties (a shearing working condition and a cross stretching working condition);
s23: building a riveting simulation model;
s24: calibrating riveting mechanical property and failure parameters;
s3: testing the mechanical property of the glue-riveted composite connection and correcting simulation parameters;
s31: processing a glue-rivet composite connection sample (a shearing working condition and a cross stretching working condition);
s32: testing the mechanical properties (shearing condition and cross stretching condition) of the glue-riveting composite connection;
s33: constructing a glue-riveting composite connection simulation model;
s34: correcting the mechanical property and failure parameters of the glue-riveted composite connection;
s4: and (5) correcting the riveting performance of different plate thickness combinations made of the same material.
S41: correcting the riveting mechanical property;
s42: correcting riveting failure parameters;
s5: viscose defect simulation method.
S6: and (5) carrying out simulation analysis on the collision working condition of the glue-riveting composite connection of the whole vehicle structure.
Further, in the step S11,
the experimental working conditions comprise a shearing working condition and a cross stretching working condition;
further, in the step S14,
(1) The viscose adopts a solid 6-face solid grid.
(2) The mechanical property calibration mainly comprises a force rising stage and a peak force.
(3) The main index of the failure parameter calibration is failure loading displacement.
Further, in the step S23,
(1) And the riveting simulation modeling adopts a CONSTRAINED _ SPR2 keyword carried by DYNA. The normal load and the tangential load have the following relations with the rivet displacement:
normal direction:
Figure GDA0003890089610000031
tangential direction:
Figure GDA0003890089610000032
wherein:
Figure GDA0003890089610000037
f t max respectively indicating the maximum normal load and the maximum tangential load of the joint under the stretching and shearing working conditions;
Figure GDA0003890089610000033
respectively indicating the normal displacement and the tangential displacement when the rivet fails in the normal loading process and the tangential loading process; eta max Is a dimensionless value and is used to characterize the damage degree of the joint.
Figure GDA0003890089610000034
About eta under the stretching condition and the shearing condition respectively max The dimensionless function of (2), whose expression is as follows:
Figure GDA0003890089610000035
Figure GDA0003890089610000036
ξ n 、ξ t to describe the degree of damage to the joint under maximum tensile and maximum shear loads, respectively.
Further, in the step S24, the calibration indexes include a climbing stage, a force platform, and a fracture loading displacement.
Further, in the step S33,
1. modeling the viscose area. Dividing the viscose area into three parts; a first part: no viscose entity unit is established within the range of the size diameter of the rivet; a second part: a continuous entity viscose unit is established on one ring at the outermost side; and a third part: and arranging the viscose defects in the other viscose areas according to the defect rate of 5 percent of the total area. The size of the rubber grid is controlled to be about 1 mm. 2. And taking the mechanical property and the failure parameter of the viscose calibrated in the step S14 as simulation input parameters of the viscose. 3. Taking the riveting mechanical property parameters calibrated by S24 as simulation input parameters of riveting
Further, in the step S34,
(1) The mechanical property parameters of the viscose glue are corrected upwards, and the correction range is between 1.0 and 3.0; and correcting the failure displacement parameter downwards within the range of 0.2-1.0.
(2) Downward correction is carried out on the mechanical property parameters of riveting, and the correction range is between 0.7 and 1.0; and correcting the failure displacement parameter of the rivet downwards, wherein the correction range is between 0.7 and 1.0.
Further, in the step S41,
and (3) correcting the mechanical property of the rivet:
defining the rivet strength as FN;
(1) Measuring the undercut LH of a base sheet thickness combination base
(2) Measuring the undercut LH of a target sheet thickness combination target
(3) Correcting the mechanical property of the combined riveting of the target plate thickness based on the following empirical formula:
the strength of the rivet is as follows: FN (FN) device target =FN base α LH β T
Wherein alpha is LH Correction factor, beta, for undercut amount versus rivet mechanical properties T The thickness correction coefficient is obtained.
Wherein alpha is LH The calculation steps are as follows:
the undercut amount and riveting mechanical property empirical curve is as follows:
f(x)=FN max (C+Eln(1+GX))
wherein
X is the ratio of the undercut amount to the thickness of the underlying layer, and the value range is (0, 1), and the calculation formula is as follows:
Figure GDA0003890089610000041
wherein LH is the amount of undercut; t2 is the thickness of the lower plate;
FN max the maximum value of the riveting mechanical property at the moment X = 1;
c is an initial intensity coefficient at the moment when the undercut amount is X =1, and can be modified according to engineering experience;
E. g is a curve correction coefficient used for adjusting the slope of the undercut amount curve, and the following constraint is required to be met:
when X =1, C + Eln (1 + g) =1,
correction coefficient formula alpha of undercut amount to riveting mechanical property LH The following were used:
Figure GDA0003890089610000051
if LH origin If < 0.05, f (X) origin )=C
The empirical formula of the thickness of the upper plate and the lower plate to the riveting mechanical property correction is as follows:
defining T1 as the thickness of an upper plate (rivet head direction), and T2 as the thickness of a lower plate (rivet leg direction);
wherein beta is T The calculation formula of (c) is as follows:
Figure GDA0003890089610000052
further, in the step S42,
correcting rivet failure parameters:
correcting the normal failure displacement of the rivet:
Figure GDA0003890089610000053
wherein eta is a correction coefficient, correction is required according to actual conditions, and the numeric area is (0.8-1.2).
Correcting the tangential failure displacement of the rivet:
Figure GDA0003890089610000054
v is a correction coefficient, correction is required according to actual conditions, and the value range of v is (0.8-1.2).
Further, in the step S5,
the viscose area is divided into three parts;
a first part: no viscose unit is established within the diameter range of the rivet;
a second part: the viscose defect part is mainly divided into three rings from inside to outside, the first ring is a viscose defect-free area, the second ring is a viscose defect-less area, the viscose defect rate is 20 percent, and the viscose defect rate is adjusted according to the actual situation; the third ring is a region with more defects, the viscose defect rate is about 30-40%, and the adjustment is carried out according to the actual situation;
and a third part: and the non-adhesive area accounts for 5-10% of the total design area in the adhesive design range due to the process and other reasons and needs to be adjusted according to actual conditions.
Compared with the prior art, the simulation method for the collision working condition of the glue-riveting composite connection, provided by the invention, has the advantages that reliable mechanical properties such as gluing and riveting and failure parameters are obtained through simulation calibration, and parameter correction is carried out through glue-riveting mechanical property experiments and failure experiments; in the application of a whole vehicle structure, the mechanical property and failure parameters of the glue-riveting composite connection combined by different thicknesses of the same material are empirically corrected, and the glue defect simulation method is invented, so that the precision of a glue-riveting simulation model under a collision working condition can be greatly improved, and the accurate prediction of the deformation of the glue-riveting composite connection structure in the collision process is realized.
Drawings
FIG. 1 is a block diagram of the system architecture of the present invention;
FIG. 2 is a schematic view of a lap joint (shear condition) in an embodiment of the invention;
FIG. 3 is a schematic diagram of cross stretching (pull-out condition) in an embodiment of the present invention;
FIG. 4 is a simulation calibration chart of the mechanical properties and failure parameters of the glue riveting lap joint in the embodiment of the invention;
FIG. 5 is a force displacement curve for a viscose lap tensile test;
FIG. 6 is a viscose simulation model;
FIG. 7 is a simulated force displacement curve of the viscose mechanical property simulated calibration;
FIG. 8 is a force displacement curve for a rivet lap joint test;
FIG. 9 is a calibration curve of the mechanical properties of the rivet;
FIG. 10 is a diagram showing a force-displacement curve of the glue-rivet composite connection;
FIG. 11 is a schematic view of a modeling plane for the glue area;
FIG. 12 shows the calibration results of the glue-rivet composite connection;
FIG. 13 is a schematic view of a glue defect plane (full car modeling);
Detailed Description
Unless defined otherwise, technical terms used in the following examples have the same meanings as commonly understood by one of ordinary skill in the art to which the present invention belongs. The test reagents used in the following examples, unless otherwise specified, are all conventional biochemical reagents; the experimental methods are conventional methods unless otherwise specified.
The present invention will be described in detail below with reference to the following examples and the accompanying drawings.
Taking 6082 aluminum plates for glue-riveting composite connection as an example, the method for simulating the collision working condition of the glue-riveting composite connection provided by the invention comprises the following steps:
1. pure viscose glue mechanical property test and simulation parameter calibration
1.1 prepare a bonding template, the area of the glued area is consistent with the average area of the glue of the actual whole vehicle (neglecting the rivet and the area of the non-glued area around the rivet).
The adhesive sheet manner is divided into lap joint and cross form, as shown in fig. 2 and 3.
The thickness of the upper and lower adhesive aluminum plates is 2.7mm.
And performing lap joint stretching and cross stretching tests on the viscose board to obtain a test force displacement curve.
The force displacement curve for the lap joint test is shown in figure 5 below:
1.2, building a finite element model according to the sample size and the boundary conditions in the step 1.1, wherein the finite element model mainly comprises the following contents.
Drawing the middle surface of the CAD data of the upper aluminum plate and the lower aluminum plate, drawing grids, wherein the average size of the grids is 2mm, the material is 6082, and the thickness is 2.7mm.
And (3) carrying out solid grid modeling on the gluing area, wherein the gluing area adopts a solid 6-face solid grid, the material adopts MAT-SPOTWLED cards, the initial gluing strength and gluing failure parameters are set, and a gluing simulation model is shown in figure 6.
1.3, simulating the mechanical property of the viscose by referring to the experimental working condition, comparing the simulation with the experimental force displacement curve, and carrying out simulation calibration by adjusting the elastic modulus, the strength parameter and the failure parameter of the viscose to obtain a force displacement curve with better consistency with the experiment, wherein the calibrated force displacement curve is shown in fig. 7.
2. Riveting mechanical property test and simulation parameter calibration:
2.1 rivet mechanical Property test
(1) Determining that the size of the rivet is consistent with that of the actual vehicle;
(2) The thickness of the riveted upper and lower aluminum plates is 2.7mm, and the riveting form is lap joint and cross, as shown in figures 2 and 3.
(3) And performing lap tensile and cross tensile tests on the riveting plate to obtain a test force displacement curve, as shown in fig. 8.
2.2 build SPR riveting model to carry out emulation and mark, mainly include:
(1) The position of the center of the welding spot is determined.
(2) The CONSTRAIN _ SPR2 card is created and mainly comprises the following parameters:
MID, the number of the main board part; SID: numbering slave plates Part; NSID: a set of nodes positioned by rivets; THICK: the total thickness of the master and slave plates; d: the diameter of the rivet; FN: rivet pull-out strength; FT: pure shear rivet strength; DN: normal destructive displacement; DT: tangential destructive displacement; XIN: fraction of destructive displacement (proportion) at maximum normal force; XIT: fraction of destructive displacement at maximum tangential force (proportion);
(3) 6082 is adopted as the material of the riveting plate.
(4) And defining the boundary condition of the simulation model according to the experimental working condition of the riveting mechanical property.
(5) By adjusting the peak load and the peak load displacement under the pure stretching and pure shearing working conditions, a force displacement curve consistent with the experiment is obtained, and a calibrated simulation force displacement curve is shown in fig. 9.
3. Mechanical property test and simulation parameter correction of glue-riveting composite connection
3.1 testing mechanical properties of glue-riveted composite connection
(1) Preparing a glue riveting test piece, wherein the requirements of the glue riveting test piece are as follows:
and the adhesive area of the glue riveting test piece is consistent with the average adhesive area of the real vehicle. The size of the rivet is consistent with that of a real vehicle. The glue riveting plate is 2.7mm and is made of 6082;
(2) The glue riveting plates are in lap joint and cross form.
(3) And (3) carrying out lap joint stretching and cross stretching working condition experiments on the glue riveting test sample to obtain a force displacement curve, wherein a force displacement curve of the glue riveting composite connection experiment is shown in figure 10.
3.2 glue-riveted composite connection simulation parameter correction
(1) Rubber riveting simulation model building
1) Modeling the viscose area.
Dividing the viscose area into three parts;
a first part: no viscose entity unit is established within the diameter range of the rivet size;
a second part: a continuous entity viscose unit is established on one ring at the outermost side;
and a third part: and arranging the viscose defects in the other viscose areas according to the defect rate of 5 percent of the total area.
The size of the viscose mesh is controlled to be about 1 mm.
The plan view of the adhesive area is shown in fig. 11 below:
2) Riveting modeling reference step 2.2
(2) And (3) taking the viscose glue and riveting mechanical property parameters and failure parameters calibrated in the steps 1 and 2 as initial input of a glue riveting simulation model.
(3) And (5) correcting the mechanical property and the failure parameter of the adhesive.
Comparing the force displacement curve with the experimental force displacement curve, and calibrating the force displacement curve before the viscose fails by adjusting the strength parameter and the failure parameter of the viscose; the correction process was based on the following experience:
the mechanical property parameters of the viscose glue are corrected upwards, and the correction range is between 1.0 and 3.0; and correcting the failure displacement parameter downwards within the range of 0.2-1.0.
(4) And (4) correcting riveting mechanical property and failure parameters.
And comparing the force displacement curve with the experimental force displacement curve, and correcting the riveting mechanical property and failure parameters by adjusting the strength parameters and failure parameters of the rivet.
Downward correction is carried out on the mechanical property parameters of riveting, and the correction range is between 0.7 and 1.0; the failure displacement parameter of the rivet is corrected downwards, the correction range is between 0.7 and 1.0, and the calibrated force displacement curve is shown in a 12-figure.
(5) After correction, the force displacement curve in the simulation and the force displacement curve of the experiment keep better consistency, and the attached figure 4 shows that the calibration result of the force displacement curve of the glue riveting lap joint under the lap joint working condition can be used as a reference. And obtaining mechanical property parameters and failure parameters of the glue riveting composite connection available for engineering.
4. And (5) correcting the riveting performance of the same material and different plate thickness combinations.
In the application of a whole vehicle structure, riveting thickness combinations of the same materials may be different, and the glue-riveting composite connection performance of the same materials with different thicknesses needs to be corrected.
4.1 correcting the mechanical property of the rivet:
defining rivet strength FN; measuring the undercut LH of the base sheet thickness combination base (ii) a Measuring the undercut LH of a target sheet thickness combination target
Correcting the mechanical property of the combined riveting of the target plate thickness based on the following empirical formula:
rivet strength: fN target =FN base α LH β T
Wherein alpha is LH Correction factor, beta, for undercut amount versus riveting performance T The thickness correction coefficient is obtained.
(1) Wherein alpha is LH The calculation steps are as follows:
the undercut amount and riveting mechanical property empirical curve is as follows:
f(x)=FN max (C+Eln(1+GX))
wherein X is the ratio of the undercut amount to the thickness of the bottom layer plate, the value range is (0, 1), and the calculation formula is as follows:
Figure GDA0003890089610000091
wherein LH is the amount of undercut; t2 is the thickness of the lower plate;
FN max the maximum value of the riveting mechanical property at the moment X = 1;
c is an initial intensity coefficient at the moment when the undercut amount is X =1, which can be modified according to engineering experience;
E. g is a curve correction coefficient used for adjusting the slope of the undercut amount curve, and the following constraints need to be satisfied:
when X =1, C + Eln (1 + g) =1,
correction coefficient formula alpha of undercut amount to riveting mechanical property LH The following were used:
Figure GDA0003890089610000092
if LH origin < 0.05, then f (X) origin )=C
(2) The empirical formula of the thickness of the upper plate and the lower plate to the riveting mechanical property correction is as follows:
defining T1 as the thickness of an upper plate (in the direction of a rivet head) and T2 as the thickness of a lower plate (in the direction of a rivet leg);
wherein beta is T The calculation formula of (c) is as follows:
Figure GDA0003890089610000101
4.2 rivet failure parameter modification:
(1) Correcting the normal failure displacement of the rivet:
Figure GDA0003890089610000102
wherein eta is a correction coefficient, correction is needed according to actual conditions, and the value range is (0.8-1.2).
(2) Correcting the tangential failure displacement of the rivet:
Figure GDA0003890089610000103
v is a correction coefficient, correction is needed according to actual conditions, and the value range of v is (0.8-1.2).
5. Viscose defect simulation method.
In the glue riveting composite connection of the whole vehicle structure, the glue defect exists, and the glue defect is simulated through the following steps.
The viscose area is equivalently divided into three major parts;
a first part: no viscose unit is established within the diameter range of the rivet;
a second part: the viscose defect part is mainly divided into three rings from inside to outside, the first ring is a viscose defect-free area, the second ring is a viscose defect-less area, the viscose defect rate is 20 percent, and the viscose defect rate is adjusted according to the actual situation; the third ring is a region with more defects, the viscose defect rate is about 30-40%, and the adjustment is carried out according to the actual situation;
and a third part: and the non-adhesive area accounts for 5-10% of the total design area in the adhesive design range due to the process and other reasons and needs to be adjusted according to actual conditions. In the entire vehicle modeling, the schematic diagram of the viscose defect plane is shown in fig. 13:
6. and (5) carrying out simulation analysis on the collision working condition of the glue-riveting composite connection of the whole vehicle structure.
And according to the actual condition of the whole vehicle, establishing a glue-riveting composite connection model of the whole vehicle structure by referring to the steps.
With the increasing of the light-weight requirement of automobiles, the glue-rivet composite connection is more and more widely applied to the automobile body. The accurate glue riveting simulation method can help realize reasonable arrangement of the positions and the number of glue riveting composite connections, realize accurate prediction of collision performance, and ensure that light-weight work can be effectively unfolded.
The invention provides a simulation method for a collision working condition of glue-riveting composite connection, which is characterized in that simulation calibration is respectively carried out on the mechanical properties and failure parameters of glue and riveting, and then the mechanical properties and failure parameters of the glue-riveting composite connection are corrected and simulated calibration is carried out on the basis. In the application of whole vehicle collision simulation analysis, a riveting performance correction empirical formula and a viscose defect simulation method of the same material with different thicknesses are defined, the precision of a glue riveting simulation model under a collision working condition can be greatly improved, and the accurate prediction of deformation of a glue riveting composite connection structure in a collision process is realized. And according to the actual condition of the whole vehicle, establishing a glue-riveting composite connection model of the whole vehicle structure by referring to the steps.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and should not be taken as limiting the scope of the present invention, which is intended to cover any modifications, equivalents, improvements, etc. within the spirit and scope of the present invention.

Claims (9)

1. A glue-riveting composite connection collision working condition simulation method is characterized by comprising the following steps: the method comprises the following steps:
s1: testing the mechanical property of pure viscose glue and calibrating simulation parameters, testing the mechanical property of riveting and calibrating simulation parameters, testing the mechanical property of rubber-riveting composite connection and correcting the simulation parameters;
s2: correcting the combined riveting performance of different plate thicknesses of the same material;
s3: a viscose defect simulation method, which is used for simulating and analyzing the collision working condition of the glue-rivet composite connection of the whole vehicle structure;
in the correction of the mechanical properties of the rivet:
correcting the normal failure displacement of the rivet:
Figure FDA0003956333150000011
wherein eta is a correction coefficient, correction is needed according to actual conditions, and the value range of eta is 0.8-1.2;
DN target is a riveting target normal failure displacement;
T2 target target lower plate thickness;
T2 origin is the initial lower plate thickness;
DN origin the initial riveting normal failure displacement is realized;
correcting the tangential failure displacement of the rivet:
Figure FDA0003956333150000012
wherein
DT target A target rivet tangential failure displacement is obtained;
T2 target a target lower plate thickness;
T2 origin is the initial lower plate thickness;
DN origin the initial rivet tangential failure displacement;
v is a correction coefficient, correction is needed according to actual conditions, and the value range of v is 0.8-1.2;
and (3) taking the pure viscose glue, the riveted mechanical property parameters and the failure parameters calibrated in the step (S1) as initial inputs of a glue riveting simulation model.
2. The glue-riveted composite connection collision working condition simulation method according to claim 1, characterized in that: s1, processing a pure viscose test sample, and performing a mechanical property experiment on the pure viscose test sample; building a pure viscose simulation model; calibrating the mechanical property and failure parameters of the pure viscose simulation model; carrying out experimental working condition processing on the riveting sample; carrying out mechanical property test on the processed riveting sample, building a riveting simulation model, and calibrating the riveting mechanical property and failure parameters; processing the experimental working condition of the glue-riveting composite connection sample, testing the mechanical property of the processed glue-riveting composite connection, building a simulation model of the glue-riveting composite connection, and correcting the mechanical property and the failure parameter of the glue-riveting composite connection.
3. The glue-riveted composite connection collision working condition simulation method according to claim 2, characterized in that: the sample adhesive area in the pure adhesive test sample processing should refer to the range covered by the actual single adhesive riveting structure adhesive, and the area of the rivet area is ignored.
4. The glue riveting composite connection collision working condition simulation method according to claim 2, characterized in that: the experimental working conditions in the pure viscose test sample processing include a shearing working condition and a cross stretching working condition, a bonding sample plate is prepared, the area of a viscose area is consistent with the average area of actual complete viscose, the mode of the viscose plate is divided into a lap joint mode and a cross form, lap joint stretching and cross stretching tests are carried out on the viscose plate to obtain a test force displacement curve, solid grid modeling is carried out on the viscose area in the calibration of the mechanical property and failure parameters of a pure viscose simulation model, the viscose adopts a solid 6-face solid grid, initial viscose strength and viscose failure parameters are set, the calibration of the mechanical property and the mechanical property of the pure viscose simulation model comprises a force rising stage and a peak force in the calibration of the failure parameters, the calibration index of the failure parameters is failure loading displacement, the simulation of the mechanical property of the viscose is carried out according to the experimental working conditions, the simulation calibration is compared with the experimental force displacement curve, the simulation calibration is carried out by adjusting the elastic modulus, the strength parameters and the viscose failure parameters to obtain a force displacement curve which is consistent with the experiment, the riveting simulation modeling adopts a DYNA self-carried CONSTRAINEED _ 2 keyword, and the normal load and the rivet displacement have the following relationship:
normal direction:
Figure FDA0003956333150000031
tangential direction:
Figure FDA0003956333150000032
wherein:
Figure FDA0003956333150000033
f t max respectively indicating the maximum normal load and the maximum tangential load of the joint under the stretching and shearing working conditions;
Figure FDA0003956333150000034
respectively indicating the normal displacement and the tangential displacement when the rivet fails in the normal loading process and the tangential loading process; eta max Is a dimensionless value used to characterize the damage degree of the joint;
Figure FDA0003956333150000035
about eta under the stretching condition and the shearing condition respectively max The dimensionless function of (2), whose expression is as follows:
Figure FDA0003956333150000036
Figure FDA0003956333150000037
ξ n 、ξ t respectively describing the damage degree of the joint under the maximum tensile load and the maximum shearing load;
the calibration indexes in the calibration of riveting mechanical property and failure parameters are divided into a climbing stage, a force platform and fracture loading displacement.
5. The glue-riveted composite connection collision working condition simulation method according to claim 2, characterized in that: the method comprises the following steps of establishing a glue area region in the building of a glue-riveting composite connection simulation model, dividing the glue area into three parts, wherein the first part is a part which does not establish a glue entity unit within the diameter range of the rivet size, and the second part is a part which establishes a continuous entity glue unit for one ring at the outermost side; the third part is the rest viscose areas, viscose defects are arranged according to the defect rate of 5% of the total area, the size of a viscose grid is controlled to be 0.8mm-1.2mm, viscose mechanical property and failure parameters calibrated by a pure viscose simulation model are used as simulation input parameters of viscose, and riveting mechanical property parameters calibrated by riveting mechanical property are used as simulation input parameters of riveting.
6. The glue riveting composite connection collision working condition simulation method according to claim 2, characterized in that: the mechanical property and the failure parameter of the glue-riveting composite connection are corrected to be upward corrected according to the mechanical property parameter of the viscose glue, and the correction range is between 1.0 and 3.0; downward correction is carried out on the failure displacement parameter, the correction range is between 0.2 and 1.0, and downward correction is carried out on the mechanical property parameter of riveting, and the correction range is between 0.7 and 1.0; and correcting the failure displacement parameter of the rivet downwards, wherein the correction range is between 0.7 and 1.0.
7. The glue riveting composite connection collision working condition simulation method according to claim 1, characterized in that: and S2, correcting riveting mechanical properties and riveting failure parameters of the combined riveting of different plate thicknesses of the same material.
8. The glue-riveted composite connection collision working condition simulation method according to claim 7, characterized in that: during the correction of the mechanical property of the rivet:
defining rivet strength FN:
(1) Measuring the undercut LH of a base sheet thickness combination origin
(2) Measuring the undercut LH of a target sheet thickness combination target
(3) Correcting the mechanical property of the combined riveting of the target plate thickness based on the following empirical formula:
rivet strength: FN (FN) target =FN origin α LH β T
FN target Target rivet strength;
FN origin the strength of the basic rivet;
wherein alpha is LH Correction factor, beta, for undercut amount versus riveting performance T The thickness correction coefficient is used;
wherein alpha is LH The calculation steps are as follows:
the undercut amount versus riveting mechanics empirical curve is as follows:
f(x)=FN max (C+Eln(1+GX))
wherein X is the ratio of the undercut amount to the thickness of the underlying layer, the value range is (0, 1), and the calculation formula is as follows:
Figure FDA0003956333150000051
wherein LH is the amount of undercut; t2 is the thickness of the lower plate; FN (FN) max The maximum value of the riveting mechanical property at the moment of X = 1; c is an initial intensity coefficient at the moment of X =1, and is modified according to engineering experience;
E. g is a curve correction coefficient for adjusting the slope of the undercut amount curve,
the following constraints need to be satisfied:
when X =1, C + Eln (1 + g) =1,
correction coefficient formula alpha of undercut amount to riveting mechanical property LH The following were used:
Figure FDA0003956333150000052
if LH target < 0.05, then f (X) origin )=C
X origin The ratio of the initial undercut amount to the thickness of the underlayer sheet;
X target is the ratio of the target undercut amount to the underlayer plate thickness;
LH target the amount of undercut combined for measuring the target sheet thickness.
The empirical formula for correcting the riveting mechanical property by the thickness of the upper plate and the lower plate is as follows: defining T1 as the thickness of the upper plate in the direction of the rivet head, and T2 asThe thickness of the lower plate in the direction of the rivet leg; wherein beta is T The calculation formula of (a) is as follows:
Figure FDA0003956333150000053
T1 target target upper plate thickness;
T2 target a target lower plate thickness;
T1 origin is the initial upper plate thickness;
T2 origin is the initial lower plate thickness.
9. The glue riveting composite connection collision working condition simulation method according to claim 1, characterized in that: the viscose area in the S3 is equivalently divided into three parts; a first part: no viscose unit is established within the diameter range of the rivet; a second part: the viscose defect part is divided into three rings from inside to outside, the first ring is a viscose defect-free area, the second ring is a viscose defect-less area, the viscose defect rate is 20 percent, and the viscose defect rate is adjusted according to the actual situation; the third ring is a region with more defects, the viscose defect rate is 30-40%, and the third ring is adjusted according to the actual situation; and a third part: and the non-adhesive area accounts for 5% -10% of the total design area in the adhesive design range due to process reasons and needs to be adjusted according to actual conditions.
CN202110120476.9A 2021-01-28 2021-01-28 Simulation method for rubber riveting composite connection collision working condition Active CN112966407B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110120476.9A CN112966407B (en) 2021-01-28 2021-01-28 Simulation method for rubber riveting composite connection collision working condition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110120476.9A CN112966407B (en) 2021-01-28 2021-01-28 Simulation method for rubber riveting composite connection collision working condition

Publications (2)

Publication Number Publication Date
CN112966407A CN112966407A (en) 2021-06-15
CN112966407B true CN112966407B (en) 2023-02-28

Family

ID=76271755

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110120476.9A Active CN112966407B (en) 2021-01-28 2021-01-28 Simulation method for rubber riveting composite connection collision working condition

Country Status (1)

Country Link
CN (1) CN112966407B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107202756A (en) * 2017-05-12 2017-09-26 简式国际汽车设计(北京)有限公司 A kind of simulating analysis of glass viscose glue failure
CN110909431A (en) * 2019-12-02 2020-03-24 合肥工业大学智能制造技术研究院 Design method of dissimilar material connecting structure for automobile and connecting structure thereof
CN111651881A (en) * 2020-06-01 2020-09-11 中国第一汽车股份有限公司 Method for simplifying lock riveting simulation failure parameters
CN111680358A (en) * 2020-05-09 2020-09-18 中国第一汽车股份有限公司 Collision simulation method for automobile aluminum alloy section parts

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107415345A (en) * 2017-07-11 2017-12-01 昆明理工大学 A kind of foamed metal sandwich board self-piercing riveting manufacturing process

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107202756A (en) * 2017-05-12 2017-09-26 简式国际汽车设计(北京)有限公司 A kind of simulating analysis of glass viscose glue failure
CN110909431A (en) * 2019-12-02 2020-03-24 合肥工业大学智能制造技术研究院 Design method of dissimilar material connecting structure for automobile and connecting structure thereof
CN111680358A (en) * 2020-05-09 2020-09-18 中国第一汽车股份有限公司 Collision simulation method for automobile aluminum alloy section parts
CN111651881A (en) * 2020-06-01 2020-09-11 中国第一汽车股份有限公司 Method for simplifying lock riveting simulation failure parameters

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
碳/玻璃纤维混编复合材料和钢板的单搭接接头的力学性能分析与预测;陈洋等;《重庆大学学报》;20200215(第02期);第68-80页 *
自冲铆连接性能分析及在车身前纵梁碰撞仿真中的应用;解宇;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》;20200115(第01期);第42-44页 *
陈洋等.碳/玻璃纤维混编复合材料和钢板的单搭接接头的力学性能分析与预测.《重庆大学学报》.2020,(第02期), *

Also Published As

Publication number Publication date
CN112966407A (en) 2021-06-15

Similar Documents

Publication Publication Date Title
Nguyen et al. Simulation of impact on sandwich structures
He et al. Numerical and experimental investigations of self-piercing riveting
CN110287550A (en) White body solder joint optimization method based on density variable method and analysis of Fatigue-life
Moraes et al. Effect of process deformation history on mechanical performance of AM60B to AA6082 self-pierce riveted joints
CN108984909B (en) Mar-Lin model-based method for analyzing residual strength of composite material structure of airplane containing large damage
Baker et al. An affordable methodology for replacing metallic aircraft panels with advanced composites
CN112966407B (en) Simulation method for rubber riveting composite connection collision working condition
Rahman et al. Fatigue life prediction of spot-welded structures: a finite element analysis approach
CN116882232A (en) Method for predicting damage state of carbon fiber composite material laminated plate under low-speed impact
Correia et al. Parametric study of composite curved adhesive joints
Rans The role of rivet installaton on the fatigue performance of riveted lap joints
Okafor et al. Design and analysis of adhesively bonded thick composite patch repair of corrosion grind-out and cracks on 2024 T3 clad aluminum aging aircraft structures
Boscolo et al. A modelling technique for calculating stress intensity factors for structures reinforced by bonded straps. Part I: Mechanisms and formulation
Wagner et al. A simplified and semi-analytical bolted joint model for crash and impact simulations of composite structures
Grujicic et al. Self-piercing riveting process and joint modeling and simulations
Boscolo et al. Design and modelling of selective reinforcements for integral aircraft structures
Fleuret et al. Complex wing spar design in carbon fiber reinforced composite for a light aerobatic aircraft
Kabir et al. Slow growth damage assessment of bonded metal-composite joint using a numerical approach
Folgar et al. Joining TWIP-Steel simulation models
Viganò et al. Investigation on the fatigue life of an Al 8090 helicopter riveted T-joint
Taki Optimization of flat Z-stiffened panel subjected to compression
Klimmek et al. Aircraft Loads–A Wide Range of Disciplinary and Process-Related Issues in Simulation and Experiment
Trego et al. Evaluation of damage tolerance analysis tools for lap joints
Heng et al. Prediction of Damage Tolerance in Metallic Structure Repaired with a Co-Cured Composite Patch
Ekladious et al. Numerical and experimental studies on the joint performance of fastened, bonded, and hybrid thin metal joints used in aircraft structures

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