CN110532702A - A kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method - Google Patents
A kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method Download PDFInfo
- Publication number
- CN110532702A CN110532702A CN201910821548.5A CN201910821548A CN110532702A CN 110532702 A CN110532702 A CN 110532702A CN 201910821548 A CN201910821548 A CN 201910821548A CN 110532702 A CN110532702 A CN 110532702A
- Authority
- CN
- China
- Prior art keywords
- fiber
- interface
- mechanical performance
- composite material
- volume
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Moulding By Coating Moulds (AREA)
Abstract
A kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method disclosed by the invention, belongs to Materials Science and Engineering application field.Implementation method are as follows: the three-dimensional representation elementary volume, volume element RVE geometrical model of fiber random distribution is established based on random sequence additive process, carry out FEM meshing and to fiber element, base unit and boundary element assign respective material parameter, by the deterioration of modulus of debonding interface unit with modeling interface debonding defect, applying periodic boundary condition and required displacement load containing interfacial detachment defect finite element model boundary, composite material Equivalent Mechanical performance is sought based on the method for homogenizing, the Equivalent Mechanical performance prediction under different interfacial detachment defect conditions is realized by changing interfacial detachment ratio, it compares the Equivalent Mechanical performance prediction result containing the random unsticking in interface of engineer application and prepares the mechanical property measured result of finished product, improve long fiber reinforcement composite material preparation process, sufficiently and safely play composite material bearing structure power Learn performance.
Description
Technical field
There are the long fiber reinforcements of random debonding defect at the interface that the present invention relates to a kind of between fiber and resin matrix
The prediction technique of type composite material Equivalent Mechanical performance, belongs to Materials Science and Engineering application field.
Background technique
Long-fiber-reinforced resin composite material is since the high the features such as comparing strong, Gao Bigang and low-density of itself, is by increasingly
Mostly it is applied to the various fields such as auto industry, sport and aerospace.Long-fiber-reinforced resin composite wood is a kind of design
Strong composite material, so that fiber and resin have complementary advantages, it is insufficient to improve one-component material property by rationally designing
Disadvantage, so that synergistic effect is generated between the two, so that composite material overall performance be made to meet engineering better than its composition material
In different demands.
The unsticking process that the current research for long-fiber-reinforced resin interface is concentrated mainly on fiber/resin interface is ground
Study carefully, especially the gradually unsticking failure procedure of the carbon fiber/epoxy resin interface in the case where bearing external load function, and ignores compound
The material influence of the existing random debonding defect in fiber/resin interface to composite material Equivalent Mechanical performance after preparation molding is ground
Study carefully.Since most pulp freeness are small, surface-active is low, wetability and poor with the two-phase binding performance of resin, even if
It is designed and manufactured under stringent quality control conditions, is formed by fiber/resin interface in this course and also unavoidably deposits
In a certain proportion of random unsticking, this seriously affects this kind of composite inner stress transfer and deformation under external load function, leads
It causes its true mechanical property to be far below the desired value at design initial stage, and then may result in bearing structure premature failure.In addition, fine
The factors such as the fiber random distribution under random unsticking position and different carbon fibrous body fractions at dimension/resin boundary surface also increase
The complexity that the research random unsticking in fiber/resin interface influences corresponding composite material Equivalent Mechanical performance.Therefore it is studying
When the Equivalent Mechanical performance of long fiber reinforcement composite material, the random debonding defect at fiber/resin interface is taken into account can be with
Load-carrying properties assessment accuracy when structure design is greatly improved, both can be unable to give full play material to avoid overly conservative
Can, and safety accident can be led to avoid excessive use.Based on this, the random debonding defect in fiber/resin interface is to multiple
The influence research of condensation material Equivalent Mechanical performance just has very major and immediate significance in engineer application.
Summary of the invention
For fibre/random debonding defect in base interface deficiency is not considered in existing composite material Equivalent Mechanical performance prediction,
A kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method technical problems to be solved disclosed by the invention are: being based on
The composite material different component method that homogenizes using finite element method establishes three-dimensional representation elementary volume, volume element model, by interfacial detachment
Defective proportion reflection consider fibre/random debonding defect in base interface long fiber reinforcement composite inner there are fibre/base interface with
Machine debonding defect situation is degenerated by the modulus to unsticking unit collection with the fixed debonding defect ratio in simulated implementation interface
In the case of Equivalent Mechanical performance prediction, pass through change interfacial detachment defective proportion realize interface difference debonding defect ratio condition
Under Equivalent Mechanical performance prediction；Establishing interfacial detachment defective proportion according to prediction result influences relationship song to Equivalent Mechanical performance
Line, that is, the random unsticking for establishing interface influences relation curve to Equivalent Mechanical performance, on this basis, by the practical work of engineer application
The Equivalent Mechanical performance prediction result of the random unsticking in interface and the actual measurement mechanical property for preparing finished product compare under condition, further
Guidance and improvement long fiber reinforcement composite material prepare end properties, and then solve long fiber reinforcement composite material application field
Correlation engineering problem.
The purpose of the present invention is what is be achieved through the following technical solutions.
A kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method disclosed by the invention, is added based on random sequence
Addition establishes the three-dimensional representation elementary volume, volume element RVE geometrical model of fiber random distribution, carries out FEM meshing and to fiber
Unit, base unit and boundary element assign respective material parameter, and wherein the deterioration of modulus of debonding interface unit is with modeling interface
Debonding defect is applying periodic boundary condition and required displacement load, base containing interfacial detachment defect finite element model boundary
Composite material Equivalent Mechanical performance is sought in the method for homogenizing, realizes different interfacial detachment defect items by changing interfacial detachment ratio
Equivalent Mechanical performance prediction under part, compare engineer application the Equivalent Mechanical performance prediction result containing the random unsticking in interface and
The mechanical property measured result of finished product is prepared, and then instructs and improve long fiber reinforcement composite material preparation process, sufficiently and is pacified
Composite material bearing structure mechanical property is played entirely.
A kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method disclosed by the invention, includes the following steps:
Step 1 establishes the three-dimensional representation elementary volume, volume element RVE geometry of fiber random distribution based on random sequence adding method
Model.
Preferably, establishing fiber random distribution using Python based on random sequence adding method in step 1
Three-dimensional representation elementary volume, volume element RVE geometrical model.
Step 1.1: three-dimensional representation volume list is determined according to the calculating of long fibre radius, fiber volume fraction and fibre number
First RVE geometrical model size, and the three-dimensional cubic structure according to identified size generation resin matrix.
Step 1.2: random selection is a little used as the center of circle, and binding fiber diameter in the one side of the three-dimensional cube of generation
Determine first position of fibers.
Step 1.3: on first fiber boundary of generation, being enclosed in first fiber in conjunction with determining interfacial thickness generation
Fibre/base interface of surrounding.
Step 1.4: in the case where having generated first fiber, its cofibre and correspondence are gradually generated according to similar approach
Fibre/base interface, the subsequent fiber gradually generated need to be met and existing fiber is non-intersecting, until the fiber volume for reaching required contains
Until amount, that is, realize the three-dimensional representation elementary volume, volume element RVE geometry mould that fiber random distribution is established based on random sequence adding method
Type.
Step 2, the three-dimensional representation elementary volume, volume element RVE model established to step 1 carries out FEM meshing, in institute
Node is extracted at the interface location of some fiber peripherals, and interface node set and consistent with its is generated according to extracted node
Node serial number list；Boundary element set and element number list are generated using extracted interface node.
Step 2 implementation method is as follows:
Step 2.1: grid dividing is carried out to three-dimensional representation elementary volume, volume element RVE model.
Step 2.2: extracting node at the interface location of each fiber peripheral, generate interface node set and and one
The node serial number list of cause.
Step 2.3: utilizing extracted interface node, generate boundary element set and arranged with its consistent element number
Table.
Step 3: extracting all fibres node in finite element model, generates fiber element collection using extracted fiber nodes
Close and with its consistent element number list, assign fibrous material parameter to the unit of the fiber element set of generation；Extract institute
There is matrix node, matrix node set and consistent element number list therewith is generated using extracted matrix node, to generation
Base unit set unit assign basis material parameter；Unit concentration in generated interface randomly selects interface list
The boundary element modulus randomly selected degenerate as random debonding defect existing for interface by member, and remaining boundary element is made
Corresponding non-debonding interface material parameter, i.e. the assignment processing of realization each component material parameter are assigned for non-debonding interface unit.
Each component material parameter described in step 3 refers to fibrous material parameter in step 3, basis material parameter, debonding interface
Material parameter and non-debonding interface material parameter.
Unit concentration in step 3 in generated interface randomly selects boundary element, preferably, having generated
Interface unit focus utilization Python in the generation of random series method based on Mason's Rotation Algorithm randomly select boundary
Face unit is as unsticking unit.
Step 4: in the generated fibre containing fixed proportion/random unsticking in base interface three-dimensional representation elementary volume, volume element RVE
Finite element model boundary applies periodic boundary condition and required displacement load.
Step 5: based on the three-dimensional representation elementary volume, volume element RVE finite element model in step 4, by interfacial detachment defective proportion
Reflection consider the long fiber reinforcement composite inner of the random debonding defect in interface there are the random debonding defect situation in interface, according to
The each component material parameter of step 3 assignment processing solves composite material Equivalent Mechanical performance based on Homogenization, that is, realizes
Equivalent Mechanical performance prediction under the fixed debonding defect ratio situation in interface.
Step 6: the Equivalent Mechanical performance under the fixed debonding defect ratio situation in the interface based on step 1 to step 5 is pre-
Survey method, by changing interfacial detachment defective proportion to realize the Equivalent Mechanical performance under the conditions of different interfacial detachment defective proportions
Prediction；According to prediction result establish interfacial detachment defective proportion on Equivalent Mechanical performance influence relation curve, that is, establish interface with
Machine unsticking influences relation curve to Equivalent Mechanical performance.
Step 7: being realized based on the random unsticking in interface that step 6 is established influences relation curve, root to Equivalent Mechanical performance
The Equivalent Mechanical performance prediction under the conditions of actual condition interfacial detachment defective proportion is realized according to step 1 to step 6, and engineering is answered
It is carried out pair with the Equivalent Mechanical performance prediction result of the random unsticking in interface under actual condition and the actual measurement mechanical property for preparing finished product
Than further guidance and improvement long fiber reinforcement composite material prepare end properties, and then solve long fiber reinforcement composite wood
Expect application field correlation engineering problem, and then sufficiently and safely plays composite material bearing structure mechanical property.
Further, fiber random distribution is established using Python based on random sequence additive process in the step 1
Composite three dimensional represent elementary volume, volume element RVE geometrical model, be to generate column long fibre by root in the base, and any two
Fiber cannot be overlapped or intersect, and include the following steps:
Step a: according to fiber radius r, fiber volume fraction VfThe three-dimensional representation elementary volume, volume element RVE established with fiber count n
Moulded dimension calculation formula are as follows:
Wherein: L is RVE model side length.
Step b: Python generating random number initial fiber F is utilized1Central coordinate of circle (x1,y1), generate second fiber F2
Central coordinate of circle (x2,y2), calculate F1, F2Center away from D:
Step c: if D >=2r+2d, wherein d is interfacial thickness, then second fiber is met the requirements, and then generates third root
Fiber and judging whether, which is all satisfied, is equal to 2r+2d away from being all larger than with preceding two fibrillar centers, if with any one fibrillar center away from
Above-mentioned requirements are unsatisfactory for, then regenerate third root fiber until meeting with the centers of preceding two fibers away from requiring.And so on,
Its cofibre and corresponding interface are gradually generated according to similar approach, need to meet the subsequent fiber gradually generated with existing fiber
Center away from be more than or equal to 2r+2d, until reaching required fiber volume fraction, that is, realize be based on random sequence addition side
Method establishes the three-dimensional representation elementary volume, volume element RVE geometrical model of fiber random distribution.
Further, apply periodic boundary condition in the step 4 with the required load that is displaced in the hope of different equivalent
Mechanical property, periodic boundary condition formula:
Wherein: u (x, y, z) indicates cyclic shift field, and subscript i represents x, y, the direction z, subscript+and-it respectively indicates along i
The positive direction and negative direction of axis, c'iIndicate constant.
Further, composite material Equivalent Mechanical performance is solved based on Homogenization in the step 5, obtained compound
Material basic material performance parameter equivalence value, homogenize formula:
Wherein, E is the equivalent elastic modulus of material, and v is the equivalent Poisson's ratio of material, and G is the equivalent shear modulus of material；
section nodeini-direction forcesini-directionIndicate institute in a certain boundary face of selected three-dimensional representation elementary volume, volume element RVE
There is node along the load in the direction i, sectional area indicates the boundary face area, axial strainini-directionIt indicates
The boundary face is acted on the lower strain along the direction i generated by the direction i load；straininj-directionIndicate three-dimensional representation volume
The strain that unit R VE boundary face is generated along the direction j, strainini-directionIndicate the boundary three-dimensional representation elementary volume, volume element RVE
The strain that face is generated along the direction i；Tensors of shear strain indicates a certain boundary face of three-dimensional representation elementary volume, volume element RVE
The tangential strain generated when outer carry is born, " i, j " respectively represent x, y, the direction z, and i ≠ j to subscript.Plane 0-x-y is each to same
Property face.
The utility model has the advantages that
1, in the equivalent modulus method of existing prediction long fiber reinforcement composite material, fibre/Ji Jie of physical presence is not considered
The influence of the random unsticking in face, causes predicted equivalent elastic modulus bigger than normal.A kind of long fiber reinforcement disclosed by the invention is compound
Material Equivalent Mechanical performance prediction method using finite element method establishes three based on the composite material different component method that homogenizes
Dimension represents elementary volume, volume element model, is reacted by interfacial detachment defective proportion and considers that fibre/random debonding defect in base interface long fibre increases
Strong composite inner is degenerated there are the random debonding defect situation in interface by the modulus to unsticking unit collection to simulate
It realizes the Equivalent Mechanical performance prediction under immobile interface debonding defect ratio situation, obtains being difficult to accurately measure interface by testing
Determine the composite material Equivalent Mechanical performance under unsticking ratio.
2, a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method disclosed by the invention, by changing interface
Debonding defect ratio realizes the Equivalent Mechanical performance prediction under the conditions of different interfacial detachment defective proportions；It is established according to prediction result
Interfacial detachment defective proportion influences relationship to Equivalent Mechanical performance, that is, establish the random unsticking in interface influences to close on Equivalent Mechanical performance
System by the Equivalent Mechanical performance prediction result of the random unsticking in interface under engineer application actual condition and is prepared on this basis
The actual measurement mechanical property of product compares, and further guidance and improvement long fiber reinforcement composite material prepare end properties, into
And solve the problems, such as long fiber reinforcement composite material application field correlation engineering.
3, a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method disclosed by the invention is based on mesomechanics
The finite element method of model, according to the basic mechanical property parameters of fiber, the material property of matrix, determination the random unsticking ratio in interface
Example can calculate corresponding Equivalent Mechanical performance parameter after composite molding, simple and easy.
Detailed description of the invention
Fig. 1 is a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method flow chart disclosed by the invention；
Fig. 2 is the RVE finite element model generated；
Fig. 3 is that there are the schematic diagrames of 6% random debonding defect at fibre/base interface；
Fig. 4 is the edge load for solving composite material Equivalent Mechanical performance parameter and applying to RVE model, wherein figure a is asked
Solve Ez；Scheme b and solves ExAnd vxy；Scheme c and solves Gzx；Scheme d and solves Gxy。
Fig. 5 is the concrete moduli change curve for considering the random unsticking in interface and influencing, wherein the figure equivalent transverse elasticity mould of a
Amount and corresponding decline percentage；Scheme b equivalent shear modulus and corresponding decline percentage.
Specific embodiment
Objects and advantages in order to better illustrate the present invention with reference to the accompanying drawing make furtherly the content of present invention
It is bright.
As shown in Figure 1, a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method disclosed in the present embodiment, tool
Body realizes that steps are as follows:
Step 1 establishes the three-dimensional representation elementary volume, volume element RVE model of fiber random distribution.
Step 1.1: utilizing r=3.55 μm of long fibre radius, fiber volume fraction Vf=60% and fibre number n=30 is logical
It crosses calculating and determines RVE side length L, and generate the three-dimensional cubic structure of resin matrix according to identified size.RVE moulded dimension
Calculation formula are as follows:
It is obtained after calculating and corresponds to VfL=44.49 μm of RVE side length when=60%.
Step 1.2: random selection is a little used as the center of circle in the one side of the cube of generation, and binding fiber diameter determines
First position of fibers.
Step 1.3: with the fiber center of circle and combining determining interfacial thickness d=0.068 μ on first fiber boundary of generation
M generates the fibre/base interface for being enclosed in first fiber peripheral.
Step 1.4: in the case where having generated first fiber, its cofibre and correspondence are gradually generated according to similar approach
Fibre/base interface, the subsequent fiber gradually generated need to be met and existing fiber is non-intersecting, until the fiber volume for reaching required contains
Until amount, that is, realize the three-dimensional representation elementary volume, volume element RVE geometry mould that fiber random distribution is established based on random sequence adding method
Type.Generation of random series method using Python based on Mason's Rotation Algorithm generates initial fiber F at random1Central coordinate of circle
(x1,y1), generate second fiber F2Central coordinate of circle (x2,y2), calculate F1, F2Center away from D:
If D >=7.236 μm, second fiber is met the requirements, and is then generated third root fiber and is judged whether to be all satisfied
With preceding two fibrillar centers away from be all larger than be equal to 7.236 μm, if above-mentioned away from being unsatisfactory for any generated fibrillar center
It is required that then regenerating third root fiber until meeting with the centers of preceding two fibers away from requiring.And so on, according to similar side
Method gradually generates its cofibre and corresponding interface, need to meet center of the subsequent fiber gradually generated with existing fiber away from big
In being equal to 7.236 μm, until reaching required fiber volume fraction, that is, realizes and fibre is established based on random sequence adding method
Tie up the three-dimensional representation elementary volume, volume element RVE geometrical model of random distribution.
Step 2, the three-dimensional representation elementary volume, volume element RVE model established to step 1 carries out FEM meshing, in institute
Node is extracted at the interface location of some fiber peripherals, and interface node set and consistent with its is generated according to extracted node
Node serial number list；Boundary element set and element number list are generated using extracted interface node.Including walking as follows
It is rapid:
Step 2.1: grid dividing is carried out to three-dimensional representation elementary volume, volume element RVE model.
Step 2.2: in fibre/base interface extraction node of each fiber peripheral, generating interface node set and and one
The node serial number list of cause.
Step 2.3: utilizing extracted interface node, generate boundary element set and element number list.
Step 3: extracting all fibres node, utilizes extracted fiber nodes generation fiber element set and consistent with its
Element number list, assign fibrous material parameter to the unit of the fiber element set of generation；Extract all substrates node, benefit
Matrix node set and consistent element number list therewith are generated with extracted matrix node, to the base unit set of generation
Unit assign basis material parameter；Unit concentration in generated interface randomly selects boundary element, will randomly select
Boundary element modulus degenerate as random debonding defect existing for interface, remaining boundary element is used as non-debonding interface list
Member assigns corresponding non-debonding interface material parameter, i.e. realization fiber element, base unit, non-debonding interface unit, unsticking circle
The assignment of the material parameter of face unit is handled.Wherein the fibrous material parameter of transverse isotropy is set as Ef1=184GPa, Ef2=
15GPa, vf1=0.2, vf2=0.2, Gf1=15GPa, Gf2=7GPa；Isotropic reisn base material parameter is set as Em=
3.06GPa vm=0.38.
Unit concentration in step 3 in generated interface randomly selects boundary element, preferably, having generated
Interface unit focus utilization Python in the generation of random series method based on Mason's Rotation Algorithm randomly select boundary
Face unit simultaneously establishes set as unsticking unit, and it is 1 that its modulus, which is arranged, to unsticking unit collection, and corresponding Poisson's ratio is set as 0.1, remains
Remaining non-debonding interface assigns above-mentioned basis material parameter.Fig. 2 is the RVE finite element model generated, and Fig. 3 is corresponding fibre/Ji Jie
Contain 6% random unsticking schematic diagram in face.
Step 4: in the generated fibre containing fixed proportion/random unsticking in base interface three-dimensional representation elementary volume, volume element RVE
Finite element model boundary applies periodic boundary condition and required displacement load.Periodic boundary condition formula:
Wherein: u (x, y, z) indicates cyclic shift field, and subscript i represents x, y, the direction z, subscript+and-it respectively indicates along i
The positive direction and negative direction of axis.
For seeking different equivalent mechanical property parameters, the displacement load being necessary to apply different, corresponding schematic diagram such as Fig. 4
It is shown.
Step 5: based on the three-dimensional representation elementary volume, volume element RVE finite element model in step 4, by interfacial detachment defective proportion
Reflection consider the long fiber reinforcement composite inner of the random debonding defect in interface there are the random debonding defect situation in interface, according to
The each component material parameter of step 3 assignment processing solves composite material Equivalent Mechanical performance based on Homogenization, that is, realizes
Equivalent Mechanical performance prediction under the fixed debonding defect ratio situation in interface.Homogenize formula:
Wherein, E is the equivalent elastic modulus of material, and v is the equivalent Poisson's ratio of material, and G is the equivalent shear modulus of material；
section nodeini-direction forcesini-directionIndicate institute in a certain boundary face of selected three-dimensional representation elementary volume, volume element RVE
There is node along the load in the direction i, sectional area indicates the boundary face area, axial strainini-directionIt indicates
The boundary face is acted on the lower strain along the direction i generated by the direction i load；straininj-directionIndicate three-dimensional representation volume
The strain that unit R VE boundary face is generated along the direction j, strainini-directionIndicate the boundary three-dimensional representation elementary volume, volume element RVE
The strain that face is generated along the direction i；Tensors of shear strain indicates a certain boundary face of three-dimensional representation elementary volume, volume element RVE
The tangential strain generated when outer carry is born, " i, j " respectively represent x, y, the direction z, and i ≠ j to subscript.Plane 0-x-y is each to same
Property face.
Step 6: the Equivalent Mechanical performance under the fixed debonding defect ratio situation in the interface based on step 1 to step 5 is pre-
Survey method, by changing interfacial detachment defective proportion to realize the Equivalent Mechanical performance under the conditions of different interfacial detachment defective proportions
Prediction；According to prediction result establish interfacial detachment defective proportion on Equivalent Mechanical performance influence relation curve, that is, establish interface with
Machine unsticking, which is realized, influences relation curve to Equivalent Mechanical performance.
Fig. 5 is the concrete moduli change curve for considering the random unsticking in interface and influencing.From figure it can be seen that in corpus fibrosum
Product content VfConsider that fibre/random unsticking in base interface concrete moduli calculated value is less than in composite material under=60% not considering
Fibre/random unsticking in base interface predicted value.Further from figure curve it can be found that concrete moduli with the random unsticking ratio in interface
Increase and non-linear downward trend is presented, and decrease speed is gradually increased.When unsticking non-compared to interface, when interface exists at random
When debonding defect ratio reaches 24%, equivalent transverse modulus of elasticity and equivalent shear modulus down ratio respectively reach 27% He
24%, this illustrates the influence of the consideration random debonding defect in interface that the method for the present invention can be more accurate to mechanical property, in advance
Composite material after machine-shaping is effectively predicted contains a certain proportion of fibre/influence of the random unsticking in base interface to mechanical property,
When being used as composite material in practical engineering applications bearing structure, it can avoid can not effectively playing composite material excellent
Mechanical property, while can also be to avoid the premature failure of bearing structure.
Step 7: being established random debonding interface based on step 6 and realized influences relationship to Equivalent Mechanical performance, according to step
Equivalent Mechanical performance prediction under the conditions of one to step 6 realization actual condition debonding interface defective proportion, by engineer application reality
The Equivalent Mechanical performance prediction result of random debonding interface and the actual measurement mechanical property for preparing finished product compare under operating condition, into one
Step guidance and improvement long fiber reinforcement composite material prepare end properties, and then solve long fiber reinforcement composite material application neck
Domain correlation engineering problem, and then sufficiently and safely play composite material bearing structure mechanical property.
According to prediction result establish interfacial detachment defective proportion on Equivalent Mechanical performance influence relationship it is as follows, fiber radius is
3.55 μm, interfacial thickness is 0.068 μm, fiber volume fraction 40%-60%, is divided into 5%, the random debonding defect ratio in interface
Example is 0-24%, is divided into 6%, and the random debonding defect ratio in median surface is 0% to indicate the bonding of interface perfection at this time, is not taken off
Viscous defect.The interfacial detachment defect and the equivalent transverse modulus of elasticity of long fiber composites and equivalent shear modulus relationship of prediction are such as
Shown in following table, wherein machine direction is set as in the z-direction:
Above-described specific descriptions have carried out further specifically the purpose of invention, technical scheme and beneficial effects
It is bright, it should be understood that the above is only a specific embodiment of the present invention, the protection model being not intended to limit the present invention
It encloses, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in the present invention
Protection scope within.
Claims (10)
1. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method, it is characterised in that: include the following steps,
Step 1 establishes the three-dimensional representation elementary volume, volume element RVE geometrical model of fiber random distribution based on random sequence adding method；
Step 2, the three-dimensional representation elementary volume, volume element RVE model established to step 1 carries out FEM meshing, all
At the interface location of fiber peripheral extract node, and according to extracted node generate interface node set and with its consistent section
Point numbered list；Boundary element set and element number list are generated using extracted interface node；
Step 3: extract finite element model in all fibres node, using extracted fiber nodes generation fiber element set and
With its consistent element number list, fibrous material parameter is assigned to the unit of the fiber element set of generation；Extract all bases
Body node generates matrix node set and consistent element number list therewith using extracted matrix node, to the base of generation
The unit of body unit set assigns basis material parameter；Unit concentration in generated interface randomly selects boundary element,
The boundary element modulus randomly selected degenerate as random debonding defect existing for interface, remaining boundary element is not as
Debonding interface unit assigns corresponding non-debonding interface material parameter, i.e. the assignment processing of realization each component material parameter；
Each component material parameter described in step 3 refers to fibrous material parameter in step 3, basis material parameter, debonding interface material
Parameter and non-debonding interface material parameter；
Step 4: limited in the generated fibre containing fixed proportion/random unsticking in base interface three-dimensional representation elementary volume, volume element RVE
Meta-model boundary applies periodic boundary condition and required displacement load；
Step 5: based on the three-dimensional representation elementary volume, volume element RVE finite element model in step 4, reflected by interfacial detachment defective proportion
There are the random debonding defect situations in interface for the long fiber reinforcement composite inner of the consideration random debonding defect in interface, according to step
The each component material parameter of three assignment processing solves composite material Equivalent Mechanical performance based on Homogenization, i.e. realization interface
Equivalent Mechanical performance prediction under fixed debonding defect ratio situation.
2. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as described in claim 1, it is characterised in that:
It further include step 6, the Equivalent Mechanical performance prediction under the fixed debonding defect ratio situation in the interface based on step 1 to step 5
Method, by changing interfacial detachment defective proportion to realize that the Equivalent Mechanical performance under the conditions of different interfacial detachment defective proportions is pre-
It surveys；Establishing interfacial detachment defective proportion according to prediction result influences relation curve to Equivalent Mechanical performance, that is, it is random to establish interface
Unsticking influences relation curve to Equivalent Mechanical performance.
3. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 2, it is characterised in that:
It further include step 7, being realized based on the random unsticking in interface that step 6 is established influences relation curve to Equivalent Mechanical performance, according to
Step 1 to step 6 realizes the Equivalent Mechanical performance prediction under the conditions of actual condition interfacial detachment defective proportion, by engineer application
The Equivalent Mechanical performance prediction result of the random unsticking in interface and the actual measurement mechanical property for preparing finished product compare under actual condition,
Further guidance and improvement long fiber reinforcement composite material prepare end properties, and then solve long fiber reinforcement composite material and answer
With field correlation engineering problem, and then sufficiently and safely play composite material bearing structure mechanical property.
4. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 1,2 or 3, feature
Be: step 1 implementation method is as follows,
Step 1.1: three-dimensional representation elementary volume, volume element is determined according to the calculating of long fibre radius, fiber volume fraction and fibre number
RVE geometrical model size, and the three-dimensional cubic structure according to identified size generation resin matrix；
Step 1.2: random selection is a little used as the center of circle in the one side of the three-dimensional cube of generation, and binding fiber diameter determines
First position of fibers；
Step 1.3: on first fiber boundary of generation, being enclosed in first fiber peripheral in conjunction with determining interfacial thickness generation
Fibre/base interface；
Step 1.4: in the case where having generated first fiber, its cofibre and corresponding is gradually generated according to similar approach
Fibre/base interface, need to meet the subsequent fiber gradually generated and existing fiber is non-intersecting, until reaching required fiber volume fraction
Until, that is, realize the three-dimensional representation elementary volume, volume element RVE geometrical model that fiber random distribution is established based on random sequence adding method.
5. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 4, it is characterised in that:
Step 2 implementation method is as follows,
Step 2.1: grid dividing is carried out to three-dimensional representation elementary volume, volume element RVE model；
Step 2.2: extracting node at the interface location of each fiber peripheral, generate interface node set and consistent with it
Node serial number list；
Step 2.3: utilize extracted interface node, generate boundary element set and with its consistent element number list.
6. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 5, it is characterised in that:
Establish the three-dimensional representation elementary volume, volume element of fiber random distribution in step 1 using Python based on random sequence adding method
RVE geometrical model.
7. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 6, it is characterised in that:
In step 3, the random sequence based on Mason's Rotation Algorithm in the unit focus utilization Python of generated interface is raw
Boundary element is randomly selected as unsticking unit at method.
8. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 7, it is characterised in that:
Establish the composite three dimensional generation of fiber random distribution in the step 1 using Python based on random sequence additive process
Table elementary volume, volume element RVE geometrical model is to generate column long fibre by root in the base, and any two fibers cannot be overlapped or phase
It hands over, includes the following steps,
Step a: according to fiber radius r, fiber volume fraction VfThe three-dimensional representation elementary volume, volume element RVE model scale established with fiber count n
Very little calculation formula are as follows:
Wherein: L is RVE model side length.
Step b: Python generating random number initial fiber F is utilized1Central coordinate of circle (x1,y1), generate second fiber F2Circle
Heart coordinate (x2,y2), calculate F1, F2Center away from D:
Step c: if D >=2r+2d, wherein d is interfacial thickness, then second fiber is met the requirements, and then generates third root fiber
And judge whether to be all satisfied with preceding two fibrillar centers away from being all larger than equal to 2r+2d, if with any one fibrillar center away from discontented
Sufficient above-mentioned requirements then regenerate third root fiber until meeting with the centers of preceding two fibers away from requiring；And so on, according to
Similar approach gradually generates its cofibre and corresponding interface, need to meet the subsequent fiber gradually generated and in existing fiber
The heart until reaching required fiber volume fraction, that is, is realized and is built based on random sequence adding method away from 2r+2d is more than or equal to
The three-dimensional representation elementary volume, volume element RVE geometrical model of vertical fiber random distribution.
9. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 8, it is characterised in that:
Apply periodic boundary condition in the step 4 from required displacement load in the hope of different Equivalent Mechanical performances, periodical side
Boundary's condition formula:
Wherein: u (x, y, z) indicates cyclic shift field, and subscript i represents x, y, the direction z, subscript+and-it respectively indicates along i axis
Positive direction and negative direction, c 'iIndicate constant.
10. a kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method as claimed in claim 9, feature exist
In: composite material Equivalent Mechanical performance is solved based on Homogenization in the step 5, obtains the basic material of composite material
Energy parameter equivalent value, homogenize formula:
Wherein, E is the equivalent elastic modulus of material, and v is the equivalent Poisson's ratio of material, and G is the equivalent shear modulus of material；
section nodein i-direction forcesin i-directionIt indicates in a certain boundary face of selected three-dimensional representation elementary volume, volume element RVE
For all nodes along the load in the direction i, sectional area indicates the boundary face area, axial strainin i-directionTable
Show that the boundary face is acted on the lower strain along the direction i generated by the direction i load；strainin j-directionIndicate three-dimensional representation body
The strain that product unit RVE boundary face is generated along the direction j, strainin i-directionIndicate the side three-dimensional representation elementary volume, volume element RVE
The strain that interface is generated along the direction i；Tensors of shear strain indicates a certain boundary three-dimensional representation elementary volume, volume element RVE
The tangential strain generated when outer carry is born in face, and " i, j " respectively represent x, y, the direction z, and i ≠ j to subscript；Plane 0-x-y be it is each to
Same sex face.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910821548.5A CN110532702B (en) | 2019-09-02 | 2019-09-02 | Long fiber reinforced composite material equivalent mechanical property prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910821548.5A CN110532702B (en) | 2019-09-02 | 2019-09-02 | Long fiber reinforced composite material equivalent mechanical property prediction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110532702A true CN110532702A (en) | 2019-12-03 |
| CN110532702B CN110532702B (en) | 2021-04-27 |
Family
ID=68666011
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910821548.5A Active CN110532702B (en) | 2019-09-02 | 2019-09-02 | Long fiber reinforced composite material equivalent mechanical property prediction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110532702B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111639451A (en) * | 2020-06-08 | 2020-09-08 | 吉林大学 | Refined modeling simulation method of two-dimensional plain woven fiber reinforced composite material |
| CN112632819A (en) * | 2020-12-18 | 2021-04-09 | 沈阳航空航天大学 | Method for predicting basic mechanical property parameters of continuous fiber reinforced composite material |
| CN112666051A (en) * | 2020-12-11 | 2021-04-16 | 中航复合材料有限责任公司 | Method for evaluating interface performance of fiber reinforced composite material and application thereof |
| CN112733408A (en) * | 2021-02-23 | 2021-04-30 | 江西省科学院应用物理研究所 | Method for generating high volume fraction two-dimensional numerical model by combining fiber micromotion and hard filling |
| CN113312824A (en) * | 2021-06-16 | 2021-08-27 | 西北工业大学 | Mesomechanics-based unidirectional fiber composite material mechanical property prediction method |
| CN113570544A (en) * | 2021-05-31 | 2021-10-29 | 清华大学 | Cancellous bone modeling method and device, storage medium and electronic equipment |
| CN114023400A (en) * | 2021-10-19 | 2022-02-08 | 上海索辰信息科技股份有限公司 | Method for rapidly predicting equivalent characteristics of composite material under different volume fractions |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108090963A (en) * | 2017-12-25 | 2018-05-29 | 大连理工大学 | A kind of numerical computation method of fibre reinforced composites thermal residual strain at low temperature |
| CN110175419A (en) * | 2019-05-30 | 2019-08-27 | 新疆大学 | Fan blade composite material mesomechanics damage development analysis method |
-
2019
- 2019-09-02 CN CN201910821548.5A patent/CN110532702B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108090963A (en) * | 2017-12-25 | 2018-05-29 | 大连理工大学 | A kind of numerical computation method of fibre reinforced composites thermal residual strain at low temperature |
| CN110175419A (en) * | 2019-05-30 | 2019-08-27 | 新疆大学 | Fan blade composite material mesomechanics damage development analysis method |
Non-Patent Citations (4)
| Title |
|---|
| A. CAPORALE等: "Micromechanical analysis of interfacial debonding in unidirectional fiber-reinforced composites", 《COMPUTERS AND STRUCTURES》 * |
| LINA RIAÑO等: "Validation of a Representative Volume Element for unidirectional fiber reinforced composites: Case of a monotonic traction in its cross section", 《COMPOSITE STRUCTURES》 * |
| SADIK L.OMAIREY等: "Development of an ABAQUS plugin tool for periodic RVE homogenisation", 《ENGINEERING WITH COMPUTERS》 * |
| 金杰: "基于细观力学理论的颗粒增强复合材料脱粘损伤研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111639451A (en) * | 2020-06-08 | 2020-09-08 | 吉林大学 | Refined modeling simulation method of two-dimensional plain woven fiber reinforced composite material |
| CN112666051A (en) * | 2020-12-11 | 2021-04-16 | 中航复合材料有限责任公司 | Method for evaluating interface performance of fiber reinforced composite material and application thereof |
| CN112666051B (en) * | 2020-12-11 | 2022-06-10 | 中航复合材料有限责任公司 | Method for evaluating interface performance of fiber reinforced composite material and application thereof |
| CN112632819A (en) * | 2020-12-18 | 2021-04-09 | 沈阳航空航天大学 | Method for predicting basic mechanical property parameters of continuous fiber reinforced composite material |
| CN112632819B (en) * | 2020-12-18 | 2024-01-02 | 沈阳航空航天大学 | Continuous fiber reinforced composite material basic mechanical property parameter prediction method |
| CN112733408A (en) * | 2021-02-23 | 2021-04-30 | 江西省科学院应用物理研究所 | Method for generating high volume fraction two-dimensional numerical model by combining fiber micromotion and hard filling |
| CN112733408B (en) * | 2021-02-23 | 2023-03-31 | 江西省科学院应用物理研究所 | Method for generating high volume fraction two-dimensional numerical model by combining fiber micromotion and hard filling |
| CN113570544A (en) * | 2021-05-31 | 2021-10-29 | 清华大学 | Cancellous bone modeling method and device, storage medium and electronic equipment |
| CN113312824A (en) * | 2021-06-16 | 2021-08-27 | 西北工业大学 | Mesomechanics-based unidirectional fiber composite material mechanical property prediction method |
| CN113312824B (en) * | 2021-06-16 | 2024-03-19 | 西北工业大学 | Method for predicting mechanical properties of unidirectional fiber composite material based on mesomechanics |
| CN114023400A (en) * | 2021-10-19 | 2022-02-08 | 上海索辰信息科技股份有限公司 | Method for rapidly predicting equivalent characteristics of composite material under different volume fractions |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110532702B (en) | 2021-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110532702A (en) | A kind of long fiber reinforcement composite material Equivalent Mechanical performance prediction method | |
| Devireddy et al. | Physical and mechanical behavior of unidirectional banana/jute fiber reinforced epoxy based hybrid composites | |
| CN105022898B (en) | A kind of the crack performance measurement and optimization method of adhesive bonding of composites structure | |
| Rachchh et al. | Mechanical characterisation of rattan fibre polyester composite | |
| Sapuan et al. | Tensile and flexural strengths of coconut spathe-fibre reinforced epoxy composites | |
| Gohil et al. | Analytical investigation and comparative assessment of interphase influence on elastic behavior of fiber reinforced composites | |
| CN105653775B (en) | One kind being mechanically fixed engaging portion modeling method | |
| Khaji et al. | Strain energy release rate in combination with reinforcement isotropic solid model (SERIS): A new mixed-mode I/II criterion to investigate fracture behavior of orthotropic materials | |
| Adamczak et al. | Estimating the approximation uncertainty for digital materials subjected to stress relaxation tests | |
| CN105468826B (en) | The design method of composite material | |
| Nasir et al. | Experimental study towards determination of critical feed for minimization of delamination damage in drilling flax natural fibre composites | |
| CN105046076B (en) | Three layers of single cell structure micro mechanical property calculation method of carbon fibre composite are laminated | |
| Rajulu et al. | Void content, density and weight reduction studies on short bamboo fiber–epoxy composites | |
| Mohammadabadi et al. | Influence of a biaxially corrugated core geometry on flexural stiffness of wood-strand composite sandwich panels | |
| Chokshi et al. | Experimental investigation and mathematical modeling of longitudinally placed natural fiber reinforced polymeric composites including interphase volume fraction | |
| Pradhan et al. | Study of mechanical and abrasive wear properties of lantana camara particulate reinforced epoxy composite | |
| Kalita et al. | Mechanical characterization and finite element investigation on properties of PLA-jute composite | |
| US6333092B1 (en) | Fractal interfacial enhancement of composite delamination resistance | |
| Sathiyamurthy et al. | Modelling and optimization of mechanical behaviors of Al 2 O 3-coir-polyester composites using response surface methodology | |
| Bere et al. | Fabrication and mechanical characterization of short fiber-glass epoxy composites | |
| Monteiro et al. | Interfacial shear strength in lignocellulosic fibers incorporated polymeric composites | |
| CN110348166A (en) | A kind of virtual materials parameter visualization recognition methods of basalt fibre resin concrete joint surface | |
| Prabhakaran et al. | Modelling and simulation of natural fibre/epoxy composites-prediction of stress state and deformations | |
| Pereira et al. | Modeling and analysis for wear performance of coconut shell powder filled glass fiber composite using Taguchi approach | |
| Kumar et al. | Optimization of process parameters in drilling of short fiber (Kenaf) composite |
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 |