CN108168997A - D braided composites method for testing performance based on space group P* symmetry - Google Patents
D braided composites method for testing performance based on space group P* symmetry Download PDFInfo
- Publication number
- CN108168997A CN108168997A CN201810026755.7A CN201810026755A CN108168997A CN 108168997 A CN108168997 A CN 108168997A CN 201810026755 A CN201810026755 A CN 201810026755A CN 108168997 A CN108168997 A CN 108168997A
- Authority
- CN
- China
- Prior art keywords
- stress
- space group
- symmetry
- braided composites
- braided
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- 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]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0216—Finite elements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses based on space groupThe D braided composites method for testing performance of symmetry, is related to field of compound material field, can detect based on space groupWhether the D braided composites of symmetry fail.The present invention is based on space groupsSymmetrical structure, it is first determined the structure snd size of the minimum representative unit cell of this three-dimensional knitting material；Then model is established using Abaqus finite element softwares, simulate the axial compression test of D braided composites hexagonal prisms pipe and obtains result of calculation；According to Hashin criterion and Von Mises stress criterions, prediction is based on space groupThe maximum stress that can be born before the woven composite failure of symmetry.The present invention is solved because of this non-mass production of novel three-dimensional woven composite, the problem of mechanical property is difficult to test, different models can be established according to different braided materials, can also different models be established according to different braiding conditions, the optimal case to make this composite material braiding provides reliable theoretical foundation.
Description
Technical field
The present invention relates to field of compound material more particularly to based on space groupThe D braided composites of symmetry
Method for testing performance.
Background technology
It is low and be easily layered that D braided composites overcome poor performance between conventional composite materials laminate, impact flexibility
Shortcoming has many advantages, such as good tensile property, shock resistance and cutting performance, is widely used in the high-techs such as Aeronautics and Astronautics
Skill field.Based on space groupThe D braided composites of symmetry are according to crystal space dot matrix and space symmetr group's structure
The novel weaved composite material found out, due to the also non-mass production of the material of this structure, so to the prediction of its performance just
Seem particularly significant.
However still lack be directed to based on space group in the prior artThe property of the D braided composites of symmetry
It can detection method.
Invention content
The present invention provides based on space groupThe D braided composites method for testing performance of symmetry, can detect
Based on space groupWhether the D braided composites of symmetry fail.
In order to achieve the above objectives, the present invention adopts the following technical scheme that, based on space groupThe 3 D weaving of symmetry is answered
Condensation material method for testing performance, including：
S1, under the premise of periodical principle and integrality principle is met, according toThe structure size of space group marks off
The minimum representative volume unit cell of three-dimensional knitting material, and establish unit cell geometrical model；
S2, minimum representative volume unit cell is replicated, and is combined by way of array, obtain standarized component, standarized component
With matrix is compound obtains stressed member；
The material parameter of S3, acquisition D braided composites to be measured；
S4, material parameter is given to stressed member；
S5, basis are compressed axially boundary condition, load and the contact conditions of experimental conditions setting stressed member, and to component
Grid division；
S6, quasi-static and dynamic axial compression stress FEM calculation is carried out to stressed member, obtains stress-strain song
Line, stress distribution cloud atlas, Strain Distribution cloud atlas；
S7, the direct stress, shear stress are brought into failure criteria, obtains compressive strength/compression of three-dimensional knitting material
Modulus and failure conditions.When the plastic strain of material is equal to or more than the equivalent plastic strain specified, it is believed that material occurs
It destroys, obtains the greatest compressive strength of D braided composites.
Further, according toThe structure of hexagonal prisms in the single cell structure of the structure determination three-dimensional knitting material of space group,
Its size is hexagon bottom surface length of side 4mm, and a height of 12mm, the section radius of cylindrical fiber is 1.3mm.
Further, the shape of minimum representative volume unit cell is hexagonal prisms tubulose, by a hexagonal prism fiber and six
The cylinder fiber of root different angle is formed.
Further, in S3, material parameter includes：Density, elasticity modulus, Poisson's ratio, yield stress, tangent modulus,
Hardening parameter, failure strain, true stress and plastic strain parameter.
Further, failure criteria includes Hashin criterion and Von-Mises stress criterions.
Beneficial effects of the present invention are：
It is imitated and adapted to using Finite Element Method and is complicated based on space groupThe 3 D weaving composite wood of symmetry
Material establishes accurately stressed member space structure and size, then different material property parameters is assigned to stressed member, with
And boundary condition, contact and load under the conditions of actual loaded, it goes to approach unlimited unknown quantity using the unknown quantity of limited quantity
Real system, measurement result are accurate；And then it solves to obtain the power of the D braided composites by finite element analysis software
Performance is learned, compared with laboratory facilities, Finite Element Method can deeply probe into the D braided composites in dynamic at a high speed
The mechanical response and failure mechanism in loading process are compressed, cumbersome production, experimental procedure is avoided, has saved time and economy
Cost, confidence level as a result is also very high, meets the requirement that digital simulation technology is widely used in engineer application.
Description of the drawings
It to describe the technical solutions in the embodiments of the present invention more clearly, below will be to needed in the embodiment
Attached drawing is briefly described, it should be apparent that, the accompanying drawings in the following description is only some embodiments of the present invention, for ability
For the those of ordinary skill of domain, without creative efforts, it can also be obtained according to these attached drawings other attached
Figure.
Fig. 1 is space groupThe layout viewing of hexagonal prisms and cylindrical fiber in structure；
Fig. 2 is based on space groupThe structure chart of the woven composite unit cell of symmetry；
Fig. 3 is the fibre bundle geometrical model of analog compression experiment；
Fig. 4 is the fiber of analog compression experiment and the compound geometrical model vertical view of matrix；
Fig. 5 is the fiber of analog compression experiment and the compound geometrical model of matrix.
Specific embodiment
For those skilled in the art is made to more fully understand technical scheme of the present invention, With reference to embodiment to this
Invention is described in further detail.
An embodiment of the present invention provides based on space groupThe D braided composites method for testing performance of symmetry, this
Embodiment is described in more detail technical scheme of the present invention by taking carbon fiber and epoxy resin as an example.
Based on space groupThe D braided composites method for testing performance of symmetry, including：
S1, assume carbon fiber and epoxy resin complete wetting, between the two without interface sliding, closed according to Elastic-plastic Constitutive
System, establishes the relevant plasticity moving model of a strain rate.
It is mutually squeezed between carbon fiber, matrix, intermediate main carbon fiber section is hexagonal shape, and main carbon fiber is surrounded on side
Carbon fiber section be circular shape.Carbon fiber section is not damaged during remaining unchanged and weaving along its length.Carbon is fine
Peacekeeping matrix is according to space group in Fig. 1Structure determination individual carbon fibers braid angle, element length and arrangement mode, one
A unit cell is made of the fiber of hexagonal prism and the cylinder fiber of six roots of sensation different angle.
It is established in Abaqus based on space groupThe representative unit cell of minimum of the woven composite of symmetry, Fig. 2 are
The structure chart of unit cell.The hexagonal prisms that several unit cells and an inside dig out respective shapes are assembled up, such as Fig. 3, Fig. 4, Fig. 5
It is shown.
The material parameter of S3, acquisition D braided composites to be measured, material parameter include：Density, elasticity modulus, Poisson
Than, yield stress, tangent modulus, hardening parameter, failure strain parameter.
S4, to fibre bundle assign carbon fiber attribute, to cuboid matrix assign epoxy resin attribute, including density,
Elasticity modulus, Poisson's ratio, yield stress, tangent modulus, hardening parameter, failure strain parameter.
The way of contact of carbon fiber and epoxy resin is defined as plane-plane contact, sets epoxy resin as interarea, carbon fiber
Shu Weicong faces.
S5, using the highest C3D8R solid elements of computational accuracy come to fibre bundle grid division, using C3D4 solid elements
To matrix grid division.
Z-direction is defined as by machine direction by local coordinate, the bottom surface of cuboid is fixed, is applied in the upper end of cuboid
The axial compression displacement load for adding a direction upward, sets several groups of different compression speeds as the case may be and displacement is long
Degree loads certain time.
S6, for carbon fiber, the damage type of fibre bundle is judged using Hash i n criterion.
The fibre bundle compressive damage criterion in L directions is：
Fibre bundle compression and shear property model criterion on TZ directions are：
Wherein,Represent fibre bundle L respectively, T, the compressive strength in Z-direction,Table respectively
Show fibre bundle in LT, TZ, the shear strength on ZL directions.σL, σT, σZ, σLT,σZL,σTZIt represents fibre bundle part direct stress and cuts
Shearing stress.In this model, L, T, Z-direction distinguish the X of representative model, and Y, Z-direction, LT, TZ, ZL directions represent XY, YZ respectively,
ZX directions.
For epoxy resin-base, matrix damage process is judged using Von-Mises stress criterions：
Wherein, σmFor the intensity of basis material, σ1、σ2、σ3Represent first, second, third principal stress.τ12, τ23, τ31It represents
The shear stress of aggregate material.
Stress distribution, misalignment in result of calculation is come the failure conditions of judgement material.
Beneficial effects of the present invention are：
It is imitated and adapted to using Finite Element Method and is complicated based on space groupThe 3 D weaving composite wood of symmetry
Material establishes accurately stressed member space structure and size, then different material property parameters is assigned to stressed member, with
And boundary condition, contact and load under the conditions of actual loaded, it goes to approach unlimited unknown quantity using the unknown quantity of limited quantity
Real system, measurement result are accurate；And then it solves to obtain the power of the D braided composites by finite element analysis software
Performance is learned, compared with laboratory facilities, Finite Element Method can deeply probe into the D braided composites in dynamic at a high speed
The mechanical response and failure mechanism in loading process are compressed, cumbersome production, experimental procedure is avoided, has saved time and economy
Cost, confidence level as a result is also very high, meets the requirement that digital simulation technology is widely used in engineer application.
More than, only specific embodiment of the invention, but protection scope of the present invention is not limited thereto is any to be familiar with
In the technical scope disclosed by the present invention, the change or replacement that can be readily occurred in should all be covered those skilled in the art
Within protection scope of the present invention.Therefore, protection scope of the present invention should be subject to the protection scope in claims.
Claims (5)
1. based on space groupThe D braided composites method for testing performance of symmetry, which is characterized in that including：
S1, under the premise of periodical principle and integrality principle is met, according toThe structure size of space group, marks off three-dimensional
The minimum representative volume unit cell of braided material, and establish a unit cell geometrical model；
S2, the minimum representative volume unit cell is replicated, and is combined by way of array, obtain standarized component, the standard
Component and matrix is compound obtains stressed member；
The material parameter of S3, acquisition D braided composites to be measured；
S4, the material parameter is given to the stressed member；
S5, basis are compressed axially boundary condition, load and contact conditions that experimental conditions set the stressed member, and division
Part grid；
S6, quasi-static and dynamic axial compression stress FEM calculation is carried out to the stressed member, obtains stress-strain song
Line, stress distribution cloud atlas, Strain Distribution cloud atlas, so as to obtain direct stress and shear stress；
S7, the direct stress, the shear stress are brought into failure criteria, obtains compressive strength/compression of three-dimensional knitting material
Modulus and failure conditions, when the plastic strain of material is equal to or more than the equivalent plastic strain specified, it is believed that material occurs
It destroys, obtains the greatest compressive strength of D braided composites.
2. according to the method described in claim 1, it is characterized in that, according toThe structure determination three-dimensional knitting material of space group
The structure of hexagonal prisms in single cell structure, size be hexagon bottom surface length of side 4mm, a height of 12mm, the section radius of cylindrical fiber
For 1.3mm.
3. according to the method described in claim 1, it is characterized in that, the shape of the minimum representative volume unit cell is hexagonal prisms
Tubulose.
4. according to the method described in claim 1, it is characterized in that, in the S3, the material parameter includes：Density, bullet
Property modulus, Poisson's ratio, yield stress, tangent modulus, hardening parameter, failure strain, true stress and plastic strain parameter.
5. according to the method described in claim 1, it is characterized in that, the failure criteria includes Hashin criterion and Von-
Mises stress criterions.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810026755.7A CN108168997A (en) | 2018-01-11 | 2018-01-11 | D braided composites method for testing performance based on space group P* symmetry |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810026755.7A CN108168997A (en) | 2018-01-11 | 2018-01-11 | D braided composites method for testing performance based on space group P* symmetry |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108168997A true CN108168997A (en) | 2018-06-15 |
Family
ID=62517949
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810026755.7A Pending CN108168997A (en) | 2018-01-11 | 2018-01-11 | D braided composites method for testing performance based on space group P* symmetry |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108168997A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114139308A (en) * | 2021-11-23 | 2022-03-04 | 南京航空航天大学 | Resin-based composite material damage prediction method considering nonlinear characteristics of matrix |
| CN114184497A (en) * | 2020-09-14 | 2022-03-15 | 广州汽车集团股份有限公司 | Method and device for testing fabric shearing performance parameters |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005250767A (en) * | 2004-03-03 | 2005-09-15 | Yokohama Rubber Co Ltd:The | Dynamic characteristic simulation method for composite material and dynamic behavior simulation device for composite material |
| CN201908202U (en) * | 2010-12-24 | 2011-07-27 | 河南科技大学 | Three-dimensional knitting material |
| WO2011108468A1 (en) * | 2010-03-02 | 2011-09-09 | 日本電気株式会社 | Material constant estimation system and method of estimating material constant |
| CN102358990A (en) * | 2011-09-21 | 2012-02-22 | 河南科技大学 | Space group P* symmetry based three-dimensional braiding process and technology equipment thereof |
| CN103366085A (en) * | 2013-06-08 | 2013-10-23 | 上海交通大学 | Multiscale prediction method for mechanical property of woven composite material |
| CN103455712A (en) * | 2013-08-20 | 2013-12-18 | 北京航空航天大学 | Method for predicting mechanical modulus of triaxial weaving type composite material grid |
| CN105354390A (en) * | 2015-12-01 | 2016-02-24 | 北京航空航天大学 | New method for predicting two-way tensile modulus and strength of planar oblique woven composite material |
-
2018
- 2018-01-11 CN CN201810026755.7A patent/CN108168997A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005250767A (en) * | 2004-03-03 | 2005-09-15 | Yokohama Rubber Co Ltd:The | Dynamic characteristic simulation method for composite material and dynamic behavior simulation device for composite material |
| WO2011108468A1 (en) * | 2010-03-02 | 2011-09-09 | 日本電気株式会社 | Material constant estimation system and method of estimating material constant |
| CN201908202U (en) * | 2010-12-24 | 2011-07-27 | 河南科技大学 | Three-dimensional knitting material |
| CN102358990A (en) * | 2011-09-21 | 2012-02-22 | 河南科技大学 | Space group P* symmetry based three-dimensional braiding process and technology equipment thereof |
| CN103366085A (en) * | 2013-06-08 | 2013-10-23 | 上海交通大学 | Multiscale prediction method for mechanical property of woven composite material |
| CN103455712A (en) * | 2013-08-20 | 2013-12-18 | 北京航空航天大学 | Method for predicting mechanical modulus of triaxial weaving type composite material grid |
| CN105354390A (en) * | 2015-12-01 | 2016-02-24 | 北京航空航天大学 | New method for predicting two-way tensile modulus and strength of planar oblique woven composite material |
Non-Patent Citations (3)
| Title |
|---|
| WENSUO MA ET AL: "Processing and properties of plate 3D braided material based on space croup P3 symmetry", 《APPLIED MECHANICS AND MATERIALS》 * |
| 张芳芳: "编织复合材料力学性能及热物理性能预报研究", 《中国博士学位论文全文数据库 工程技术I辑》 * |
| 张超等: "三维编织复合材料渐进损伤模拟及强度预测", 《复合材料学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114184497A (en) * | 2020-09-14 | 2022-03-15 | 广州汽车集团股份有限公司 | Method and device for testing fabric shearing performance parameters |
| CN114184497B (en) * | 2020-09-14 | 2023-12-22 | 广州汽车集团股份有限公司 | Method and device for testing shear performance parameters of fabric |
| CN114139308A (en) * | 2021-11-23 | 2022-03-04 | 南京航空航天大学 | Resin-based composite material damage prediction method considering nonlinear characteristics of matrix |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhao et al. | Modeling the transverse tensile and compressive failure behavior of triaxially braided composites | |
| CN103942441B (en) | Carbon fibre composite estimating method for fatigue life based on stress ratio impact | |
| Zhu et al. | Micromechanical model for simulating the fracture process of rock | |
| CN106777769B (en) | Predict the finite element method of the progressive failure of composite material by multilayer slab under low velocity impact | |
| Penava et al. | Determination of the elastic constants of plain woven fabrics by a tensile test in various directions | |
| CN110348110B (en) | Self-adaptive identification method for rigidity of bolt-to-connector | |
| CN113420376B (en) | Multi-scale-based impact mechanical property simulation method for carbon fiber composite material | |
| CN111832209B (en) | Dynamic progressive failure analysis method for multi-scale model of composite material | |
| CN110175419A (en) | Fan blade composite material mesomechanics damage development analysis method | |
| Gao et al. | Failure analysis of C/SiC composites plate with a hole by the PFA and DIC method | |
| CN111666699B (en) | Rock mass engineering trans-scale simulation calculation method based on REV full-area coverage | |
| CN112100765A (en) | High-low cycle composite fatigue test piece of turbine disc cold extrusion reinforced hole structure and design method | |
| CN106874547B (en) | A method of the prediction full wound composite gas cylinder burst pressure of fiber | |
| CN108709856A (en) | A kind of concrete structure crack monitoring and warning system and method for early warning | |
| CN108168997A (en) | D braided composites method for testing performance based on space group P* symmetry | |
| Anıl et al. | Low velocity impact behavior of shear deficient RC beam strengthened with cfrp strips | |
| Zheng et al. | Arch-dam crack deformation monitoring hybrid model based on XFEM | |
| CN110502793A (en) | A kind of unidirectional ceramic matric composite off-axis tension testpieces optimum design method | |
| CN108197398A (en) | A kind of finite element method of D braided composites failure predicted based on space group P4 | |
| Bozkurt et al. | Simulation of drop-weight impact test on composite laminates using finite element method | |
| CN105046076A (en) | Method for calculating micro-mechanical properties of three-layer single-cell structure of laminated carbon fibre composite material | |
| Marinescu et al. | Comparative Study of Tensile Behaviour for Layered Polymeric Composite Materials | |
| Kikuchi | Study on multiple surface crack growth and coalescence behaviors | |
| CN112632737B (en) | Method and device for determining neutral axis of embedded part | |
| De Luycker et al. | Numerical analysis of 3D interlock composite preforming |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180615 |