CN103605865B - A kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method - Google Patents
A kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method Download PDFInfo
- Publication number
- CN103605865B CN103605865B CN201310627964.4A CN201310627964A CN103605865B CN 103605865 B CN103605865 B CN 103605865B CN 201310627964 A CN201310627964 A CN 201310627964A CN 103605865 B CN103605865 B CN 103605865B
- Authority
- CN
- China
- Prior art keywords
- model
- multidirectional
- composite
- load
- cohesive zone
- 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.)
- Expired - Fee Related
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method, comprise the following steps: (1), according to the structural parameters of composite multidirectional paving layer test coupon, sets up geometrical model;(2) calculate the key parameter of reflection interface behavior, be respectively provided with material properties;(3) geometrical model of paving layer test coupon multidirectional to composite carries out stress and strain model, composition model respectively, sets up three-dimensional finite element model;(4) according to actual conditions and the loaded-up condition of composite structure, the load of FEM model, boundary condition are determined;(5) FEM model based on cohesive zone model is carried out computational analysis, extract the load displacement curve of load(ing) point, the delamination behavior of analog composite material multidirectional laying plate, it is thus achieved that load displacement curve maximum load value, it was predicted that the derogatory behaviour of multidirectional laying plate.
Description
Technical field
The present invention relates to a kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method, it is adaptable to
The research of the widely used composite in the fields such as Aero-Space multidirectional laying plate delamination behavior and prediction.
Background technology
Advanced composite material is widely used at present in aerospace structure, and with aluminium alloy, titanium alloy, alloy
Steel becomes aerospace four big structural materials together.Wherein, carbon fiber resin matrix composite has the same of relatively high specific strength
Time, have specific stiffness concurrently high, anticorrosive, fatigue behaviour good and the plurality of advantages such as performance designability and become Modern aircraft structure
In the main a kind of composite used.Along with batch production and the application of toughened resin of middle modulus high-strength carbon fiber, this kind of
Composite application in aircaft configuration is also gradually developed into, by secondary load-carrying construction, the main force support structure that force-bearing situation is complicated.And answer
The safety of structure of condensation material is the most increasingly valued by the people.
Different from metal material, the failure mode of composite and failure mechanism are increasingly complex various.Owing to component material is fine
Dimension and the greatest differences of matrix mechanical property, the interlayer performance of composite laminated plate is far below performance in its face, at different loads
In the case of, the layering that low velocity impact or geometric discontinuity cause extends under fatigue load, finally causes the flexing of sublayer plate,
It it is a kind of main damage form of practice of composite laminate structures.Delamination behavior can cause strength and stiffness to incur loss,
Integrality and safety period for structure are all extremely disadvantageous, in some instances it may even be possible to the calamitous structure causing outside to be difficult to discover is broken
Bad.Therefore, during design and analyzing, the prediction to composite delamination behavior is a highly 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, accordingly
Composite unidirectional laminate test method (ASTM D5528-01 and ASTM D6671/D6671M-06) compared
Ripe.In generation and the build phase of delamination, cohesive zone model (CZM) is widely applied.Have passed through substantial amounts of
Research work, use the vertical cohesive zone model of this structures linear can well analog composite material unidirectional laminate crack initiation with
Growth, gained predicts the outcome and existing result of the test has good uniformity.
Compared with laying plate unidirectional with composite, the application in engineering practice of the composite multidirectional laying plate is more extensive.Compound
On the laying interface that material multidirectional laying plate machine direction is different, it is easy to form crackle and the delamination of crackle occurs.But,
The work that composite interlayer performance is carried out at present mostly also focuses in the research of composite single-ply plate, for composite wood
Expect that the delamination rule at multidirectional laying interface is rare to relate to.Due in composite multidirectional laying flaggy and the failure mechanism of interlayer
Extremely complex, there is fiber bridge joint, wear the various of failure such as layer extension, and can document for reference extremely limited, for multiple
The numerical simulation of condensation material multidirectional laying plate delamination, remains one at present and has much challenging work.In order in this respect
Carry out new research, to predict the damage propatagtion in engineering structure, while playing material property advantage, ensure structure peace
Entirely, need to introduce the cohesive zone model improved, preferably reflect the failure procedure of composite multidirectional laying plate.
Summary of the invention
The technology of the present invention solves problem: overcome the deficiencies in the prior art, it is provided that a kind of composite based on cohesive zone model is many
To laying plate delamination Behavior modeling method, it is adaptable to engineer applied, multidirectional laying plate ply stacking-sequence can be assisted to design, fall
Low experimentation cost, laying plate inefficacy damage mode multidirectional to composite is effectively analyzed and predicts simultaneously, preferably ensures
Safety of structure.
The technical solution of the present invention, the multidirectional laying plate delamination Behavior modeling of a kind of composite based on cohesive zone model
Method, it is achieved step is as follows:
Step A, utilizes the d solid modeling function of finite element software Abaqus, according to the knot of composite multidirectional laying plate
Structure parameter, sets up three geometrical models that composite multidirectional paving layer test coupon is corresponding;
Step B, select the linear constitutive relation of cohesive zone model, arrange cohesive zone related interfaces parameter and test record critical should
Become energy release rate, model is respectively provided with laying direction and material properties;
Three geometrical models are respectively divided grid, composition model, set up three-dimensional finite element model by step C;
Step D, according to actual conditions and the loaded-up condition of composite structure, force analysis situation, determine test specimen load,
Boundary condition, and be applied on FEM model;
Step E, carries out stress analysis to FEM model based on cohesive zone model, it is thus achieved that result of calculation, extracts load(ing) point
Load displacement curve, simulates composite based on cohesive zone model multidirectional laying plate delamination behavior, it is thus achieved that load displacement is bent
Line maximum load value, it was predicted that the derogatory behaviour of multidirectional laying plate.
Described step A utilizes the d solid modeling function of finite element software Abaqus, and the detailed process setting up geometrical model is:
(A1) open Abaqus software, utilize d solid modeling function, according to composite multidirectional paving layer test coupon structure
Parameter value, draws three geometric profile figures that test specimen upper, middle, and lower part is corresponding in two dimensional surface;
(A2) by entity stretching mode, multidirectional 3-D geometric model, the most upper and lower two parts through-thickness are set up respectively
Being stretched to the half of actual (real) thickness value, the mid portion through-thickness at cohesive zone place is stretched to 0.01mm.
Described step B determines the key parameter of reflection interface behavior, and the process that implements arranging material properties is:
(B1) definition local material direction, divides upper and lower two parts geometrical model according to each overlay thickness value, creates composite wood
Material laying also defines laying angle, needs when defining laying angle to set up corresponding local coordinate system according to actual conditions;
(B2) according to material property parameter, the material properties of upper and lower two parts geometrical model is defined;
(B3) select the linear constitutive relation of cohesive zone model, introduce, by subprogram, the facing with crack length change that test obtains
Boundary's strain energy release rate, calculates initial interface rigidity, boundary strength, and viscosity etc. reflects the interface parameter of interface behavior,
The material properties of definition cohesive zone place mid portion geometrical model.
Described step C carries out stress and strain model, composition model respectively to the geometrical model of test specimen, sets up the tool of three-dimensional finite element model
Body process is:
(C1) to upper and lower two parts along three direction Even seedings of length, and grid division, this grid need not ten
Divide fine;
(C2) calculate cohesive zone length, determine cohesive force unit size, quantity, according to cohesive force unit size and number of unit,
To mid portion Even seeding, divide the finest grid;
(C3) three geometrical models are assembled successively according to from the order of upper, middle and lower, arrange at two assembly and connection interfaces and tie up
Conclude a contract or treaty bundle, in addition for preventing laminate upper and lower part from running through, set up between the two sides, upper and lower part at precrack and connect
Touch, set up three-dimensional finite element model.
Described step D determines that geometrical model load, boundary condition realize process and be:
(D1) judge that I type or I/II mixed type loaded-up condition, paving layer test coupon multidirectional to composite carry out force analysis, really
Determine the restraint condition of test specimen and suffered load situation;
(D2) restraint condition above is arranged on the FEM model of load as boundary condition, simultaneously in test specimen finite element
The load of force analysis gained is applied on model.
Described step E carries out computational analysis to cohesive zone model, it is thus achieved that result of calculation, extracts the load displacement curve of load(ing) point,
Laying plate delamination Behavior modeling multidirectional to composite implements process and is:
(E1) FEM model is carried out stress analysis, calculate the strain energy rate of cohesive force unit, by reading cohesive force
Unit coordinate information, further according to critical strain energy release rate-crack length relation curve, determines the critical strain of cohesive force unit
Energy release rate, compares currency with critical value, it is judged that cohesive force element failure and the spread scenarios of crackle;
(E2) obtaining result of finite element, utilize the post-processing function of finite element software, the load displacement extracting load(ing) point is bent
Line, analog composite material multidirectional laying plate delamination behavior, it is thus achieved that load displacement curve maximum load value, when real load reaches
During to maximum, the extension of crackle needs not continue to increase load, thus realizes the prediction to multidirectional laying dash-board injury behavior.
Present invention advantage compared with prior art is:
(1) present invention composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method, overcomes pin
The deficiency that unidirectional laying plate is studied, it is contemplated that the different fiber direction shadows to multidirectional laying plate delamination behavior
Ring, the delamination behavior of analog composite material multidirectional laying plate, thus effectively analyze and predict the derogatory behaviour of multidirectional laying plate.
(2) present invention composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method, can be effective
Analyze and simulate the delamination behavior of the composite multidirectional laying plate at multiple interface, according to Numerical Simulation Results, assist multidirectional
The design in laying plate laying direction.
(3) present invention composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method, only needs to set up
Cohesive zone model, carries out finite element analysis, laying plate delamination behavior multidirectional to composite can carry out numerical simulation, suitable
The prediction of the composite multidirectional laying plate delamination failure behaviour in engineer applied and analysis, can significantly shorten development
In the cycle, reduce experimentation cost.
Accompanying drawing explanation
Fig. 1 is the flowchart of the present invention;
Fig. 2 is three geometrical model schematic diagrames that composite multidirectional paving layer test coupon is corresponding;
Fig. 3 is to assemble geometrical model schematic diagram;
Fig. 4 is finite element method (fem) analysis flow chart based on cohesive zone model.
Detailed description of the invention
As it is shown in figure 1, the inventive method is implemented as:
1, according to composite multidirectional paving layer test coupon values of the structural parameters, Abaqus software is set up the geometrical model of test specimen.
Open Abaqus software, utilize d solid modeling function, according to composite multidirectional paving layer test coupon values of the structural parameters,
In two dimensional surface, draw the geometric profile figure that test specimen is corresponding respectively, by entity stretching mode, set up 3-D geometric model respectively,
As in figure 2 it is shown, wherein a, c two parts through-thickness is stretched to the half of actual (real) thickness value, b part through-thickness stretches
0.01mm。
2, determine the key parameter of reflection interface behavior, laying angle and material properties are set.
Definition local material direction, divides the geometrical model of a and c according to each overlay thickness value, creates composite plys and determines
Justice laying angle, needs when defining laying angle to set up corresponding local coordinate system according to actual conditions, joins according to material property
Number, definition a and c geometrical model material properties, select the linear constitutive relation of cohesive zone model, introduce test obtain with splitting
The critical strain energy release rate of line length change, calculates the reflection interface row such as initial interface rigidity, boundary strength, viscosity
For parameter, definition b model is the material properties of corresponding cohesive zone.
3, geometrical model is carried out respectively stress and strain model, composition model, set up three-dimensional finite element model.
To a, c two parts along three direction Even seedings of length, and grid division, this grid need not very fine,
Calculate cohesive zone length, determine cohesive force unit size, quantity, according to cohesive force unit size and number of unit, to b part
Even seeding, divides the finest grid, and by geometrical model according to from a, the order of b, c assembles successively, at geometrical model
The assembly and connection interface of a/b and b/c arranges binding constraints, sets up contact, have as it is shown on figure 3, set up three-dimensional at precrack
Limit meta-model.
4, geometrical model load, boundary condition are determined.
Judge that I type or I/II mixed type loaded-up condition, paving layer test coupon multidirectional to composite carry out force analysis, determine test specimen
Restraint condition and suffered load situation, restraint condition above is arranged on the FEM model of load as boundary condition,
On test specimen FEM model, apply the load of force analysis gained simultaneously.
5, analog composite material multidirectional laying plate delamination behavior.
FEM model is carried out stress analysis, calculates the strain energy rate of correspondence position, it is thus achieved that crack length, further according to facing
Boundary's strain energy release rate-crack length relation curve, determines the critical energy release rate of this position, by currency with critical
Value compares, it is judged that cohesive force element failure and the spread scenarios of crackle, idiographic flow is as shown in Figure 4, it is thus achieved that FEM calculation
As a result, utilize the post-processing function of finite element software, extract the load displacement curve of load(ing) point, simulate based on cohesive zone model
Composite multidirectional laying plate delamination behavior, it is thus achieved that load displacement curve maximum load value, it is judged that when real load reaches
When being worth greatly, the extension of crackle needs not continue to increase load, thus realizes the prediction to multidirectional laying dash-board injury behavior.
Non-elaborated part of the present invention belongs to the known technology of those skilled in the art.
Claims (3)
1. composite based on a cohesive zone model multidirectional laying plate delamination Behavior modeling method, it is characterised in that include
Following steps:
Step A, utilizes the d solid modeling function of finite element software Abaqus, according to the knot of composite multidirectional laying plate
Structure parameter, sets up three geometrical models that composite multidirectional paving layer test coupon is corresponding;
In described step A, structural parameters according to composite multidirectional paving layer test coupon are set up the geometrical model of test specimen and are realized process
For:
(A1) open Abaqus software, utilize d solid modeling function, according to composite multidirectional paving layer test coupon structure
Parameter value, draws three geometric profile figures that test specimen upper, middle, and lower part is corresponding in two dimensional surface;
(A2) by entity stretching mode, multidirectional 3-D geometric model, the most upper and lower two parts through-thickness are set up respectively
Being stretched to the half of actual (real) thickness value, the mid portion through-thickness at cohesive zone place is stretched to 0.01mm;
Step B, select the linear constitutive relation of cohesive zone model, arrange cohesive zone related interfaces parameter and test record critical should
Become energy release rate, model is respectively provided with laying direction and material properties;
Described step B determines the key parameter of reflection interface behavior, and arrange laying direction and material properties implements process
For:
(B1) definition local material direction, divides upper and lower two parts geometrical model according to each overlay thickness value, creates composite wood
Material laying also defines laying angle, needs when defining laying angle to set up corresponding local coordinate system according to actual conditions;
(B2) according to material property parameter, the material properties of upper and lower two parts geometrical model is defined;
(B3) select the linear constitutive relation of cohesive zone model, introduce, by subprogram, the facing with crack length change that test obtains
Boundary's strain energy release rate, calculates initial interface rigidity, boundary strength, and viscosity etc. reflects the interface parameter of interface behavior,
The material properties of definition cohesive zone place mid portion geometrical model;
Three geometrical models are respectively divided grid, composition model, set up three-dimensional finite element model by step C;
Described step C carries out stress and strain model, composition model respectively to the geometrical model of test specimen, sets up the tool of three-dimensional finite element model
Body process is:
(C1) to upper and lower two parts along three direction Even seedings of length, and grid division, this grid need not ten
Divide fine;
(C2) calculate cohesive zone length, determine cohesive force unit size, quantity, according to cohesive force unit size and number of unit,
To mid portion Even seeding, divide the finest grid;
(C3) three geometrical models are assembled successively according to from the order of upper, middle and lower, arrange at two assembly and connection interfaces and tie up
Conclude a contract or treaty bundle, in addition for preventing laminate upper and lower part from running through, set up between the two sides, upper and lower part at precrack and connect
Touch, set up three-dimensional finite element model;
Step D, according to actual conditions and the loaded-up condition of composite structure, force analysis situation, determine test specimen load,
Boundary condition, and be applied on FEM model;
Step E, carries out stress analysis to FEM model based on cohesive zone model, it is thus achieved that result of calculation, after utilizing finite element
Process function, extract the load displacement curve of load(ing) point, simulate the multidirectional paving Laminates With Delamination of composite based on cohesive zone model and expand
Exhibition behavior, it is thus achieved that load displacement curve maximum load value, it was predicted that multidirectional laying dash-board injury behavior.
A kind of composite based on cohesive zone model the most according to claim 1 multidirectional laying plate delamination behavior mould
Plan method, it is characterised in that: described step D determines that geometrical model load, boundary condition realize process and be:
(D1) judge that I type or I/II mixed type loaded-up condition, paving layer test coupon multidirectional to composite carry out force analysis, really
The restraint condition determining test specimen is loaded with institute;
(D2) restraint condition above is arranged on the FEM model of load as boundary condition, simultaneously in test specimen finite element
The load of force analysis gained is applied on model.
A kind of composite based on cohesive zone model the most according to claim 1 multidirectional laying plate delamination behavior mould
Plan method, it is characterised in that: described step E carries out computational analysis to cohesive zone model, it is thus achieved that result of calculation, extracts load(ing) point
Load displacement curve, laying plate delamination Behavior modeling multidirectional to composite implements process and is:
(E1) FEM model is carried out stress analysis, calculate the strain energy rate of cohesive force unit, by reading cohesive force
Unit coordinate information, further according to critical strain energy release rate-crack length relation curve, determines the critical strain of cohesive force unit
Energy release rate, compares currency with critical value, it is judged that cohesive force element failure and the spread scenarios of crackle;
(E2) obtaining result of finite element, utilize the post-processing function of finite element software, the load displacement extracting load(ing) point is bent
Line, analog composite material multidirectional laying plate delamination behavior, it is thus achieved that load displacement curve maximum load value, when real load reaches
During to maximum, the extension of crackle needs not continue to increase load, thus realizes the prediction to multidirectional laying dash-board injury behavior.
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 CN103605865A (en) | 2014-02-26 |
CN103605865B true 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 (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391783A (en) * | 2017-06-02 | 2017-11-24 | 西安交通大学 | A kind of asymmetric Combined Cycle Unit shafting vibration monitoring method based on finite element theory |
CN107832492A (en) * | 2017-10-11 | 2018-03-23 | 大连理工大学 | A kind of corrosion of steel structure Fatigue Damage Calculation method based on cohesive zone model |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6297955B2 (en) * | 2014-09-16 | 2018-03-20 | 株式会社Ihi | Design support apparatus, design support method, and program |
CN105335575B (en) * | 2015-11-25 | 2019-06-28 | 中国航空工业集团公司沈阳飞机设计研究所 | A kind of finite element modeling method of the structure-integrated laminate of electronic antenna |
US10775286B2 (en) * | 2016-01-28 | 2020-09-15 | The Boeing Company | Method for test data reduction of composite intralaminar failure mode |
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 |
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 |
CN107330202B (en) * | 2017-07-05 | 2020-06-30 | 北京航空航天大学 | Method for determining I-type layered bridging rule of composite material |
CN107356478B (en) * | 2017-07-07 | 2019-10-15 | 吉林大学 | Long fiber reinforced thermoplastic composites Multiscale Simulation Method |
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 |
CN107451360B (en) * | 2017-07-28 | 2021-03-09 | 中国航空工业集团公司西安飞机设计研究所 | Numerical simulation method for composite material quarter turn continuous layering |
CN108268692A (en) * | 2017-11-21 | 2018-07-10 | 北京航空航天大学 | Automatic fiber placement composite material curing deformation forecasting procedure |
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 |
CN108256281B (en) * | 2018-03-26 | 2020-03-10 | 中国矿业大学 | Strength prediction method considering lapping interface morphology and lapping object gradient property |
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 |
CN109190258B (en) * | 2018-09-06 | 2023-01-13 | 中国航发北京航空材料研究院 | Composite material laminate I-type layering simulation method based on cohesion unit superposition |
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 |
CN110032765B (en) * | 2019-03-11 | 2022-06-21 | 武汉科技大学 | Method, system, device and medium for optimizing refractory material fracture parameters |
CN110047563B (en) * | 2019-03-19 | 2021-10-26 | 首钢集团有限公司 | Composite material lightning damage simulation method and device |
CN110457734B (en) * | 2019-06-05 | 2023-02-28 | 东南大学 | Fiber layer laying method for bolt connection FRP (fiber reinforced Plastic) sectional material component and FRP sectional material component |
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 |
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 |
CN111553104B (en) * | 2020-05-08 | 2023-12-05 | 东南大学 | Bridge deck pavement structure simulation analysis method |
CN112487557B (en) * | 2020-12-03 | 2022-10-14 | 上海交通大学 | Method for predicting interface failure and microscopic crack propagation of composite material under hydraulic permeation load |
CN112733324B (en) * | 2020-12-17 | 2023-06-16 | 重庆大学 | Composite material laminate fatigue layered expansion behavior prediction method considering stress ratio and fiber bridging influence |
CN112595593B (en) * | 2020-12-29 | 2021-10-29 | 哈尔滨工业大学 | 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 |
CN113971354A (en) * | 2021-10-15 | 2022-01-25 | 厦门大学 | Finite element simulation method for nano composite material bending fracture damage research |
CN114201814B (en) * | 2021-11-02 | 2023-05-30 | 南京航空航天大学 | Civil conveyor 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 |
CN114896763B (en) * | 2022-04-18 | 2024-07-12 | 重庆大学 | Prediction method of composite material I/II mixed fracture response surface |
CN117524361B (en) * | 2023-10-24 | 2024-04-26 | 天津大学 | Fatigue crack initiation prediction method for cycle cohesive force model with variable accumulation length |
Citations (2)
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 |
CN103345545A (en) * | 2013-06-14 | 2013-10-09 | 北京航空航天大学 | Method for predicating composite material Pi-shaped non-planar glue joint strength based on triangular envelopes |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201013227D0 (en) * | 2010-08-06 | 2010-09-22 | Rolls Royce Plc | A composite material and method |
-
2013
- 2013-11-28 CN CN201310627964.4A patent/CN103605865B/en not_active Expired - Fee Related
Patent Citations (2)
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 |
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)
Title |
---|
Fatigue delamination growth rates and thresholds of composite laminates under mixed mode loading;Jianyu Zhang等;《International Journal of Fatigue》;20120730;第40卷;第7-15页 * |
Simulation of multi-directional crack growth in braided composite T-piece specimens using cohesive models;J.CHEN;《Fatigue & Fracture of Engineering Materials & Structures》;20110228;第34卷(第2期);第123-130页 * |
碳纤维复合材料层合板损伤脱层过程的有限元分析;苏罗丹等;《第17届全国复合材料学术会议论文(复合材料力学分论坛)集》;20121001;正文第303页右栏第2段,第304页左栏第2-4段,图1,3-6,9-10 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391783A (en) * | 2017-06-02 | 2017-11-24 | 西安交通大学 | A kind of asymmetric Combined Cycle Unit shafting vibration monitoring method based on finite element theory |
CN107391783B (en) * | 2017-06-02 | 2020-04-28 | 西安交通大学 | Finite element theory-based asymmetric combined cycle unit shafting vibration monitoring method |
CN107832492A (en) * | 2017-10-11 | 2018-03-23 | 大连理工大学 | A kind of corrosion of steel structure Fatigue Damage Calculation method based on cohesive zone model |
Also Published As
Publication number | Publication date |
---|---|
CN103605865A (en) | 2014-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103605865B (en) | A kind of composite based on cohesive zone model multidirectional laying plate delamination Behavior modeling method | |
Kärger et al. | Development and validation of a CAE chain for unidirectional fibre reinforced composite components | |
CN107451307B (en) | Method for multi-scale calculation of equivalent stiffness matrix of complex composite structure | |
Wang et al. | Bolted glulam beam-column connections under different combinations of shear and bending | |
CN107451308B (en) | Multi-scale calculation method for equivalent heat conduction coefficient of complex composite material structure | |
Xu et al. | Mechanical behavior and modeling of dowelled steel-to-timber moment-resisting connections | |
Ammendolea et al. | Crack propagation modeling in functionally graded materials using Moving Mesh technique and interaction integral approach | |
Feih et al. | Adhesive and composite failure prediction of single-L joint structures under tensile loading | |
Khelifa et al. | Finite element analysis of flexural strengthening of timber beams with Carbon Fibre-Reinforced Polymers | |
Wang et al. | Finite element simulation of the failure process of single fiber composites considering interface properties | |
CN109190258B (en) | Composite material laminate I-type layering simulation method based on cohesion unit superposition | |
CN107451309B (en) | Method for multi-scale calculation of equivalent thermal expansion coefficient of complex composite material structure | |
Khalili et al. | Numerical analysis of adhesively bonded T-joints with structural sandwiches and study of design parameters | |
Cao et al. | Computational investigation of both intra-laminar matrix cracking and inter-laminar delamination of curved composite components with cohesive elements | |
Li et al. | Mechanical properties of L-joint with composite sandwich structure | |
Buka-Vaivade et al. | Experimental verification of design procedure for elements from cross-laminated timber | |
Wu et al. | Damage assessment of braided composite tube subjected to repeated transverse impact | |
Li et al. | Advances in the characterization of residual stress in composite structures | |
Qian et al. | Effect of gap and overlap fiber placement defects on the delamination behavior of L-shaped composite laminates | |
Borrelli et al. | On the use of global–local kinematic coupling approaches for delamination growth simulation in stiffened composite panels | |
Carello et al. | Building block approach’for structural analysis of thermoplastic composite components for automotive applications | |
Turco et al. | Influence of width and thickness of composite laminates on the flexural behavior of reinforced concrete beams and slabs | |
Littell et al. | Effect of microscopic damage events on static and ballistic impact strength of triaxial braid composites | |
Mishra et al. | Failure initiation in composite structures under low-velocity impact: Analytical studies | |
Wang et al. | A FEA simulation model for thin-walled C-section composite beam assembling with R-angle deviation |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20161128 |