CN109190225A - A kind of finite element method for simulating hard coat underbead crack and Interface Crack - Google Patents
A kind of finite element method for simulating hard coat underbead crack and Interface Crack Download PDFInfo
- Publication number
- CN109190225A CN109190225A CN201810971796.3A CN201810971796A CN109190225A CN 109190225 A CN109190225 A CN 109190225A CN 201810971796 A CN201810971796 A CN 201810971796A CN 109190225 A CN109190225 A CN 109190225A
- Authority
- CN
- China
- Prior art keywords
- hard coat
- model
- geometrical model
- finite element
- pressure head
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention provides a kind of finite element methods for simulating hard coat underbead crack and Interface Crack, method provided by the invention has combined cohesive force elements method and extended finite element method, the Interface Cracking behavior of hard coat can be intended by mulch very well by adding the cohesive force unit between hard coat and substrate, and the crack propagation region of the three-dimensional finite element model defined in hard coat can the effective film internal fissure that extends along free routing of simulation.Method provided by the invention can simultaneously simulate hard coat interfacial failure and the extension of film internal fissure, according to Numerical Simulation Results, the experimental designs of the standby hard coat with strong binding performance of guidance system and improvement.Method provided by the invention need to only be established and finite element model similar in experiment condition, finite element analysis is carried out, the appearance to hard coat crackle can be realized and propagation behavior carries out numerical simulation, the analysis suitable for hard coat failure behaviour, facilitate guiding experiment design, reduces experimental cost.
Description
Technical field
The present invention relates to hard coat defeat technique fields, and in particular to a kind of simulation hard coat underbead crack and interface
The finite element method of crackle.
Background technique
Hard coat proposes the performance of improvement workpiece as a kind of protection materials haveing excellent performance, application and popularization
The life and reliability of high workpiece, economize on resources and protect environment etc. have great importance.Wherein, diamond-film-like
(Diamond-like carbon, DLC) due to having many advantages, such as high rigidity, low-friction coefficient and good corrosion resistance,
The fields such as machining, mold manufacture, aerospace are widely used.However, hard coat such as DLC film is in many substrates
The problem of binding force on especially metal base material surface is weak, is easy to cause film layer to crack and leads to coating failure causes people
Extensive concern.The study found that failure mode of the DLC film under load effect mainly includes film internal fissure and Interface Crack,
The generation and extension of both crackles can promote film layer finally to fall off from substrate surface, influence the further of DLC film excellent properties
Using.Therefore, accurately the generation of analysis crackle and propagation behavior, exploration influence the factor of crack propagation, have for experiment preparation
The hard coats such as the DLC film of strong binding performance, which provide directive function, just seems necessary.
Both at home and abroad to coating because the Problem of Failure caused by crack propagation has carried out a large amount of theoretical and numerical analysis.Mesh
Before, the finite element methods such as virtual crack closure techniques, cohesive force unit, extension finite element are usually used in analyzing coating crack expansion
Exhibition, but virtual crack closure techniques are unable to simulating crack germinating, cohesive force unit can only simulate splitting along specified path extension
Line.For hard coat, hardness and brittleness with higher, propagation direction is difficult to determine after crackle generates, and extension has
Limiting first method has the advantages that precrack position and the crackle extended along free routing can not can be simulated, and is very suitable for
Simulate the failure of the hard coats such as DLC film under external load function.In recent years, extension finite element method has been applied to coating material mistake
In the research of effect, but mostly can only one of independent analysis film internal fissure or Interface Crack, can analyze simultaneously along any
The film internal fissure of Directional Extension and the finite element method of the Interface Crack propagation behavior extended along specific direction are seldom.
Summary of the invention
The purpose of the present invention is to provide a kind of finite element analysis sides for simulating hard coat underbead crack and Interface Crack
Method, method provided by the invention can instruct the experiment of indentation method evaluated for film/base binding performance, have strong combine for experiment preparation
The hard coat of performance provides guidance, reduces experimental cost, while carrying out finite element analysis to the factor for influencing crack propagation, more
Crack propagation bring is inhibited to endanger well.
In order to achieve the above-mentioned object of the invention, the present invention the following technical schemes are provided:
The present invention provides a kind of finite element methods for simulating hard coat underbead crack and Interface Crack, including with
Lower step:
(1) utilize finite element software d solid modeling function, according to simulation ballast effect push-down head, hard coat and
The values of the structural parameters of substrate establishes pressure head geometrical model, hard coat geometrical model and substrate geometry model respectively；
(2) pressure head geometrical model, hard coat geometrical model and substrate geometry model in the step (1) are carried out respectively
Grid dividing generates pressure head stand-alone entity, hard coat stand-alone entity and substrate stand-alone entity；Based on the hard coat geometry
The bottom surface of model generates cohesive force unit；
(3) material properties of hard coat geometrical model and substrate geometry model in the step (2) are defined, are selected
Select cohesive force constitutive relation and failure parameter；
(4) group is carried out to pressure head geometrical model, hard coat geometrical model and substrate geometry model in the step (2)
Dress, obtains three-dimensional finite element model, the crack propagation area of the three-dimensional finite element model is defined in hard coat geometrical model
Domain；
(5) according to the actual conditions of hard coat and loaded-up condition, analysis mode, boundary condition and loading method are determined,
And load substep is applied in the step (4) in three-dimensional finite element model；
(6) stress analysis is carried out to three-dimensional finite element model in the step (5), utilizes the post-processing function of finite element software
The relationship of energy assumed (specified) load and pressure depth, generation and propagation behavior of the analysis hard coat Interface Crack with film internal fissure, probes into shadow
Ring the factor of crack propagation.
Preferably, the shape of pressure head geometrical model is half cone-shaped, hard coat geometrical model and base in the step (1)
The shape of bottom geometrical model is semi-cylindrical.
Preferably, in the step (1) pressure head geometrical model, hard coat geometrical model and substrate geometry model foundation
Method, comprising the following steps:
Using the d solid modeling function of finite element software, according to simulation ballast effect push-down head, hard coat and base
It is corresponding with pressure head, hard coat and substrate respectively to draw the part of upper, middle and lower three in a two-dimensional plane for the values of the structural parameters at bottom
Geometric profile；In such a way that rotation generates entity, by three geometric profiles of gained respectively along 180 ° of center axis rotation, obtain
Pressure head geometrical model, hard coat geometrical model and substrate geometry model.
Preferably, the step (2) specifically:
By the way of diametrically inclined in two-way respectively to hard coat geometrical model and substrate geometry model arrangement seed,
It, to grid dividing is carried out to pressure head geometrical model arrangement seed by the way of cone bottom inclined in two-way, is generated isolated real using along the vertex of a cone
Body；It is then based on the bottom surface of hard coat geometrical model, cohesive force unit is generated using deflection method, selects cohesive force cell type
And definition unit is deleted.
Preferably, the step (3) specifically:
According to the physical parameter of hard coat, the material properties of hard coat geometrical model, failure mode and critical are defined
Parameter defines the material properties of substrate geometry model；Select the linear constitutive relation of cohesive force, definition interfaces rigidity, boundary strength
And viscosity, definition interfaces failure mode and critical parameters, define the material properties of cohesive force unit.
Preferably, the method for three-dimensional finite element model is assembled in the step (4), comprising the following steps:
The successively group in accordance with the order from top to bottom by pressure head geometrical model, hard coat geometrical model and substrate geometry model
Dress, establishes surface-to-surface contact between pressure head geometrical model and the upper surface of hard coat geometrical model, in hard coat geometry mould
Binding constraints is established between the lower surface and substrate geometry model of type.
Preferably, the step (5) specifically:
According to the actual conditions and loaded-up condition of hard coat, using static Gneral analysis mode, and incremental step is determined
Quantity and size；Staff cultivation is applied to substrate geometry model, constrains pressure head geometrical model along the movement and rotation in other directions, only
Leave movement along the vertical direction；According to experiment condition, load substep is applied on pressure head geometrical model.
Preferably, the step (6) specifically:
Stress analysis is carried out to three-dimensional finite element model, it is deep using the post-processing function assumed (specified) load and pressure of finite element software
Relationship: by read stress distribution situation, judge the generation and spread scenarios of film internal fissure；It is rigid by reading cohesive force unit
Amount of degradation is spent, judges the failure of cohesive force unit and the generation of Interface Crack and spread scenarios；When institute is extended according to film internal fissure
The length for stating film internal fissure changes with time situation, calculates the variation relation of length between Interface Crack and film internal fissure, adjusts
Hard coat thickness and elasticity modulus are saved, to further determine that the factor for influencing hard coat crack propagation.
The present invention provides a kind of finite element methods for simulating hard coat underbead crack and Interface Crack, mainly have
Following advantages:
(1) method provided by the invention has combined cohesive force elements method and extended finite element method, addition hard coat with
The fine mulch of cohesive force unit energy between substrate intends the Interface Cracking behavior of hard coat, utilizes extension limited in hard coat
First method (being realized by defining enrichment unit) can the effective film internal fissure that extends along free routing of simulation.
(2) method provided by the invention can simultaneously simulate hard coat interfacial failure and the extension of film internal fissure,
According to Numerical Simulation Results, the experimental designs of the standby hard coat with strong binding performance of guidance system and improvement.
(3) method provided by the invention need to only establish with finite element model similar in experiment condition, carry out finite element analysis,
The appearance to hard coat crackle can be realized and propagation behavior carries out numerical simulation, point suitable for hard coat failure behaviour
Analysis facilitates guiding experiment design, reduces experimental cost.
Therefore, method provided by the invention can instruct the experiment of indentation method evaluated for film/base binding performance, for experiment preparation
Hard coat with strong binding performance provides guidance, reduces experimental cost, while having to the factor for influencing crack propagation
The analysis of limit, the preferably harm of inhibition crack propagation bring.
Detailed description of the invention
Fig. 1 is the process of the finite element method of simulation hard coat underbead crack provided by the invention and Interface Crack
Figure；
Fig. 2 is the schematic diagram for simulating corresponding three geometrical models of hard coat crack propagation；
Fig. 3 is three geometrical model grid divisions and the schematic diagram for generating stand-alone entity；
Fig. 4 is assembling and the geometrical model schematic diagram for establishing contact relation；
Fig. 5 is the relational graph for analyzing load and pressure depth in crack propagation process；
Schematic diagram when Fig. 6 is film internal fissure and Interface Crack starting；
Fig. 7 is the schematic diagram after the completion of crack propagation.
Specific embodiment
The present invention provides a kind of finite element methods for simulating hard coat underbead crack and Interface Crack, including with
Lower step:
(1) utilize finite element software d solid modeling function, according to simulation ballast effect push-down head, hard coat and
The values of the structural parameters of substrate establishes pressure head geometrical model, hard coat geometrical model and substrate geometry model respectively；
(2) pressure head geometrical model, hard coat geometrical model and substrate geometry model in the step (1) are carried out respectively
Grid dividing generates pressure head stand-alone entity, hard coat stand-alone entity and substrate stand-alone entity；Based on the hard coat geometry
The bottom surface of model generates cohesive force unit；
(3) material properties of hard coat geometrical model and substrate geometry model in the step (2) are defined, are selected
Select cohesive force constitutive relation and failure parameter；
(4) group is carried out to pressure head geometrical model, hard coat geometrical model and substrate geometry model in the step (2)
Dress, obtains three-dimensional finite element model, the crack propagation area of the three-dimensional finite element model is defined in hard coat geometrical model
Domain；
(5) according to the actual conditions of hard coat and loaded-up condition, analysis mode, boundary condition and loading method are determined,
And load substep is applied in the step (4) in three-dimensional finite element model；
(6) stress analysis is carried out to three-dimensional finite element model in the step (5), utilizes the post-processing function of finite element software
The relationship of energy assumed (specified) load and pressure depth, generation and propagation behavior of the analysis hard coat Interface Crack with film internal fissure, probes into shadow
Ring the factor of crack propagation.
The present invention utilizes the d solid modeling function of finite element software, is applied according to simulation ballast effect push-down head, hard
The values of the structural parameters of layer and substrate, establishes pressure head geometrical model, hard coat geometrical model and substrate geometry model respectively.This hair
It is bright there is no special restriction for the finite element software, using finite element software well known to those skilled in the art, tool
Body such as ABAQUS software.In the present invention, pressure head geometrical model is preferably shaped to half cone-shaped, hard coat geometrical model with
The shape of substrate geometry model is preferably all semi-cylindrical.The pressure head geometrical model, hard coat geometrical model and substrate are several
The method for building up of what model, preferably includes following steps:
Using the d solid modeling function of finite element software, according to simulation ballast effect push-down head, hard coat and base
It is corresponding with pressure head, hard coat and substrate respectively to draw the part of upper, middle and lower three in a two-dimensional plane for the values of the structural parameters at bottom
Geometric profile；In such a way that rotation generates entity, by three geometric profiles of gained respectively along 180 ° of center axis rotation, obtain
Pressure head geometrical model, hard coat geometrical model and substrate geometry model.
After establishing pressure head geometrical model, hard coat geometrical model and substrate geometry model, the present invention is respectively to the pressure
Head geometrical model, hard coat geometrical model and substrate geometry model carry out grid dividing, generate pressure head stand-alone entity, hard applies
Layer stand-alone entity and substrate stand-alone entity；Based on the bottom surface of the hard coat geometrical model, cohesive force unit is generated.The present invention
It is preferred that respectively to hard coat geometrical model and substrate geometry model arrangement seed, use by the way of diametrically inclined in two-way
Grid dividing is carried out to pressure head geometrical model arrangement seed to the mode of cone bottom inclined in two-way along the vertex of a cone, generates stand-alone entity；So
Bottom surface afterwards based on hard coat geometrical model generates cohesive force unit using deflection method, selects cohesive force cell type and determines
Adopted element deletion.Present invention preferably employs the modes of diametrically inclined in two-way respectively to hard coat geometrical model and substrate geometry
Model arranges that seed carries out grid dividing, i.e., among geometrical model around, size of mesh opening is gradually increased, wherein grid ruler
Very little size can be calculated and compared to choose by the geometrical model to different size of mesh opening, hard coat geometrical model institute
There is grid of the grid preferably than substrate geometry model corresponding portion finer.Present invention preferably employs along the vertex of a cone, to boring, bottom is two-way to incline
Oblique mode carries out grid dividing to pressure head geometrical model arrangement seed, i.e., from the vertex of a cone of geometrical model to cone bottom, size of mesh opening
It is gradually increased.Grid dividing is carried out to the pressure head geometrical model, hard coat geometrical model and substrate geometry model, is generated lonely
After vertical entity, the present invention is preferably based on the bottom surface of hard coat geometrical model, generates cohesive force unit using deflection method, interface is viscous
Connecing layer unit type is cohesive force unit (COH3D8) and definition unit deletion, i.e. the Stiffness Deterioration amount (SDEG) of cohesive force unit
Reach 1, indicates that unit is entirely ineffective, form crackle.
After generation pressure head stand-alone entity, hard coat stand-alone entity and substrate stand-alone entity and cohesive force unit, this hair
The bright material properties to the pressure head geometrical model, hard coat geometrical model and substrate geometry model are defined, in selection
Poly- power constitutive relation and failure parameter.The present invention preferably according to the physical parameter of hard coat, defines hard coat geometrical model
Material properties, failure mode and critical parameters, define substrate geometry model material properties；This linear structure of cohesive force is selected to close
System, definition interfaces rigidity, boundary strength and viscosity, definition interfaces failure mode and critical parameters define cohesive force unit
Material properties.
After generating stand-alone entity and cohesive force unit, the present invention is to the pressure head geometrical model, hard coat geometrical model
It is assembled with substrate geometry model, obtains three-dimensional finite element model, the three-dimensional is defined in hard coat geometrical model to be had
Limit the crack propagation region of meta-model.In the present invention, the method for assembling the three-dimensional finite element model preferably includes following step
It is rapid:
The successively group in accordance with the order from top to bottom by pressure head geometrical model, hard coat geometrical model and substrate geometry model
Dress, establishes surface-to-surface contact between pressure head geometrical model and the upper surface of hard coat geometrical model, in hard coat geometry mould
Binding constraints is established between the lower surface and substrate geometry model of type.
The present invention is preferably for entire hard coat geometrical model to be defined as the crack propagation of the three-dimensional finite element model
Region.
Assembling three-dimensional finite element model and after defining the crack propagation region of the three-dimensional finite element model, the present invention according to
The actual conditions and loaded-up condition of hard coat determine analysis mode, boundary condition and loading method, and load substep are applied
Onto the three-dimensional finite element model.The present invention is logical using static state preferably according to the actual conditions of hard coat and loaded-up condition
With analysis mode, and determine the quantity and size of incremental step；Staff cultivation is applied to substrate geometry model, constrains pressure head geometrical model
Movement and rotation along other directions, leave behind movement along the vertical direction；According to experiment condition, load substep is applied to pressure
On head geometrical model.
After load substep is applied in the three-dimensional finite element model, the present invention carries out the three-dimensional finite element model
Stress analysis analyzes hard coat Interface Crack using the relationship that the post-processing function assumed (specified) load and pressure of finite element software are deep
The factor for influencing crack propagation is probed into generation and propagation behavior with film internal fissure.The present invention is preferably to three-dimensional finite element model
Stress analysis is carried out, utilizes the relationship that the post-processing function assumed (specified) load and pressure of finite element software are deep: by reading stress distribution
Situation judges the generation and spread scenarios of film internal fissure；By reading cohesive force element stiffness amount of degradation, cohesive force unit is judged
Failure and Interface Crack generation and spread scenarios；The length of the film internal fissure is at any time when being extended according to film internal fissure
Situation of change calculates the variation relation of length between Interface Crack and film internal fissure, adjusts hard coat thickness and springform
Amount, to further determine that the factor for influencing hard coat crack propagation.
Below in conjunction with the embodiment in the present invention, the technical solution in the present invention is clearly and completely described.It is aobvious
So, described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based on the reality in the present invention
Example is applied, every other embodiment obtained by those of ordinary skill in the art without making creative efforts all belongs to
In the scope of protection of the invention.
Embodiment 1
Fig. 1 is the process of the finite element method of simulation hard coat underbead crack provided by the invention and Interface Crack
Figure simulates the generation and extension of DLC film Interface Crack and film crackle referring to Fig.1, comprising the following steps:
(1) the d solid modeling function of utilizing ABAQUS software, according to simulation ballast effect push-down head, DLC film and base
The values of the structural parameters at bottom establishes pressure head geometrical model, DLC film geometrical model and substrate geometry model respectively, specific as follows:
ABAQUS software is opened, using d solid modeling function, according to simulation ballast effect push-down head, DLC film and base
It is corresponding with pressure head, DLC film and substrate several respectively to draw the part of upper, middle and lower three in a two-dimensional plane for the values of the structural parameters at bottom
What profile；In such a way that rotation generates entity, by three geometric profiles of gained respectively along 180 ° of center axis rotation, pressure head is obtained
Geometrical model, DLC film geometrical model and substrate geometry model, wherein pressure head geometrical model shape is half cone-shaped, and DLC film is several
What model and substrate geometry model are semi-cylindrical, as shown in Figure 2.
(2) grid is carried out to pressure head geometrical model, DLC film geometrical model and substrate geometry model in the step (1) to draw
Point, generate pressure head stand-alone entity, DLC film stand-alone entity and substrate stand-alone entity；Based on the bottom surface of the DLC film geometrical model,
Cohesive force unit is generated, specific as follows:
Seed is arranged to DLC film geometrical model and substrate geometry model respectively by the way of diametrically inclined in two-way, i.e.,
Among geometrical model around, size of mesh opening is gradually increased, wherein all grids of DLC film geometrical model all compare substrate geometry
The grid of model corresponding portion is finer；Using along the vertex of a cone to cone bottom inclined in two-way by the way of to pressure head geometrical model arrange kind
Son is gradually increased from the vertex of a cone to cone bottom size of mesh opening；After the completion of grid dividing, stand-alone entity is generated, it is several to be then based on DLC film
The bottom surface of what model generates cohesive force unit using deflection method, selects cohesive force cell type (COH3D8) and definition unit is deleted
It removes, i.e. the Stiffness Deterioration amount (SDEG) of cohesive force unit reaches 1, indicates that unit is entirely ineffective, forms crackle.Each section geometry mould
Corresponding stand-alone entity is as shown in Figure 3 to generation for the grid dividing of type.
(3) to the material properties of pressure head geometrical model, DLC film geometrical model and substrate geometry model in the step (2)
It is defined, selects cohesive force constitutive relation and failure parameter, specific as follows:
According to the physical parameter of DLC film, material properties, failure mode and the critical parameters of DLC film geometrical model are defined, it is fixed
The material properties of adopted substrate geometry model；Select the linear constitutive relation of cohesive force, definition interfaces rigidity, boundary strength and viscosity system
Number, definition interfaces failure mode and critical parameters define the material properties of cohesive force unit.
(4) pressure head geometrical model, DLC film geometrical model and substrate geometry model in the step (2) are assembled, is obtained
To three-dimensional finite element model, the crack propagation region of the three-dimensional finite element model is defined, specific as follows:
Pressure head geometrical model, DLC film geometrical model and substrate geometry model are successively assembled in accordance with the order from top to bottom,
Surface-to-surface contact is established between pressure head geometrical model and the upper surface of DLC film geometrical model, in the following table of DLC film geometrical model
Binding constraints is established between face and substrate geometry model, realizes the foundation of three-dimensional finite element model；By entire DLC film geometrical model
It is defined as crack propagation region, as shown in Figure 4.
(5) according to the actual conditions of DLC film and loaded-up condition, analysis mode, boundary condition and loading method are determined, and will
Load substep is applied in the step (4) in three-dimensional finite element model, specific as follows:
According to the actual conditions and loaded-up condition of DLC film, using static Gneral analysis mode, and the quantity of incremental step is determined
And size；Staff cultivation is applied to substrate geometry model, pressure head geometrical model is constrained along the movement and rotation in other directions, leaves behind
Movement along the vertical direction；According to experiment condition, load substep is applied on pressure head geometrical model.
(6) stress analysis is carried out to three-dimensional finite element model in the step (5), utilizes the post-processing function of finite element software
The relationship of energy assumed (specified) load and pressure depth, generation and propagation behavior of the simulation DLC film Interface Crack with film internal fissure, probes into influence and splits
The factor of line extension, specific as follows:
Stress analysis is carried out to the three-dimensional finite element model, by the post-processing function assumed (specified) load of ABAQUS software and
Deep relationship is pressed, as shown in Figure 5；By reading stress distribution situation, the initiating events of film internal fissure are judged, such as Fig. 6 (a) institute
Show；By reading cohesive force element stiffness amount of degradation, failure and the Interface Crack initiating events of cohesive force unit, such as Fig. 6 are determined
(b) shown in；After the completion of analysis, Interface Crack and film internal fissure are completed to extend, and form biggish crackle with interface in DLC film, such as
Shown in Fig. 7；The length of the film internal fissure changes with time situation when being extended according to film internal fissure, seeks Interface Crack and film
The variation relation of length between internal fissure adjusts the parameters such as DLC film thickness and elasticity modulus, influences DLC to further determine that
The factor of film crack propagation.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (8)
1. a kind of finite element method for simulating hard coat underbead crack and Interface Crack, comprising the following steps:
(1) the d solid modeling function of utilizing finite element software, according to simulation ballast effect push-down head, hard coat and substrate
Values of the structural parameters, establish pressure head geometrical model, hard coat geometrical model and substrate geometry model respectively；
(2) grid is carried out to pressure head geometrical model, hard coat geometrical model and substrate geometry model in the step (1) respectively
It divides, generates pressure head stand-alone entity, hard coat stand-alone entity and substrate stand-alone entity；Based on the hard coat geometrical model
Bottom surface, generate cohesive force unit；
(3) material properties of hard coat geometrical model and substrate geometry model in the step (2) are defined, in selection
Poly- power constitutive relation and failure parameter；
(4) pressure head geometrical model, hard coat geometrical model and substrate geometry model in the step (2) are assembled, is obtained
To three-dimensional finite element model, the crack propagation region of the three-dimensional finite element model is defined in hard coat geometrical model；
(5) according to the actual conditions of hard coat and loaded-up condition, analysis mode, boundary condition and loading method are determined, and will
Load substep is applied in the step (4) in three-dimensional finite element model；
(6) stress analysis is carried out to three-dimensional finite element model in the step (5), utilizes the post-processing function meter of finite element software
The relationship of load and pressure depth, the generation of analysis hard coat Interface Crack and film internal fissure and propagation behavior are calculated, influence is probed into and splits
The factor of line extension.
2. the method according to claim 1, wherein the shape of pressure head geometrical model is half in the step (1)
The shape of cone, hard coat geometrical model and substrate geometry model is semi-cylindrical.
3. method according to claim 1 or 2, which is characterized in that pressure head geometrical model, hard apply in the step (1)
The method for building up of layer geometrical model and substrate geometry model, comprising the following steps:
Using the d solid modeling function of finite element software, according to simulation ballast effect push-down head, hard coat and substrate
It is corresponding with pressure head, hard coat and substrate several respectively to draw the part of upper, middle and lower three in a two-dimensional plane for values of the structural parameters
What profile；In such a way that rotation generates entity, by three geometric profiles of gained respectively along 180 ° of center axis rotation, pressure head is obtained
Geometrical model, hard coat geometrical model and substrate geometry model.
4. the method according to claim 1, wherein the step (2) specifically:
Respectively to hard coat geometrical model and substrate geometry model arrangement seed, use by the way of diametrically inclined in two-way
Grid dividing is carried out to pressure head geometrical model arrangement seed to the mode of cone bottom inclined in two-way along the vertex of a cone, generates stand-alone entity；So
Bottom surface afterwards based on hard coat geometrical model generates cohesive force unit using deflection method, selects cohesive force cell type and determines
Adopted element deletion.
5. the method according to claim 1, wherein the step (3) specifically:
According to the physical parameter of hard coat, material properties, failure mode and the critical parameters of hard coat geometrical model are defined,
Define the material properties of substrate geometry model；Select the linear constitutive relation of cohesive force, definition interfaces rigidity, boundary strength and viscosity
Coefficient, definition interfaces failure mode and critical parameters define the material properties of cohesive force unit.
6. the method according to claim 1, wherein assembling the side of three-dimensional finite element model in the step (4)
Method, comprising the following steps:
Pressure head geometrical model, hard coat geometrical model and substrate geometry model are successively assembled in accordance with the order from top to bottom,
Surface-to-surface contact is established between pressure head geometrical model and the upper surface of hard coat geometrical model, in hard coat geometrical model
Lower surface and substrate geometry model between establish binding constraints.
7. the method according to claim 1, wherein the step (5) specifically:
According to the actual conditions and loaded-up condition of hard coat, using static Gneral analysis mode, and the quantity of incremental step is determined
And size；Staff cultivation is applied to substrate geometry model, pressure head geometrical model is constrained along the movement and rotation in other directions, leaves behind
Movement along the vertical direction；According to experiment condition, load substep is applied on pressure head geometrical model.
8. the method according to claim 1, wherein the step (6) specifically:
Stress analysis is carried out to three-dimensional finite element model, utilizes the pass that the post-processing function assumed (specified) load and pressure of finite element software are deep
System: by reading stress distribution situation, judge the generation and spread scenarios of film internal fissure；It is moved back by reading cohesive force element stiffness
Change amount judges the failure of cohesive force unit and the generation of Interface Crack and spread scenarios；Film when being extended according to film internal fissure
The length of internal fissure changes with time situation, calculates the variation relation of length between Interface Crack and film internal fissure, adjusts hard
Matter coating layer thickness and elasticity modulus, to further determine that the factor for influencing hard coat crack propagation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810971796.3A CN109190225B (en) | 2018-08-24 | 2018-08-24 | Finite element analysis method for simulating internal cracks and interface cracks of hard coating |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810971796.3A CN109190225B (en) | 2018-08-24 | 2018-08-24 | Finite element analysis method for simulating internal cracks and interface cracks of hard coating |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109190225A true CN109190225A (en) | 2019-01-11 |
| CN109190225B CN109190225B (en) | 2022-11-18 |
Family
ID=64919477
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810971796.3A Active CN109190225B (en) | 2018-08-24 | 2018-08-24 | Finite element analysis method for simulating internal cracks and interface cracks of hard coating |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109190225B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111859663A (en) * | 2020-07-20 | 2020-10-30 | 哈尔滨工业大学 | Design method of anti-icing and deicing coating on large-area coating surface |
| CN112184857A (en) * | 2020-08-14 | 2021-01-05 | 杭州群核信息技术有限公司 | Data generation system |
| CN112651150A (en) * | 2020-11-30 | 2021-04-13 | 成都大学 | Al (aluminum)2O3Propagation finite element simulation analysis method for multiple crack behaviors on surface of 316L stainless steel tritium resistance system |
| CN112651154A (en) * | 2020-12-11 | 2021-04-13 | 成都大学 | Al with rough substrate2O3Propagation finite element simulation method for multiple cracks on surface of/316L stainless steel tritium resistance system |
| CN113094956A (en) * | 2021-04-12 | 2021-07-09 | 大连理工大学 | Multilayer hard coating damping loss factor calculation method considering interlayer dissipation |
| CN113971354A (en) * | 2021-10-15 | 2022-01-25 | 厦门大学 | Finite element simulation method for nano composite material bending fracture damage research |
| CN117524361A (en) * | 2023-10-24 | 2024-02-06 | 天津大学 | Fatigue crack initiation prediction method for cycle cohesive force model with variable accumulation length |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030024323A1 (en) * | 2001-06-25 | 2003-02-06 | Jy-An Wang | Fracture toughness determination using spiral-grooved cylindrical specimen and pure torsional loading |
| CN103345545A (en) * | 2013-06-14 | 2013-10-09 | 北京航空航天大学 | Method for predicating composite material Pi-shaped non-planar glue joint strength based on triangular envelopes |
| CN103605865A (en) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | Delamination propagation behavior simulation method of composite material multidirectional laminated plate on basis of cohesion model |
-
2018
- 2018-08-24 CN CN201810971796.3A patent/CN109190225B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030024323A1 (en) * | 2001-06-25 | 2003-02-06 | Jy-An Wang | Fracture toughness determination using spiral-grooved cylindrical specimen and pure torsional loading |
| CN103345545A (en) * | 2013-06-14 | 2013-10-09 | 北京航空航天大学 | Method for predicating composite material Pi-shaped non-planar glue joint strength based on triangular envelopes |
| CN103605865A (en) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | Delamination propagation behavior simulation method of composite material multidirectional laminated plate on basis of cohesion model |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111859663A (en) * | 2020-07-20 | 2020-10-30 | 哈尔滨工业大学 | Design method of anti-icing and deicing coating on large-area coating surface |
| CN112184857A (en) * | 2020-08-14 | 2021-01-05 | 杭州群核信息技术有限公司 | Data generation system |
| CN112184857B (en) * | 2020-08-14 | 2023-09-29 | 杭州群核信息技术有限公司 | Data generation system |
| CN112651150A (en) * | 2020-11-30 | 2021-04-13 | 成都大学 | Al (aluminum)2O3Propagation finite element simulation analysis method for multiple crack behaviors on surface of 316L stainless steel tritium resistance system |
| CN112651150B (en) * | 2020-11-30 | 2022-09-02 | 成都大学 | Al (aluminum) 2 O 3 Propagation finite element simulation analysis method for multiple crack behaviors on surface of 316L stainless steel tritium resistance system |
| CN112651154A (en) * | 2020-12-11 | 2021-04-13 | 成都大学 | Al with rough substrate2O3Propagation finite element simulation method for multiple cracks on surface of/316L stainless steel tritium resistance system |
| CN113094956A (en) * | 2021-04-12 | 2021-07-09 | 大连理工大学 | Multilayer hard coating damping loss factor calculation method considering interlayer dissipation |
| CN113094956B (en) * | 2021-04-12 | 2024-05-17 | 大连理工大学 | Multilayer hard coating damping loss factor calculation method considering interlayer dissipation |
| CN113971354A (en) * | 2021-10-15 | 2022-01-25 | 厦门大学 | Finite element simulation method for nano composite material bending fracture damage research |
| 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 |
|---|---|
| CN109190225B (en) | 2022-11-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109190225A (en) | A kind of finite element method for simulating hard coat underbead crack and Interface Crack | |
| CN104200039B (en) | A kind of diaclase occurrence quantitative forecasting technique | |
| CN110555229B (en) | Meshless solid mechanics simulation method, electronic equipment and storage medium | |
| CN104965947B (en) | A kind of bridge three-dimensional design method based on BIM technology | |
| Albaric et al. | Crustal rheology and depth distribution of earthquakes: Insights from the central and southern East African Rift System | |
| CN100549660C (en) | Three-dimensional cracks expansion analogy method based on embedded type finite element model | |
| CN111488703B (en) | Method for simulating fracture toughness of layered rock based on expanded finite element | |
| CN106644868B (en) | A kind of measuring method of two dimension not convex Random Aggregate ambient interfaces concentration | |
| CN108090292A (en) | A kind of width string fan blade two dimensional finite element modeling method | |
| CN106980735B (en) | Numerical simulation method for thermal fracture of brittle material | |
| CN113223158A (en) | Blast hole charging amount calculation method based on blasting rock three-dimensional solid model | |
| CN111159794A (en) | Geometric damage rheological analysis method for mechanical properties of multi-fracture rock sample | |
| CN109064560A (en) | Geological conditions modeling method during a kind of constructing metro tunnel | |
| CN111563333B (en) | Landslide numerical analysis method considering influence of soil body strain on strength | |
| CN115034161B (en) | Self-adaptive time step calculation method for stable three-dimensional hydraulic fracture expansion calculation and acceleration | |
| Wang et al. | Study on the influence of hole defects with different shapes on the mechanical behavior and damage law of coal and rock | |
| CN112711889B (en) | Dynamic heat-fluid-solid simulation method based on LS-DYNA software | |
| CN109472046A (en) | Complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method | |
| CN110321569A (en) | It is a kind of to plug the method for numerical simulation influenced on Adjacent Pile Foundation suitable for shoe | |
| CN111894545B (en) | Method for determining proppant pumping scheme | |
| CN114491730A (en) | Dynamic stability analysis iteration method and device for high-speed railway roadbed structure | |
| CN109359358A (en) | It is a kind of for the side channel geometry optimization design method with side channel compact tension specimen | |
| Mau et al. | Arch dam system identification using vibration test data | |
| CN115235928A (en) | Method for approximately obtaining metal elastic-plastic parameters based on Knoop indenter indentation method | |
| CN110532669A (en) | A method of it is modeled for Machine Joint Surfaces contact stiffness |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |