CN103605865A - Delamination propagation behavior simulation method of composite material multidirectional laminated plate on basis of cohesion model - Google Patents

Delamination propagation behavior simulation method of composite material multidirectional laminated plate on basis of cohesion model Download PDF

Info

Publication number
CN103605865A
CN103605865A CN201310627964.4A CN201310627964A CN103605865A CN 103605865 A CN103605865 A CN 103605865A CN 201310627964 A CN201310627964 A CN 201310627964A CN 103605865 A CN103605865 A CN 103605865A
Authority
CN
China
Prior art keywords
model
multidirectional
load
compound substance
finite element
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
Application number
CN201310627964.4A
Other languages
Chinese (zh)
Other versions
CN103605865B (en
Inventor
赵丽滨
龚愉
辛安
山美娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201310627964.4A priority Critical patent/CN103605865B/en
Publication of CN103605865A publication Critical patent/CN103605865A/en
Application granted granted Critical
Publication of CN103605865B publication Critical patent/CN103605865B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a delamination propagation behavior simulation method of a composite material multidirectional laminated plate on the basis of a cohesion model. The method includes the following steps that (1) a geometric model is built according to structural parameters of a composite material multidirectional laminated plate test piece, (2) key parameters capable of reflecting interfacial behaviors are calculated, and material attributes are respectively set, (3) meshing is carried out on the geometric model of the composite material multidirectional laminated plate test piece, a module is assembled, and a three-dimensional finite element model is built, (4) a load and boundary conditions of the finite element model are determined according to practical situations and the load state of a composite material structure, (5) the finite element model based on the cohesion model is calculated and analyzed, a load displacement curve of a loading point is extracted, delamination propagation behaviors of the composite material multidirectional laminated plate are simulated to obtain the largest load value of the load displacement curve, and damage behaviors of the multidirectional laminated plate are predicted.

Description

The multidirectional laying plate of a kind of compound substance based on cohesive zone model delamination Behavior modeling method
Technical field
The present invention relates to the multidirectional laying plate of a kind of compound substance based on cohesive zone model delamination Behavior modeling method, be applicable to research and the prediction of the behavior of the multidirectional laying plate of the widely used compound substance in the fields such as Aero-Space delamination.
Background technology
Advanced composite material is widely used at present in aerospace structure, and has become aerospace four macrostructure materials together with aluminium alloy, titanium alloy, alloy steel.Wherein, carbon fiber resin matrix composite, when having compared with high specific strength, has specific stiffness concurrently high, anticorrosive, fatigue behaviour good and the plurality of advantages such as performance designability and become the main a kind of compound substance using in aircaft configuration in modern age.Along with the batch production of middle modulus high-strength carbon fiber and the application of toughened resin, this class compound substance application in aircaft configuration is also developed into the main force support structure of force-bearing situation complexity gradually by inferior load-carrying construction.And the safety of structure of compound substance is also more and more valued by the people.
Different from metal material, the failure mode of compound substance and failure mechanism are more complicated various.Due to component material---the greatest differences of fiber and matrix mechanical property, the interlayer performance of composite laminated plate is far below performance in its face, in different loads situation, the layering that low velocity impact or geometric discontinuity cause is expanded under fatigue load, finally causing the flexing of sublayer plate, is a kind of main damage form of practice of composite laminate structures.Delamination behavior meeting causes strength and stiffness to incur loss, for integrality and safety period of structure be all extremely disadvantageous, the calamitous structural failure that even may cause outside to be difficult for discovering.Therefore, in design and analysis process, to the prediction of compound substance delamination behavior, be, a very important job.
DCB test and MMB test are often used in the fracture toughness under research I type and I/II mixed type loaded-up condition, and corresponding compound substance unidirectional laminate test method (ASTM D5528-01 and ASTM D6671/D6671M-06) is comparative maturity.In generation and the build phase of delamination, cohesive zone model (CZM) is widely applied.Passed through a large amount of research work, adopting linear this to build vertical cohesive zone model can well analog composite material unidirectional laminate crack initiation and growth, and gained predicts the outcome and has test findings and have good consistance.
Compare with the unidirectional laying plate of compound substance, the application of the multidirectional laying plate of compound substance in engineering practice is more extensive.On the different laying interface of the multidirectional laying plate of compound substance machine direction, be easy to the delamination that forms crackle and crackle occurs.Yet the work that composite layer intersexuality the subject of knowledge and the object of knowledge carries out at present mostly also concentrates in the research of compound substance single-ply plate, for rare the relating to of delamination rule at the multidirectional laying of compound substance interface.Due to very complicated with the failure mechanism of interlayer in the multidirectional laying flaggy of compound substance, there is fiber bridge joint, wear the various of failure such as layer expansion, and document that can be for reference is very limited, for the numerical simulation of the multidirectional laying plate of compound substance delamination, remain at present one and have much challenging work.In order to carry out in this respect new research, to predict the damage expansion in engineering structure, in performance material property advantage, ensure structural safety, need to introduce improved cohesive zone model, reflect better the failure procedure of the multidirectional laying plate of compound substance.
Summary of the invention
Technology of the present invention is dealt with problems: overcome the deficiencies in the prior art, the multidirectional laying plate of a kind of compound substance based on cohesive zone model delamination Behavior modeling method is provided, be applicable to engineering application, can assist multidirectional laying plate laying sequence design, reduce experimentation cost, the multidirectional laying plate of compound substance inefficacy damage mode is carried out to effective analysis and prediction simultaneously, ensure better safety of structure.
Technical solution of the present invention, the multidirectional laying plate of a kind of compound substance based on cohesive zone model delamination Behavior modeling method, performing step is as follows:
Steps A, utilizes the d solid modeling function of finite element software Abaqus, according to the structural parameters of the multidirectional laying plate of compound substance, sets up three geometric models that the multidirectional paving layer test coupon of compound substance is corresponding;
Step B, selects the linear constitutive relation of cohesive zone model, and the critical strain energy release rate that cohesive zone related interfaces parameter and test record is set, and model is arranged respectively to laying direction and material properties;
Step C, to three geometric models difference grid divisions, composition model, sets up three-dimensional finite element model;
Step D, according to the actual conditions of composite structure and loaded-up condition, force analysis situation, determines load, the boundary condition of test specimen, and is applied on finite element model;
Step e, finite element model based on cohesive zone model is carried out to stress analysis, obtain result of calculation, extract the load displacement curve of load(ing) point, the behavior of the multidirectional laying plate of the compound substance of simulation based on cohesive zone model delamination, obtain load displacement curve maximum load value, predict the derogatory behaviour of multidirectional laying plate.
Described steps A is utilized the d solid modeling function of finite element software Abaqus, and the detailed process of setting up geometric model is:
(A1) open Abaqus software, utilize d solid modeling function, according to the multidirectional paving layer test coupon of compound substance values of the structural parameters, in two dimensional surface, draw three geometric profile figure corresponding to test specimen upper, middle, and lower part;
(A2) by entity stretching mode, set up respectively multidirectional 3-D geometric model, wherein upper and lower two parts through-thickness is stretched to half of actual (real) thickness value, and the center section through-thickness at cohesive zone place is stretched to 0.01mm.
Described step B determines the key parameter of reflection interface behavior, and the specific implementation process that material properties is set is:
(B1) definition local material direction, divides upper and lower two parts geometric model according to each laying one-tenth-value thickness 1/10, creates compound substance laying and defines laying angle, when definition laying angle, need to set up corresponding local coordinate system according to actual conditions;
(B2), according to material property parameter, define the material properties of upper and lower two parts geometric model;
(B3) select the linear constitutive relation of cohesive zone model, by subroutine, introduce the critical strain energy release rate changing with crack length that test obtains, calculate initial interface rigidity, boundary strength, the interface parameter of the reflection such as coefficient of viscosity interface behavior, the material properties of definition cohesive zone place center section geometric model.
Described step C carries out respectively grid division to the geometric model of test specimen, composition model, and the detailed process of setting up three-dimensional finite element model is:
(C1) to upper and lower two parts along three direction Even seedings of length, and grid division, this grid does not need very meticulous;
(C2) calculate cohesive zone length, determine cohesion unit size, quantity, according to cohesion unit size and number of unit, to center section Even seeding, divide comparatively meticulous grid;
(C3) three geometric models are assembled successively according to the order from upper, middle and lower, at two assembly and connection interfaces, binding constraints is set, in addition for preventing that laminate upper and lower part from running through, between the two sides, upper and lower part at precrack place, set up contact, set up three-dimensional finite element model.
Described step D determines that geometric model load, boundary condition implementation procedure are:
(D1) judgement I type or I/II mixed type loaded-up condition, carry out force analysis to the multidirectional paving layer test coupon of compound substance, determines restraint condition and the suffered load situation of test specimen;
(D2) using restraint condition above as boundary condition, be arranged on the finite element model of load, on test specimen finite element model, apply the load of force analysis gained simultaneously.
Described step e is carried out computational analysis to cohesive zone model, obtains result of calculation, extracts the load displacement curve of load(ing) point, to the specific implementation process of the multidirectional laying plate of compound substance delamination Behavior modeling, is:
(E1) finite element model is carried out to stress analysis, calculate the strain energy rate of cohesion unit, by reading cohesion unit coordinate information again according to critical strain energy release rate-crack length relation curve, determine the critical strain energy release rate of cohesion unit, by currency and critical value comparison, the spread scenarios of judgement cohesion element failure and crackle;
(E2) obtain result of finite element, utilize the post-processing function of finite element software, extract the load displacement curve of load(ing) point, the behavior of the multidirectional laying plate of analog composite material delamination, obtain load displacement curve maximum load value, when real load reaches maximal value, the expansion of crackle does not need to continue to increase load, thereby realizes the prediction to the behavior of multidirectional laying dash-board injury.
The present invention's advantage is compared with prior art:
(1) the present invention is based on the multidirectional laying plate of the compound substance delamination Behavior modeling method of cohesive zone model, overcome the deficiency of studying for unidirectional laying plate, consider the impact of different fiber directions on the behavior of multidirectional laying plate delamination, the delamination behavior of the multidirectional laying plate of analog composite material, thereby the effectively derogatory behaviour of the multidirectional laying plate of analysis and prediction.
(2) the present invention is based on the multidirectional laying plate of the compound substance delamination Behavior modeling method of cohesive zone model, can effectively analyze and simulate the delamination behavior of the multidirectional laying plate of compound substance at multiple interface, according to Numerical Simulation Results, the design of auxiliary multidirectional laying plate laying direction.
(3) the present invention is based on the multidirectional laying plate of the compound substance delamination Behavior modeling method of cohesive zone model, only need set up cohesive zone model, carry out finite element analysis, can carry out numerical simulation to the behavior of the multidirectional laying plate of compound substance delamination, be applicable to prediction and the analysis of the multidirectional laying plate of the compound substance delamination failure behaviour in engineering application, can significantly shorten the lead time, reduce experimentation cost.
Accompanying drawing explanation
Fig. 1 is realization flow figure of the present invention;
Fig. 2 is three geometric model schematic diagram that the multidirectional paving layer test coupon of compound substance is corresponding;
Fig. 3 is assembling geometric model schematic diagram;
Fig. 4 is the finite element method (fem) analysis process flow diagram based on cohesive zone model.
Embodiment
As shown in Figure 1, the inventive method specific implementation is:
1,, according to the multidirectional paving layer test coupon of compound substance values of the structural parameters, in Abaqus software, set up the geometric model of test specimen.
Open Abaqus software, utilize d solid modeling function, according to the multidirectional paving layer test coupon of compound substance values of the structural parameters, in two dimensional surface, draw respectively the geometric profile figure that test specimen is corresponding, by entity stretching mode, set up respectively 3-D geometric model, as shown in Figure 2, a wherein, c two parts through-thickness is stretched to half of actual (real) thickness value, b part through-thickness stretching 0.01mm.
2, the key parameter of determining reflection interface behavior, arranges laying angle and material properties.
Definition local material direction, according to each laying one-tenth-value thickness 1/10, divide the geometric model of a and c, create compound substance laying and define laying angle, when definition laying angle, need to set up corresponding local coordinate system according to actual conditions, according to material property parameter, the material properties of definition a and c geometric model, select the linear constitutive relation of cohesive zone model, introduce the critical strain energy release rate changing with crack length that test obtains, calculate initial interface rigidity, boundary strength, the parameter of the reflection such as coefficient of viscosity interface behavior, the material properties that definition b model is corresponding cohesive zone.
3, geometric model is carried out respectively to grid division, composition model, sets up three-dimensional finite element model.
To a, c two parts are along three direction Even seedings of length, and grid division, this grid does not need very meticulous, calculate cohesive zone length, determine cohesion unit size, quantity, according to cohesion unit size and number of unit, to b part Even seeding, divide comparatively meticulous grid, by geometric model according to from a, b, the order of c is assembled successively, assembly and connection interface at geometric model a/b and b/c arranges binding constraints, at precrack place, set up contact, as shown in Figure 3, set up three-dimensional finite element model.
4, determine geometric model load, boundary condition.
Judgement I type or I/II mixed type loaded-up condition, the multidirectional paving layer test coupon of compound substance is carried out to force analysis, determine restraint condition and the suffered load situation of test specimen, using restraint condition above as boundary condition, be arranged on the finite element model of load, on test specimen finite element model, apply the load of force analysis gained simultaneously.
5, the multidirectional laying plate of analog composite material delamination behavior.
Finite element model is carried out to stress analysis, calculate the strain energy rate of correspondence position, obtain crack length, again according to critical strain energy release rate-crack length relation curve, determine the critical energy release rate of this position, by currency and critical value comparison, the spread scenarios of judgement cohesion element failure and crackle, idiographic flow as shown in Figure 4, obtain result of finite element, utilize the post-processing function of finite element software, extract the load displacement curve of load(ing) point, the behavior of the multidirectional laying plate of the compound substance of simulation based on cohesive zone model delamination, obtain load displacement curve maximum load value, judgement is when real load reaches maximal value, the expansion of crackle does not need to continue to increase load, thereby realize the prediction to the behavior of multidirectional laying dash-board injury.
Non-elaborated part of the present invention belongs to those skilled in the art's known technology.

Claims (7)

1. the multidirectional laying plate of the compound substance based on a cohesive zone model delamination Behavior modeling method, is characterized in that comprising the following steps:
Steps A, utilizes the d solid modeling function of finite element software Abaqus, according to the structural parameters of the multidirectional laying plate of compound substance, sets up three geometric models that the multidirectional paving layer test coupon of compound substance is corresponding;
Step B, selects the linear constitutive relation of cohesive zone model, and the critical strain energy release rate that cohesive zone related interfaces parameter and test record is set, and model is arranged respectively to laying direction and material properties;
Step C, to three geometric models difference grid divisions, composition model, sets up three-dimensional finite element model;
Step D, according to the actual conditions of composite structure and loaded-up condition, force analysis situation, determines load, the boundary condition of test specimen, and is applied on finite element model;
Step e, finite element model based on cohesive zone model is carried out to stress analysis, obtain result of calculation, utilize finite element post-processing function, extract the load displacement curve of load(ing) point, the behavior of the multidirectional laying plate of the compound substance of simulation based on cohesive zone model delamination, obtains load displacement curve maximum load value, predicts the behavior of multidirectional laying dash-board injury.
2. the multidirectional laying plate of a kind of compound substance based on cohesive zone model according to claim 1 delamination Behavior modeling method, is characterized in that: the geometric model implementation procedure of setting up test specimen according to the structural parameters of the multidirectional paving layer test coupon of compound substance in described steps A is:
(A1) open Abaqus software, utilize d solid modeling function, according to the multidirectional paving layer test coupon of compound substance values of the structural parameters, in two dimensional surface, draw three geometric profile figure corresponding to test specimen upper, middle, and lower part;
(A2) by entity stretching mode, set up respectively multidirectional 3-D geometric model, wherein upper and lower two parts through-thickness is stretched to half of actual (real) thickness value, and the center section through-thickness at cohesive zone place is stretched to 0.01mm.
3. the multidirectional laying plate of a kind of compound substance based on cohesive zone model according to claim 1 delamination Behavior modeling method, it is characterized in that: described step B determines the key parameter of reflection interface behavior, and the specific implementation process that laying direction and material properties are set is:
(B1) definition local material direction, divides upper and lower two parts geometric model according to each laying one-tenth-value thickness 1/10, creates compound substance laying and defines laying angle, when definition laying angle, need to set up corresponding local coordinate system according to actual conditions;
(B2), according to material property parameter, define the material properties of upper and lower two parts geometric model;
(B3) select the linear constitutive relation of cohesive zone model, by subroutine, introduce the critical strain energy release rate changing with crack length that test obtains, calculate initial interface rigidity, boundary strength, the interface parameter of the reflection such as coefficient of viscosity interface behavior, the material properties of definition cohesive zone place center section geometric model.
4. the multidirectional laying plate of a kind of compound substance based on cohesive zone model according to claim 1 delamination Behavior modeling method, it is characterized in that: described step C carries out respectively grid division to the geometric model of test specimen, composition model, the detailed process of setting up three-dimensional finite element model is:
(C1) to upper and lower two parts along three direction Even seedings of length, and grid division, this grid does not need very meticulous;
(C2) calculate cohesive zone length, determine cohesion unit size, quantity, according to cohesion unit size and number of unit, to center section Even seeding, divide comparatively meticulous grid;
(C3) three geometric models are assembled successively according to the order from upper, middle and lower, at two assembly and connection interfaces, binding constraints is set, in addition for preventing that laminate upper and lower part from running through, between the two sides, upper and lower part at precrack place, set up contact, set up three-dimensional finite element model.
5. the multidirectional laying plate of a kind of compound substance based on cohesive zone model according to claim 1 delamination Behavior modeling method, is characterized in that: described step D determines that geometric model load, boundary condition implementation procedure are:
(D1) judgement I type or I/II mixed type loaded-up condition, carry out force analysis to the multidirectional paving layer test coupon of compound substance, determines restraint condition and the suffered load of test specimen;
(D2) using restraint condition above as boundary condition, be arranged on the finite element model of load, on test specimen finite element model, apply the load of force analysis gained simultaneously.
6. the multidirectional laying plate of a kind of compound substance based on cohesive zone model according to claim 1 delamination Behavior modeling method, it is characterized in that: described step e is carried out computational analysis to cohesive zone model, obtain result of calculation, the load displacement curve that extracts load(ing) point, to the specific implementation process of the multidirectional laying plate of compound substance delamination Behavior modeling is:
(E1) finite element model is carried out to stress analysis, calculate the strain energy rate of cohesion unit, by reading cohesion unit coordinate information again according to critical strain energy release rate-crack length relation curve, determine the critical strain energy release rate of cohesion unit, by currency and critical value comparison, the spread scenarios of judgement cohesion element failure and crackle;
(E2) obtain result of finite element, utilize the post-processing function of finite element software, extract the load displacement curve of load(ing) point, the behavior of the multidirectional laying plate of analog composite material delamination, obtain load displacement curve maximum load value, when real load reaches maximal value, the expansion of crackle does not need to continue to increase load, thereby realizes the prediction to the behavior of multidirectional laying dash-board injury.
7. the multidirectional paving layer test coupon of compound substance, it is characterized in that: the multidirectional paving layer test coupon of this compound substance is according to obtaining after the multidirectional laying plate of a kind of compound substance based on the cohesive zone model delamination Behavior modeling methods analyst described in claim 1 to 6 any one, overcome the deficiency of studying for unidirectional laying plate, consider the impact of different fiber directions on the behavior of multidirectional laying plate delamination, obtain load displacement curve maximum load value, thereby effectively predict the derogatory behaviour of multidirectional laying plate, according to Numerical results, can also assist the design of multidirectional laying plate laying order.
CN201310627964.4A 2013-11-28 2013-11-28 A kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method Expired - Fee Related CN103605865B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310627964.4A CN103605865B (en) 2013-11-28 2013-11-28 A kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310627964.4A CN103605865B (en) 2013-11-28 2013-11-28 A kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method

Publications (2)

Publication Number Publication Date
CN103605865A true CN103605865A (en) 2014-02-26
CN103605865B CN103605865B (en) 2016-08-17

Family

ID=50124086

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310627964.4A Expired - Fee Related CN103605865B (en) 2013-11-28 2013-11-28 A kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method

Country Status (1)

Country Link
CN (1) CN103605865B (en)

Cited By (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105335575A (en) * 2015-11-25 2016-02-17 中国航空工业集团公司沈阳飞机设计研究所 Finite element modeling method of electronic antenna structural integrated laminated board
CN106897483A (en) * 2017-01-08 2017-06-27 浙江大学 A kind of method that hail model is set up based on insertion cohesive force between any Unit two
CN106934127A (en) * 2017-02-28 2017-07-07 江铃汽车股份有限公司 The stress acquisition methods of the metallic framework of powertrain mounting system
CN107016156A (en) * 2016-01-28 2017-08-04 波音公司 The method that test data for failure mode in composite layer simplifies
CN107077517A (en) * 2014-09-16 2017-08-18 株式会社Ihi Design assistant device, design aiding method and program
CN107256320A (en) * 2017-07-28 2017-10-17 中科机器人科技有限公司 A kind of small-sized hand based on composite plys throws the finite element method of unmanned plane
CN107330202A (en) * 2017-07-05 2017-11-07 北京航空航天大学 A kind of composite I types layering bridging rule based on virtual hierarchy expanding test determines method
CN107357973A (en) * 2017-06-26 2017-11-17 北京航空航天大学 A kind of Forecasting Methodology of any interface layer I types delamination fractures toughness stationary value of CFRP layer plate
CN107356478A (en) * 2017-07-07 2017-11-17 吉林大学 Long fiber reinforced thermoplastic composites Multiscale Simulation Method
CN107451360A (en) * 2017-07-28 2017-12-08 中国航空工业集团公司西安飞机设计研究所 The method for numerical simulation of the continuous laying of composite right-angled bend
CN108256281A (en) * 2018-03-26 2018-07-06 中国矿业大学 A kind of intensity prediction method for considering overlap joint interface topography and overlapping object graded properties
CN108268692A (en) * 2017-11-21 2018-07-10 北京航空航天大学 Automatic fiber placement composite material curing deformation forecasting procedure
CN108595754A (en) * 2018-03-20 2018-09-28 南京航空航天大学 The emulation mode of interlayer toughened composite laminated plate
CN108595724A (en) * 2017-12-21 2018-09-28 刘长喜 Composite material revolving meber design method
CN108681650A (en) * 2018-05-29 2018-10-19 内蒙古工业大学 A kind of warp-knitted spacer fabric composite material stress analysis method
CN108804787A (en) * 2018-05-26 2018-11-13 南京理工大学 The method for being inserted into cohesive force unit simulation Rock Bridge Failure based on batch
CN108984809A (en) * 2017-05-31 2018-12-11 天津工业大学 A kind of method of the new raising tensile capacity of composite laminated plate containing eccentric orfice
CN109190258A (en) * 2018-09-06 2019-01-11 中国航发北京航空材料研究院 A kind of composite-material laminates I type method for layering simulation based on the superposition of cohesion unit
CN109190225A (en) * 2018-08-24 2019-01-11 湘潭大学 A kind of finite element method for simulating hard coat underbead crack and Interface Crack
CN110032765A (en) * 2019-03-11 2019-07-19 武汉科技大学 A kind of optimization method, system, device and the medium of refractory material fragmentation parameters
CN110047563A (en) * 2019-03-19 2019-07-23 首钢集团有限公司 A kind of composite material lightning damage analogy method and device
CN110222420A (en) * 2019-06-06 2019-09-10 吉林大学 A kind of failure judgment method of viscous riveting connection structure
CN110376055A (en) * 2019-06-20 2019-10-25 重庆大学 A kind of CFRP layer plate layering failure behaviour prediction technique based on novel cohesive force constitutive relation
CN110457734A (en) * 2019-06-05 2019-11-15 东南大学 One kind being bolted FRP section components fiber laying laying method and FRP section components
CN110826284A (en) * 2019-11-21 2020-02-21 北京航空航天大学 Modeling and analyzing method for composite material laminated board of interweaved and laminated mixed laying layers
CN110968916A (en) * 2018-09-28 2020-04-07 天津大学 Seabed metallurgy composite pipe die generator containing circumferential surface crack V-shaped groove welding seam and application thereof
CN111104751A (en) * 2019-12-24 2020-05-05 肇庆市海特复合材料技术研究院 Strength analysis method of composite material connecting rod
CN111553104A (en) * 2020-05-08 2020-08-18 东南大学 Simulation analysis method for bridge deck pavement structure
CN112487557A (en) * 2020-12-03 2021-03-12 上海交通大学 Method for predicting interface failure and microscopic crack propagation of composite material under hydraulic permeation load
CN112595593A (en) * 2020-12-29 2021-04-02 哈尔滨工业大学 Method for obtaining interface strength of three-dimensional woven composite material
CN112699585A (en) * 2020-12-29 2021-04-23 中国航空工业集团公司西安飞机设计研究所 Finite element modeling method for composite thick-layer laminated plate joint
CN112733324A (en) * 2020-12-17 2021-04-30 重庆大学 Composite material laminate fatigue layering expansion behavior prediction method considering stress ratio and fiber bridging influence
CN113049361A (en) * 2019-12-27 2021-06-29 中国航空工业集团公司西安飞机设计研究所 Laminated board III type crack delamination failure test method and test piece
CN113128097A (en) * 2021-04-29 2021-07-16 浙江理工大学 Method for simulating and predicting heat transfer performance of porous nanofiber medium
CN113158502A (en) * 2021-01-19 2021-07-23 西北工业大学 Reliability analysis method and device, storage medium and electronic equipment
CN113971354A (en) * 2021-10-15 2022-01-25 厦门大学 Finite element simulation method for nano composite material bending fracture damage research
CN114201814A (en) * 2021-11-02 2022-03-18 南京航空航天大学 Civil transport aircraft windshield collision simulation modeling method based on cohesion model
CN114357823A (en) * 2021-12-08 2022-04-15 湖北工业大学 Curved surface moving load analysis method based on ABAQUS
CN114896763A (en) * 2022-04-18 2022-08-12 重庆大学 Method for predicting composite material I/II mixed fracture response surface
CN117524361A (en) * 2023-10-24 2024-02-06 天津大学 Fatigue crack initiation prediction method for cycle cohesive force model with variable accumulation length

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107391783B (en) * 2017-06-02 2020-04-28 西安交通大学 Finite element theory-based asymmetric combined cycle unit shafting vibration monitoring method
CN107832492B (en) * 2017-10-11 2021-07-06 大连理工大学 Steel structure corrosion fatigue damage calculation method based on cohesion model

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101613944A (en) * 2009-07-10 2009-12-30 天津工业大学 A kind of method that is used to improve shock resistance of resin-based carbon fiber composite material
US20120034089A1 (en) * 2010-08-06 2012-02-09 Rolls-Royce Plc Composite material and method
CN103345545A (en) * 2013-06-14 2013-10-09 北京航空航天大学 Method for predicating composite material Pi-shaped non-planar glue joint strength based on triangular envelopes

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101613944A (en) * 2009-07-10 2009-12-30 天津工业大学 A kind of method that is used to improve shock resistance of resin-based carbon fiber composite material
US20120034089A1 (en) * 2010-08-06 2012-02-09 Rolls-Royce Plc Composite material and method
CN103345545A (en) * 2013-06-14 2013-10-09 北京航空航天大学 Method for predicating composite material Pi-shaped non-planar glue joint strength based on triangular envelopes

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
J.CHEN: "Simulation of multi-directional crack growth in braided composite T-piece specimens using cohesive models", 《FATIGUE & FRACTURE OF ENGINEERING MATERIALS & STRUCTURES》 *
JIANYU ZHANG等: "Fatigue delamination growth rates and thresholds of composite laminates under mixed mode loading", 《INTERNATIONAL JOURNAL OF FATIGUE》 *
苏罗丹等: "碳纤维复合材料层合板损伤脱层过程的有限元分析", 《第17届全国复合材料学术会议论文(复合材料力学分论坛)集》 *

Cited By (62)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107077517A (en) * 2014-09-16 2017-08-18 株式会社Ihi Design assistant device, design aiding method and program
CN107077517B (en) * 2014-09-16 2020-09-01 株式会社Ihi Design support device, design support method, and program
CN105335575A (en) * 2015-11-25 2016-02-17 中国航空工业集团公司沈阳飞机设计研究所 Finite element modeling method of electronic antenna structural integrated laminated board
CN105335575B (en) * 2015-11-25 2019-06-28 中国航空工业集团公司沈阳飞机设计研究所 A kind of finite element modeling method of the structure-integrated laminate of electronic antenna
CN107016156A (en) * 2016-01-28 2017-08-04 波音公司 The method that test data for failure mode in composite layer simplifies
CN106897483A (en) * 2017-01-08 2017-06-27 浙江大学 A kind of method that hail model is set up based on insertion cohesive force between any Unit two
CN106897483B (en) * 2017-01-08 2019-05-17 浙江大学 A method of hail model is established based on cohesive force is inserted between any Unit two
CN106934127A (en) * 2017-02-28 2017-07-07 江铃汽车股份有限公司 The stress acquisition methods of the metallic framework of powertrain mounting system
CN106934127B (en) * 2017-02-28 2020-06-23 江铃汽车股份有限公司 Method for acquiring stress state of metal framework of power assembly suspension system
CN108984809A (en) * 2017-05-31 2018-12-11 天津工业大学 A kind of method of the new raising tensile capacity of composite laminated plate containing eccentric orfice
CN107357973B (en) * 2017-06-26 2020-08-04 北京航空航天大学 Prediction method for I-type layered fracture toughness stable value of any interlayer interface of CFRP (carbon fiber reinforced polymer) laminate
CN107357973A (en) * 2017-06-26 2017-11-17 北京航空航天大学 A kind of Forecasting Methodology of any interface layer I types delamination fractures toughness stationary value of CFRP layer plate
CN107330202A (en) * 2017-07-05 2017-11-07 北京航空航天大学 A kind of composite I types layering bridging rule based on virtual hierarchy expanding test determines method
CN107330202B (en) * 2017-07-05 2020-06-30 北京航空航天大学 Method for determining I-type layered bridging rule of composite material
CN107356478A (en) * 2017-07-07 2017-11-17 吉林大学 Long fiber reinforced thermoplastic composites Multiscale Simulation Method
CN107356478B (en) * 2017-07-07 2019-10-15 吉林大学 Long fiber reinforced thermoplastic composites Multiscale Simulation Method
CN107451360B (en) * 2017-07-28 2021-03-09 中国航空工业集团公司西安飞机设计研究所 Numerical simulation method for composite material quarter turn continuous layering
CN107451360A (en) * 2017-07-28 2017-12-08 中国航空工业集团公司西安飞机设计研究所 The method for numerical simulation of the continuous laying of composite right-angled bend
CN107256320A (en) * 2017-07-28 2017-10-17 中科机器人科技有限公司 A kind of small-sized hand based on composite plys throws the finite element method of unmanned plane
CN108268692A (en) * 2017-11-21 2018-07-10 北京航空航天大学 Automatic fiber placement composite material curing deformation forecasting procedure
CN108595724A (en) * 2017-12-21 2018-09-28 刘长喜 Composite material revolving meber design method
CN108595724B (en) * 2017-12-21 2022-01-14 刘长喜 Composite material rotating part design method
CN108595754A (en) * 2018-03-20 2018-09-28 南京航空航天大学 The emulation mode of interlayer toughened composite laminated plate
CN108256281A (en) * 2018-03-26 2018-07-06 中国矿业大学 A kind of intensity prediction method for considering overlap joint interface topography and overlapping object graded properties
CN108804787A (en) * 2018-05-26 2018-11-13 南京理工大学 The method for being inserted into cohesive force unit simulation Rock Bridge Failure based on batch
CN108804787B (en) * 2018-05-26 2022-06-28 南京理工大学 Method for simulating rock bridge through based on mass insertion of cohesive force units
CN108681650A (en) * 2018-05-29 2018-10-19 内蒙古工业大学 A kind of warp-knitted spacer fabric composite material stress analysis method
CN109190225B (en) * 2018-08-24 2022-11-18 湘潭大学 Finite element analysis method for simulating internal cracks and interface cracks of hard coating
CN109190225A (en) * 2018-08-24 2019-01-11 湘潭大学 A kind of finite element method for simulating hard coat underbead crack and Interface Crack
CN109190258B (en) * 2018-09-06 2023-01-13 中国航发北京航空材料研究院 Composite material laminate I-type layering simulation method based on cohesion unit superposition
CN109190258A (en) * 2018-09-06 2019-01-11 中国航发北京航空材料研究院 A kind of composite-material laminates I type method for layering simulation based on the superposition of cohesion unit
CN110968916B (en) * 2018-09-28 2023-11-17 天津大学 Submarine metallurgical composite pipe die generator containing circumferential surface crack V-shaped groove weld and application thereof
CN110968916A (en) * 2018-09-28 2020-04-07 天津大学 Seabed metallurgy composite pipe die generator containing circumferential surface crack V-shaped groove welding seam and application thereof
CN110032765B (en) * 2019-03-11 2022-06-21 武汉科技大学 Method, system, device and medium for optimizing refractory material fracture parameters
CN110032765A (en) * 2019-03-11 2019-07-19 武汉科技大学 A kind of optimization method, system, device and the medium of refractory material fragmentation parameters
CN110047563B (en) * 2019-03-19 2021-10-26 首钢集团有限公司 Composite material lightning damage simulation method and device
CN110047563A (en) * 2019-03-19 2019-07-23 首钢集团有限公司 A kind of composite material lightning damage analogy method and device
CN110457734A (en) * 2019-06-05 2019-11-15 东南大学 One kind being bolted FRP section components fiber laying laying method and FRP section components
CN110222420A (en) * 2019-06-06 2019-09-10 吉林大学 A kind of failure judgment method of viscous riveting connection structure
CN110376055A (en) * 2019-06-20 2019-10-25 重庆大学 A kind of CFRP layer plate layering failure behaviour prediction technique based on novel cohesive force constitutive relation
CN110826284B (en) * 2019-11-21 2021-04-16 北京航空航天大学 Modeling and analyzing method for composite material laminated board of interweaved and laminated mixed laying layers
CN110826284A (en) * 2019-11-21 2020-02-21 北京航空航天大学 Modeling and analyzing method for composite material laminated board of interweaved and laminated mixed laying layers
CN111104751A (en) * 2019-12-24 2020-05-05 肇庆市海特复合材料技术研究院 Strength analysis method of composite material connecting rod
CN111104751B (en) * 2019-12-24 2023-04-07 肇庆市海特复合材料技术研究院 Strength analysis method of composite material connecting rod
CN113049361A (en) * 2019-12-27 2021-06-29 中国航空工业集团公司西安飞机设计研究所 Laminated board III type crack delamination failure test method and test piece
CN111553104A (en) * 2020-05-08 2020-08-18 东南大学 Simulation analysis method for bridge deck pavement structure
CN111553104B (en) * 2020-05-08 2023-12-05 东南大学 Bridge deck pavement structure simulation analysis method
CN112487557A (en) * 2020-12-03 2021-03-12 上海交通大学 Method for predicting interface failure and microscopic crack propagation of composite material under hydraulic permeation load
CN112733324A (en) * 2020-12-17 2021-04-30 重庆大学 Composite material laminate fatigue layering expansion behavior prediction method considering stress ratio and fiber bridging influence
CN112733324B (en) * 2020-12-17 2023-06-16 重庆大学 Composite material laminate fatigue layered expansion behavior prediction method considering stress ratio and fiber bridging influence
CN112699585A (en) * 2020-12-29 2021-04-23 中国航空工业集团公司西安飞机设计研究所 Finite element modeling method for composite thick-layer laminated plate joint
CN112595593A (en) * 2020-12-29 2021-04-02 哈尔滨工业大学 Method for obtaining interface strength of three-dimensional woven composite material
CN112699585B (en) * 2020-12-29 2024-04-09 中国航空工业集团公司西安飞机设计研究所 Finite element modeling method for thick-layer laminated board joint of composite material
CN113158502A (en) * 2021-01-19 2021-07-23 西北工业大学 Reliability analysis method and device, storage medium and electronic equipment
CN113128097B (en) * 2021-04-29 2023-11-17 浙江理工大学 Method for simulating and predicting heat transfer performance of porous nanofiber medium
CN113128097A (en) * 2021-04-29 2021-07-16 浙江理工大学 Method for simulating and predicting heat transfer performance of porous nanofiber medium
CN113971354A (en) * 2021-10-15 2022-01-25 厦门大学 Finite element simulation method for nano composite material bending fracture damage research
CN114201814A (en) * 2021-11-02 2022-03-18 南京航空航天大学 Civil transport aircraft windshield collision simulation modeling method based on cohesion model
CN114357823A (en) * 2021-12-08 2022-04-15 湖北工业大学 Curved surface moving load analysis method based on ABAQUS
CN114896763A (en) * 2022-04-18 2022-08-12 重庆大学 Method for predicting composite material I/II mixed fracture response surface
CN117524361A (en) * 2023-10-24 2024-02-06 天津大学 Fatigue crack initiation prediction method for cycle cohesive force model with variable accumulation length
CN117524361B (en) * 2023-10-24 2024-04-26 天津大学 Fatigue crack initiation prediction method for cycle cohesive force model with variable accumulation length

Also Published As

Publication number Publication date
CN103605865B (en) 2016-08-17

Similar Documents

Publication Publication Date Title
CN103605865A (en) Delamination propagation behavior simulation method of composite material multidirectional laminated plate on basis of cohesion model
Bergara et al. Fatigue crack propagation in complex stress fields: Experiments and numerical simulations using the Extended Finite Element Method (XFEM)
Vassilopoulos et al. Fatigue of fiber-reinforced composites
Qi et al. An approach to predict the mechanical properties of CFRP based on cross-scale simulation
CN107451307B (en) Method for multi-scale calculation of equivalent stiffness matrix of complex composite structure
CN109190258B (en) Composite material laminate I-type layering simulation method based on cohesion unit superposition
Qiao et al. Cohesive fracture simulation and failure modes of FRP–concrete bonded interfaces
Xiaohui et al. A new zig-zag theory and C 0 plate bending element for composite and sandwich plates
Wahab Simulating mode I fatigue crack propagation in adhesively-bonded composite joints
Yu et al. Experimental study and DEM modelling of bolted composite lap joints subjected to tension
Lampeas et al. Interlaminar stresses calculation using a stacked-shell finite element modeling approach
Geng et al. An interface damage model for high-cycle fatigue
Lin et al. Theoretical models and reliability assessment of steel-UHPC-steel composite beams in offshore structures
Lu et al. Evaluation of debonding failure of reinforced concrete girders strengthened in flexure with FRP laminates using finite element modeling
CN115391947A (en) Virtual analysis method and system for composite material structure of rail transit vehicle
Žmindák et al. Finite Element Implementation of Failure and Damage Simulation in Composite Plates
Guida et al. NiTi SMA wires coupled with Kevlar fabric: A real application of an innovative aircraft LE slat system in SMAHC material
Petrov Numerical simulation of fatigue failure of composite materials under compression
Zhuang et al. A progressive damage model for carbon fiber-reinforced polymer laminates subjected to fastener pull-through failure
Rao et al. Degradation model based on Tsai-Hill factors to model the progressive failure of fiber metal laminates
Dong et al. Reinforcement of fatigue damaged steel structures using CFRP lamellas-Part 3: Numerical simulation
Zhu et al. Experimental and numerical characterization of pultruded GFRP structural elements with stiffened Web-Flange junctions
Vassilopoulos et al. Introduction to the fatigue of fiber-reinforced polymer composites
Burke-Veliz et al. Fatigue crack shielding and deflection in plain bearings under large-scale yielding
Fang et al. X-FEM co-simulation of delamination and matrix cracking in fiber metal laminated structures under fatigue loading

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20161128

CF01 Termination of patent right due to non-payment of annual fee