Discrete Element Modelling of Damage Evolution of Concrete Considering Meso-Structure of ITZ

The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.

Stress-corrosion coupled damage localization induced by secondary phases in bio-degradable Mg alloys:phase-field modeling

In this study,a phase-field scheme that rigorously obeys conservation laws and irreversible thermodynamics is developed for modeling stress-corrosion coupled damage(SCCD).The coupling constitutive relationships of the deformation,phase-field damage,mass transfer,and electrostatic field are derived from the entropy inequality.The SCCD localization induced by secondary phases in Mg is numerically simulated using the implicit iterative algorithm of the self-defined finite elements.The quantitative evaluation of the SCCD of a C-ring is in good agreement with the experimental results.To capture the damage localization,a micro-galvanic corrosion domain is defined,and the buffering effect on charge migration is explored.Three cases are investigated to reveal the effect of localization on corrosion acceleration and provide guidance for the design for resistance to SCCD at the crystal scale.

Revised damage evolution equation for high cycle fatigue life prediction of aluminum alloy LC4 under uniaxial loading

The fatigue life prediction for components is a difficult task since many factors can affect the final fatigue life. Based on the damage evolution equation of Lemaitre and Desmorat, a revised two-scale damage evolution equation for high cycle fatigue is presented according to the experimental data, in which factors such as the stress amplitude and mean stress are taken into account. Then, a method is proposed to obtain the material parameters of the revised equation from the present fatigue experimental data. Finally, with the utilization of the ANSYS parametric design language （APDL） on the ANSYS platform, the coupling effect between the fatigue damage of materials and the stress distribution in structures is taken into account, and the fatigue life of specimens is predicted. The outcome shows that the numerical prediction is in accord with the experimental results, indicating that the revised two-scale damage evolution model can be well applied for the high cycle fatigue life prediction under uniaxial loading.

Numerical study of fatigue damage of asphalt concrete using cohesive zone model

In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model （CZM） is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.

Aero-engine Blade Fatigue Analysis Based on Nonlinear Continuum Damage Model Using Neural Networks

Fatigue life and reliability of aero-engine blade are always of important significance to flight safety.The establishment of damage model is one of the key factors in blade fatigue research.Conventional linear Miner＇s sum method is not suitable for aero-engine because of its low accuracy.A back propagation neutral network（BPNN） based on the combination of Levenberg-Marquardt（LM） and finite element method（FEM） is used to describe process of nonlinear damage accumulation behavior in material and predict fatigue life of the blade.Fatigue tests of standard specimen made from TC4 are carried out to obtain material fatigue parameters and S-N curve.A nonlinear continuum damage model（CDM）,based on the BPNN with one hidden layer and ten neurons,is built to investigate the nonlinear damage accumulation behavior,in which the results from the tests are used as training set.Comparing with linear models and previous nonlinear models,BPNN has the lowest calculation error in full load range.It has significant accuracy when the load is below 500 MPa.Especially,when the load is 350 MPa,the calculation error of the BPNN is only 0.4%.The accurate model of the blade is built by using 3D coordinate measurement technology.The loading cycle in fatigue analysis is defined from takeoff to cruise in 10 min,and the load history is obtained from finite element analysis（FEA）.Then the fatigue life of the compressor blade is predicted by using the BPNN model.The final fatigue life of the aero-engine blade is 6.55 104 cycles（10 916 h） based on the BPNN model,which is effective for the virtual design of aero-engine blade.

Numerical modeling of fracture propagation of supercritical CO_(2)compound fracturing

The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.

Analysis of the damage mechanism of strainbursts by a global-local modeling approach

Strainburst is the most common type of rockbursts.The research of strainburst damage mechanisms is helpful to improve and optimize the rock support design in the burst-prone ground.In this study,an improved global-local modeling approach was first adopted to study strainburst damage mechanisms.The extracted stresses induced by multiple excavations from a three-dimensional(3D)global model established by fast Lagrangian analysis of continua in 3 dimensions(FLAC3D)are used as boundary conditions for a two-dimensional(2D)local model of a deep roadway built by universal distinct element code(UDEC)to simulate realistic stress loading paths and conduct a detailed analysis of rockburst damage from both micro and macro perspectives.The results suggest that the deformation and damage level of the roadway gradually increase with the growth of surrounding rock stress caused by the superposition of mining-or excavation-induced stresses of the panel and nearby roadways.The significant increase of surrounding rock stresses will result in more accumulated strain energy in two sidewalls,providing a necessary condition for the strainburst occurrence in the dynamic stage.The strainburst damage mechanism for the study site combines three types of damage:rock ejection,rock bulking,and rockfall.During the strainburst,initiation,propagation,and development of tensile cracks play a crucial role in controlling macroscopic failure of surrounding rock masses,although the shear crack always accounts for the main proportion of damage levels.The deformation and damage level of the roadway during a strainburst positively correlate with the increasing peak particle velocities(PPVs).The yielding steel arch might not dissipate kinetic energy and mitigate strainburst damage effectively due to the limited energy absorption capacity.The principles to control and mitigate strainburst damage are proposed in this paper.This study presents a systematic framework to investigate strainburst damage mechanisms using the global-local modeling approach.

Modeling asymmetric fracture mechanics of Mg alloy wire in drawing process

In this study,a numerical analysis was conducted on the ductile fracture of a 2-mm diameter Mg-1Zn-0.5Mn-0.5Sr-0.1Ca alloy wire during drawing.The hexagonally close-packed crystal structure of Mg alloys causes asymmetric fracture behavior,especially in the compression region.The aim of this study is to develop a comprehensive damage model for Mg alloy wire that accurately predicts ductile fracture,with a focus on the compression region.A novel experimental method was introduced to measure the ductile fracture of Mg alloy wires under different stress states.The wire drawing process was simulated using the Generalized Incremental Stress-State dependent damage(GISSMO)Model and the Semi-Analytical Model for Polymers(SAMP)model.The damage model's prediction and the experimental results were found to be in excellent agreement,especially in determining crack initiation.Computational analysis established a safe zone diagram for die angle and reduction ratio,and experimental validation confirmed the feasibility of this approach.The proposed damage model can provide a practical and reliable analysis for optimizing the drawing process of Mg alloy wire.

Damage Detection in Reinforced Concrete Berthing Jetty Using a Plasticity Model Approach

A conventional method of damage modeling by a reduction in stiffness is insufficient to model the complex non-linear damage characteristics of concrete material accurately.In this research,the concrete damage plasticity constitutive model is used to develop the numerical model of a deck beam on a berthing jetty in the Abaqus finite element package.The model constitutes a solid section of 3D hexahedral brick elements for concrete material embedded with 2D quadrilateral surface elements as reinforcements.The model was validated against experimental results of a beam of comparable dimensions in a cited literature.The validated beam model is then used in a three-point load test configuration to demonstrate its applicability for preliminary numerical evaluation of damage detection strategy in marine concrete structural health monitoring.The natural frequency was identified to detect the presence of damage and mode shape curvature was found sensitive to the location of damage.

A Coupling Model of the Discontinuous Deformation Analysis Method and the Finite Element Method

Neither the finite element method nor the discontinuous deformation analysis method can solve problems very well in rock mechanics and engineering due to their extreme complexities. A coupling method combining both of them should have wider applicability. Such a model coupling the discontinuous deforma- tion analysis method and the finite element method is proposed in this paper. In the model, so-called line blocks are introduced to deal with the interaction via the common interfacial boundary of the discontinuous deformation analysis domain with the finite element domain. The interfacial conditions during the incre- mental iteration process are satisfied by means of the line blocks. The requirement of gradual small dis- placements in each incremental step of this coupling method is met through a displacement control proce- dure. The model is simple in concept and is easy in numerical implementation. A numerical example is given. The displacement obtained by the coupling method agrees well with those obtained by the finite ele- ment method, which shows the rationality of this model and the validity of the implementation scheme.

Numerical simulation of aircraft arresting process with an efficient full-scale rigid-flexible coupling dynamic model

The arresting process of carrier-based aircraft is widely recognized as a challenging task,characterized by the highest accident rate among all carrier-based aircraft operations.Dynamic simulation plays a crucial role in assessing the intricate responses of the arresting process,favoring the design of carrier-based aircraft.An efficient and accurate rigid-flexible coupling model for analyzing the dynamic response of the arresting process is proposed.By combining the dynamic characteristics of airframe,landing gear,arresting hook and arresting gear system,the rigid-flexible coupling dynamic model is established to reflect the relative motion of the coupling parts and arresting load.The dynamic model is verified through simulations of landing gear landing drops and by comparing the arresting simulation results with corresponding data in the US military standard.Additionally,simulations of the arresting process under different off-center distance and aircraft yaw angle are conducted to obtain the dynamic response of the aircraft during the arresting process.The result indicates that the rigid-flexible coupling dynamic model proposed is effective for analyzing the arresting dynamics response of carrier-based aircraft.The axial force of the arresting cable on both sides of the hook engagement point,pitch and yaw angle of aircraft are inconsistent under yaw and off-center arresting.The analysis method and obtained results provide valuable references for assessing the dynamic responses of carrier-based aircraft during arresting process and offer valuable in-sights in the design of carrier-based aircraft.

Numerical modeling for the coupled thermo-mechanical processes and spalling phenomena in sp Pillar Stability Experiment (APSE)

In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both continuum and discontinuum based numerical methods. Two-dimensional （2D） and three- dimensional （3D） finite element method （FEM） and 2D distinct element method （DEM） with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone （EDZ） during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.

基于SBFEM的剪切型裂纹动态开裂模拟

剪切型断裂是岩土工程中常见的破坏模式,了解剪切破坏机理并准确预测剪切型裂纹的萌生、扩展过程对保障工程结构的安全性与稳定性具有重要意义.文章建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)和非局部宏-微观损伤模型的剪切型裂纹动态开裂模拟方法,定义了基于偏应变概念的物质点对的正伸长量,可作为预测剪切型裂纹扩展行为的动态开裂准则,一点的损伤定义为该点影响域范围内连接的物质键损伤的加权平均值,而物质键的损伤则与基于偏应变概念的物质点对的正伸长量相关联,并引入能量退化函数建立结构域几何拓扑损伤与能量损失之间的关系,将拓扑损伤与应力应变联系起来,通过能量退化函数修正了SBFEM的刚度系数矩阵,得到了子域在损伤状态下的刚度矩阵,推导了考虑结构损伤的SBFEM动力控制方程,采用Newmark隐式算法对控制方程进行时间离散.最后,通过3个典型算例验证了建议的模型可较好地模拟剪切型断裂问题,能够很好地捕捉剪切型裂纹的扩展路径,并得到较为准确的载荷-位移曲线.

考虑气体的渗流应力耦合模型在黄壁庄水库中的应用

目前建立的多孔介质耦合模型通常认为气体是不可压缩的,无法描述气体的输运过程,无法反映气体对水体入渗和固体变形的影响。为了描述多孔介质中固液气三相的渗流应力耦合过程,基于多孔介质理论和有效应力原理,提出了一种耦合了改进守恒水平集(ICLS)方法并考虑气体可压缩性的固液气三相渗流应力耦合模型。将该方法应用于黄壁庄水库坝前冒泡现象的渗流应力耦合分析中,模拟了水位骤升后,固体变形、水体入渗和气体运移过程,及其耦合效应。计算结果表明,该方法能够简单直观的描述水气两相的运移过程及渗流应力耦合过程,且模拟结果与实际观测结果基本一致。

A novel two-way method for dynamically coupling a hydrodynamic model with a discrete element model(DEM)

The effect of floating objects has so far been little considered for hazard risk assessment and structure design, despite being an important factor causing structural damage in flood-prone and coastal areas. In this work, a novel two-way method is proposed to fully couple a shock-capturing hydrodynamic model with a discrete element model （DEM） for simulation of complex debris-enriched flow hydrodynamics. After being validated against an idealized analytical test, the new coupled model is used to reproduce flume experiments of floating debris driven by dam-break waves. The numerical results agree satisfactorily with the experimental measurements, demonstrating the model＇s capability and efficiency in simulating complex fluid-debris interactions induced by violent shallow flows.

Multi-physics coupling field finite element analysis on giant magnetostrictive materials smart component

