Numerical modeling and dynamic analysis of the 2017 Xinmo landslide in Maoxian County, China

A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eighty-three people were killed or went missing and more than 103 houses were destroyed. In this paper, the geological conditions of the landslide are analyzed via field investigation and high-resolution imagery. The dynamic process and runout characteristics of the landslide are numerically analyzed using a depth-integrated continuum method and Mac Cormack-TVD finite difference algorithm.Computational results show that the evaluated area of the danger zone matchs well with the results of field investigation. It is worth noting that soil sprayed by the high-speed blast needs to be taken into account for such kind of large high-locality landslide. The maximum velocity is about 55 m/s, which is consistent with most cases. In addition, the potential danger zone of an unstable block is evaluated. The potential risk area evaluated by the efficient depthintegrated continuum method could play a significant role in disaster prevention and secondary hazard avoidance during rescue operations.

Numerical analysis of the effects of rock bolts on stress redistribution around a roadway

Besides opening geometry, in situ stress and material properties, opening support also has significant effects on stress redistribution around a roadway. To investigate these effects of rock bolts on the stress redistribution around a roadway, a series of numerical studies were carried out using the finite difference method. Since the stress changes around a roadway caused by rock bolting is small relative to the in situ stress, they cannot obviously be observed in stress contour plots. To overcome this difficulty, a new result processing methodology was developed using the contouring program Surfer. With this methodology, the effects of rock bolts on stress redistribution can obviously be analyzed. Numerical results show that in the three patterns of rock bolts installed in the roof, in the roof and the two lateral sides, and in all the four sides of the rectangular roadway, the maximum stress magnitude of the increase is 0.931 MPa, 2.46 MPa,and 6.5 MPa, respectively; the bolt number of 5 can form an integrated ground arch; the appropriate length and pre-tensioned force of the rock bolt is 2.0 m and 60 k N, respectively. What is more, the ground arch action under the function of rock bolting is able to be effectively examined. The rock bolts dramatically increase the minor principal stress around a roadway which results in significant increase in material strength. Consequently, the major principal stress that the material can carry will greatly increase.With adequate supports, an integrated ground arch which is critical for the stability of roadway will be formed around the roadway.

Numerical simulation study of the failure evolution process and failure mode of surrounding rock in deep soft rock roadways

Based on the safety coefficient method,which assigns rock failure criteria to calculate the rock mass unit,the safety coefficient contour of surrounding rock is plotted to judge the distribution form of the fractured zone in the roadway.This will provide the basis numerical simulation to calculate the surrounding rock fractured zone in a roadway.Using the single factor and multi-factor orthogonal test method,the evolution law of roadway surrounding rock displacements,plastic zone and stress distribution under different conditions is studied.It reveals the roadway surrounding rock burst evolution process,and obtains five kinds of failure modes in deep soft rock roadway.Using the fuzzy mathematics clustering analysis method,the deep soft surrounding rock failure model in Zhujixi mine can be classified and patterns recognized.Compared to the identification results and the results detected by geological radar of surrounding rock loose circle,the reliability of the results of the pattern recognition is verified and lays the foundations for the support design of deep soft rock roadways.

A Study on Numerical Simulation of Core-Shooting Process

In order to evaluate the main factors influencing the core-shooting process and to optimize the design of core boxes, the fluid-particle (air-sand) model has been built based on the two-phase flow theory. The fluid phase, air, and the particulate phase, sand granules, have been treated as a continuum. By using this model, it is possible to simulate the flow and compaction behavior of sand particles during the core-shooting process. To benchmark the calculated results, the shooting process has been recorded with a digital high speed camera, and the inlet condition of sand particles has also been achieved by using the camera. The preliminary results have showed that the calculation is in agreement with the testing results.

Study on Numerical Simulation of Mold Filling and Heat Transfer in Die Casting Process

A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow of mold filling process of die casting. In addition, the mathematical model for simulating the heat transfer in die casting process has also been established. The computation program has been developed by the authors with the finite difference method (FDM) recently. As verification, the mold filling process of a S-shaped die casting has been simulated and the simulation results coincide with that of the benchmark test. Finally, as a practical application, the gating design of a motorcycle component was modified by the mold filling simulation and the dies design of another motorcycle component was optimized by the heat transfer simulation. All the optimized designs were verified by the production practice.

