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.

Improvements of cohesive zone model on artificial compliance and discontinuous force

The cohesive zone model(CZM)has been used widely and successfully in fracture propagation,but some basic problems are still to be solved.In this paper,artificial compliance and discontinuous force in CZM are investigated.First,theories about the cohesive element(local coordinate system,stiffness matrix,and internal nodal force)are presented.The local coordinate system is defined to obtain local separation;the stiffness matrix for an eight-node cohesive element is derived from the calculation of strain energy;internal nodal force between the cohesive element and bulk element is obtained from the principle of virtual work.Second,the reason for artificial compliance is explained by the effective stiffnesses of zero-thickness and finite-thickness cohesive elements.Based on the effective stiffness,artificial compliance can be completely removed by adjusting the stiffness of the finite-thickness cohesive element.This conclusion is verified from 1D and 3D simulations.Third,three damage evolution methods(monotonically increasing effective separation,damage factor,and both effective separation and damage factor)are analyzed.Under constant unloading and reloading conditions,the monotonically increasing damage factor method without discontinuous force and healing effect is a better choice than the other two methods.The proposed improvements are coded in LS-DYNA user-defined material,and a drop weight tear test verifies the improvements.

A Coastal Zone Segmentation Variational Model and Its Accelerated ADMM Method

Effective and efficient SAR image segmentation has a significant role in coastal zone interpretation. In this paper, a coastal zone segmentation model is proposed based on Potts model. By introducing edge self-adaption parameter and modifying noisy data term, the proposed variational model provides a good solution for the coastal zone SAR image with common characteristics of inherent speckle noise and complicated geometrical details. However, the proposed model is difficult to solve due to to its nonlinear, non-convex and non-smooth characteristics. Followed by curve evolution theory and operator splitting method, the minimization problem is reformulated as a constrained minimization problem. A fast alternating minimization iterative scheme is designed to implement coastal zone segmentation. Finally, various two-stage and multiphase experimental results illustrate the advantage of the proposed segmentation model, and indicate the high computation efficiency of designed numerical approximation algorithm.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

The efficiency and accuracy of probability diagram, spatial statistic and fractal methods in the identification of shear zone gold mineralization: a case study of the Saqqez gold ore district,NW Iran

In this study,geochemical anomaly separation was carried out with methods based on the distribution model,which includes probability diagram(MPD),fractal(concentration-area technique),and U-statistic methods.The main objective is to evaluate the efficiency and accuracy of the methods in separation of anomalies on the shear zone gold mineralization.For this purpose,samples were taken from the secondary lithogeochemical environment(stream sediment samples)on the gold mineralization in Saqqez,NW of Iran.Interpretation of the histograms and diagrams showed that the MPD is capable of identifying two phases of mineralization.The fractal method could separate only one phase of change based on the fractal dimension with high concentration areas of the Au element.The spatial analysis showed two mixed subpopulations after U=0 and another subpopulation with very high U values.The MPD analysis followed spatial analysis,which shows the detail of the variations.Six mineralized zones detected from local geochemical exploration results were used for validating the methods mentioned above.The MPD method was able to identify the anomalous areas higher than 90%,whereas the two other methods identified 60%(maximum)of the anomalous areas.The raw data without any estimation for the concentration was used by the MPD method using aminimum of calculations to determine the threshold values.Therefore,the MPD method is more robust than the other methods.The spatial analysis identified the detail soft hegeological and mineralization events that were affected in the study area.MPD is recommended as the best,and the spatial U-analysis is the next reliable method to be used.The fractal method could show more detail of the events and variations in the area with asymmetrical grid net and a higher density of sampling or at the detailed exploration stage.

Sintering zone prediction in direct metal laser sintering by finite element method

A three-dimensional finite element thermal model in direct metal laser sintering(DMLS) including the effect of powder-to-solid transition were established to predict sintering zone, which benefited the determination of suitable process parameters in DMLS. The nonlinear transient model of the metals thermal conductivity for powder-to-solid transition was developed. The model uses solid thermal properties of material in liquid-phase zone, transitional ones in sintering or sintered zone and powder ones in unsintered zones of powder bed to predict, respectively. Sintering zone boundary was estimated by maximum temperature history profile. Experiments were carried out using multi-component Cu-based metal powder. Compared experimental and predicted results, the mean error of sintering depth and width are 7.8% and 14.4%, respectively, which confirms the accuracy of the FEM prediction.

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.

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.

Towards fully automatic modelling of the fracture process in quasi-brittle and ductile materials:a unified crack growth criterion