This study presents a new method to solve the difficult problem of precise machining a non-cylinder pinhole of a piston using embedded giant magnetostrictive material (GMM) in the component. We propose the finite element model of GMM smart component in electric, magnetic, and mechanical fields by step computation to optimize the design of GMM smart com-ponent. The proposed model is implemented by using COMSOL multi-physics V3.2a. The effects of the smart component on the deformation and the system resonance frequencies are studied. The results calculated by the model are in excellent agreement (relative errors are below 10%) with the experimental values.

Finite element implementation of coupled hydro-mechanical modeling of transversely isotropic porous media in DEAL.II

In geotechnical engineering,modeling geo-structures is challenging,particularly in cases where the interaction between the structures and soil or rock is complex.Most wellknown commercial modeling software is based on homogenous and isotropic materials.However,soil and rock are often modeled in heterogeneous and anisotropic media because of the inherent anisotropy of sedimentary rock masses and their stratified structure.In recent decades,coupled hydro-mechanical(HM)interactions in isotropic porous media have been studied;however,the behavior of transversely isotropic porous media is rarely considered.In addition,it is difficult for commercial software such as Plaxis and Flac3D to express complex rock formation where the anisotropy of the material and the associated cracks and fractures could be assembled into a single model.In this study,a finite element implementation using Differential Equation Analysis Library(DEAL.II),an open-source library of finite element codes,was developed to model the fully coupled HM behavior of transversely isotropic porous media.The proposed implementation can be applied to both isotropic and transversely isotropic porous media based on Biot’s theory.The developed code can be used to model poroelastic media with(1)equations of linear elasticity for the solid matrix and(2)diffusion equations for fluid flow based on mass and linear-momentum conservation laws.We verified the performance and accuracy of the code through two examples,i.e.,Mandel’s problem with a compared analytical solution and a tunnel excavation process with the Flac3D software.On the basis of these numerical applications,we present the code to model the behavior of various geo-structures such as tunnels and pile–soil interactions with anisotropic materials.

钛合金室温受压蠕变损伤本构模型

为了解决工程上在评估深潜器耐压壳体安全可靠性时不考虑拉压蠕变破坏机理产生的差异,会导致较大计算误差问题,本文提出了适用于室温蠕变的蠕变损伤本构模型。通过分析钛合金材料蠕变特性,考虑了室温环境下位错堆积造成的蠕变阻力,并且位错堆积的现象与时间有关。基于微分自洽法提出了适用于室温环境下的钛合金蠕变本构模型,得到受压情况下TC4ELI的相关材料参数,并将该模型用USDFLD和CREEP子程序进行定义。研究表明:对于受压结构而言,使用通过拉伸蠕变实验得到的本构模型的计算结果会过于保守,对于含V型缺口平板受压时等效蠕变应变的相对误差为168.20%。本文所提出的模型适用于环肋耐压壳结构的蠕变损伤分析。

考虑材料参数空间变异性的SBFEM开裂模拟

论文建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)框架的非局部宏微观损伤模型,考虑材料细观物理参数的空间变异性,探讨了材料参数的空间变异性对结构开裂过程的影响。结果表明:考虑材料参数空间变异性后,裂纹扩展路径具有不确定性,建议的模型能够很好地反应材料内在的随机性;随着结构受力情况的复杂化和结构本体缺陷的增多,裂纹开裂模式的变异性也会增大。自相关长度和参数变异系数对结构开裂分析结果有重要影响。

基于损伤模型的混凝土冻融循环过程数值模拟研究

冻融损伤是寒区混凝土主要破坏方式之一,冻融过程是混凝土材料逐渐破坏演变的过程。引入损伤理论建立了混凝土水热力损伤耦合模型,该模型可以模拟混凝土冻融过程中水的相变、孔隙水压等变化过程,通过累积损伤真实反映破坏过程。通过与相关文献中的试验数据对比,验证了模型的准确性,并研究了中心最低温度、降温速率以及初始孔隙率对混凝土试件冻融损伤的影响规律。结果表明:随着试件中心最低温度的降低、降温速率的增大和初始孔隙率的增大,混凝土的累积损伤值逐渐增大;损伤区虽然从混凝土试件四周逐渐向内发展,但在降温速率和初始孔隙率过大时,最大损伤区和裂缝却出现在棱角附近的试样内部,并向外部扩展,该结果解释了棱柱体混凝土冻融破坏后棱角最先掉落的现象。