A numerical model of coastline deformation for sandy beach at downstream of a jetty

A reformed numerical model based on the 'one-line theory' for beach deformation is presented. In this model, the change of beach slope during coastline procession is considered.A wave numerical model combined with wave refraction, diffraction and reflection is used to simulate wave climate to increase numerical accuracy.The results show that the numerical model has a good precision based on the adequate field data. The results can be applied to practical engineering.

Numerical simulation of non-gaussian process of wind waves

In this paper further mathematical analysis on 'correlation transfer technique' by Polge el al. is carried out, the tenable conditions and the extent of suitability for the said method are proved as well. In consideration of the influence of skewness of the sea surface elevation on spectral shape, a 'quasi-correlation transfer techique' is developed by the modification of the simulated target spectrum. Meanwhile, the numerical simulation of the non-Gaussian process of wind waves is carried out in view of the two conditions of the surface elevation probability distribution and the spectrum. By using its simulated results, the influence of skewness of the sea surface elevation on two parameters in the distribution of wave heights (which had been fitted by using the Weibull distribution) is analysed. The 'quasi- correlation tranfer technique' is verified and compared with the selection wave data observed in the Jiaozhou Bay in the period of 1980 to 1981. Results make clear (hat, as far as the statistical distribution of the wave heights and the distribution of the maximum (minimum) values of the sea surface elevation are concerned, the said method is obviously superior to the conventional mothed of the linear wave superposition, and that the simulated results are closer to the observation data.

Numerical simulation of cutting process of the slice components cutting machine

By the use of ANSYS/LS-DYNA FEA software,numerical simulation on the cutting process of cutting plates with a reamer was carried out in the paper. The logical improvement was brought forward and the phenomenon of stress concentration was deceased by weighted analysis. The effects of different cut velocities and cutting thickness on life-spans of reamers were investigated, and the cutting parameters were optimized to satisfy the cutting precision and cutting efficiency. The study will provide a guide for the practical production.

The Design of a Graphical User Environment for Numerical Simulation of Powder Forming Processes

As computer simulation increasingly supports engine er ing design and manufacture, the requirement for a computer software environment providing an integration platform for computational engineering software increas es. A key component of an integrated environment is the use of computational eng ineering to assist and support solutions for complex design. Computer methods fo r structural, flow and thermal analysis are well developed and have been used in design for many years. Many software packages are now available which provi de an advanced capability. However, they are not designed for modelling of powde r forming processes. This paper describes the powder compaction software (PCS_SU T), which is designed for pre- and post-processing for computational simulatio n of the process compaction of powder. In the PCS_SUT software, the adaptive analysis of transient metal powder forming process is simulated by the finite element method based on deformation theories . The error estimates and adaptive remeshing schemes are applied for updated co -ordinate analysis. A generalized Newmark scheme is used for the time domain di scretization and the final nonlinear equations are solved by a Newton-Raphson p rocedure. An incremental elasto-plastic material model is used to simulate the compaction process. To describe the constitutive model of nonlinear behaviour of powder materials, a combination of Mohr-Coulomb and elliptical yield cap model is applied. This model reflects the yielding, frictional and densification char acteristics of powder along with strain and geometrical hardening which occurs d uring the compaction process. A hardening rule is used to define the dependence of the yield surface on the degree of plastic straining. A plasticity theory for friction is employed in the treatment of the powder-tooling interface. The inv olvement of two different materials, which have contact and relative movement in relation to each other, must be considered. A special formulation for friction modelling is coupled with a material formulation. The interface behaviour betwee n the die and the powder is modelled by using an interface element mesh. In the present paper, we have demonstrated pre- and post-processor finite elem ent software, written in Visual Basic, to generate the graphical model and visua lly display the computed results. The software consist of three main part: · Pre-processor: It is used to create the model, generate an app ropriate finite element grid, apply the appropriate boundary conditions, and vie w the total model. The geometric model can be used to associate the mesh with th e physical attributes such as element properties, material properties, or loads and boundary conditions. · Analysis: It can deal with two-dimensional and axi-symmetric applications for linear and non-linear behaviour of material in static and dyna mic analyses. Both triangular and quadrilateral elements are available in the e lement library, including 3-noded, 6-noded and 7-noded (T6B1) triangles and 4 -noded, 8-noded and 9-noded quadrilaterals. The direct implicit algorithm bas ed on the generalized Newmark scheme is used for the time integration and an aut omatic time step control facility is provided. For non-linear iteration, choice s among fully or modified Newton-Raphson method and quasi-Newton method, using the initial stiffness method, Davidon inverse method or BFGS inverse method, ar e possible. · Post-processor: It provides visualization of the computed resu lts, when the finite element model and analysis have been completed. Post-proce ssing is vital to allow the appropriate interpretation of the completed results of the finite element analysis. It provides the visual means to interpret the va st amounts of computed results generated. Finally, the powder behaviour during the compaction of a multi-level component is numerically simulated by the PCS_SUT software, as shown in Fig.1. The predict ive compaction forces at different displacements are computed and compared with the available experimental

Numerical Study on Dry Deposition Processes in Canopy Layer

A coupling model between the canopy layer(CL) and atmospheric boundary layer (ABL) for the study of dry deposition velocity is developed. The model consists of six parts: chemical species conservation equation including absorptive factor; the species uptake action including detailed vertical variation of absorptive element in CL; momentum exchange in CL which is represented by a first-order closure momentum equation with an additional larger-scale diffusive term; momentum exchange in ABL which is described by a complete set of the ABL turbulent statistic parameters; absorptivity (or solubility or reflection) at the surface including effects of the physical and chemical characters of the species, land type, seasonal and diurnal variations of the meteorological variables; and deposition velocity derived by distributions of the species with height in CL. Variational rules of the concentration and deposition velocity with both height and time are simulated with the model for both corn and forest canopies. Results predicted with the bulk deposition velocity derived in the paper consist well with experimental data.

On Numerical Modelling of Industrial Powder Compaction Processes for Large Deformation of Endochronic Plasticity at Finite Strains

Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the compaction processes are to obtain a compact with the geometrical requirements, without cracks, and with a uniform distribution of density. Design of such proc esses consist, essentially, in determine the sequence and relative displacements of die and punches in order to achieve such goals. A.B. Khoei presented a gener al framework for the finite element simulation of powder forming processes based on the following aspects; a large displacement formulation, centred on a total and updated Lagrangian formulation; an adaptive finite element strategy based on error estimates and automatic remeshing techniques; a cap model based on a hard ening rule in modelling of the highly non-linear behaviour of material; and the use of an efficient contact algorithm in the context of an interface element fo rmulation. In these references, the non-linear behaviour of powder was adequately desc ribed by the cap plasticity model. However, it suffers from a serious deficiency when the stress-point reaches a yield surface. In the flow theory of plasticit y, the transition from an elastic state to an elasto-plastic state appears more or less abruptly. For powder material it is very difficult to define the locati on of yield surface, because there is no distinct transition from elastic to ela stic-plastic behaviour. Results of experimental test on some hard met al powder show that the plastic effects were begun immediately upon loading. In such mater ials the domain of the yield surface would collapse to a point, so making the di rection of plastic increment indeterminate, because all directions are normal to a point. Thus, the classical plasticity theory cannot deal with such materials and an advanced constitutive theory is necessary. In the present paper, the constitutive equations of powder materials will be discussed via an endochronic theory of plasticity. This theory provides a unifi ed point of view to describe the elastic-plastic behaviour of material since it places no requirement for a yield surface and a ’loading function’ to disting uish between loading an unloading. Endochronic theory of plasticity has been app lied to a number of metallic materials, concrete and sand, but to the knowledge of authors, no numerical scheme of the model has been applied to powder material . In the present paper, a new approach is developed based on an endochronic rate independent, density-dependent plasticity model for describing the isothermal deformation behavior of metal powder at low homologous temperature. Although the concept of yield surface has not been explicitly assumed in endochronic theory, it is shown that the cone-cap plasticity yield surface (Fig.1), which is the m ost commonly used plasticity models for describing the behavior of powder materi al can be easily derived as a special case of the proposed endochronic theory. Fig.1 Trace of cone-cap yield function on the meridian pl ane for different relative density As large deformation is observed in powder compaction process, a hypoelastic-pl astic formulation is developed in the context of finite deformation plasticity. Constitutive equations are stated in unrotated frame of reference that greatly s implifies endochronic constitutive relation in finite plasticity. Constitutive e quations of the endochronic theory and their numerical integration are establish ed and procedures for determining material parameters of the model are demonstra ted. Finally, the numerical schemes are examined for efficiency in the model ling of a tip shaped component, as shown in Fig.2. Fig.2 A shaped tip component. a) Geometry, boundary conditio n and finite element mesh; b) density distribution at final stage of

