CN107506595A - Material damage prediction method based on critical condition dynamic mobile - Google Patents
Material damage prediction method based on critical condition dynamic mobile Download PDFInfo
- Publication number
- CN107506595A CN107506595A CN201710765954.5A CN201710765954A CN107506595A CN 107506595 A CN107506595 A CN 107506595A CN 201710765954 A CN201710765954 A CN 201710765954A CN 107506595 A CN107506595 A CN 107506595A
- Authority
- CN
- China
- Prior art keywords
- stress
- equation
- research object
- area
- proportion
- 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
- 230000006378 damage Effects 0.000 title claims abstract description 67
- 239000000463 material Substances 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012544 monitoring process Methods 0.000 claims abstract description 8
- 230000006835 compression Effects 0.000 claims description 8
- 238000007906 compression Methods 0.000 claims description 8
- 238000010008 shearing Methods 0.000 claims description 6
- 230000003542 behavioural effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011439 discrete element method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention provides a kind of material damage prediction method based on critical condition dynamic mobile, comprises the following steps：Based on step 1. is judged by theoretical research, numerical analysis and scene, the critical state point or line of present situation material damage are determined, and determines destruction area and non-demolition area, by measuring study, establishes geometric description equation；The proportion distribution characteristics of step 2. analysis and research object, establishes proportion equation；Step 3. analyzes the boundary condition stress characteristics in destruction area and non-demolition area, establishes boundary condition stress equation；The different zones that research object is divided into rule by step 4. are studied, and analyze both continuously and discontinuously Physical Mechanical Quantities, establish corresponding Physical Mechanical Quantities method for solving；Step 5. chooses stress and represents equation, and calculates each constant coefficient；Step 6. bond material peak strength criterion, constitutive equation and field monitoring situation, according to critical condition dynamic mobile rule, implement the prediction of material damage.
Description
Technical field
The invention belongs to mechanics, civil engineering and Geological Engineering etc. and power and deform the technical fields of Related Computational Methods,
More particularly to a kind of material damage prediction method based on critical condition dynamic mobile.
Technical background
Natural calamity, the mankind build structures destruction etc. and cause property and human loss in various degree, thus to nature
The generation of the destructions such as disaster, it is very important from implementation prediction on the basis of power and deformation.Existing stress and deformation are asked
Solution is often established on the numerical computation methods such as finite elements, and Finite Element is for two-dimensional structure using the calculating fanned out from point to area
Mode, three-dimensional structure is used to put the calculation for body, therefore for different size of unit, result of calculation is often not
Identical；When carrying out numerical computations, generally use linear method solves nonlinear problem, i.e. initial stress method, due to using not
Same convergence, the result for causing to calculate also can be different.
The content of the invention
The present invention is, and it is an object of the present invention to provide a kind of be based on critical condition dynamic mobile in order to solving the above problems and carry out
Material damage prediction method, force boundary condition, dynamic balance condition and critical condition breakdown point stress condition etc. can be based on,
Implement the prediction of material damage.
The present invention to achieve these goals, employs following scheme：
The present invention provides a kind of material damage prediction method based on critical condition dynamic mobile, it is characterised in that bag
Include following steps：Based on step 1. is judged by theoretical research, numerical analysis and scene, the critical shape of present situation material damage is determined
State point or line (respectively for two dimension or three-dimensional), and destruction area and non-demolition area are determined, while measuring study object (contains destruction area
With non-demolition area) macroscopical geometric properties, establish the geometric description equation corresponding with macroscopical geometric properties；Step 2. analysis is ground
Study carefully the proportion distribution characteristics of object, establish proportion equation of the research object in survey region；Step 3. analysis and research object destroys
Area and the boundary condition stress characteristics in non-demolition area, establish the boundary condition stress side corresponding with boundary condition stress characteristics
Journey；Research object is divided into regular different zones and studied (such as by step 4. according to the feature of research object：For two
Dimension：Triangle, quadrangle etc., for three-dimensional：Tetrahedron, hexahedron etc.), particularly destruction area and non-demolition area are divided
Open, and analyze both continuously and discontinuously that Physical Mechanical Quantities are (such as：Back zone is destroyed for the thrust load caused landslide based on being destroyed with compression shear, its
Shearing stress and driving both downslide stress are discontinuous etc., but the compression of vertical sliding surface and back-pressure stress are continuous), foundation
Discontinuous Physical Mechanical Quantities method for solving is (such as：For landslide, due to shearing stress and drive both downslide stress discontinuous, then
Shearing stress can be according to patent：A kind of computational methods of progressive disruption of slope potential water use, the patent No.：
201510658880.6 solve)；In Physical Mechanical Quantities continuum, corresponding stress etc. must be continuous；Step 5. is chosen should
Power represents equation, and stress represents that equation meets corresponding power, stress and (or) the moment of flexure of research object destruction area and non-demolition area
Equilibrium equation (see：A kind of stress solving method based on force boundary and equilibrium condition, the patent No.：201611034900.3)；Facing
Boundary's state point or line (respectively for two dimension or three-dimensional), meet the corresponding peak strength criterion of material, and calculate each constant coefficient；
Step 6. bond material peak strength criterion is (such as：Mole coulomb criterion, Griffth criterions etc.), in that case it can be decided that after research object
Continuous Damage course, that is, determine Latent destruction face；It is right with reference to constitutive equation on the basis of to mechanical characteristic labor
The deformation behaviour of research object is studied, and with reference to field monitoring, is calculated deformation and is analyzed with field monitoring, it is determined that grinding
Study carefully the behavioural characteristic of object, according to critical condition dynamic mobile rule, implement the prediction of material damage.The present invention is in vacation
If research object stress meets border and equilibrium condition, breakdown point meets criterion of strength, so as to the destruction for research object
Journey implements prediction, can be with stress of the prediction research object in random geometry, deformation and failure distribution characteristics.The party
The it is proposed of method by dam, bridge, side slope, roadbed, house, tunnel, be mingled with, tunnel, the various construction structures such as culvert move
Static state plus unloading fracture implementation Process prediction.
Material damage prediction method provided by the invention based on critical condition dynamic mobile, can also have such
Feature：In step 1, in the case where research object is D geometric modeling, studied with Two-dimensional Theoretical, numerical analysis and scene
Based on judgement, the critical state point of present situation material damage is determined, and determines destruction area and non-demolition area, while measuring study pair
Macroscopical geometric properties of elephant, establish the geometric description equation corresponding with macroscopical geometric properties.
Material damage prediction method provided by the invention based on critical condition dynamic mobile, can also have such
Feature：In step 2, in the case where research object is D geometric modeling, the proportion of survey region includes the ratio of X-direction
Weight γw,xWith the proportion γ of Y directionw,y；In the case where research object is three-dimensional geometry, the proportion of survey region includes
The proportion γ of X-directionw,x, Y direction proportion γw,yWith the proportion γ of Z-directionw,z。
Material damage prediction method provided by the invention based on critical condition dynamic mobile, can also have such
Feature：In step 3, if destruction area crushes, according to broken boundary condition；Otherwise destruction area and non-demolition area are all abided by
According to boundary condition feature.
Material damage prediction method provided by the invention based on critical condition dynamic mobile, can also have such
Feature：In step 4, back zone is destroyed for the thrust load caused landslide based on being destroyed with compression shear, its shearing stress glides with driving should
Both power is discontinuous, but the compression of vertical sliding surface and back-pressure stress are continuous.
The effect of invention and effect
1st, dynamic mobile prediction methods of the invention, its stress represent that equation meets research object destruction area and non-broken
The corresponding power of bad block, stress and (or) moment of flexure equilibrium equation；In critical state point or line, meet the corresponding peak strength of material
Criterion；
2nd, prediction methods of the invention：The main stress for solving the discontinuous problems such as destruction area stress, strain reason
By solution, and according to the prediction of critical condition dynamic mobile implementation material damage；
3rd, prediction method of the invention to dam, bridge, side slope, roadbed, house, tunnel, be mingled with, tunnel, culvert etc.
The sound state destructive process of various construction structures implements prediction, has impetus to correlative study and application.
Brief description of the drawings
Fig. 1 is the critical condition movement dynamically prediction schematic diagram of the material damage involved by the embodiment of the present invention.
Embodiment
Below in conjunction with accompanying drawing to the material damage prediction method of the present invention based on critical condition dynamic mobile
Specific embodiment is described in detail.
It is pre- to disclose a kind of material damage prediction based on critical condition dynamic mobile so that two dimension comes down as an example for the present embodiment
Report method, specifically comprises the following steps：
Step 1. is studied (see Chinese invention patent with Two-dimensional Theoretical：It is a kind of to be made every effort to based on force boundary and equilibrium condition
Solution, the patent No.：201611034900.3), numerical analysis (such as：Existing Finite Element, discrete element method etc.) and scene sentence
It is set to basis, determines the critical state point (D points in such as Fig. 1) of present situation material damage, and determines destruction area (such as Tu1Zhong ABDE areas)
With non-demolition area (such as Tu1Zhong BCD areas), while macroscopical geometric properties of measuring study object (containing destruction area and non-demolition area),
Establish the geometric description equation corresponding with macroscopical geometric properties (such as：In Fig. 1：AB, BC, EA be linear equation (such as：For straight
Line BC：Y=KBCx+bBC, KBC,bBCRespectively slope and intercept), DE is curvilinear equation)；
Step 2. research object is D geometric modeling, then the proportion of survey region includes the proportion γ of X-directionw,xAnd Y
The proportion γ of direction of principal axisw,yFollowing (the note of equation：Other representations can also be taken) be：
γw,x=γ0,x+a4,1x+a4,2y+a4,3x2+a4,4xy+a4,5y2+a4,6x3+a4,7x2y+a4,8xy2+…(1)
γw,y=γ0,y+a5,1x+a5,2y+a5,3x2+a5,4xy+a5,5y2+a5,6x3+a5,7x2y+a5,8xy2+…(2)
Step 3. is according to claim 3, the feature in object destruction area of analyzing and researching and non-demolition area, if destruction area is broken
It is broken, then according to broken boundary condition；Otherwise destruction area and non-demolition area are all in accordance with boundary condition feature, such as：Meet holy text
Southern principle, and establish the boundary condition stress equation corresponding with boundary condition stress characteristics.
If destruction area (ABDE in such as Fig. 1) is broken,：AB boundary conditions stress can be expressed as：
σxx|AB=0, σyy|AB=0, τxy|AB=0 (3)
If destruction area (ABDE in such as Fig. 1) is completely,：
AB boundary conditions stress may be：
σxx|AB=0, σyy|AB=0, τxy|AB=0 (4)
Or
Or
AE boundary conditions stress may be：
σxx|AE=0, σyy|AE=0, τxy|AE=0 (7)
Or
Or
Rather than destruction area (BCD in such as Fig. 1), then：
BC boundary conditions stress may be：
Wherein：
σN BC=l2σxx BC+m2σyy BC+2lmτxy BC (11)
In formula：L, m be BC faces outside normal orientation cosine, σxx BC,σyy BC,τxy BCFor the boundary stress on BC sides,
σxx,σyy,τxy：Respectively X, Y-direction direct stress and shear stress；η：Along straight line BC unit length, i.e.,
Step 4. inquires into ABCDE feature, ABCDE is divided into two research objects of ABDE and BCD, wherein ABDE is edge
Sliding surface destruction area, BCD are non-demolition area, discontinuous along both the shear stress of DE sliding surfaces and friction resistance to sliding stress in destruction area, when
Right displacement is also discontinuous, and other Physical Mechanical Quantities are continuous.
In analysis, both continuously and discontinuously Physical Mechanical Quantities are (such as：After being destroyed for the thrust load caused landslide based on being destroyed with compression shear
Area, its shearing stress and driving both downslide stress are discontinuous etc., but the compression of vertical sliding surface and back-pressure stress are continuous)
On the basis of, discontinuous frictional resistance shear stress is according to Chinese invention patent ----((number of patent application：201610860012.0 hair
Bright name of patent application：A kind of progressive disruption of slope overall process New calculating method or number of patent application：2014100250810, hair
Bright name of patent application：The slide face boundary that a kind of slope stability calculates) calculated.
Step 5. chooses stress and represents equation, and the stress represents that equation meets that research object destruction area and non-demolition area are relative
The power answered, stress and (or) moment of flexure equilibrium equation (see：A kind of stress solving method based on force boundary and equilibrium condition, the patent No.：
201611034900.3)；In critical state point or line (respectively for two dimension or three-dimensional), meet that the corresponding peak strength of material is accurate
Then, and each constant coefficient is calculated；Its step is as follows：
For ABDE research objects, by coordinate XOY in Fig. 1, the power expression formula of stress 3 is chosen：
σxx=a1,0+a1,1x+a1,2y+a1,3x2+a1,4xy+a1,5y2+a1,6x3+a1,7x2y+a1,8xy2+a1,9y3 (13)
σyy=a2,0+a2,1x+a2,2y+a2,3x2+a2,4xy+a2,5y2+a2,6x3+a2,7x2y+a2,8xy2+a2,9y3 (14)
τxy=a3,0+a3,1x+a3,2y+a3,3x2+a3,4xy+a3,5y2+a3,6x3+a3,7x2y+a3,8xy2+a3,9y3 (15)
Above-mentioned stress must is fulfilled for equilibrium equation, in the situation that proportion is constant (i.e.：γw,x=0, γw,y=γ) under, from
And eliminate 15 constant coefficients.Remaining 15 constant determineds.
Using equation (3), 9 constant coefficients, remaining 6 constant coefficients can be eliminated.
UtilizeWith research object ABDE in X and Y
Direction of principal axis dynamic balance, and meet existing peak stress criterion (such as in the stress of D points：Mole coulomb criterion, Griffth are accurate
Then), then research object ABDE stress and strain solution can be obtained, certainly, for research object ABDE, its torque should also be put down
Weighing apparatus.
For research object BCD, by X'O'Y' coordinates in Fig. 1, the power expression formula of stress 3 is chosen：
σxx'=b1,0+b1,1x'+b1,2y'+b1,3x'2+b1,4x'y'+b1,5y'2+b1,6x'3+b1,7x'2y'+b1,8x'y'2+
b1,9y'3 (16)
σyy'=b2,0+b2,1x'+b2,2y'+b2,3x'2+b2,4x'y'+b2,5y'2+b2,6x'3+b2,7x'2y'+b2,8x'y'2+
b2,9y'3 (17)
τxy'=b3,0+b3,1x'+b3,2y'+b3,3x'2+b3,4x'y'+b3,5y'2+b3,6x'3+b3,7x'2y'+b3,8x'y'2+
b3,9y'3 (18)
Stress using research object ABDE and BCD on BD sides is equal, in that case it can be decided that 12 constant coefficients, remaining 6 are not
Know constant coefficient.Utilize y'=0, σ 'y=0 or y'=0, σ 'x=0 may decide that 3 constant coefficients, meet existing peak value using D points
Stress criterion is (such as：Mole coulomb criterion, Griffth criterions), while using in formula (10)：
In two, research object can be obtained
BCD stress solutions.
It was found from above-mentioned solution procedure, it, which is solved, has more solution features.
Step 6. bond material peak strength criterion is (such as：Mole coulomb criterion, Griffth criterions etc.), in that case it can be decided that grind
The follow-up destruction angle of object is studied carefully, with reference to displacement monitoring, in that case it can be decided that subsequent burst path, that is, determine Latent destruction face；
On the basis of mechanical characteristic labor, with reference to constitutive equation, the deformation behaviour of research object can be studied, with reference to
Field monitoring, calculate deformation and be analyzed with field monitoring, determine the behavioural characteristic of research object, moved according to critical condition
State movement law, implement the prediction of material damage.
Above example is only to put to the proof, and is not the restriction to embodiments of the present invention.In addition to the implementation, originally
Invention also has other embodiment.All technical schemes formed using equivalent substitution or equivalent transformation, all fall within application claims
Protection domain.
Claims (5)
1. a kind of material damage prediction method based on critical condition dynamic mobile, it is characterised in that comprise the following steps：
Step 1. by theoretical research, numerical analysis and scene judge based on, determine present situation material damage critical state point or
Line, and destruction area and non-demolition area are determined, while macroscopical geometric properties of measuring study object, establish and macroscopical geometric properties phase
Corresponding geometric description equation；
The proportion distribution characteristics of step 2. analysis and research object, establishes proportion equation of the research object in survey region；
Step 3. analysis and research object destruction area and the boundary condition stress characteristics in non-demolition area, establish special with boundary condition stress
Levy corresponding boundary condition stress equation；
Research object is divided into regular different zones and studied, particularly will by step 4. according to the feature of research object
Destruction area and non-demolition area are separated, and analyze both continuously and discontinuously Physical Mechanical Quantities, are established discontinuous Physical Mechanical Quantities and are asked
Solution method；In Physical Mechanical Quantities continuum, corresponding stress must be continuous；
Step 5. chooses stress and represents equation, and the stress represents that equation meets that research object destruction area and non-demolition area are corresponding
Power, stress or moment of flexure equilibrium equation；In critical state point or line, meet the corresponding peak strength criterion of material, and calculate
Each constant coefficient；
Step 6. bond material peak strength criterion, in that case it can be decided that the follow-up Damage course of research object, that is, determine potential broken
Bad face；On the basis of to mechanical characteristic labor, with reference to constitutive equation, the deformation behaviour of research object is studied,
With reference to field monitoring, calculate deformation and be analyzed with field monitoring, the behavioural characteristic of research object is determined, according to critical shape
State dynamic mobile rule, implement the prediction of material damage.
2. the material damage prediction method according to claim 1 based on critical condition dynamic mobile, it is characterised in that：
In the step 1, research object be D geometric modeling in the case of, studied with Two-dimensional Theoretical, numerical analysis and
Based on scene judges, the critical state point of present situation material damage is determined, and determines destruction area and non-demolition area, while measurement is ground
Study carefully macroscopical geometric properties of object, establish the geometric description equation corresponding with macroscopical geometric properties.
3. the material damage prediction method according to claim 1 based on critical condition dynamic mobile, it is characterised in that：
In the step 2, in the case where research object is D geometric modeling, the proportion of survey region includes X-direction
Proportion γw,xWith the proportion γ of Y directionw,y；In the case where research object is three-dimensional geometry, the proportion of survey region
Proportion γ including X-directionw,x, Y direction proportion γw,yWith the proportion γ of Z-directionw,z。
4. the material damage prediction method according to claim 1 based on critical condition dynamic mobile, it is characterised in that：
In the step 3, if destruction area crushes, according to broken boundary condition；Otherwise destruction area and non-demolition area are complete
Portion is in accordance with boundary condition feature.
5. the material damage prediction method according to claim 1 based on critical condition dynamic mobile, it is characterised in that：
In the step 4, back zone is destroyed for the thrust load caused landslide based on being destroyed with compression shear, its shearing stress glides with driving
Both stress is discontinuous, but the compression of vertical sliding surface and back-pressure stress are continuous.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710765954.5A CN107506595B (en) | 2017-08-30 | 2017-08-30 | material damage prediction method based on critical state dynamic movement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710765954.5A CN107506595B (en) | 2017-08-30 | 2017-08-30 | material damage prediction method based on critical state dynamic movement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107506595A true CN107506595A (en) | 2017-12-22 |
| CN107506595B CN107506595B (en) | 2019-12-10 |
Family
ID=60694447
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710765954.5A Expired - Fee Related CN107506595B (en) | 2017-08-30 | 2017-08-30 | material damage prediction method based on critical state dynamic movement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107506595B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101487846B1 (en) * | 2014-01-24 | 2015-02-02 | 서울대학교산학협력단 | Method for analyzing 2D shallow water flow through assignment of wetting/drying condition using critical drying depth method |
| CN105335607A (en) * | 2015-10-12 | 2016-02-17 | 湖北工业大学 | Slope progressive failure latent slide surface calculating method |
| CN106503421A (en) * | 2016-09-28 | 2017-03-15 | 西南交通大学 | A kind of progressive disruption of slope overall process computational methods |
| CN106599382A (en) * | 2016-11-23 | 2017-04-26 | 湖北工业大学 | Stress solution method based on force boundary and balance conditions |
-
2017
- 2017-08-30 CN CN201710765954.5A patent/CN107506595B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101487846B1 (en) * | 2014-01-24 | 2015-02-02 | 서울대학교산학협력단 | Method for analyzing 2D shallow water flow through assignment of wetting/drying condition using critical drying depth method |
| CN105335607A (en) * | 2015-10-12 | 2016-02-17 | 湖北工业大学 | Slope progressive failure latent slide surface calculating method |
| CN106503421A (en) * | 2016-09-28 | 2017-03-15 | 西南交通大学 | A kind of progressive disruption of slope overall process computational methods |
| CN106599382A (en) * | 2016-11-23 | 2017-04-26 | 湖北工业大学 | Stress solution method based on force boundary and balance conditions |
Non-Patent Citations (4)
| Title |
|---|
| GAOHUI WANG.ET.: "Damage prediction of concrete gravity dams subjected to underwater explosion shock loading", 《ENGINEERING FAILURE ANALYSIS》 * |
| YINGFA LU: "Deformation and failure mechanism of slope in three dimensions", 《JOURNAL OF ROCK MECHANICS AND GEOTECHNICAL ENGINEERING》 * |
| 卢应发等: "边坡力的分布特征和稳定性分析", 《岩土工程学报》 * |
| 徐颖: "强降雨作用下类土质滑坡演化过程及破坏机理研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107506595B (en) | 2019-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106599382B (en) | Stress solution method based on force boundary and balance condition | |
| Morris et al. | Simulations of fracture and fragmentation of geologic materials using combined FEM/DEM analysis | |
| Zhu et al. | High rock slope stability analysis using the enriched meshless Shepard and least squares method | |
| Gupta et al. | Dynamic analysis of axially loaded end-bearing pile in a homogeneous viscoelastic soil | |
| Ramezani et al. | Analysis of micropolar elastic beams | |
| Dehghanpoor et al. | Soil-pile-superstructure interaction effects in seismically isolated bridges under combined vertical and horizontal strong ground motions | |
| Pham et al. | Limit load space of rigid footing under eccentrically inclined load | |
| Latha et al. | Elasto-plastic analysis of jointed rocks using discrete continuum and equivalent continuum approaches | |
| He et al. | Numerical estimation and prediction of stress-dependent permeability tensor for fractured rock masses | |
| Mohamed | Three-dimensional fully parameterized triangular plate element based on the absolute nodal coordinate formulation | |
| Zhang et al. | Verification of a DDA-based hydro-mechanical model and its application to dam foundation stability analysis | |
| Leonel et al. | Non-linear boundary element formulation with tangent operator to analyse crack propagation in quasi-brittle materials | |
| Vo-Minh et al. | Upper bound limit analysis of circular tunnel in cohesive-frictional soils using isogeometric analysis based on Bézier extraction | |
| Karim et al. | Numerical analysis of piled raft foundation on clayey soil | |
| Shu et al. | A surface-to-surface scheme for 3D contact problems by boundary face method | |
| An et al. | Arbitrary discontinuities in the numerical manifold method | |
| Tao et al. | Failure and stability analysis of Jinping-I arch dam based on geomechanical model test and nonlinear numerical analysis | |
| Arslan et al. | Finite element simulation of localization in granular materials by micropolar continuum approach | |
| CN107506595A (en) | Material damage prediction method based on critical condition dynamic mobile | |
| Junior et al. | Brittle fracture and hydroelastic simulations based on moving particle simulation | |
| Tang et al. | Low-order mixed finite element analysis of progressive failure in pressure-dependent materials within the framework of the Cosserat continuum | |
| Bakuckas Jr | Comparison of boundary correction factor solutions for two symmetric cracks in a straight-shank hole | |
| de Figueiredo et al. | Analysis of bookshelf mechanisms using the mechanics of Cosserat generalized continua | |
| Haghani et al. | Comparative study of smeared crack and extended finite element method for predicting the crack propagation in concrete gravity dams | |
| Xiao et al. | Absolute nodal coordinate formulation for dynamic analysis of reinforced concrete structures |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191210 |