Fully automatic finite element(FE) modelling of the fracture process in quasi-brittle materials such as concrete and rocks and ductile materials such as metals and alloys,is of great significance in assessing structural integrity and presents tre-mendous challenges to the engineering community. One challenge lies in the adoption of an objective and effective crack propagation criterion. This paper proposes a crack propagation criterion based on the principle of energy conservation and the cohesive zone model(CZM) . The virtual crack extension technique is used to calculate the differential terms in the criterion. A fully-automatic discrete crack modelling methodology,integrating the developed criterion,the CZM to model the crack,a simple remeshing procedure to accommodate crack propagation,the J2 flow theory implemented within the incremental plasticity framework to model the ductile materials,and a local arc-length solver to the nonlinear equation system,is developed and im-plemented in an in-house program. Three examples,i.e.,a plain concrete beam with a single shear crack,a reinforced concrete(RC) beam with multiple cracks and a compact-tension steel specimen,are simulated. Good agreement between numerical predictions and experimental data is found,which demonstrates the applicability of the criterion to both quasi-brittle and ductile materials.

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers:A Case Study for the Tonga Subducting Slab

For the purpose of investigating the influence of metastable olivine（MO） phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near～410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab：（i） slab pull force exerts above a depth of～230 km;（ii） MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;（iii） as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;（iv） a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.

Large Eddy Simulation for Wave Breaking in the Surf Zone

In this paper, the large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads to the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, Such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale model is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this paper has been used to study the propagation of the solitary wave in constant water depth and the shoaling of the non-breaking solitary wave on a beach. The model is based on large eddy simulation, and to track free-surface movements, the Tokyo University Modified Marker and Cell (TUMMAC) method is employed. In order to ensure the accuracy of each component of this wave mathematical model, several steps have been taken to verify calculated solutions; with either analytical solutions or experimental data. For non-breaking waves, very accurate results are obtained for a solitary wave propagating over a constant depth and on a beach. Application of the model to cnoidal wave breaking in the surf zone shows that the model results are in good agreement with analytical solution and experimental data. From the present model results, it can be seen that the turbulent eddy viscosity increases from the bottom to the water surface in surf zone. In the eddy viscosity curve, there is a turn-point obviously, dividing water depth into two parts, in the upper part, the eddy viscosity becomes very large near the wave breaking position.

Fracture development around deep underground excavations: Insights from FDEM modelling

Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method （FDEM） code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone （EDZ） around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Translational - torsional coupled model based nonlinear dynamic analysis of an NGW planetary gear train

This paper aims to investigate the nonlinear dynamic behaviors of an NGW planetary gear train with multi-clearances and manufacturing/assembling errors. For this purpose, an analytical translational- torsional coupled dynamic model is developed considering the effects of time-varying stiffness, gear backlashes and component errors. Based on the proposed model, the nonlinear differential equations of motion are derived and solved iteratively by the Runge-Kutta method. An NGW planetary gear reducer with three planets is taken as an example to analyze the effects of nonlinear factors. The results indicate that the backlashes induce complicated nonlinear dynamic behaviors in the gear train. With the increment of the backlashes, the gear system has experienced periodic responses, quasi-periodic response and chaos responses in sequence. When the planetary gear system is in a chaotic motion state, the vibration amplitude increases sharply, causing severe vibration and noise. The present study provides a fundamental basis for design and parameter optimization of NGW planetary gear trains.

Stability of Stochastic Logistic Model with Ornstein-Uhlenbeck Process for Cell Growth of Microorganism in Fermentation Process

Current research is concerned with the stability of stochastic logistic equation with Ornstein-Uhlenbeck process. First, this research proves that the stochastic logistic model with Ornstein-Uhlenbeck process has a positive solution. After that, it also introduces the sufficient conditions for stochastically stability of stochastic logistic model for cell growth of microorganism in fermentation process for positive equilibrium point by using Lyapunov function. In addition, this research establishes the sufficient conditions for zero solution as mentioned in Appendix A due to the cell growth of microorganism μmax, which cannot be negative in fermentation process. Furthermore, for numerical simulation, current research uses the 4-stage stochastic Runge-Kutta (SRK4) method to show the reality of the results.

New Fuzzy Fractional Epidemic Model Involving Death Population

In this research,we propose a new change in classical epidemic models by including the change in the rate of death in the overall population.The existing models like Susceptible-Infected-Recovered(SIR)and Susceptible-Infected-Recovered-Susceptible(SIRS)include the death rate as one of the parameters to estimate the change in susceptible,infected and recovered populations.Actually,because of the deficiencies in immunity,even the ordinary flu could cause death.If people’s disease resistance is strong,then serious diseases may not result in mortalities.The classical model always assumes a closed system where there is no new birth or death,no immigration or emigration,while in reality,such assumptions are not realistic.Moreover,the classical epidemic model does not report the change in population due to death caused by a disease.With this study,we try to incorporate the rate of change in the population of death caused by a disease,where the model is framed to reduce the curve of death along with the susceptible and infected populations.Since the rate of change turned out to be very small,we have tried to estimate it fractionally.Thus,the model is defined using fuzzy logic and is solved by two different methods:a Laplace Adomian decomposition method(LADM)and a differential transform method(DTM)for an arbitrary order α.To test its accuracy,we compared the results of both DTM and LADM with the fourth-order Runge-Kutta method(RKM-4)at α=1.

Numerical Simulation of Two-Dimensional Dendritic Growth Using Phase-Field Model

In this article, we study the phase-field model of solidification for numerical simulation of dendritic crystal growth that occurs during the casting of metals and alloys. Phase-field model of solidification describes the physics of dendritic growth in any material during the process of under cooling. The numerical procedure in this work is based on finite difference scheme for space and the 4th-order Runge-Kutta method for time discretization. The effect of each physical parameter on the shape and growth of dendritic crystal is studied and visualized in detail.

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers

Kinematic models compute the temperature distribution by prescribing a constant convergent velocity for the subducting slab,resulting in an artificial velocity discontinuity,which may accelerate the heating of subducting slabs.For the purpose to moderate the influence of such artificial discontinuity, we construct a 2D thermal model for subduction zones with a velocity boundary layer,within which the velocities decrease linearly with the distance from the interfaces of slabs.Temperatures are calculated

Acoustic Monitoring of Anomaly Stressed Zones, Determination Their Positions, Surfaces, Evaluation of Catastrophic Risk

Abstract Problem statement. Self-organization is not a universal property of matter, it exists under certain internal and external conditions and this is not associated with a special class of substances. The study of the morphology and dynamics of migration of anomalous zones associated with increased stresses is of particular importance in the development of deep deposits, complicated by dynamic phenomena in the form of rock impacts. Applied method and design: An important tool for this study is geophysical exploration. To describe the geological environment in the form of an array of rocks with its natural and technogenic heterogeneity, one should use its more adequate description, which is a discrete model of the medium in the form of a piecewise heterogeneous block medium with embedded heterogeneities of a lower rank than the block size. This nesting can be traced several times, i.e. changing the scale of the study;we see that heterogeneities of a lower rank now appear in the form of blocks for heterogeneities of the next rank. A simple averaging of the measured geophysical parameters can lead to distorted ideas about the structure of the medium and its evolution. Typical results: We have analyzed the morphology of the structural features of disintegration zones before a strong dynamic phenomenon. The introduction of the proposed integrated passive and active geophysical monitoring into the mining system, aimed at studying the transient processes of the redistribution of stress-strain and phase states, can help prevent catastrophic dynamic manifestations during the development of deep-located deposits. Concluding note (Practical value/implications): Active geophysical monitoring methods should be tuned to a model of a hierarchical heterogeneous environment. Iterative algorithms for 2-D modeling and interpretation for sound diffraction and a linearly polarized transverse elastic wave on the inclusion with a hierarchical elastic structure located in the J-th layer of the N-layer elastic medium are constructed.

Studying the Impact of Vaccination Strategy and Key Parameters on Infectious Disease Models

In the current work, we study two infectious disease models and we use nonlinear optimization and optimal control theory which helps to find strategies towards transmission control and to forecast the international spread of the infectious diseases. The relationship between epidemiology, mathematical modeling and computational tools lets us to build and test theories on the development and fighting with a disease. This study is motivated by the study of epidemiological models applied to infectious diseases in an optimal control perspective. We use the numerical methods to display the solutions of the optimal control problems to find the effect of vaccination on these models. Finally, global sensitivity analysis LHS Monte Carlo method using Partial Rank Correlation Coefficient (PRCC) has been performed to investigate the key parameters in model equations. This present work will advance the understanding about the spread of infectious diseases and lead to novel conceptual understanding for spread of them.

基于Roe-Siegmund循环内聚力模型焊趾疲劳裂纹萌生仿真

为建立焊趾疲劳裂纹萌生行为仿真方法,基于Roe-Siegmund循环内聚力模型对ABAQUS进行二次开发,形成用于反映疲劳累积损伤的VUMAT子程序,依据文献方法及试验数据获取Q345焊接区域材料的内聚力参数,通过Voronoi图法和内聚力单元法生成具有疲劳累积损伤特性和晶粒特征的微观模型,并与宏观对接焊缝模型合并,进行多尺度疲劳裂纹萌生仿真模拟,获得裂纹萌生路径、临界循环次数及晶粒力学响应特征.结果表明,此方法能够自发地选择符合实际情况的裂纹萌生位置以及短裂纹扩展路径,完成对焊接区域材料微观断裂过程的模拟,不同仿真组获取的临界循环次数存在于一定的分布范围内,模型中的累积内聚力长度需通过试验数据进行拟合.