An iterative wavefront sensing algorithm for high-contrast imaging systems

Wavefront sensing from multiple focal plane images is a promising technique for high-contrast imaging systems.However,the wavefront error of an optics system can be properly reconstructed only when it is very small.This paper presents an iterative optimization algorithm for the direct measurement of large static wavefront errors from only one focal plane image.We first measure the intensity of the pupil image to get the pupil function of the system and acquire the aberrated image on the focal plane with a phase error that will be measured.Then we induce a dynamic phase on the tested pupil function and calculate the associated intensity of the reconstructed image on the focal plane.The algorithm will then try to minimize the intensity difference between the reconstructed image and the aberrated test image in the focal plane,where the induced phase is a variable of the optimization algorithm.The simulation shows that the wavefront of an optical system can theoretically be reconstructed with high precision,which indicates that such an iterative algorithm may be an effective way to perform wavefront sensing for high-contrast imaging systems.

Numerical magnitude processing in abacus-trained children with superior mathematical ability: an EEG study

Distance effect has been regarded as the best established marker of basic numerical magnitude processes and is related to individual mathematical abilities. A larger behavioral distance effect is suggested to be concomitant with lower mathematical achievement in children. However, the relationship between distance effect and superior mathematical abilities is unclear. One could get superior mathematical abilities by acquiring the skill of abacus-based mental calculation （AMC）, which can be used to solve calculation problems with exceptional speed and high accuracy. In the current study, we explore the relationship between distance effect and superior mathematical abilities by examining whether and how the AMC training modifies numerical magnitude processing. Thus, mathematical competencies were tested in 18 abacus-trained children （who accepted the AMC training） and 18 non-trained children. Electroencephalography （EEG） waveforms were recorded when these children executed numerical comparison tasks in both Arabic digit and dot array forms. We found that： （a） the abacus-trained group had superior mathematical abilities than their peers; （b） distance effects were found both in behavioral results and on EEG wave- forms; （c） the distance effect size of the average amplitude on the late negative-going component was different between groups in the digit task, with a larger effect size for abacus-trained children; （d） both the behavioral and EEG distance effects were modulated by the notation. These results revealed that the neural substrates of magnitude processing were modified by AMC training, and suggested that the mechanism of the representation of numerical magnitude for children with superior mathematical abilities was different from their peers. In addition, the results provide evidence for a view of non-abstract numerical representation.

Fast compression and reconstruction of astronomical images based on compressed sensing

With the fast increase in the resolution of astronomical images, the question of how to process and transfer such large images has become a key issue in astronomy. We propose a new real-time compression and fast reconstruction algorithm for astronomical images based on compressive sensing techniques. We first reconstruct tile Original signal with fewer measurements, according to its compressibility. Then, based on the characteristics of astronomical images, we apply Daubechies orthogonal wavelets to obtain a sparse representation. A matrix representing a random Fourier ensembleis used to obtain a sparse representation in a lower dimensional space. For reconstructing the image, we propose a novel minimum total variation with block addptive sensing to balance the accuracy and eomputation time. Our experimental results show that the proposed algorithm can efficiently reconstruct colorful astronomicai images with high resolution and improve the applicability of compressed sensing.

Simulation Research on the Effect of Spreading Process Parameters on the Quality of Lunar Regolith Powder Bed in Additive Manufacturing

Lunar surface additive manufacturing with lunar regolith is a key step in in-situ resource utilization.The powder spreading process is the key process,which has a major impact on the quality of the powder bed and the precision of molded parts.In this study,the discrete element method(DEM)was adopted to simulate the powder spreading process with a roller.The three powder bed quality indicators,including the molding layer offset,voidage fraction,and surface roughness,were established.Besides,the influence of the three process parameters,which are roller’s translational speed,rotational speed,and powder spreading layer thickness on the powder bed quality indicators was also analyzed.The results show that with the reduction of the powder spreading layer thickness and the increase of the rotational speed,the offset increased significantly;when the translational speed increased,the offset first increased and then decreased,which resulted in an extreme value;with the increase of the layer thickness and the decrease of the translational speed,the values for voidage fraction and surface roughness significantly reduced.The powder bed quality indicators were adopted as the optimization objective,and the multi-objective parameter optimization was carried out.The predicted optimal powder spreading parameters and powder bed quality indicators were then obtained.Moreover,the optimal values were then verified.This study can provide informative guidance for in-situ manufacturing at the moon in future deep space exploration missions.

NUMERICAL EXPERIMENT WITH PROCESSES FOR EFFECT OF AUSTRALIAN COLD AIR ACTIVITY ON EAST-ASIAN SUMMER MONSOON

Based on diagnostic results,a numerical study is made of the processes of Australian cold air activity affecting East Asian summer monsoon by using Kuo-Qian P-σ incorporated coordinate five-layer primitive equation spherical band model.Analysis is done of the response to the Southern Hemisphere circulation with and without cold air activity in Australia of the flow,rainfall and diabatic heating fields in the monsoon area of Asia,especially,East Asia,with special attention to the intensification and northward march of the monsoon due to the activity.It is found that the processes for the effect transmission are very analogous to the meridional propagation of quasi-40-day oscillation,together with the meridional wind disturbance showing south-north travel and the flow/rainfall fields exhibitirg corresponding movement in this direction,only with a 12-day lag.

NUMERICAL SIMULATION OF MICROPHYSICAL PROCESSES IN CUMULONIMBUS——PART Ⅱ CASE STUDIES OF SHOWER,HAILSTORM AND TORRENTIAL RAIN

A shower cloud observed in Jiangxi,a hailstorm observed in Hebei and“75.8”torrential rain in Henan are simulated with our microphysical model in a one-dimensional framework.The model,using the radio- sonde data as input,gets its output which shows agreement in many aspects as compared with observations in each case.The glaciation of small cumulus cloud,low precipitation efficiency of hailstorm and the per- sistence of torrential rain are demonstrated.It is also shown that the Bergeron process has little influence, but the warm-rain process plays an important role in the formation of precipitation in cumulonimbus with a warm cloud base.

A NUMERICAL SIMULATION OF THE GEOSTROPHIC ADJUSTMENT PROCESS

In this paper,a numerical simulation of the geostrophic adjustment process with C-grid network is illustrated.A difference scheme which has the energy and potential vorticity conserving relation consistent with the differential equations is given,and the effect of some time difference schemes on dispersion of the gravity-inertia wave is discussed.An improved forward-backward time integration scheme is proposed for keeping the computational stability.The effect of various boundary conditions for a finite region model On the gravity-inertia wave is shown by some calculated results.

Transient simulation of regression rate on thrust regulation process in hybrid rocket motor

The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB（hydroxyl-terminated polybutadiene） hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.

Microstructure Evolution of Different Forging Processes for12%Cr Steel During Hot Deformation

Five forging experiments were designed and conducted to investigate the effect of process parameters on microstructure evolution during hot deformation for X12CrMoWVNbN10-1-1 steel.The experimental results indicated that average grain size became finer with the increasing number of upsetting and stretching.Especially,the size of stretching three times with upsetting twice had the most remarkable effect on refinement,and the size was only 27.36%of the original one.Moreover,the stress model was integrated into the software and finite element models were established.Simulation results demonstrated that the strain at center point of workpiece was far larger than critical strain value in each process,so that dynamic recrystallization(DRX) occurred in each workpiece,which implied DRX could occur for several times with the increasing number of upsetting and stretching,and uniform finer microstructure would be obtained.However,the results also showed that higher temperature was an unfavorable factor for grain refinement,so the times of heating should be limited for workpiece,and as many forging processes as possible should be finished in once heating.