CN105046737B - A kind of method for numerical simulation generating Three-dimensional damage spot figure - Google Patents
A kind of method for numerical simulation generating Three-dimensional damage spot figure Download PDFInfo
- Publication number
- CN105046737B CN105046737B CN201510222376.1A CN201510222376A CN105046737B CN 105046737 B CN105046737 B CN 105046737B CN 201510222376 A CN201510222376 A CN 201510222376A CN 105046737 B CN105046737 B CN 105046737B
- Authority
- CN
- China
- Prior art keywords
- spot
- stress
- dimensional
- crit
- numerical simulation
- 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.)
- Active
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The method for numerical simulation provided by the present invention for generating Three-dimensional damage spot figure, 3-D geometric model is being established to object described in the Three-dimensional damage spot figure, subdivision carefully sees grid, set the boundary condition of the numerical simulation, on the basis of primary condition and material parameter, choose the yield surface with angle point, consider the constitutive relation of strength softening, non-associated flow law, using finite element dynamic relaxation, numerical simulation is carried out to geometrical model and is solved, export Three-dimensional damage spot figure, to directly from mechanical equation, it does not need first to assume any physical mechanism, produce Three-dimensional damage spot figure, it is a kind of method for directly generating Three-dimensional damage spot figure, actual conditions of the spot figure generated closer to material damage, to the mechanical property of more intuitive and true reflection material.
Description
Technical field
The invention belongs to material numerical analogue technique fields, and in particular to a kind of Numerical-Mode for generating Three-dimensional damage spot figure
Quasi- method.
Background technique
For material, for load-bearing solid-state material, load damage or fatigue damage are to determine the material longevity
An important factor for life.Spot figure, be it is a kind of space or on the time with certain regularity non-homogeneous macrostructure figure, and
The non-linear macrograph of solution and simulation material damage Emergence and Development, end.Research material damages the process of spot figure growth, for
Understanding Important Project dynamics catastrophe is bred, the mechanism that occurs is with especially important basic meaning.
When loading to material, solid material is damaged, and is initially loaded, the damage spot figure for describing material damage is substantially
Uniformly, with the repetition of load or temporal extension or the increase of load, the damage of material is continued to develop, correspondingly, damage
Blemish figure gradually forms the high damage zone in part, and damage from substantially evenly changing to localized failure in damage spot figure
It is more likely to develop in high damage zone.From the development process of pattern dynamics angle analysis material damage, localized failure goes out
Now derived from the unstability of uniformly damage spot figure, various damage spot figures are ultimately formed, the unstability of spot figure, the pole of respective material intensity are damaged
Limit value, thus the damage of forming material.
Damaging spot map generalization is a complicated non-linear process.In the prior art, the spot figure commonly side of description is damaged
Method is to use reaction-diffusion equation, but this method needs first to assume various physical mechanisms, is a kind of round-about way, gives birth to
At damage spot figure and actual damage situation have certain gap.There is presently no the direct sides for generating material three-dimensional damage spot figure
Method.
Summary of the invention
The purpose of the embodiment of the present invention is that a kind of method for numerical simulation for generating Three-dimensional damage spot figure is proposed, by mechanics side
Journey directly generates the Three-dimensional damage spot figure of material damage, so that the damage of more objective, true simulation material occurs, development, knot
The process of beam instructs the implementation of engineering and the maintenance of material.
According to an aspect of the present invention, it proposes a kind of method for numerical simulation for generating Three-dimensional damage spot figure, the sides
Method includes:
The yield surface with angle point, the constitutive relation for considering strength softening, non-associated flow law are chosen, using finite element
Dynamic relaxation carries out numerical simulation and solves output Three-dimensional damage spot figure.
In above scheme, the method also includes:
3-D geometric model is established to object described in the Three-dimensional damage spot figure, subdivision carefully sees grid；
Set the boundary condition of the numerical simulation；
Set the primary condition of the numerical simulation；
Set material parameter.
In above scheme, the progress numerical simulation simultaneously solves output Three-dimensional damage spot figure, further are as follows: to the geometry
Model carries out numerical simulation and solves.
In above scheme, the boundary condition includes at least stress boundary condition, displacement boundary conditions.
In above scheme, the primary condition is included at least: stress initial value, the calculating for being displaced initial value, numerical simulation
Step number initial value, timing initial value.
In above scheme, the setting material parameter, further are as follows: setting internal friction angle, tensile strength, is cut at cohesive strength
Swollen angle.
The method for numerical simulation provided by the present invention for generating Three-dimensional damage spot figure, is retouched to the Three-dimensional damage spot figure
The object stated establishes boundary condition, primary condition and the material that 3-D geometric model, subdivision carefully see grid, the setting numerical simulation
On the basis of expecting parameter, the yield surface with angle point, the constitutive relation for considering strength softening, non-associated flow law are chosen, is adopted
With finite element dynamic relaxation, numerical simulation is carried out to geometrical model and is solved, Three-dimensional damage spot figure is exported, thus directly from power
It learns equation to set out, does not need first to assume any physical mechanism, that is, produce Three-dimensional damage spot figure, be that one kind directly generates three-dimensional damage
The method of blemish figure, spot figure generated closer to material damage actual conditions, thus more intuitive and true reflection material
The mechanical property of material.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill of field, without any creative labor, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is the method for numerical simulation flow diagram of the generation Three-dimensional damage spot figure of first embodiment of the invention；
Fig. 2 is the method for numerical simulation flow diagram of the generation Three-dimensional damage spot figure of second embodiment of the invention；
Fig. 3 is the geometrical model damage spot figure grid chart that second embodiment of the invention is established；
Fig. 4 is that second embodiment of the invention uses the softening curve graph for considering the constitutive relation of strength softening；
Fig. 5 is damage spot figure the first calculated result figure of second embodiment of the invention；
Fig. 6 is damage spot figure the second calculated result figure of second embodiment of the invention.
Specific embodiment
Those skilled in the art of the present technique are appreciated that unless expressly stated, singular " one " used herein, " one
It is a ", " described " and "the" may also comprise plural form.It is to be further understood that being arranged used in specification of the invention
Diction " comprising " refer to that there are the feature, integer, step, operation, element and/or component, but it is not excluded that in the presence of or addition
Other one or more features, integer, step, operation, element, component and/or their group.It should be understood that when we claim member
Part is " connected " or when " coupled " to another element, it can be directly connected or coupled to other elements, or there may also be
Intermediary element.In addition, " connection " used herein or " coupling " may include being wirelessly connected or coupling.Wording used herein
"and/or" includes one or more associated any cells for listing item and all combinations.
Those skilled in the art of the present technique are appreciated that unless otherwise defined, all terms used herein (including technology art
Language and scientific term) there is meaning identical with the general understanding of those of ordinary skill in fields of the present invention.Should also
Understand, those terms such as defined in the general dictionary, which should be understood that, to be had and the meaning in the context of the prior art
The consistent meaning of justice, and unless defined as here, it will not be explained in an idealized or overly formal meaning.
In order to facilitate understanding of embodiments of the present invention, embodiments of the present invention are described below in detail, by reference to attached drawing
The embodiment of description is exemplary, and for explaining only the invention, and is not construed as limiting the claims.
Fig. 1 is the method for numerical simulation flow diagram of the generation Three-dimensional damage spot figure of first embodiment of the invention.
As shown in Figure 1, the method for numerical simulation of the generation Three-dimensional damage spot figure of the present embodiment, includes the following steps:
Step 101,3-D geometric model is established to object described in the Three-dimensional damage spot figure, subdivision carefully sees grid.
In this step, object described in the Three-dimensional damage spot figure refers to solid material under normal conditions.It is general solid
Body material in use, can bear pressure or be damaged because of fatigue, the map generalization of Three-dimensional damage spot and change here
Change process describes the Emergence and Development process of material damage.Subdivision carefully sees grid, is the refinement to the geometrical model.
Step 102, the boundary condition of the numerical simulation is set.
Boundary condition in this step is boundary condition when carrying out numerical simulation solution to the 3-D geometric model.
The boundary condition includes at least stress boundary condition, displacement boundary conditions.Preferably, boundary condition here can be by setting
The mode of three-dimensional system of coordinate is determined to realize.It can be to set boundary stress as 0, side with conditions setting in three-dimensional system of coordinate
Boundary's displacement is 0；In the case where applying pure shear stress, set practical pure shear stress initial value can also be used to answer for boundary
Power.
Step 103, the primary condition of the numerical simulation is set.
Primary condition in this step is primary condition when carrying out numerical simulation solution to the 3-D geometric model,
Be it is corresponding with boundary condition, the boundary condition is for boundary, and primary condition is for inside.Therefore, this
In primary condition be all internal relevant parameter.The primary condition include at least internal stress initial value, displacement initial value,
Calculating step number initial value, the timing initial value of numerical simulation.For example, set internal stress as 0, initial displacement is set as 0, numerical value
The calculating step number of simulation is initially 0, and timing initial value is 0, at this point, corresponding internal strain and strain rate may be 0.This
In value be only an example, in practice, need according to circumstances to be set.Initial value be only and subsequent meter
The relative value compared.
It should be noted that the boundary condition and the primary condition cannot be 0 simultaneously.For example, when setting the side
When boundary's stress is 0, it is not possible to set the internal stress initial value as 0.
Further, according to the primary condition of above-mentioned setting, primary Calculation mass matrix M.
Step 104, the yield surface with angle point is chosen.
Step 105, the constitutive relation for considering strength softening is chosen.
Step 106, non-associated flow law is chosen.
Step 107, material parameter is set.
The material parameter in this step includes: internal friction angleCohesive strength c, tensile strength sigmat, dilative angle ψ.It carries out
When setting,For internal friction angle initial value, c0For cohesive strength initial value.
Above-mentioned steps 101 to 107, out-of-order requirement are the condition selection and setting to further numerical simulation and calculating,
After completing the setting of selection and parameter of above-mentioned corresponding conditions, step 108 is executed.
Step 108, within the scope of the boundary condition, using the primary condition as starting point, using the yield surface, this structure
Relationship, non-associated flow law are carried out numerical simulation to the geometrical model and are solved using finite element dynamic relaxation.
Step 109, Three-dimensional damage spot figure is exported according to the result of the solution.
The method for numerical simulation of the generation Three-dimensional damage spot figure of the present embodiment is strong using the yield surface with angle point, consideration
The constitutive relation and non-associated flow law for spending softening, are solved by finite element dynamic slave method, thus directly from power
It learns equation to set out, does not need first to assume any physical mechanism, that is, produce Three-dimensional damage spot figure, be that one kind directly generates three-dimensional
The method for damaging spot figure, spot figure generated closer to material damage actual conditions, thus more intuitive and true reflection
The mechanical property of material.
Fig. 2 is the method for numerical simulation flow diagram of the generation Three-dimensional damage spot figure of second embodiment of the invention.
As shown in Fig. 2, the method for numerical simulation of the generation Three-dimensional damage spot figure of the present embodiment, includes the following steps:
Step 201, the cylindrical geometry model that radius is 2cm, subdivision unit are established.
This step is that geometrical model is established to actual circular cylindrical solid material, and Fig. 3 is the geometry that the present embodiment is established
MODEL DAMAGE spot figure grid chart.As shown in figure 3, first choice establishes three-dimensional system of coordinate, it is 2cm's that radius is established in the coordinate system
Cylindrical geometry model sets 0 point of x-axis direction as one of boundary condition, and 0 point of y-axis direction is another perimeter strip
Part.
Step 202, the initial value of internal stress, internal displacement, strain rate, acceleration is set as 0.
When carrying out numerical simulation, makes material and applying external force to circular cylindrical solid material deformation occurs, then not
Apply external force when, the stress point shape of material becomes 0, therefore, in this step set initial displacement be 0, material strain rate with add
Speed is also 0.This step is to set primary condition, and primary condition here refers to the interior condition of material, with the side
Boundary's condition is opposite.
Step 203, using the Mohr-Coulomb yield surface with angle point:
ft=σ3+σt (2)
Wherein, the fsFor shear strength function, ftFor tensile strength function, τ is shear stress, and c is cohesiveness, and ψ is swollen to cut
Angle, σnFor unit direct stress, σ3For minor principal stress, σtFor tensile strength.
Step 204, using non-associated flow law:
gs=τ-tan (ψ) σn (3)
gt=σ3 (4)
Wherein, gsTo shear potential function, gtTo stretch potential function,τFor shear stress, ψ is dilative angle, σnFor unit direct stress,
σ3For minor principal stress.
Step 205, using consider strength softening constitutive relation, with the inclined strain controlling c of plasticity,Softening process, modeling
The inclined strain increment of property is defined as follows:
Wherein, Δ ε1,Δε2,Δε3It is three principal strain increments.ΔγpFor the inclined strain increment of plasticity.
Fig. 4 is by the present embodiment using the softening curve graph for the constitutive relation for considering strength softening.As shown in figure 4, with
Plasticity shearing strain (γP) increase, the internal friction angle of materialTheir initial value compared with cohesive strength (c)By
It is decrescence few.
Step 206, definition material parameterc0, σt, ψ.
Preferably, the material real parameter in the present embodimentValue range be 0~80 ° ,=40 °, c0Value range be 0.5 ×
10e3Pa~10 × 10e3Pa, σtValue range is 0.5 × 10e3Pa~10 × 10e3Pa, ψ value range are 0~90 °.
Further, a special case in the present embodiment, set material parameter asc0=10e3Pa, σt=
10e3Pa, ψ=10 °.
Step 207, apply pure shear boundary condition in the cylindrical cross-section direction:
X in formula, y are boundary coordinates, and r is cylindrical radius, v0For start node speed.Preferably, v0=10e-5m/step。
Step 208, force calculation (fn,Δtcrit), wherein n is to calculate step number.
Further, this step includes:
Step A, initialization: fn=0, Δ tcrit=∞；
Step B calculates external node power fext,n；
Step C recycles all unit e, specifically includes:
Step C1, aggregation units modal displacement and speed；
Step C2, internal node power
Step C3, to integral dot cycle ξQ, it specifically includes:
Step C31 executes step 210 if n=0；
Step C32 calculates deformation: Dn-1/2(ξQ),Fn(ξQ),En(ξQ)；
Step C33, using the yield surface with angle point；
Step C34, using non-associated flow law；
Step C35, using the constitutive relation for considering strength softening；
Step C36 calculates stress σn(ξQ)；
Step C37, regeneration function power
Step C38 terminates integral dot cycle.
Step C4, computing unit external node power, fe ext,n；
Step C5, fe n=fe ext,n-fe int,n；
Step C6 calculates Δ te crit, if Δ te crit< Δ tcritSo Δ tcrit=Δ te crit；
Step C7, by fe nAssemble entirety fn；
Step D, end unit e circulation.
Step 209, acceleration a is calculatedn=M-1(fn-Cdamp v0)。
Step 210, renewal time: tn+1=tn+Δtn+1/2,tn+1/2=(tn+tn+1)/2。
Step 211, node speed: v is updated for the first timen+1/2=vn+(tn+1/2-tn)an。
Step 212, apply velocity boundary conditions: if node I belongs to
Step 213, modal displacement: d is updatedn+1=dn+Δtn+1/2vn+1/2。
Step 214, n=n+1 repeats step 208, executes step 215 after step 208.
Step 215, acceleration a is calculatedn+1。
Step 216, node speed: v is updated secondn+1=vn+1/2+(tn+1-tn+1/2)an+1。
Step 218, it calculates kinetic energy and checks energy balance.
Step 219, it updates and calculates step: n ← n+1.
Step 220, by stress σnCalculate damage.
Step 221, judge whether to complete simulation, if completing simulation, be transferred to step 222；If not completing simulation, step is returned
Rapid 208.
It, can be in such a way that active be set, as run by setting to whether the judgement simulated is completed in this step
Calculating step number is 100 steps, and after calculating 100 steps, simulation is completed.It can also such as be terminated by setting by way of passively calculating
The condition of simulation judges that simulation is completed when the condition is satisfied.Here condition can be set according to actual needs.
Step 222, Three-dimensional damage spot figure is exported.
In above-mentioned steps, σnFor direct stress,For angle of friction,For initial friction angle, C0For original cohesion, σtTo draw
Intensity is stretched, a is node acceleration, and f is nodal force, CdampFor damping matrix, tcritFor crash time step, D is deformation-rate tensor, F
For deformation gradient, E is Green strain, and B is strain transition matrix, and Q is point, fintFor internal node power, fextFor outside segments
Point power, J is Jacobian,For elastic potential energy,For boundary node rate, xIFor boundary node I coordinate.
Fig. 5 is damage spot figure the first calculated result figure of second embodiment of the invention；Fig. 6 is second embodiment of the invention
Damage the second calculated result of spot figure figure.During actual loaded, it is normally set up calculating step number, it, will according to set condition
The loading procedure is divided into several subprocess or stage, while presetting the calculating step number that each stage needs to be implemented.For example, total
5000 steps are calculated, the 1st to the 1000th step is the first stage；1001st step to 5000 steps is second stage.In the present embodiment, by institute
It states calculating process and is divided into two stages, the first stage executes presetting calculating step number, and obtains the first damage spot figure；Second-order
Duan Zhihang presetting calculating step number, and obtain the second damage spot figure.As shown in figure 5, the first damage spot figure, i.e. damage spot
Figure the first calculated result figure, the rounded spot figure shape arrangement of damage spot figure, the maximum value of damage are located at the middle part of round spot figure.It is all
Adjacent spot figure is all by reversed shear-deformable generation, and this spot figure is referred to as super lattice-like spot figure.As shown in fig. 6,
The second damage spot figure, i.e. damage the second calculated result of spot figure figure describe the continuation with load, some circle spot figures
It is elongated, rotates, aggregate into ribbon, and still maintain circular spot figure between band, referred to as superlattices turn ribbon damage
Spot figure.
The method for numerical simulation of the generation Three-dimensional damage spot figure of the present embodiment is strong using the yield surface with angle point, consideration
The constitutive relation and non-associated flow law for spending softening, are solved by finite element dynamic slave method, thus directly from power
It learns equation to set out, does not need first to assume any physical mechanism, that is, produce Three-dimensional damage spot figure, be that one kind directly generates three-dimensional damage
The method of blemish figure, spot figure generated closer to material damage actual conditions, thus more intuitive and true reflection material
The mechanical property of material.
Those of ordinary skill in the art will appreciate that: attached drawing is the schematic diagram of one embodiment, module in attached drawing or
Process is not necessarily implemented necessary to the present invention.
As seen through the above description of the embodiments, those skilled in the art can be understood that the present invention can
It realizes by means of software and necessary general hardware platform.Based on this understanding, technical solution of the present invention essence
On in other words the part that contributes to existing technology can be embodied in the form of software products, the computer software product
It can store in storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used so that a computer equipment
(can be personal computer, server or the network equipment etc.) executes the certain of each embodiment or embodiment of the invention
Method described in part.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for device or
For system embodiment, since it is substantially similar to the method embodiment, so describing fairly simple, related place is referring to method
The part of embodiment illustrates.Apparatus and system embodiment described above is only schematical, wherein the conduct
The unit of separate part description may or may not be physically separated, component shown as a unit can be or
Person may not be physical unit, it can and it is in one place, or may be distributed over multiple network units.It can root
According to actual need that some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.Ordinary skill
Personnel can understand and implement without creative efforts.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
In the technical scope disclosed by the present invention, any changes or substitutions that can be easily thought of by anyone skilled in the art,
It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with scope of protection of the claims
Subject to.
Claims (4)
1. a kind of method for numerical simulation for directly generating Three-dimensional damage spot figure, which is characterized in that the described method includes:
Step S1, establishes three-dimensional system of coordinate, establishes 3-D geometric model, subdivision to object described in the Three-dimensional damage spot figure
It is thin to see grid；
Step S2 sets 0 point of x-axis direction as one of boundary condition, and 0 point of y-axis direction is another boundary condition；
Step S3 sets the initial value of internal stress, internal displacement, strain rate, acceleration as 0；
Step S4, using the Mohr-Coulomb yield surface with angle point:
ft=σ3+σt (2)
Wherein, the fsFor shear strength function, ftFor tensile strength function, τ is shear stress, and c is cohesiveness,For interior friction
Angle, σnFor direct stress, σ3For minor principal stress, σtFor tensile strength；
Step S5, using non-associated flow law:
gs=τ-tan (ψ) σn (3)
gt=σ3 (4)
Wherein, gsTo shear potential function, gtTo stretch potential function, τ is shear stress, and ψ is dilative angle, σnFor direct stress, σ3For small master
Stress；
Step S6, using consider strength softening constitutive relation, with the inclined strain controlling c of plasticity,Softening process, plasticity answers partially
It is as follows to become increment definition:
Wherein, Δ ε1,Δε2,Δε3It is three principal strain increments, Δ γpFor the inclined strain increment of plasticity；
Step S7, definition material parameterc0, σt, ψ；
Step S8 selectes one cross-wise direction of the 3-D geometric model, and applies pure shear boundary condition:
X in formula, y are boundary coordinate, v0For start node speed；
Step S9, force calculation (fn,Δtcrit), wherein n is to calculate step number；fnFor the nodal force of the n-th step, Δ tcritIt is critical
Time step is poor；
Step S10 calculates acceleration an=M-1(fn-Cdampv0), M is mass matrix, v0For start node speed；
Step S11, renewal time: tn+1=tn+Δtn+1/2,tn+1/2=(tn+tn+1)/2；Δtn+1/2For the time step of setting；
Step S12 updates node speed: v for the first timen+1/2=vn+(tn+1/2-tn)an；
Step S13 applies velocity boundary conditions: if node I belongs to
Step S14 updates modal displacement: dn+1=dn+Δtn+1/2vn+1/2；
Step S15, n=n+1 repeat step S9, execute step S16 after step S9；
Step S16 calculates acceleration an+1；
Step S17, second of update node speed: vn+1=vn+1/2+(tn+1-tn+1/2)an+1；
Step S18 calculates kinetic energy and checks energy balance；
Step S19 updates and calculates step: n ← n+1；
Step S20, by stress σnCalculate damage；
Step S21 judges whether to complete simulation, if completing simulation, is transferred to step S22；If not completing simulation, return step
S9；
Step S22 exports Three-dimensional damage spot figure；
Wherein, in above-mentioned steps formula, σnFor direct stress,For angle of friction,For initial friction angle, c0For original cohesion, σt
For tensile strength, a is node acceleration, and f is nodal force, CdampFor damping matrix, tcritIt is walked for the crash time,For boundary section
Point velocity, xIFor boundary node I coordinate.
2. the method for numerical simulation according to claim 1 for directly generating Three-dimensional damage spot figure, which is characterized in that described
Step S9 further comprises:
Step A, initialization: fn=0, Δ tcrit=∞；
Step B calculates external node power fext,n；
Step C recycles all unit e, specifically includes:
Step C1, aggregation units modal displacement and speed are transposed matrix and rate matrices；
Step C2, internal node power fe int,n=0；
Step C3, to point ξQCirculation, specifically includes:
Step C31 executes step S11 if n=0；
Step C32 calculates deformation: Dn-1/2(ξQ),Fn(ξQ),En(ξQ)；Wherein, " Dn-1/2(ξQ) " indicate point ξQIn (n-1/
2) deformation-rate tensor of step, Fn(ξQ) " indicate point ξQDeformation gradient, En(ξQ) " indicate point ξQGreen strain；
Step C33, using the yield surface with angle point；
Step C34, using non-associated flow law；
Step C35, using the constitutive relation for considering strength softening；
Step C36 calculates stress σn(ξQ)；
Step C37, regeneration function power
Step C38 terminates integral dot cycle；
Step C4, computing unit external node power fe ext,n；
Step C5, fe n=fe ext,n-fe int,n；
Step C6 calculates Δ te critIf Δ te crit< Δ tcritSo Δ tcrit=Δ te crit；
Step C7, by fe nAssemble entirety fn；
Step D, end unit e circulation；
Wherein, D is deformation-rate tensor, and F is deformation gradient, and E is Green strain, and B is strain transition matrix, and Q is point, fint
For internal node power, fextFor external node power, J is Jacobian,For elastic potential energy.
3. the method for numerical simulation according to claim 1 or 2 for directly generating Three-dimensional damage spot figure, which is characterized in that institute
It states in step S7,Value range is 0~80 °, c0Value range is 0.5 × 10e3Pa~10 × 10e3Pa, σtValue range is
0.5×10e3Pa~10 × 10e3Pa, ψ value range are 0~90 °.
4. the method for numerical simulation according to claim 1 or 2 for directly generating Three-dimensional damage spot figure, which is characterized in that institute
State in step S21 to whether complete simulation judgement:
By way of actively setting N value, when the calculating step number of operation reaches N value, simulation is completed；
Or,
By passively setting the condition of end simulation, when the condition is satisfied, judge that simulation is completed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510222376.1A CN105046737B (en) | 2015-05-04 | 2015-05-04 | A kind of method for numerical simulation generating Three-dimensional damage spot figure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510222376.1A CN105046737B (en) | 2015-05-04 | 2015-05-04 | A kind of method for numerical simulation generating Three-dimensional damage spot figure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105046737A CN105046737A (en) | 2015-11-11 |
| CN105046737B true CN105046737B (en) | 2019-02-15 |
Family
ID=54453253
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510222376.1A Active CN105046737B (en) | 2015-05-04 | 2015-05-04 | A kind of method for numerical simulation generating Three-dimensional damage spot figure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105046737B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106649995A (en) * | 2016-11-18 | 2017-05-10 | 重庆邮电大学 | Numerical simulation research method for microscopic damage in PBX explosive particle compression moulding process |
| CN106919770A (en) * | 2017-03-21 | 2017-07-04 | 安徽理工大学 | A kind of damage variable based on numerical simulation determines method |
| CN109187843A (en) * | 2018-08-14 | 2019-01-11 | 安徽师范大学 | Rotation speed regulating and controlling prepares the method and three-dimensional spiral spot figure of three-dimensional spiral figure spirit spot figure |
| CN109443992A (en) * | 2018-09-29 | 2019-03-08 | 安徽师范大学 | The method and three-dimensional concentric ring spot figure of three-dimensional concentric ring spot figure are prepared in simple system |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102819633A (en) * | 2012-07-27 | 2012-12-12 | 哈尔滨工业大学 | Method for establishing constitutive relation between welding thermal cycle temperature and thermal deformation history material and senior nonlinear finite element analysis software (MSC. MARC) secondary development |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2471073A1 (en) * | 2003-06-17 | 2004-12-17 | Zale Tabakman | System, method and computer for communicating using interactive three-dimensional renderings |
-
2015
- 2015-05-04 CN CN201510222376.1A patent/CN105046737B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102819633A (en) * | 2012-07-27 | 2012-12-12 | 哈尔滨工业大学 | Method for establishing constitutive relation between welding thermal cycle temperature and thermal deformation history material and senior nonlinear finite element analysis software (MSC. MARC) secondary development |
Non-Patent Citations (3)
| Title |
|---|
| Three-dimensional finite-element simulation of the dynamic Brazilian tests on concrete cylinders;Gonzalo Ruiz等;《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING》;20001231;第963-994页 |
| 准脆性材料损伤斑图生长的数值模拟;王利 等;《工程力学》;20110430;第28卷(第4期);第238-244,256页 |
| 岩土工程极限分析有限元法及其应用;郑颖人 等;《土木工程学报》;20050131;第38卷(第1期);第91-98,104页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105046737A (en) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105046737B (en) | A kind of method for numerical simulation generating Three-dimensional damage spot figure | |
| Mehrabian et al. | A novel technique for optimal placement of piezoelectric actuators on smart structures | |
| Samir et al. | Damage detection in CFRP composite beams based on vibration analysis using proper orthogonal decomposition method with radial basis functions and cuckoo search algorithm | |
| Nariman-Zadeh et al. | Evolutionary design of generalized polynomial neural networks for modelling and prediction of explosive forming process | |
| Veuve et al. | Active control for mid-span connection of a deployable tensegrity footbridge | |
| Luo et al. | Topology optimization of hyperelastic structures with frictionless contact supports | |
| Talha et al. | Stochastic vibration characteristics of finite element modelled functionally gradient plates | |
| Hutson et al. | Combining laser microsurgery and finite element modeling to assess cell-level epithelial mechanics | |
| Lum et al. | Structural optimization for flexure-based parallel mechanisms–Towards achieving optimal dynamic and stiffness properties | |
| CN105488328A (en) | Fatigue crack growth rate prediction method based on artificial neuron network | |
| CN105243185B (en) | A kind of tension integral structure based on Monte Carlo method looks for shape method | |
| Raja et al. | Active control of tensegrity structures under random excitation | |
| Du et al. | Force analytic method for rolling gaits of tensegrity robots | |
| CN112329290A (en) | Finite element discrete element multi-scale coupling calculation method for construction process simulation | |
| Groh et al. | Snaking and laddering in axially compressed cylinders | |
| Chen et al. | Mechanical models of the cellular cytoskeletal network for the analysis of intracellular mechanical properties and force distributions: A review | |
| Nguyen et al. | Fast analysis and prediction approach for geometrically nonlinear bending analysis of plates and shells using artificial neural networks | |
| Srinivas et al. | Computational methodologies for vibration-based damage assessment of structures | |
| Wang et al. | Damage detection utilising the artificial neural network methods to a benchmark structure | |
| Alagappan et al. | Wave propagation due to impact through layered polymer composites | |
| Ojha et al. | Adaptive search space decomposition method for pre-and post-buckling analyses of space truss structures | |
| CN110020485A (en) | One kind being based on bolted suspension type thin-walled column shell inherent characteristic analysis method | |
| CN109918838A (en) | A kind of method and system of prediction large and complex structure military service performance containing splicing | |
| Daş et al. | Free vibration analysis of isotropic plates using regressive ensemble learning | |
| Zhang | Local bifurcations of nonlinear viscoelastic panel in supersonic flow |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |