CN108664731A - A kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller - Google Patents
A kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller Download PDFInfo
- Publication number
- CN108664731A CN108664731A CN201810448171.9A CN201810448171A CN108664731A CN 108664731 A CN108664731 A CN 108664731A CN 201810448171 A CN201810448171 A CN 201810448171A CN 108664731 A CN108664731 A CN 108664731A
- Authority
- CN
- China
- Prior art keywords
- model
- curing
- composite material
- resin
- residual stress
- 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
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000004088 simulation Methods 0.000 title claims abstract description 35
- 239000011347 resin Substances 0.000 claims abstract description 95
- 229920005989 resin Polymers 0.000 claims abstract description 91
- 238000004364 calculation method Methods 0.000 claims abstract description 65
- 239000000835 fiber Substances 0.000 claims abstract description 54
- 230000008878 coupling Effects 0.000 claims abstract description 26
- 238000010168 coupling process Methods 0.000 claims abstract description 26
- 238000005859 coupling reaction Methods 0.000 claims abstract description 26
- 230000008569 process Effects 0.000 claims abstract description 25
- 239000000463 material Substances 0.000 claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 10
- 230000009477 glass transition Effects 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 238000007711 solidification Methods 0.000 claims description 3
- 230000008023 solidification Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 2
- 238000010008 shearing Methods 0.000 claims 1
- 238000004458 analytical method Methods 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 4
- 230000020169 heat generation Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000004913 activation Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000001808 coupling effect Effects 0.000 description 2
- 238000000465 moulding Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Reinforced Plastic Materials (AREA)
- Moulding By Coating Moulds (AREA)
Abstract
The invention discloses a kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller, specifically include following process：Step 1, composite material macroscopic view multi- scenarios method computation model is established；Step 2, macroscopical multi- scenarios method computation model initial value and calculation of boundary conditions are assigned；Step 3, mesh generation is carried out to the model established in step 1, and finite element solving is carried out to model；Step 4, the thin thermodynamic parameter for seeing resin material constitutive model in Googol motion controller model is calculated according to step 3 acquired results；Step 5, the RVE models for including single fiber and resin matrix are established using representative volume unit, using the thermodynamic parameter obtained by step 4 as input parameter, the thin FEM calculation for seeing Googol motion controller model is carried out, the multi-scale Simulation of Googol motion controller is completed.Method of the present invention by the way that multi-field coupling theory and multi-scale Simulation to be combined so that the prediction that composite material carefully sees residual stress is more accurate.
Description
Technical Field
The invention belongs to the technical field of composite material design and manufacture, and relates to a multi-scale numerical simulation method for curing residual stress of a composite material.
Background
The composite material has the advantages of high specific strength, large specific modulus, corrosion resistance, fatigue resistance, strong designability and the like, and is widely applied to the fields of aviation, aerospace, ships, weapons, buildings, medical treatment and the like. At present, domestic enterprises mainly adopt an autoclave molding process to prepare composite material components. In the autoclave molding process, a composite material member is subjected to high temperature, high pressure and other processes, and due to factors such as anisotropy of the material, chemical shrinkage reaction of a resin matrix, mold action and the like, curing residual stress can be generated inside the composite material, so that subsequent use and assembly connection of the member are seriously influenced.
The curing residual stress can be classified into macroscopic residual stress and microscopic residual stress. The macroscopic residual stress is the residual stress formed between the layers of the composite material and can cause the composite material to generate rebound deformation after demoulding; the microscopic residual stress refers to the residual stress formed among the resin, the fiber and the interface between the resin and the fiber, which easily causes microcracks on a matrix or an interface phase, further causes initial damage and seriously affects the mechanical property of the composite material.
The composite material is prepared by combining two or more materials with different physical and chemical properties through different structural scales and layers such as microscopic, mesoscopic and macroscopic structures and curing in a high-temperature and high-pressure mode. Due to the characteristics of the multi-scale (microscopic-macroscopic), multi-phase (reinforcing phase, matrix phase and interlayer phase) material and the forming characteristics of multi-field coupling (parameters such as temperature, curing degree, viscosity, fiber volume fraction and stress), the solution of microscopic curing residual stress becomes very complicated. The existing research on the microscopic curing residual stress mostly adopts a single fiber to carry out analytical modeling calculation or selects a plurality of fibers to establish an RVE model to carry out finite element simulation calculation, and the influence of the whole size of a composite material component and resin flow is ignored. Therefore, the mesoscopic residual stress calculated by the existing numerical simulation method has a large error with an actual result, and the reliability of the result is not high.
In 2016, a microscopic residual stress calculation model considering a fiber-matrix interface phase is established by utilizing a classical laminate theory and an energy method, and the model analyzes the influence of fiber volume fraction on curing residual stress, but does not give the influence on how to calculate the fiber volume fraction, and does not consider the influence of the overall structure size of a laminate in the model.
In a paper of A composite of curing process-induced stress and current shrinkage structure with micro-scale composite structure with a thickness consistent law in 2017, Dongna Li studied the influence of the composite constitutive model on the fine curing residual stress, but the influence of resin flow and the thickness dimension of the laminated plate is not considered in the model.
The microscopic curing residual stress easily causes the problems of microcracks and the like in the composite material member, the forming quality and the subsequent use performance of the composite material member are seriously influenced, and the method is very important for accurately predicting the development rule of the microscopic curing residual stress. Therefore, when calculating the microscopic curing residual stress by the numerical simulation method, the influence of the macrostructure size of the laminate and the resin flow needs to be considered.
Disclosure of Invention
The invention aims to provide a multi-scale numerical simulation method for curing residual stress of a composite material, which is based on the idea of combining macroscopical and microscopic scales and enables the prediction of the microscopic residual stress of the composite material to be more accurate by combining a multi-field coupling theory and multi-scale simulation.
The technical scheme adopted by the invention is that the composite material curing residual stress multi-scale numerical simulation method specifically comprises the following processes:
step 1, establishing a macroscopic multi-field coupling calculation model of the composite material, wherein the macroscopic multi-field coupling calculation model comprises a thermochemical calculation model, a resin viscosity calculation model and a resin flow calculation model;
step 2, giving a macroscopic multi-field coupling calculation model initial value and calculating the boundary condition of the model;
step 3, carrying out mesh division on the model established in the step 1 by adopting entity units in finite element software, and carrying out finite element solution on the model;
step 4, calculating thermodynamic parameters of a constitutive model of the microscopic cured residual stress resin material according to the result obtained in the step 3, wherein the thermodynamic parameters comprise elastic modulus, shear modulus and Poisson ratio;
and 5, establishing an RVE model containing single fibers and a resin matrix by adopting a representative volume unit, taking the thermodynamic parameters obtained in the step 4 as input parameters, and performing finite element calculation on a microscopic curing residual stress model to complete multi-scale simulation of the curing residual stress.
The present invention is also characterized in that,
the specific process of step 3 is as follows:
calculating the internal temperature T and the curing degree a of the composite material through the thermochemical calculation model in the step 1, obtaining the internal resin viscosity η of the composite material through the obtained temperature T and the curing degree a and by utilizing a resin viscosity calculation model, and updating the internal fiber volume fraction V of the composite material through a resin flow model by utilizing the obtained temperature T, the obtained curing degree a and the obtained resin viscosity ηfTemperature T, degree of cure a, fiber volume fraction V obtained by calculationfUpdating model input parameters, using the updated model input parameters as initial values of the model at the next moment, completing the curing of the composite material in the whole circulation process, and extracting corresponding temperature T, curing degree a and fiber volume fraction V after the curing is completedf。
The specific process of step 4 is as follows:
reading the temperature T, the solidification degree a and the fiber volume fraction V which are calculated in the step 3fAnd 3, calculating to obtain the elastic modulus E in the resin curing process by taking the temperature T and the curing degree a obtained by calculation in the step 3 as input parameters of a constitutive model of the microscopic curing residual stress resin materialmShear modulus GmAnd poisson's ratio vmAs shown in formulas (7), (8) and (9), thereby realizing information transfer from the composite material macroscopic multi-field coupling model to the microscopic residual stress model;
wherein,andrespectively, the resin at the corresponding temperature TC1、TC2The modulus of elasticity of the rubber composition;poisson's ratio, T, after complete curing of the resin∞For the current glass transition temperature T of the resingDifference from the current temperature T of the resin:
T∞＝Tg-T (10)；
Tg＝164.6a2+51a+2.67 (11)。
the numerical simulation method provided by the invention has the beneficial effects that the influence of the macroscopic structure size of the composite material and the resin flow on the microscopic curing residual stress of the composite material is considered, the data interaction between the macroscopic multi-field coupling calculation model and the microscopic residual stress model is adopted, the strong coupling effect among the temperature, the curing degree and the fiber volume fraction is considered, the problem of residual stress calculation deviation in the background technology is effectively solved, the reliability of the residual stress simulation result is improved, the simulation of the residual stress is more practical, and a foundation is laid for the subsequent strength analysis.
Drawings
FIG. 1 is a block diagram of components used in an embodiment of a multi-scale numerical simulation method of residual stress in curing of a composite material according to the present invention;
FIG. 2 is a finite element model used in an embodiment of a multi-scale numerical simulation method of residual stress in curing of a composite material according to the present invention;
FIG. 3 is a schematic diagram illustrating a temperature calculation result of a macroscopic multi-field coupling model in an embodiment of a multi-scale numerical simulation method for residual stress in curing of a composite material according to the present invention;
FIG. 4 is a schematic diagram illustrating a calculation result of a degree of cure of a macroscopic multi-field coupling model in an embodiment of a multi-scale numerical simulation method for residual stress in curing of a composite material according to the present invention;
FIG. 5 is a temperature profile for an embodiment of a multi-scale numerical simulation method of residual stress in curing of a composite material according to the present invention;
FIG. 6 is a graph of the change of the degree of cure in an embodiment of a multi-scale numerical simulation method of residual stress in curing of a composite material according to the present invention;
FIG. 7 is a graph of fiber volume fraction change in an embodiment of a multi-scale numerical simulation of residual stress for curing of a composite material according to the present invention;
FIG. 8 is an RVE model of fibers and resins in an embodiment of a composite cure residual stress multi-scale numerical simulation method of the present invention;
FIG. 9 is a schematic representation of RVE model cure residual stress calculations for fibers and resins in an embodiment of a composite cure residual stress multi-scale numerical simulation method of the present invention;
FIG. 10 is a comparative graph of residual stress curves at point A of a selected component in an embodiment of a multi-scale numerical simulation of residual stress for curing of a composite material according to the present invention, respectively through a numerical simulation of an embodiment of the present invention and a comparative example.
In the figure, 1 is fiber, and 2 is resin.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
The invention relates to a multi-scale numerical simulation method for curing residual stress of a composite material, which specifically comprises the following steps:
step 1, establishing a composite material macroscopic multi-field coupling calculation model by using Comsol finite element calculation software;
the macroscopic multi-field coupling calculation model comprises a thermochemical calculation model, a resin viscosity calculation model and a resin flow calculation model;
the thermochemical calculation model is shown in the following formula (1):
wherein λ isx、λy、λzThe thermal conductivity coefficients of the composite material along the x direction, the y direction and the z direction under a global coordinate system are respectively; t is the temperature; q is the heat generation rate; rhocIs the composite density; ccIs a composite material ratioHeat capacity; t is time;
the thermochemical equation shown in the formula (1) is mainly obtained through a Fourier law and a curing kinetic equation of resin, and uneven temperature fields in the composite material easily cause uneven curing of a component in the curing process, so that curing residual stress is caused; the heat generation rate Q in the formula (1) is related to the exothermic heat of the curing reaction of the resin, and is represented by the following formula (2):
where ρ isrIs the resin density; vfIs the fiber volume fraction; hrα is the degree of curing, which represents the degree of progress of the resin curing reaction, and d α/dt is the rate of the resin curing reaction;
the viscosity calculation model of the resin is shown in the following formula (3):
wherein eta is resin viscosity eta∞is viscosity constant, U is resin viscosity flow activation energy, ηβIs a temperature independent fit coefficient; r is an ideal gas constant;
the resin flow model of the composite material is shown in equation (4) below:
wherein k isx、ky、kzPermeability of the composite material along x, y and z directions; prIs the static pressure of the resin within the composite; m isvThe relation between the stress and the strain of the object under the constrained compression condition is described for the volume change coefficient,mvthe calculation formula of (2) is shown as the following formula (5):
wherein, PfThe stress to which the skeletal structure of the fiber network is subjected, commonly referred to as the fiber effective stress; e is the porosity per unit volume of the composite, the available fiber volume fraction VfIs expressed, i.e. e ═ 1-Vf)/VfE and effective stress P of the fiberfThe relationship between them is shown in the following equation (6):
step 2, endowing the macroscopic multi-field coupling calculation model with an initial value (initial temperature T)0Degree of curing a0And fiber volume fraction Vf) And calculating the boundary conditions (temperature and pressure curing process curves) of the model;
step 3, carrying out mesh division on the model established in the step 1 by adopting entity units in finite element software, and carrying out finite element solution on the model, wherein the concrete process is as follows:
calculating the internal temperature T and the curing degree a of the composite material through a thermochemical calculation model, obtaining the internal resin viscosity η of the composite material through the obtained temperature T and curing degree a and a resin viscosity calculation model, and updating the internal fiber volume fraction V of the composite material through a resin flow model by using the obtained temperature T, curing degree a and resin viscosity ηfTemperature T, degree of cure a, fiber volume fraction V obtained by calculationfUpdating model input parameters (density, specific heat capacity, heat conduction coefficient and fiber permeability) and taking the updated model input parameters as initial values of the model at the next moment, wherein the whole cycle process is always carried out until the composite material is cured, and after the curing is finished, corresponding temperature T, curing degree a and fiber volume fraction V are extractedfInputting the data into an Excel file;
step 4, calculating thermodynamic parameters of a constitutive model of the microscopic cured residual stress resin material according to the result obtained in the step 3, wherein the thermodynamic parameters comprise elastic modulus, shear modulus and Poisson ratio;
the specific process is as follows:
reading the Excel file in the step 3, and obtaining the temperature T, the curing degree a and the fiber volume fraction V which are calculated in the step 3fAnd 3, calculating to obtain the elastic modulus E in the resin curing process by taking the temperature T and the curing degree a obtained by calculation in the step 3 as input parameters of a constitutive model of the microscopic curing residual stress resin materialmShear modulus GmAnd poisson's ratio vmThe method comprises the following steps of (1) writing calculation results of formulas (7), (8) and (9) through a subprogram UMAT (unified modeling and planning) of finite element software ABAQUS (equivalent finite element analysis), wherein the subprogram can provide input of a material constitutive model for materials used in a finite element model;
wherein,andrespectively, the resin at the corresponding temperature TC1、TC2The modulus of elasticity of the rubber composition;poisson's ratio, T, after complete curing of the resin∞For the current glass transition temperature T of the resingDifference from the current temperature T of the resin:
T∞＝Tg-T (10)；
Tg＝164.6a2+51a+2.67 (11)；
and 5, establishing an RVE model containing single fibers and a resin matrix by adopting a representative volume unit, inputting the thermodynamic parameters obtained in the step 4 into ABAQUS, taking the temperature T obtained by calculation in the step 3 as a temperature boundary condition of the microscopic curing residual stress calculation model, carrying out grid division on the model in the step 4 by adopting an entity unit in finite element software, submitting analysis operation, editing a General option in a Job module, reading in a resin UMAT subprogram written in the step 4, carrying out finite element calculation on the microscopic curing residual stress model, obtaining the microscopic curing residual stress of the fibers and the resin, and completing the multi-scale simulation of the curing residual stress.
Examples
A composite material member shown in figure 1 is selected, the effective size is 15.26cm multiplied by 3.576cm (length multiplied by width multiplied by height), the material is carbon fiber epoxy resin, the mark is AS4/3501-6, and the initial fiber volume fraction is 40.383%. The layering is unidirectional 0 degree layering along the length direction of the member, and 228 layers are paved together. Because the composite material member and the boundary condition have symmetry, in order to save computer resources, only 1/4 of the composite material model is taken for simulation analysis in the embodiment, as shown in fig. 2;
the specific process of performing the cross-scale simulation of the curing residual stress of the composite member shown in fig. 1 is as follows:
step 1, establishing a composite material macroscopic multi-field coupling calculation model by using Comsol finite element calculation software; the macroscopic multi-field coupling calculation model comprises a thermochemical calculation model, a resin viscosity calculation model and a resin flow calculation model;
the thermochemical calculation model is shown in the following formula (1):
wherein λ isx、λy、λzThe thermal conductivity coefficients of the composite material along the x direction, the y direction and the z direction under a global coordinate system are respectively; t is the temperature; q is the heat generation rate; rhocIs the composite density; ccIs the specific heat capacity of the composite material; t is time;
the thermochemical equation shown in the formula (1) is mainly obtained through a Fourier law and a curing kinetic equation of resin, and uneven temperature fields in the composite material easily cause uneven curing of a component in the curing process, so that curing residual stress is caused; the heat generation rate Q in the formula (1) is related to the exothermic heat of the curing reaction of the resin, and is represented by the following formula (2):
where ρ isrIs the resin density; vfIs the fiber volume fraction; hrα is the degree of curing, which represents the degree of progress of the resin curing reaction, and d α/dt is the rate of the resin curing reaction;
the viscosity calculation model of the resin is shown in the following formula (3):
wherein eta is resin viscosity eta∞is viscosity constant, U is resin viscosity flow activation energy, ηβIs prepared by reacting withA temperature independent fitting coefficient; r is an ideal gas constant;
the resin flow model of the composite material is shown in equation (4) below:
wherein k isx、ky、kzPermeability of the composite material along x, y and z directions; prIs the static pressure of the resin within the composite; m isvFor the volume coefficient of variation, the stress-strain relationship of an object under constrained compression is described, mvThe calculation formula of (2) is shown as the following formula (5):
wherein, PfThe stress to which the skeletal structure of the fiber network is subjected, commonly referred to as the fiber effective stress; e is the porosity per unit volume of the composite, the available fiber volume fraction VfIs expressed, i.e. e ═ 1-Vf)/VfE and effective stress P of the fiberfThe relationship between them is shown in the following equation (6):
step 3, carrying out grid division on the model established in the step 1 by adopting an entity unit in finite element software, carrying out finite element solution on the model, calculating the internal temperature T and the curing degree a of the composite material through a thermochemical calculation model, obtaining the internal resin viscosity η of the composite material by utilizing the resin viscosity calculation model according to the obtained temperature T and the curing degree a, and updating the internal fiber volume fraction V of the composite material through the resin flow model according to the obtained temperature T, the obtained curing degree a and the obtained resin viscosity ηfTemperature T, degree of cure a, fiber volume fraction obtained by calculationVfUpdating model input parameters (density, specific heat capacity, heat conduction coefficient and fiber permeability) and taking the updated model input parameters as initial values of the model at the next moment, wherein the whole cycle process is always carried out until the composite material is cured, and after the curing is finished, corresponding temperature T, curing degree a and fiber volume fraction V are extractedfInputting the data into an Excel file;
the material parameters of the resin and the fiber in the macroscopic multi-field coupling calculation model are shown in the following tables 1 and 2, and the temperature and the curing degree distribution inside the composite material laminated plate can be obtained through the macroscopic multi-field coupling calculation model, as shown in fig. 3 and 4. And (3) calculating the change rule of the temperature, the curing degree and the fiber volume fraction of the point A along with the curing time by using a macroscopic multi-field coupling model, wherein the point A is the central point of the composite material member and has coordinates of (0,0,17.88), as shown in FIGS. 5, 6 and 7.
TABLE 1 Multi-field coupling model Material parameters for resins 3501-6
TABLE 2 Multi-field coupling model Material parameters of fiber AS4
Step 4, reading the Excel file in the step 3, and obtaining the temperature T, the curing degree a and the fiber volume fraction V which are calculated in the step 3fAnd 3, calculating to obtain the elastic modulus E in the resin curing process by taking the temperature T and the curing degree a obtained by calculation in the step 3 as input parameters of a constitutive model of the microscopic curing residual stress resin materialmShear modulus GmAnd poisson's ratio vmAs shown in formulas (7), (8) and (9), thereby realizingThe information transmission from the composite material macroscopic multi-field coupling model to the microscopic residual stress model is realized by writing the calculation results of the formulas (7), (8) and (9) through a subprogram UMAT of finite element software ABAQUS, and the subprogram can provide the input of a material constitutive model for the materials used in the finite element model;
wherein,andrespectively, the resin at the corresponding temperature TC1、TC2The modulus of elasticity of the rubber composition;poisson's ratio, T, after complete curing of the resin∞For the current glass transition temperature T of the resingDifference from the current temperature T of the resin:
T∞＝Tg-T (10)；
Tg＝164.6a2+51a+2.67 (11)；
step 5, using representative volume units to create an RVE model containing individual fibers and a resin matrix, as shown in fig. 8 (1 for fiber and 2 for resin in the figure), the thermodynamic parameters obtained in step 4 were entered in the ABAQUS, as shown in table 3. Taking the temperature T of the point A calculated in the step 3 as a temperature boundary condition of the microscopic curing residual stress calculation model, performing grid division on the model in the step 4 by using an entity unit in finite element software, submitting analysis operation, editing a General option in a Job module, reading in a resin UMAT subprogram written in the step 4, and performing finite element calculation on the microscopic curing residual stress model to obtain the microscopic curing residual stress of the point A, referring to FIG. 9, FIG. 10 shows a microscopic curing residual stress curve of the point A calculated by using the method of the patent and without considering a multi-scale method, and it can be seen that the difference between the prediction results of the point A and the microscopic curing residual stress curve of the point A is very large between 3000s and 9000s when the curing process is performed.
TABLE 3 mechanical Properties parameters of resins 3501-6 and fiber AS4
In summary, for the composite material member with a relatively large thickness, the analysis result of the microscopic residual stress inside the fiber and resin matrix is not accurate because the influence of the macroscopic multi-field coupling effect and the thickness dimension of the laminated plate is not considered in the existing analysis method represented by the comparative example. Therefore, for the microscopic curing residual stress analysis of the large composite material member, the invention has the advantages that the establishment of the correlation between the macroscopic multi-field coupling model and the microscopic curing residual stress model is necessary, and the outstanding improvement effect is achieved.
Claims (3)
1. A multi-scale numerical simulation method for residual stress of composite material curing is characterized by comprising the following steps: the method specifically comprises the following steps:
step 1, establishing a macroscopic multi-field coupling calculation model of the composite material, wherein the macroscopic multi-field coupling calculation model comprises a thermochemical calculation model, a resin viscosity calculation model and a resin flow calculation model;
step 2, giving a macroscopic multi-field coupling calculation model initial value and a boundary condition of the model;
step 3, carrying out mesh division on the model established in the step 1 by adopting entity units in finite element software, and carrying out finite element solution on the model;
step 4, calculating thermodynamic parameters in a resin material constitutive model in the microscopic solidification residual stress model according to the result obtained in the step 3, wherein the thermodynamic parameters comprise elastic modulus, shear modulus and Poisson ratio;
and 5, establishing an RVE model containing single fibers and a resin matrix by adopting a representative volume unit, taking the thermodynamic parameters obtained in the step 4 as input parameters, and performing finite element calculation on a microscopic curing residual stress model to complete multi-scale simulation of the curing residual stress.
2. The multi-scale numerical simulation method for the residual stress of the cured composite material according to claim 1, wherein the method comprises the following steps: the specific process of the step 3 is as follows:
calculating the internal temperature T and the curing degree a of the composite material through the thermochemical calculation model in the step 1, obtaining the internal resin viscosity η of the composite material through the obtained temperature T and the curing degree a and by utilizing a resin viscosity calculation model, and updating the internal fiber volume fraction V of the composite material through a resin flow model by utilizing the obtained temperature T, the obtained curing degree a and the obtained resin viscosity ηfTemperature T, degree of cure a, fiber volume fraction V obtained by calculationfUpdating model input parameters, using the updated model input parameters as initial values of the model at the next moment, completing the curing of the composite material in the whole circulation process, and extracting corresponding temperature T, curing degree a and fiber volume fraction V after the curing is completedf。
3. The multi-scale numerical simulation method of the residual stress of the cured composite material according to claim 2, wherein: the specific process of the step 4 is as follows:
reading the temperature T, the solidification degree a and the fiber volume fraction V which are calculated in the step 3fAnd 3, calculating to obtain the elastic modulus E in the resin curing process by taking the temperature T and the curing degree a obtained by calculation in the step 3 as input parameters of a constitutive model of the microscopic curing residual stress resin materialmShearing and shearingModulus GmAnd poisson's ratio vmAs shown in formulas (7), (8) and (9), thereby realizing information transfer from the composite material macroscopic multi-field coupling model to the microscopic residual stress model;
wherein,andrespectively, the resin at the corresponding temperature TC1、TC2The modulus of elasticity of the rubber composition;poisson's ratio, T, after complete curing of the resin∞For the current glass transition temperature T of the resingDifference from the current temperature T of the resin:
T∞＝Tg-T (10)；
Tg＝164.6a2+51a+2.67 (11)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810448171.9A CN108664731A (en) | 2018-05-11 | 2018-05-11 | A kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810448171.9A CN108664731A (en) | 2018-05-11 | 2018-05-11 | A kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108664731A true CN108664731A (en) | 2018-10-16 |
Family
ID=63779066
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810448171.9A Pending CN108664731A (en) | 2018-05-11 | 2018-05-11 | A kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108664731A (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109857061A (en) * | 2019-01-24 | 2019-06-07 | 贵州大学 | A kind of workpiece surface residual stress regulation method based on thermal influence zone |
| CN111274667A (en) * | 2019-12-13 | 2020-06-12 | 成都材智科技有限公司 | Cross-scale material calculation software integrated calculation system and method |
| WO2020237977A1 (en) * | 2019-05-27 | 2020-12-03 | 北京工业大学 | Multi-scale simulation method for mechanical behavior of multi-phase composite material |
| CN112329297A (en) * | 2020-10-30 | 2021-02-05 | 西北工业大学 | Composite material strength detection method and system considering curing residual stress |
| CN112417547A (en) * | 2020-10-21 | 2021-02-26 | 山东大学 | Multi-scale and multi-field evolution simulation method and system for coastal sand layer grouting reinforcement strength |
| CN112906274A (en) * | 2021-02-22 | 2021-06-04 | 中国核动力研究设计院 | Visual interface and method for annealing simulation of cladding material |
| CN113033051A (en) * | 2021-03-24 | 2021-06-25 | 天津中德应用技术大学 | Compression molding process optimization method for carbon fiber reinforced composite material sheet |
| CN113221408A (en) * | 2021-05-06 | 2021-08-06 | 浙江工业大学 | Parallel NTFA theory model order reduction calculation method based on ABAQUS |
| CN113607772A (en) * | 2021-08-04 | 2021-11-05 | 西北工业大学 | Method and system for determining damage of toughened composite material during curing molding |
| CN113976920A (en) * | 2021-09-27 | 2022-01-28 | 上海交通大学 | Cross-scale control method and system for residual deformation of selective laser melting forming structure |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102229022A (en) * | 2011-06-21 | 2011-11-02 | 西安理工大学 | Ultrasonically-assisted metastable alloy rapidly-solidifying welding device and method |
| CN103473440A (en) * | 2013-08-20 | 2013-12-25 | 北京航空航天大学 | New method for predicting residual thermal stress of planar-knitted composite material |
| CN103488850A (en) * | 2013-10-15 | 2014-01-01 | 北京航空航天大学 | Method for positioning and homogenizing residual stress of thermal shock compounding |
| CN107283871A (en) * | 2017-08-24 | 2017-10-24 | 哈尔滨工业大学 | A kind of preparation method of thermoplastic resin matrix's carbon fiber titanium/titanium alloy layer plywood |
-
2018
- 2018-05-11 CN CN201810448171.9A patent/CN108664731A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102229022A (en) * | 2011-06-21 | 2011-11-02 | 西安理工大学 | Ultrasonically-assisted metastable alloy rapidly-solidifying welding device and method |
| CN103473440A (en) * | 2013-08-20 | 2013-12-25 | 北京航空航天大学 | New method for predicting residual thermal stress of planar-knitted composite material |
| CN103488850A (en) * | 2013-10-15 | 2014-01-01 | 北京航空航天大学 | Method for positioning and homogenizing residual stress of thermal shock compounding |
| CN107283871A (en) * | 2017-08-24 | 2017-10-24 | 哈尔滨工业大学 | A kind of preparation method of thermoplastic resin matrix's carbon fiber titanium/titanium alloy layer plywood |
Non-Patent Citations (5)
| Title |
|---|
| ZHENYI YUAN .ETC: ""An analytical model on through-thickness stresses and warpage of composite laminates due to toolepart interaction"", 《COMPOSITES PART B》 * |
| 元振毅等: ""热固性树脂基复合材料热隔膜成型过程数值仿真"", 《复合材料学报》 * |
| 张纪奎等: ""固化度与固化收缩对非对称复合材料层合板固化变形的影响"", 《复合材料学报》 * |
| 杨雷等: ""考虑纤维随机分布的复合材料热残余应力分析及其对横向力学性能的影响"", 《复合材料学报》 * |
| 闵荣等: ""基于黏弹性本构模型的热固性树脂基复合材料固化变形数值仿真模型"", 《复合材料学报》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109857061B (en) * | 2019-01-24 | 2022-05-17 | 贵州大学 | Workpiece surface residual stress regulation and control method based on thermal coupling model |
| CN109857061A (en) * | 2019-01-24 | 2019-06-07 | 贵州大学 | A kind of workpiece surface residual stress regulation method based on thermal influence zone |
| WO2020237977A1 (en) * | 2019-05-27 | 2020-12-03 | 北京工业大学 | Multi-scale simulation method for mechanical behavior of multi-phase composite material |
| US11798658B2 (en) | 2019-05-27 | 2023-10-24 | Beijing Universitry Of Technology | Multi-scale method for simulating mechanical behaviors of multiphase composite materials |
| CN111274667A (en) * | 2019-12-13 | 2020-06-12 | 成都材智科技有限公司 | Cross-scale material calculation software integrated calculation system and method |
| CN111274667B (en) * | 2019-12-13 | 2023-06-02 | 成都材智科技有限公司 | Cross-scale material computing software integrated computing system and method |
| CN112417547A (en) * | 2020-10-21 | 2021-02-26 | 山东大学 | Multi-scale and multi-field evolution simulation method and system for coastal sand layer grouting reinforcement strength |
| CN112329297A (en) * | 2020-10-30 | 2021-02-05 | 西北工业大学 | Composite material strength detection method and system considering curing residual stress |
| CN112906274A (en) * | 2021-02-22 | 2021-06-04 | 中国核动力研究设计院 | Visual interface and method for annealing simulation of cladding material |
| CN113033051A (en) * | 2021-03-24 | 2021-06-25 | 天津中德应用技术大学 | Compression molding process optimization method for carbon fiber reinforced composite material sheet |
| CN113033051B (en) * | 2021-03-24 | 2023-09-12 | 天津中德应用技术大学 | Compression molding process optimization method for carbon fiber reinforced composite material sheet |
| CN113221408A (en) * | 2021-05-06 | 2021-08-06 | 浙江工业大学 | Parallel NTFA theory model order reduction calculation method based on ABAQUS |
| CN113607772A (en) * | 2021-08-04 | 2021-11-05 | 西北工业大学 | Method and system for determining damage of toughened composite material during curing molding |
| CN113976920A (en) * | 2021-09-27 | 2022-01-28 | 上海交通大学 | Cross-scale control method and system for residual deformation of selective laser melting forming structure |
| CN113976920B (en) * | 2021-09-27 | 2022-08-26 | 上海交通大学 | Cross-scale control method and system for residual deformation of selective laser melting forming structure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108664731A (en) | A kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller | |
| CN106066913B (en) | Complex composite material structure equivalent material performance multi-dimension calculation method | |
| CN109241650B (en) | Carbon fiber reinforced composite material mechanical property prediction method based on cross-scale simulation | |
| CN107451307B (en) | Method for multi-scale calculation of equivalent stiffness matrix of complex composite structure | |
| CN107451308B (en) | Multi-scale calculation method for equivalent heat conduction coefficient of complex composite material structure | |
| JP6756997B2 (en) | A recording medium on which a finite element method analysis method, a finite element method analysis device, an analysis service system, and a finite element method analysis program are recorded. | |
| CN107451309B (en) | Method for multi-scale calculation of equivalent thermal expansion coefficient of complex composite material structure | |
| Altenbach | Theories for laminated and sandwich plates: A review | |
| CN106626449B (en) | Consider the design method in the composite material V-type component autoclave molding tooling type face of curing deformation | |
| CN107025366B (en) | Interactive coupling numerical simulation method for composite material autoclave molding temperature field | |
| KR101383663B1 (en) | Parametrized Material and performance properties based on virtual testing | |
| Reddy et al. | Bending analysis of laminated composite plates using finite element method | |
| CN108108578A (en) | The numerical algorithm of FG-GRC bucking of plate load factors based on gridless routing | |
| CN113033051B (en) | Compression molding process optimization method for carbon fiber reinforced composite material sheet | |
| Xing et al. | An eigenelement method of periodical composite structures | |
| CN112329297A (en) | Composite material strength detection method and system considering curing residual stress | |
| Shah et al. | Flow modeling and simulation study of vacuum assisted resin transfer molding (VARTM) process: A review | |
| Zhang et al. | Simulation of failure behavior of 2.5 D SiC/SiC variable thickness dovetail joint structures based on mesoscale model | |
| Liang et al. | Hydrodynamic performance optimization of marine propellers based on fluid-structure coupling | |
| CN112926244A (en) | Method for determining ultimate load of perforated part of composite laminated plate | |
| Amini Niaki | A three-phase integrated flow-stress framework for process modelling of composite materials | |
| Li et al. | Numerical simulation of flow and compaction during the cure of laminated composites | |
| CN111259597B (en) | Multi-scale coupling curing analysis method for resin-based composite material | |
| Xie et al. | Mechanism of curing behavior for CFRP compression molding under thermo-mechanical-chemical multi-field coupling | |
| Plaisant et al. | Finite element procedure for stress amplification factor recovering in a representative volume of composite materials |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181016 |