Numerical Method for Modeling the Constitutive Relationship of Sand under Different Stress Paths

A numerical method was used in order to establish the constitutive relationship of sands under different stress paths, Firstly, based on the numerical method modeling the constitutive law of sands, the elastoplastic constitutive relationship of sand was established for three paths： the constant proportion of principle stress path, the conventional triaxial compression （CTC） path, and the p=constant （TC） path. The yield lines of plastic volumetric strain and plastic generalized shear strain were given. Through visualization, the three dimensional surface of the stress-strain relationship in the whole stress field （p, q） obtained under the three paths was plotted. Also, by comparing the stress-strain surfaces and yield locus of the three stress paths, the differences were found to be obvious, which demonstrates that the influence of the stress paths on constitutive law was not neglected. The numerical modeling method overcame the difficulty of finding an analytical expression for plastic potential. The results simulated the experimental data with an accuracy of 90% on average, so the constitutive model established in this paper provides an effective constitutive equation for this kind of engineering, reflecting the effect of practical stress paths that occur in sands.

An Effective Numerical Method for the Solution of a Stochastic Coronavirus(2019-nCovid)Pandemic Model

Nonlinear stochastic modeling plays a significant role in disciplines such as psychology,finance,physical sciences,engineering,econometrics,and biological sciences.Dynamical consistency,positivity,and boundedness are fundamental properties of stochastic modeling.A stochastic coronavirus model is studied with techniques of transition probabilities and parametric perturbation.Well-known explicit methods such as Euler Maruyama,stochastic Euler,and stochastic Runge–Kutta are investigated for the stochastic model.Regrettably,the above essential properties are not restored by existing methods.Hence,there is a need to construct essential properties preserving the computational method.The non-standard approach of finite difference is examined to maintain the above basic features of the stochastic model.The comparison of the results of deterministic and stochastic models is also presented.Our proposed efficient computational method well preserves the essential properties of the model.Comparison and convergence analyses of the method are presented.

A baroclinic typhoon model with a moving multi-nested grid and variational adjustment initialisationI. Numerical method

A baroclinic typhoon model with a moving multi--nested grid is applied in marine environmental forecasts. This paper describes the numerical methods of the model including governing equations, finite differencing, split scheme and time integration.

A High-Order Conservative Numerical Method for Gross-Pitaevskii Equation with Time-Varying Coefficients in Modeling BEC

We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.

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.

Dynamic Damage Model of Brittle Rock and Its Application

On the basis of shock induced experiments and the ultrasonic tests of the damaged rocks, the damage evolution relation between the attenuation coefficient of sound wave and the damage dissipated energy is described. Based on the TCK and RDA models, a damage model which connects the shock compression and tensile damage is established. And then the damage model is implemented in LS DYNA3D dynamic nonlinear program. Numerical simulation of deep hole blasting of groove is studied by use of the damage model proposed. The rock damage evolution process and the distributing rules of stress field under the explosion load are described well fairly, which provides the theory basis for the engineering blasting design.

Horizontal Slow Oscillation of a Moored Tanker and Numerical Method of Hopf-Bifurcation

In this paper, a six-order nonlinear dynamic model with three degrees of freedom is presented for the study of the 'fishtail' motion of a Single Point Mooring System. The effect of parameter variations on the equilibrium state of the system is analyzed. In order to study the stability of the equilibrium state, the mooring-line length l is chosen as a bifurcation parameter, so that all eigenvalues of the Jacobian matrix under different parameters can be worked out, and then the Hopf-bifurcation point can be found. Finally, the Hopf-bifurcation periodic solution of the system is computed.

STRUCTURAL DAMAGE MODEL OF UNSATURATED EXPANSIVE SOIL AND ITS APPLICATION IN MULTI-FIELD COUPLE ANALYSIS ON EXPANSIVE SOIL SLOPE

Focused on the sensitivity to climate change and the special mechanical characteristics of undisturbed expansive soil, an elasto-plastic damage constitutive model was proposed based on the mechanics of unsaturated soil and the mechanics of damage. Undisturbed expansive soil was considered as a compound of non-damaged part and damaged part. The behavior of the non-damaged part was described using non-linear constitutive model of unsaturated soil. The property of the damaged part was described using a damage evolution equation and two yield surfaces, i.e., loading yield （LY） and shear yield （SY）. Furthermore, a consolidation model for unsaturated undisturbed expansive soil was established and a FEM program named UESEPDC was designed. Numerical analysis on solid-liquid-gas tri-phases and multi-field couple problem was conducted for four stages and fields of stress, displacement, pore water pressure, pore air pressure, water content, suction, and the damage region as well as plastic region in an expansive soil slope were obtained.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

The dynamic response and damage models of rebar reinforced polymer slabs subjected to contact and near-field explosions

Non aqueous reactive polymer materials produced by the reaction of isocyanate and polyol have been widely used in infrastructure construction,which may be subjected to explosion loads during complex service conditions.The blast response of composite materials is a crucial aspect for applications in engineering structures potentially subjected to extreme loadings.In this work,damage caused to rebar reinforced polymer slabs by surface explosive charges was studied experimentally and numerically.A total of 6 field tests were carried out to investigate the performances of the failure modes of rebar reinforced polymer slabs under contact and near-field explosions.The influence of explosive quantity(10-40 g)and stand-off distances(0-20 cm)at the damage modes were studied.The results show that the failure modes of rebar reinforced polymer slabs under near-field explosion mainly were bending and surface spalling,while under the impact of contact explosion,the failure modes were craters of the top surface,spalling of the bottom surface,and middle perforation.Furthermore,a detailed fully coupled model was developed and validated with the test data.The influences of explosive quantity and slab thickness on rebar reinforced polymer slabs under contact explosion were studied.Based on this,the calculation formula between breach diameter,explosive quantity,and slab thickness is fitted.

