CN111259597A - Multi-scale coupling curing analysis method for resin-based composite material - Google Patents
Multi-scale coupling curing analysis method for resin-based composite material Download PDFInfo
- Publication number
- CN111259597A CN111259597A CN202010148486.9A CN202010148486A CN111259597A CN 111259597 A CN111259597 A CN 111259597A CN 202010148486 A CN202010148486 A CN 202010148486A CN 111259597 A CN111259597 A CN 111259597A
- Authority
- CN
- China
- Prior art keywords
- resin
- curing
- composite material
- model
- based composite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
The invention discloses a multi-scale coupling solidification analysis method for a resin-based composite material, which comprises the following steps: establishing a microscopic analysis model of the resin-based composite material according to the periodic microstructure characteristics of the resin-based composite material; applying periodic boundary conditions to the micro-analysis model according to a homogenization theory, applying a resin curing dynamic model to obtain the curing characteristics of the micro-analysis model, and applying a structural mechanical model of a resin-based composite material to obtain the mechanical characteristics of the micro-analysis model of the composite material; establishing a macroscopic model of the resin-based composite material, and endowing the macroscopic model with curing characteristics and mechanical characteristics as attributes; and (3) applying boundary conditions to the macroscopic model of the resin-based composite material, and obtaining a curing analysis result by using a finite element analysis method. The invention can consider the curing characteristic of the resin-based composite material from a microscopic structure aiming at the macroscopic structure of the resin-based composite material, and essentially grasp the curing mechanism of the resin-based composite material, thereby improving the modeling and calculating efficiency while ensuring the accuracy.
Description
Technical Field
The invention is applied to the field of simulation of composite materials, particularly the field of curing simulation of resin-based composite materials, and particularly relates to a multi-scale coupling curing analysis method for a resin-based composite material.
Background
The composite material is obtained by mixing two or more materials by a physical or chemical mode, and comprises a matrix phase and a reinforcing phase. The forming and curing of resin-based composite materials is a very critical link in the production and manufacture of composite materials. During the curing process, the temperature field and curing degree field inside the composite material are not uniformly distributed, so that the composite material can undergo thermal expansion and curing shrinkage to different degrees, thereby causing thermal stress and curing shrinkage stress. The existence of these internal stresses has a great influence on the mechanical properties of the composite material, and particularly for thick plate members, the internal stresses may cause bending, cracking and delamination of the composite material plate, and even damage to the composite material during curing.
The current composite material molding simulation is based on a composite material macro model, and then the temperature change, the curing residual stress, the curing deformation and the like in the composite material curing deformation are obtained through resin reaction dynamics and composite material structure mechanics. However, the main problems of the above method are that the microscopic characteristics of the composite material are ignored, the essence of the curing process of the composite material is not known, corresponding phenomenological theory needs to be developed for different materials, and the application range and accuracy of the method are limited.
Disclosure of Invention
The invention provides a method for carrying out multi-scale analysis on microstructure characteristics while carrying out macroscopic solidification analysis on a composite material, aiming at the problem of structural response analysis of the composite material with periodic microstructure characteristics in the solidification process, so that the accuracy of the solidification analysis of the composite material can be greatly improved.
The invention solves the technical problems through the following technical scheme:
the invention provides a multi-scale coupling curing analysis method for a resin-based composite material, which is characterized by comprising the following steps of:
s1, establishing a microscopic analysis model of the resin-based composite material according to the periodic microstructure characteristics of the resin-based composite material;
s2, applying periodic boundary conditions to the micro-analysis model according to the homogenization theory, applying a resin curing dynamic model to obtain the curing characteristics of the micro-analysis model, and applying a structural mechanical model of the resin-based composite material to obtain the mechanical characteristics of the micro-analysis model of the composite material;
s3, establishing a macroscopic model of the resin-based composite material, and endowing the macroscopic model with curing characteristics and mechanical characteristics as attributes;
and S4, applying boundary conditions to the macroscopic model, and obtaining a solidification analysis result by using a finite element analysis method.
Preferably, in step S1, representative volume units for characterizing the periodic microstructure of the resin-based composite material are selected, a corresponding finite element analysis model is established according to the selected representative volume units, and a periodic meshing method is used for meshing to obtain a microscopic finite element analysis model.
Preferably, in step S2, a resin curing dynamic model is applied to obtain a curing degree and a curing shrinkage strain in the representative volume unit, and a structural mechanical model is used to obtain mechanical properties of the composite material micro-analysis model based on the curing degree and the curing shrinkage strain;
the resin cure kinetic model is as follows:
α is the degree of cure of the resin, T is the cure time, T is the resin temperature, AiIs a pre-exponential factor, Δ EiR is a gas constant, i is 1, 2 or 3, and m, l and n are material constants;
the resin cure shrinkage model is as follows:
the curing shrinkage strain of the resulting resin is:
preferably, the resin-based composite material includes, but is not limited to, thermoset composite materials and thermoplastic composite materials.
Preferably, the curing analysis results include, but are not limited to, curing process temperature, curing residual stress, curing deformation of the resin-based composite material.
On the basis of the common knowledge in the field, the above preferred conditions can be combined randomly to obtain the preferred embodiments of the invention.
The positive progress effects of the invention are as follows:
the invention can consider the curing characteristic of the resin-based composite material from a microscopic structure aiming at the macroscopic structure of the resin-based composite material, essentially grasp the curing mechanism of the resin-based composite material, and avoid establishing a microscopic analysis model aiming at the macroscopic structure of the whole resin-based composite material, thereby improving the modeling and calculating efficiency while ensuring the accuracy.
Drawings
FIG. 1 is a flow chart of a multi-scale coupling solidification analysis method for a resin-based composite material according to a preferred embodiment of the invention.
FIG. 2 is a schematic view of a representative volume unit of a resin-based composite material according to a preferred embodiment of the present invention.
FIG. 3 is a schematic view of a representative volume unit of a long fiber-reinforced resin-based composite material according to a preferred embodiment of the present invention.
Fig. 4 is a schematic representation of a representative volume unit of a resin-based composite material reinforced with long fibers after meshing according to a preferred embodiment of the present invention.
FIG. 5 is a schematic view of a microscopic analysis model of a long fiber-reinforced resin-based composite material with periodic boundary conditions applied in accordance with a preferred embodiment of the present invention.
FIG. 6 is a deformation diagram obtained by the curing analysis of the long fiber-reinforced resin-based composite material according to the preferred embodiment of the present invention.
Fig. 7 is a partially enlarged view of the right end portion of fig. 6.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
As shown in fig. 1, the present embodiment provides a multi-scale coupling solidification analysis method for a resin-based composite material, including the following steps:
Resin-based composites include, but are not limited to, thermoset composites and thermoplastic composites.
In step 101, a representative volume unit (see fig. 2) for characterizing periodic microstructure characteristics of the resin-based composite material is selected, a corresponding finite element analysis model is established according to the selected representative volume unit, and a periodic meshing method is adopted for meshing to obtain a microscopic finite element analysis model.
For example: taking a long fiber reinforced resin-based composite material as an example, establishing a microstructure model of the long fiber reinforced resin-based composite material according to the periodic microstructure characteristics of the long fiber reinforced resin-based composite material, selecting a representative volume unit for representing the periodic microstructure characteristics of the resin-based composite material, wherein a dark color represents fibers and a light color represents resin (see fig. 3), establishing a corresponding finite element analysis model according to the selected representative volume unit, and mapping and gridding the finite element analysis model, as shown in fig. 4.
And 102, applying periodic boundary conditions to the micro-analysis model according to a homogenization theory, obtaining the curing characteristics of the micro-analysis model by applying a resin curing dynamic model, and obtaining the mechanical characteristics of the micro-analysis model of the composite material by applying a structural mechanical model of the resin-based composite material based on the curing characteristics.
In step 102, a resin curing dynamic model is applied to obtain a curing degree and a curing shrinkage strain in a representative volume unit, and a structural mechanical model is used to obtain mechanical properties of a composite material micro-analysis model based on the curing degree and the curing shrinkage strain;
the resin cure kinetic model is as follows:
α is the degree of cure of the resin, T is the cure time, T is the resin temperature, AiIs a pre-exponential factor, Δ EiR is a gas constant, i is 1, 2 or 3, and m, l and n are material constants;
the resin cure shrinkage model is as follows:
the curing shrinkage strain of the resulting resin is:
improved curing kinetics model due to the addition of K2αm+K3αlThe model considers more factors, so that the curing degree of the composite material in the curing process can be better simulated and calculated, and the simulated and calculated curing degree is more accurate.
The improved curing shrinkage model has the advantages that the fine tuning variable E is added, so that the error is reduced, the curing shrinkage strain of the composite material in the curing process can be better simulated and calculated, and the simulated and calculated curing shrinkage strain is more accurate.
The improved model is used for calculation, so that the curing degree and curing shrinkage strain of the composite material in the curing process can be better obtained, the more accurate mechanical characteristics of the composite material micro-analysis model are further obtained, and the accuracy of the composite material curing analysis result is improved.
For example: and (3) applying periodic boundary conditions to a microscopic analysis model of the long fiber reinforced resin matrix composite material according to a homogenization theory (see figure 5), applying a resin curing dynamic model to obtain the curing characteristics of the microscopic model, and applying a structural mechanical model of the resin matrix composite material to obtain the mechanical characteristics of the composite material microscopic model.
And 103, establishing a macroscopic model of the resin-based composite material, and endowing the macroscopic model with curing characteristics and mechanical characteristics as attributes.
And 104, applying boundary conditions to the macroscopic model, and obtaining a solidification analysis result by using a finite element analysis method.
Wherein, the curing analysis result includes but is not limited to curing process temperature, curing residual stress and curing deformation of the resin-based composite material.
For example: as shown in fig. 6, a deformation diagram obtained by a curing analysis of the long-fiber-reinforced resin-based composite material, and fig. 7 is a partially enlarged view of the right end portion of the deformation diagram.
While specific embodiments of the invention have been described above, it will be appreciated by those skilled in the art that these are by way of example only, and that the scope of the invention is defined by the appended claims. Various changes and modifications to these embodiments may be made by those skilled in the art without departing from the spirit and scope of the invention, and these changes and modifications are within the scope of the invention.
Claims (5)
1. A multi-scale coupling solidification analysis method for a resin-based composite material is characterized by comprising the following steps:
s1, establishing a microscopic analysis model of the resin-based composite material according to the periodic microstructure characteristics of the resin-based composite material;
s2, applying periodic boundary conditions to the microscopic analysis model according to the homogenization theory, applying a resin curing dynamic model to obtain the curing characteristics of the microscopic analysis model, and applying a structural mechanical model of a resin-based composite material to obtain the mechanical characteristics of the composite material microscopic analysis model based on the curing characteristics;
s3, establishing a macroscopic model of the resin-based composite material, and endowing the macroscopic model with curing characteristics and mechanical characteristics as attributes;
and S4, applying boundary conditions to the macroscopic model, and obtaining a solidification analysis result by using a finite element analysis method.
2. The multi-scale coupling solidification analysis method for resin-based composite materials according to claim 1, wherein in step S1, representative volume units for characterizing periodic microstructure characteristics of the resin-based composite materials are selected, corresponding finite element analysis models are established according to the selected representative volume units, and a periodic meshing method is adopted for meshing to obtain a microscopic finite element analysis model.
3. The multi-scale coupling curing analysis method for the resin-based composite material according to claim 2, wherein in step S2, the resin curing dynamic model is applied to obtain the curing degree and the curing shrinkage strain in the representative volume unit, and the mechanical properties of the composite material micro-analysis model are obtained by using the structural mechanical model based on the curing degree and the curing shrinkage strain;
the resin cure kinetic model is as follows:
α is the degree of cure of the resin, T is the cure time, T is the resin temperature, AiIs a pre-exponential factor, Δ EiR is a gas constant, i is 1, 2 or 3, and m, l and n are material constants;
the resin cure shrinkage model is as follows:
the curing shrinkage strain of the resulting resin is:
4. the resin-based composite multi-scale coupling solidification analysis method of claim 1, wherein the resin-based composite includes but is not limited to thermosetting composite and thermoplastic composite.
5. The resin-based composite multi-scale coupling curing analysis method according to claim 1, wherein the curing analysis result includes but is not limited to curing process temperature, curing residual stress and curing deformation of the resin-based composite.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010148486.9A CN111259597B (en) | 2020-03-05 | 2020-03-05 | Multi-scale coupling curing analysis method for resin-based composite material |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010148486.9A CN111259597B (en) | 2020-03-05 | 2020-03-05 | Multi-scale coupling curing analysis method for resin-based composite material |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111259597A true CN111259597A (en) | 2020-06-09 |
| CN111259597B CN111259597B (en) | 2021-01-26 |
Family
ID=70949854
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010148486.9A Active CN111259597B (en) | 2020-03-05 | 2020-03-05 | Multi-scale coupling curing analysis method for resin-based composite material |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111259597B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113761763A (en) * | 2021-08-06 | 2021-12-07 | 上海索辰信息科技股份有限公司 | RVE multi-scale macro material property analysis method for micro and micro scale structure |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5662993A (en) * | 1995-09-08 | 1997-09-02 | General Motors Corporation | Carbon-based friction material for automotive continuous slip service |
| CN105183990A (en) * | 2015-09-07 | 2015-12-23 | 哈尔滨工业大学 | Composite material and method for multiscale response analysis of structure thereof |
| CN107451309A (en) * | 2016-05-31 | 2017-12-08 | 西北工业大学 | A kind of method of Multi-Scale Calculation complex composite material structure fiber yarn |
| CN108268691A (en) * | 2017-11-21 | 2018-07-10 | 北京航空航天大学 | The analog analysing method of automatic fiber placement composite materials property |
| CN109241650A (en) * | 2018-09-25 | 2019-01-18 | 南京航空航天大学 | Based on the Long-time Mechanical Properties of Carbon Fiber Reinforced Polymer prediction technique emulated across scale |
-
2020
- 2020-03-05 CN CN202010148486.9A patent/CN111259597B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5662993A (en) * | 1995-09-08 | 1997-09-02 | General Motors Corporation | Carbon-based friction material for automotive continuous slip service |
| CN105183990A (en) * | 2015-09-07 | 2015-12-23 | 哈尔滨工业大学 | Composite material and method for multiscale response analysis of structure thereof |
| CN107451309A (en) * | 2016-05-31 | 2017-12-08 | 西北工业大学 | A kind of method of Multi-Scale Calculation complex composite material structure fiber yarn |
| CN108268691A (en) * | 2017-11-21 | 2018-07-10 | 北京航空航天大学 | The analog analysing method of automatic fiber placement composite materials property |
| CN109241650A (en) * | 2018-09-25 | 2019-01-18 | 南京航空航天大学 | Based on the Long-time Mechanical Properties of Carbon Fiber Reinforced Polymer prediction technique emulated across scale |
Non-Patent Citations (2)
| Title |
|---|
| ZHENYI YUAN 等: "Evolution of curing residual stresses in composite using multi-scale method", 《COMPOSITES PART B》 * |
| 毕凤阳 等: "树脂基复合材料热压罐成型工艺固化变形模拟", 《兵器材料科学与工程》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113761763A (en) * | 2021-08-06 | 2021-12-07 | 上海索辰信息科技股份有限公司 | RVE multi-scale macro material property analysis method for micro and micro scale structure |
| CN113761763B (en) * | 2021-08-06 | 2023-05-30 | 上海索辰信息科技股份有限公司 | Method for analyzing properties of RVE multi-scale macroscopic materials of microscopic and microscale structures |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111259597B (en) | 2021-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gopal et al. | Optimal temperature profiles for minimum residual stress in the cure process of polymer composites | |
| Yu et al. | Deep learning-based heterogeneous strategy for customizing responses of lattice structures | |
| CN110162827A (en) | A kind of entity finite element calculation method of concrete structure time-varying effecting | |
| Alazwari et al. | Modeling and analysis of composite laminates in the presence of uncertainties | |
| CN108664731A (en) | A kind of multiple dimensioned method for numerical simulation of composite material Googol motion controller | |
| Ribeiro et al. | 3D finite element model for hybrid FRP-confined concrete in compression using modified CDPM | |
| Huang et al. | Micromechanical prediction of elastic-plastic behavior of a short fiber or particle reinforced composite | |
| CN106626449A (en) | Design method for composite material V-shaped component autoclave forming tool molded surface considering curing deformation | |
| CN107742005A (en) | A kind of fiber-reinforced composite materials structures mechanical properties prediction and control method | |
| Feng et al. | Non-uniform fiber-resin distributions of pultruded GFRP profiles | |
| Jain et al. | Mean field homogenization methods for strand composites | |
| CN111259597B (en) | Multi-scale coupling curing analysis method for resin-based composite material | |
| CN112329297A (en) | Composite material strength detection method and system considering curing residual stress | |
| CN109117573A (en) | A kind of nonlinear numerical simulation method of consideration composite material parameter three dimensional anisotropic | |
| Stabla et al. | The effect of mosaic pattern and winding angle on radially compressed filament-wound CFRP composite tubes | |
| Trofimov et al. | A review on the Representative Volume Element-based multi-scale simulation of 3D woven high performance thermoset composites manufactured using resin transfer molding process | |
| Shokrieh et al. | Modeling residual stresses in composite materials | |
| Che et al. | Warping deformation analysis of angle-ply fiber metal laminates with slippage during manufacturing | |
| CN114239157A (en) | Method for designing layering angle of composite material hydrogen storage container based on machine learning | |
| Zhang et al. | An integrated computational materials engineering method for woven carbon fiber composites preforming process | |
| Fedulov et al. | Residual stresses near the free edge of composite materials | |
| Liu et al. | Warpage optimization of the GFRP injection molding process parameters | |
| Sun et al. | Bi‐stable cured shape of CFRP variant structures with asymmetric layup: Rapid prediction and inverse design | |
| CN109501325B (en) | Method and device for predicting curing deformation of composite material member | |
| Langella et al. | Step pultrusion |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 201206 Shanghai, Pudong New Area, China (Shanghai) free trade zone, new Jinqiao Road, No. 13, building 2, floor 27 Applicant after: Shanghai suochen Information Technology Co., Ltd Address before: 201204 building 27, Lane 676, Wuxing Road, Pudong New Area, Shanghai Applicant before: SHANGHAI SUOCHEN INFORMATION TECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |