CN108388763A - A kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis - Google Patents
A kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis Download PDFInfo
- Publication number
- CN108388763A CN108388763A CN201810213287.4A CN201810213287A CN108388763A CN 108388763 A CN108388763 A CN 108388763A CN 201810213287 A CN201810213287 A CN 201810213287A CN 108388763 A CN108388763 A CN 108388763A
- Authority
- CN
- China
- Prior art keywords
- scale
- reliability
- calculation
- macro
- multiscale analysis
- 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
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention discloses a kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis.Structural reliability uses Monte Carlo method sample calculation；Couple multiscale analysis and the response of finite element method Statics of Structures, it realizes that meso-scale is transmitted to macro-scale information with homogenization theory, transmission of the macro-scale to meso-scale information is realized by regarding strain regime on macro-scale structures element integral point as thin cell element boundary condition of seeing；By ABAQUS Pathon secondary developing platforms, so that single sample calculation process and sampling process is realized in ABAQUS, achieved the effect that sampling, modeling, calculate and statistic processes automation；The present invention has differentiates more acurrate, the more accurate advantage of result of calculation than Analysis of structural reliability of the tradition based on macro-scale, failure mode.
Description
Technical field
The invention belongs to composite structure design fields, are related to multiscale analysis technology and reliability calculation method, tool
Body is related to a kind of oriented fiber reinforced composite laminate structures analysis method for reliability based on multiscale analysis.
Background technology
Composite material refers to passing through combination process group by a few class different materials such as organic polymer, inorganic non-metallic or metal
New material made of conjunction.Laminate structures are widely used a kind of structures in composite structure, there is the fiber of different directions
Laying stack composition, have stronger design.Traditional structural design and strength check are based on the deterministic design, that is, design
Variable is determining amount, however in composite structure, there is uncertainty in load, structure size, material property, make certainty
Design result and actual conditions have prodigious deviation.Structural reliability design can make up the deficiency of traditional the deterministic design, will
Uncertain factor in structure is in view of in design, making structural safety performance and more reasonably judging.On the other hand, it compares
Composite material macromechanics, composite material mesomechanics have elastic properties of materials performance, failure law more essential announcement, power
It is more acurrate to learn model, Gu, reliability results more specific to the differentiation of material damage pattern based on the structural reliability of the latter
It is more accurate；Calculation amount can be caused excessive again if establishing the analysis model of multiple material structure from meso-scale completely, or even with present
Calculated level be likely to not solve, multiscale analysis method had not only ensured macro-scale analysis efficiency but also can be to the weight in structure
Key component is wanted to be analyzed on meso-scale.In the research of existing composite structure reliability design, based on macroscopical ruler
The research of degree analysis accounts for the overwhelming majority[1-11].The shortcomings that this kind of design, is：Composite structure reliability design is set to rest on
In macro-scale, designability is confined on the macro-scales variables such as wing flapping design, overlay thickness, cannot pair and broad perspectives
The closely related thin sight characteristic of energy is analyzed and is designed；Composite material macro-scale mechanical property model and strength model are often
Not as good as the factor abundance that meso mechanical model considers, especially strength model, the former is accurate not as good as the latter so that structural reliability
Whether safe it is unable to authentic assessment structure,.For the existing composite structure fail-safe analysis based on multiscale analysis
Research[12-13], do not consider carefully to see criterion of strength, only excavates out potentiality of the multiscale analysis technology in elastic property prediction, and
There are the micro-stress distribution of larger impact, thin failure behaviour of seeing to predict not yet to consider reliability results.
Invention content
This method is a kind of oriented fiber reinforced composite laminate structures Calculation of Reliability side based on multiscale analysis
Method divides following steps：
Step 1：The present invention considers the stochastic uncertainty in directional fiber composite structure.It specifically considers following
The uncertainty of parameter：The elastic property parameter of fiber, matrix, fiber, matrix intensive parameter, macroscopical load and macroscopic view knot
The geometric parameter of structure.
Step 2：The present invention uses the Caro method sample calculation structural reliability covered, and is once to take out in step 3 to step 14
Sample calculation process.
Step 3：Establish macrostructure geometrical model.
Step 4：By the macrostructure geometrical model discretization in step 3.
Step 5：Establish material cell element model.
Step 6：By the cell element model discretization in step 5.
Step 7：Equivalent structure square of cell element is calculated with finite element method to the model definition material attribute in step 6
Battle array.
Step 8：Assign the equivalent constitutive matrix in step 7 to macroscopical discrete topology in step 4, obtaining macrostructure has
Limit meta-model.
Step 9：Macrostructure finite element in solution procedure 8.
Step 10：Obtain the strain regime of all points of macromodel in 9.
Step 11：Using the point strain regime in step 10 as the boundary condition for carefully seeing cell element model in step 7, often
One point corresponds to thin sight finite element model.
Step 12：All thin sight finite element models in solution procedure 11.
Step 13：Micro-stress field in obtaining step 12, strain field solving result.
Step 14：Input is thin to see intensive parameter, and to micro-stress field in step 13, strain field is carefully to see criterion of strength judgement
Whether occur carefully to see failure.If all points do not fail, structure safety；If there is point at least one to fail,
Then structural failure.
A kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis of the present invention
The advantages of be：
(1) present invention considers multiscale analysis in oriented composite material structural reliability design, that is, compares traditional base
In outside the structural analysis of macro-scale, it is also contemplated that the thin sight characteristic of material so that model of structural mechanics possesses multiple rulers simultaneously
The information of degree, model are more acurrate.
(2) present invention carefully sees according to cell element corresponding with point and answers in oriented composite material structural reliability design
The field of force, strain field, and failure criteria is carefully seen to judge whether material fails, compared to the macroscopic failures criterion judged roughly, this hair
Bright judgment mode is more acurrate；In addition, macroscopic failures criterion is difficult to differentiate failure mode, the present invention fails with thin see
Criterion can clearly judge which kind of failure mode that material occurs, such as fibrous fracture, matrix destroy.
Description of the drawings
Fig. 1 is a kind of oriented fiber reinforced composite laminate structures Calculation of Reliability based on multiscale analysis of the present invention
The Caro method sampling process of the arrangement flow chart of steps of method, illiteracy is omitted in flow charts, and flow chart is calculated up to single sample
Process (step 3~14 in " invention content ")；
Fig. 2 is modeling of the macrostructure (by taking tablet as an example) in ABAQUS；
Fig. 3 is the oriented fiber reinforced composite cell element model established in ABAQUS；
Fig. 4 is the finite element model of cell element；
Fig. 5 is macrostructure finite element model.
Specific implementation mode
The invention reside in provide a kind of oriented fiber reinforced composite laminate structures reliability based on multiscale analysis
Computational methods.
Monte Carlo method sample calculation structural reliability is used, coupling multiscale analysis technical computing structure statics is rung
It answers, couples numerical algorithm of the multiscale analysis based on finite element, realize that meso-scale is passed to macro-scale information by homogenization
It passs, realizes macro-scale to carefully by regarding strain regime on macro-scale structures element integral point as thin cell element boundary condition of seeing
See the transmission of dimensional information；By ABAQUS Pathon secondary developing platforms, so that entire calculating process is realized in ABAQUS,
Sampling, modeling, calculating and statistics automation are achieved the effect that；The present invention has more reliable than structure of the tradition based on macro-scale
Property analysis, failure mode differentiate more acurrate, the more accurate advantage of result of calculation.
This method divides following steps：
Step 1：According to given parameter distribution rule in MATLAB, several groups random number is generated, wherein every group random
All uncertain parameters during number is calculated comprising single sample (fiber, matrix elastic property parameter, fiber, matrix intensity
The geometric parameter of parameter, macroscopical load and macrostructure), it is stored in sample.txt files in case use of sampling.
Step 2：With the secondary development function of Python levels in ABAQUS finite element softwares, Python programs are write,
Single sample calculates the one group of random number read in sample.txt files, writes Do statement and realizes that multiple sampling calculates, from
And the frequency of failure event generation is counted, estimate failure probability.Random number is generated in advance with this mode and is stored in file, it can
So that entire sample calculation is automatically performed in finite element software and (uses its secondary development function), counted in sampling without using person
It is manually operated during calculating.Below step 3~14 introduces the embodiment of single sample calculating, the program word of entire calculating process
Speech based on ABAQUS Python.
Step 3：With the secondary development function of Python levels in ABAQUS finite element softwares, with what is read in step 2
Macrostructure geometric parameter in random number establishes macrostructure geometrical model, with slab construction for such as Fig. 2.The present invention is all
Modeling process controls the script in ABAQUS to model using Python, is easy to implement sample calculation in this way and is not necessarily to
User operate, digital independent of the present invention, calculate etc. be by ABAQUS Python secondary development function in ABAQUS
Middle realization.
Step 4：Python programs are write, oriented fiber reinforced composite cell element model such as Fig. 3 is established in ABAQUS.
The present invention establishes three-dimensional cell element model, therefore can be with analyzing three-dimensional stress state, but compares two-dimentional cell element model, number of unit,
Interstitial content increases, and calculation amount increases when solving finite element equation.
Step 5：Cellular material attribute (fiber, matrix this structure are assigned with the material parameter in the random number read in step 2
Relationship), program is write after the completion, and structural separation is obtained into cell element finite element model such as Fig. 4.
Step 6：To the equivalent three dimensional elasticity performance of discrete cell element model prediction in step 5.Predict macro equivalent performance
Based on homogenization method, Equivalent Elasticity matrixWhereinFor cell element body mean stress, size is equal to
Cell element stress field cell element volume integral is averaged, likewise,For cell element body mean strain.
Step 7：The Effective Elastic Properties predicted in step 6 are given in a manner of definition material attribute in step 3 and are built
Vertical macrostructure model, after the completion by structural separation, such as Fig. 5.Structural loads, magnitude of load are added to be derived from step 2 and read
Macroscopical load random number.
Step 8：Calculate the macrostructure finite element model in step 7.
Step 9：In read step 8 in solving result ODB files all points strain information.
Step 10：Using the strain information in step 9 as the boundary condition of cell element, apply cell element finite element in steps of 5
On model boundary, each point is finally made to correspond to thin sight finite element model.
Step 11：Calculate the thin sight finite element model in step 10.
Step 12：Result of calculation ODB files in read step 11, storage micro-stress field, strain field result of calculation.
Step 13：Thin sight intensive parameter random number in read step 3, to micro-stress field in step 12, strain field with
Thin criterion of strength of seeing judges whether to occur carefully to see failure.If all points do not fail, structure safety；If having at least one
Point fails, then structural failure.
Step 14：Next group of random number is read, return to step 3 starts sample calculation next time, until having been calculated all groups
Random number.
Step 15：The frequency to fail is recorded with program, if sampling n times (n is sufficiently large) in total, m mistake has occurred
Effect, then structural realism isStructure safe probability R=1-Pf。
Bibliography
[1] reliability optimization of Wang Xiangyang laminated composite plate structures and the Central China robust optimizing research [D] science and technology are big
It learns, 2004.
[2] composite structure stable constraint optimization design [J] the composite materials of Wu Hao, Yan Ying based on reliability
Report, 2007,24 (5):149-153.
[3] reliability design approach of Ge Rui composite laminated structures and the Central China optimization algorithm research [D] science and technology are big
It learns, 2008.
[4] Tang Yuan richnesses composite structure reliability Robust Analysis and optimization algorithm research [D] the Central China University of Science and Technology,
2012.
[5] Zhang Jiumin, Gu Rujia, Li Hongshuan wait composite material connection structure fail-safe analyses of the based on response phase method
[J] Computational Mechanics journals, 2016,33 (5):711-716.
[6] Cui Sheng modest structural reliability designs [M] National Defense Industry Press, 2014.
[7] Qin Quan, forest-road brocade, the rigid structural reliabilities STOCHASTIC FINITE ELEMENT of plum:Theoretical and engineer application [M] Tsinghua University goes out
Version society, 2006.
[8] publishing house of the Zhang Zhimin structural mechanics of complex material [M] BJ University of Aeronautics ＆ Astronautics, 1993.
[9]Motley M R,Yin L Y.Influence of uncertainties on the response and
reliability of self-adaptive composite rotors[J].Composite Structures,2012,94
(1):114-120.
[10]Gosling P D,Faimun,Polit O.Ahigh-fidelity first-order reliability
analysis for shear deformable laminated composite plates[J].Composite
Structures,2014,115(18):12-28.
[11]Lee S,Kim I G,Cho W,et al.Advanced probabilistic design and
reliability-based design optimization for composite sandwich structure[J]
.Advanced Composite Materials,2014,23(1):3-16.
[12]Zhou X Y,Gosling P D,Ullah Z,et al.Exploiting the benefits of
multi-scale analysis in reliability analysis for composite structures[J]
.Composite Structures,2016,155:197-212.
[13]Zhou X Y,Gosling P D,Ullah Z,et al.Stochastic Multi-scale Finite
Element based Reliability Analysis for Laminated Composite Structures[J]
.Applied Mathematical Modelling,2017.
Claims (5)
1. a kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis, feature exist
In：Use Monte Carlo method sample calculation structural reliability, coupling multiscale analysis technical computing structure statics response, coupling
Numerical algorithm of the multiscale analysis based on finite element is closed, realizes that meso-scale is transmitted to macro-scale information by homogenization, leads to
It crosses and realizes that macro-scale sees ruler to thin using strain regime on macro-scale structures element integral point as thin cell element boundary condition of seeing
Spend the transmission of information；By ABAQUS Pathon secondary developing platforms, so that entire calculating process is realized in ABAQUS, reach
The effect of sampling, modeling and calculating process automation.
2. a kind of oriented fiber reinforced composite laminate structures reliability meter based on multiscale analysis as described in claim 1
Calculation method, it is characterised in that：Calculating process be based on ABAQUS Pathon secondary developing platforms, so that entire calculating process is existed
Realized in ABAQUS, sample, model and calculating process automation.
3. the as described in claim a kind of oriented fiber reinforced composite laminate structures Calculation of Reliability based on multiscale analysis
Method, it is characterised in that：Cell element model boundary condition is seen using strain regime at point as thin, solves micro-stress, strain
.
4. the as described in claim a kind of oriented fiber reinforced composite laminate structures Calculation of Reliability based on multiscale analysis
Method, it is characterised in that：It is as structural response judge structure using cell element micro-stress, strain field at meso-scale upper integral point
No failure.
5. the as described in claim a kind of oriented fiber reinforced composite laminate structures Calculation of Reliability based on multiscale analysis
Method, it is characterised in that：Differentiated respectively to whether fiber and matrix fail with thin failure criteria of seeing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810213287.4A CN108388763A (en) | 2018-03-15 | 2018-03-15 | A kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810213287.4A CN108388763A (en) | 2018-03-15 | 2018-03-15 | A kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108388763A true CN108388763A (en) | 2018-08-10 |
Family
ID=63067734
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810213287.4A Pending CN108388763A (en) | 2018-03-15 | 2018-03-15 | A kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108388763A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112733341A (en) * | 2020-12-30 | 2021-04-30 | 中国海洋大学 | Progressive failure analysis method for marine fiber reinforced flexible pipe |
| CN113221408A (en) * | 2021-05-06 | 2021-08-06 | 浙江工业大学 | Parallel NTFA theory model order reduction calculation method based on ABAQUS |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105183990A (en) * | 2015-09-07 | 2015-12-23 | 哈尔滨工业大学 | Composite material and method for multiscale response analysis of structure thereof |
| CN106066913A (en) * | 2016-05-31 | 2016-11-02 | 西北工业大学 | Complex composite material structure equivalent material performance multi-dimension computational methods |
| CN107239642A (en) * | 2017-07-10 | 2017-10-10 | 南京航空航天大学 | Towards the multiple dimensioned numerical value emulation method of carbon fibre reinforced composite drilling |
| CN107451307A (en) * | 2016-05-31 | 2017-12-08 | 西北工业大学 | A kind of method of Multi-Scale Calculation complex composite material structure effective stiffness matrix |
-
2018
- 2018-03-15 CN CN201810213287.4A patent/CN108388763A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105183990A (en) * | 2015-09-07 | 2015-12-23 | 哈尔滨工业大学 | Composite material and method for multiscale response analysis of structure thereof |
| CN106066913A (en) * | 2016-05-31 | 2016-11-02 | 西北工业大学 | Complex composite material structure equivalent material performance multi-dimension computational methods |
| CN107451307A (en) * | 2016-05-31 | 2017-12-08 | 西北工业大学 | A kind of method of Multi-Scale Calculation complex composite material structure effective stiffness matrix |
| CN107239642A (en) * | 2017-07-10 | 2017-10-10 | 南京航空航天大学 | Towards the multiple dimensioned numerical value emulation method of carbon fibre reinforced composite drilling |
Non-Patent Citations (4)
| Title |
|---|
| X.-Y. ZHOU等: ""Exploiting the benefits of multi-scale analysis in reliability analysis for composite structures"", 《COMPOSITE STRUCTURES》 * |
| X.-Y. ZHOU等: ""Stochastic multi-scale finite element based reliability analysis for laminated composite structures"", 《APPLIED MATHEMATICAL MODELLING》 * |
| 张兴龙: ""纤维增强复合材料结构宏细观多尺度力学性能研究及应用"", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
| 陈滨琦: ""基于多尺度方法的复合材料层合板结构失效机理研究"", 《中国博士学位论文全文数据库 工程科技II辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112733341A (en) * | 2020-12-30 | 2021-04-30 | 中国海洋大学 | Progressive failure analysis method for marine fiber reinforced flexible pipe |
| CN112733341B (en) * | 2020-12-30 | 2022-06-07 | 中国海洋大学 | Progressive failure analysis method for marine fiber reinforced flexible pipe |
| CN113221408A (en) * | 2021-05-06 | 2021-08-06 | 浙江工业大学 | Parallel NTFA theory model order reduction calculation method based on ABAQUS |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Global optimization of reliability design for large ball mill gear transmission based on the Kriging model and genetic algorithm | |
| Chiachio et al. | Reliability in composites–A selective review and survey of current development | |
| Gao et al. | Reliability-based low-cycle fatigue damage analysis for turbine blade with thermo-structural interaction | |
| Wei et al. | Modified hybrid forecast model considering chaotic residual errors for dam deformation | |
| Li et al. | Reliability-based fatigue life investigation for a medium-scale composite hydrokinetic turbine blade | |
| CN104112069B (en) | Method for pre-estimating anisotropic heat conduction coefficient of fiber toughening composite material on basis of microstructural image recognition | |
| Bacharoudis et al. | A probabilistic approach for strength and stability evaluation of wind turbine rotor blades in ultimate loading | |
| Lynch et al. | Unstructured overset incompressible computational fluid dynamics for unsteady wind turbine simulations | |
| Liu et al. | Strength reliability analysis of aluminium–carbon fiber/epoxy composite laminates | |
| CN108388763A (en) | A kind of oriented fiber reinforced composite laminate structures reliability calculation method based on multiscale analysis | |
| CN102682175B (en) | Method for analyzing reliability of construction error of grid structure based on buckling mode combination | |
| Su et al. | Structure analysis for concrete‐faced rockfill dams based on information entropy theory and finite element method | |
| Grover et al. | Influence of parametric uncertainties on the deflection statistics of general laminated composite and sandwich plates | |
| Zhang et al. | Resistance uncertainty and structural reliability of hypar tensioned membrane structures with PVC coated polyesters | |
| Shiri et al. | Effect of strength dispersion on fatigue life prediction of composites under two-stage loading | |
| CN113836760B (en) | Turbine disk creep fatigue life reliability assessment method | |
| Jia et al. | Validation and development of trace-based approach for composite laminates | |
| Zheng et al. | Complexity of software trustworthiness and its dynamical statistical analysis methods | |
| CN105488343A (en) | Fault probability calculation method of secondary electric power equipment | |
| Bacharoudis et al. | Estimating design reliability of composite rotor blades under ultimate loading | |
| CN111985121B (en) | Method and device for converting mechanism model and pole plate structure model machine model | |
| Qiu et al. | Reliability study of fracture mechanics based non‐probabilistic interval analysis model | |
| Fanxiao et al. | Study on structural characteristics of composite smart grille based on principal component analysis | |
| Liu et al. | Analysis of time-dependent reliability of concrete gravity dam damage based on gray theory | |
| CN106022630A (en) | Operation safety assessment method for liquid oxygen filling system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180810 |