Novel damage constitutive models and new quantitative identification method for stress thresholds of rocks under uniaxial compression

Four key stress thresholds exist in the compression process of rocks,i.e.,crack closure stress(σ_(cc)),crack initiation stress(σ_(ci)),crack damage stress(σ_(cd))and compressive strength(σ_(c)).The quantitative identifications of the first three stress thresholds are of great significance for characterizing the microcrack growth and damage evolution of rocks under compression.In this paper,a new method based on damage constitutive model is proposed to quantitatively measure the stress thresholds of rocks.Firstly,two different damage constitutive models were constructed based on acoustic emission(AE)counts and Weibull distribution function considering the compaction stages of the rock and the bearing capacity of the damage element.Then,the accumulative AE counts method(ACLM),AE count rate method(CRM)and constitutive model method(CMM)were introduced to determine the stress thresholds of rocks.Finally,the stress thresholds of 9 different rocks were identified by ACLM,CRM,and CMM.The results show that the theoretical stress−strain curves obtained from the two damage constitutive models are in good agreement with that of the experimental data,and the differences between the two damage constitutive models mainly come from the evolutionary differences of the damage variables.The results of the stress thresholds identified by the CMM are in good agreement with those identified by the AE methods,i.e.,ACLM and CRM.Therefore,the proposed CMM can be used to determine the stress thresholds of rocks.

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.

NUMERICAL MODELING OF KELVIN-HELMHOLTZ INSTABILITY AND RAYLEIGH-TAYLOR INSTABILITY AT LARGE DEFORMATION STAGE——(Ⅰ) NUMERICAL METHOD AND NUMERICAL RESULTS OF RAYLEIGHTAYLOR INSTABILITY

This paper improves the discrete vortex method for modeling Kelvin-Helmholtz instability and Rayleigh-Tay- lor instability by proper choice of velocity weighted average coefficients, redistribution of markers and successive adding of computational points with the increase of interfacial deformation and gives the numerical results of Rayleigh-Taylor instability. The numerical results show that the first two techniques greatly enhance the ability of the discrete vortex method for modeling large interracial deformations and the last technique greatly reduces the computational amounts of the numerical modeling at large deformation stage. The numerical modeling of Rayleigh- Taylor instability not only reproduces some phenomena such as the roll up at the end part of the spike observed in experiments but also finds some new phenomena such as the splashes at the roll up parts which needs to be tested by experiment.

A Study on Numerical Calculation Method of Small Cluster Density in Percolation Model

Percolation theory deals with the numbers and properties of the clusters formed in the different occupation probability. In this Paper, we study the calculation method of small clusters. We calcu-lated the small cluster density of 1, 2 and 3 in the percolation model with the exact method and the numerical method. The results of the two methods are very close, which can be verified by each other. We find that the cluster density of all three kinds of small clusters reaches the highest value when the occupation probability is between 0.1 and 0.2. It is very difficult to get the analytical formula for the exact method when the cluster area is relatively large (such as the area is more than 50), so we can get the density value of the cluster by numerical method. We find that the time required calculating the cluster density is proportional to the percolation area, which is indepen-dent of the cluster size and the occupation probability.

Overhanging rock slope by design:An integrated approach using rock mass strength characterisation,large-scale numerical modelling and limit equilibrium methods

Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.

3-D numerical modelling of Domino failure of hard rock pillars in Fetr6 Chromite Mine, Iran, and comparison with empirical methods

Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.

Numerical Modeling of Response and Damage of Masonry Walls to Blast Loading

To model the damage process of masonry walls under blast loading, a dynamic continuum damage material model is constructed for brick and mortar separately. The degradation of both the stiffness and strength are governed by a damage variable. By using the proposed material model, damage and fragmentation of a typical masonry wall under blast loading at different scaled distances is calculated. The hazard level of the masonry wall to blast loading is evaluated by analyzing the numerical results.

A method to calculate design tide levels on the basis of numerical model of tidal current and its application

In order to determine the design tide levels in the areas without measured tide level data, especially in the areas where it is difficult to measure tidal levels, a calculation method based on a numerical model of tidal current is proposed. The essentials of the method are described, and its application is illustrated with an example. The results of the application show that the design tide levels calculated by the method are close to those determined by long-time measured tide level data, and its calculation precision is high, so it is feasible to use the method to determine the design tide levels in the areas.

A Finite Volume Method with Unstructured Triangular Grids for Numerical Modeling of Tidal Current

The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.

Modelling of blast-induced damage in tunnels using a hybrid finite-discrete numerical approach

This paper presents the application of a hybrid finite-discrete element method to study blast-induceddamage in circular tunnels. An extensive database of field tests of underground explosions above tunnelsis used for calibrating and validating the proposed numerical method; the numerical results areshown to be in good agreement with published data for large-scale physical experiments. The method isthen used to investigate the influence of rock strength properties on tunnel durability to withstand blastloads. The presented analysis considers blast damage in tunnels excavated through relatively weak（sandstone） and strong （granite） rock materials. It was found that higher rock strength will increase thetunnel resistance to the load on one hand, but decrease attenuation on the other hand. Thus, undercertain conditions, results for weak and strong rock masses are similar. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.