Numerical simulation of dynamic surface deformation based on DInSAR monitoring

Differential interferometric synthetic aperture radar （DInSAR） technology is a new method to monitor the dynamic surface subsidence. It can monitor the large scope of dynamic deformation process of surface subsidence basin and better reflect the surface subsidence form in different stages. But under the influence of factors such as noise and other factors, the tilt and horizontal deformation curves regularity calculated by DInSAR data are poorer and the actual deviation is larger. The tilt and horizontal deformations are the important indices for the safety of surface objects protection. Numerical simulation method was used to study the dynamic deformation of LW32 of West Cliff colliery in Australia based on the DInSAR monitoring data. The result indicates that the subsidence curves of two methods fit well and the correlation coefficient is more than 95%. The other deformations calculated by numerical simulation results are close to the theory form. Therefore, considering the influence, the surface and its subsidiary structures and buildings due to mining, the numerical simulation method based on the DInSAR data can reveal the distribution rules of the surface dynamic deformation values and supply the shortcomings of DInSAR technology. The research shows that the method has good applicability and can provide reference for similar situation.

Numerical simulation of dynamic response of operating metro tunnel induced by ground explosion

To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit dynamic nonlinear finite element software ANSYS/LS-DYNA. The blast induced wave propagation in the soil and the tunnel, and the von Mises effective stress and acceleration of the tunnel lining were presented, and the safety of the tunnel lining was evaluated based on the failure criterion. Besides, the parametric study of the soil was also carried out. The numerical results indicate that the upper part of the tunnel lining cross-section with directions ranging from 0° to 22.5° and horizontal distances 0 to 7 m away from the explosive center are the vulnerable areas, and the metro tunnel might be safe when tunnel depth is more than 7 m and TNT charge on the ground is no more than 500 kg, and the selection of soil parameters should be paid more attentions to conduct a more precise analysis.

Numerical simulation of dynamic characteristics of a cable-stayed aqueduct bridge

In this paper, a full-scale 3-D finite element model of the Jundushan cable-stayed aqueduct bridge is established with ANSYS Code. The shell, fluid, tension-only spar and beam elements are used for modeling the aqueduct deck, filled water, cables and support towers, respectively. A multi-element cable formulation is introduced to simulate the cable vibration. The dry （without water） and wet （with water） modes of the aqueduct bridge are both extracted and investigated in detail. The dry modes of the aqueduct bridge are basically similar to those of highway cable-stayed bridges. A dry mode may correspond to two types of wet modes, which are called the in-phase （with lower frequency） and out-of-phase （with higher frequency） modes. When the water-structure system vibrates in the in-phase/out-of-phase modes, the aqueduct deck moves and water sloshes in the same/opposite phase-angle, and the sloshing water may take different surface-wave modes. The wet modes of the system reflect the properties of interaction among the deck, towers, cables and water. The in-phase wet frequency generally decreases as the water depth increases, and the out-of-phase wet frequency may increase or decrease as the water depth increases.

Dynamic response and numerical simulation of Al-Sc and Al-Ti alloys under high-speed impact

Al-Sc and Al-Ti semi-infinite targets were impacted by high-speed projectiles at velocities of 0.8, 1.0, 1.2 and 1.5 km/s, respectively. It is found that the Al-Sc targets demonstrate more excellent ability to resist high-speed impact. It is concluded that different microstructures of Al-Sc and Al-Ti alloys, including different grain sizes and secondary particles precipitated in the matrix, result in their greatly different capabilities of resisting impact. Furthermore, the effect of the size range ofnanoscale A13Sc precipitate in A1-Sc alloy on the resistance of high-speed impact was investigated. In addition, computer simulations and validation of these simulations were developed which fairly accurately represented residual crater shapes/geometries. Validated computer simulations allowed representative extrapolations of impact craters well beyond the laboratory where melt and solidification occurred at the crater wall, especially for hypervelocity impact （〉5 km/s）.

Stress initialization methods for dynamic numerical simulation of rock mass with high in-situ stress

In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.

Simulation of large-scale numerical substructure in real-time dynamic hybrid testing

A solution scheme is proposed in this paper for an existing RTDHT system to simulate large-scale finite element （FE） numerical substructures. The analysis of the FE numerical substructure is split into response analysis and signal generation tasks, and executed in two different target computers in real-time. One target computer implements the response analysis task, wherein a large time-step is used to solve the FE substructure, and another target computer implements the signal generation task, wherein an interpolation program is used to generate control signals in a small time-step to meet the input demand of the controller. By using this strategy, the scale of the FE numerical substructure simulation may be increased significantly. The proposed scheme is initially verified by two FE numerical substructure models with 98 and 1240 degrees of freedom （DOFs）. Thereafter, RTDHTs of a single frame-foundation structure are implemented where the foundation, considered as the numerical substructure, is simulated by the FE model with 1240 DOFs. Good agreements between the results of the RTDHT and those from the FE analysis in ABAQUS are obtained.

Salt-Water Transport in Unsaturated Soils Under Crop Planting: Dynamics and Numerical Simulation

A laboratory salt-water dynamics experiment using unsaturated soils in packed silt loam and clay soil columns withdifferent soil texture profiles and groundwater levels under crops were conducted to study the changes of salt-waterdynamics induced by water uptake of crops and to propose the theoretical basis for the regulation and control of salt-water dynamics as well as to predict salinity levels. The HYDRUS 1D model was applied to simulate the one-dimensionalmovement of water and salt transport in the soil columns. The results showed that the salts mainly accumulated in theplow layer in the soil columns under crops. Soil water and salt both moved towards the plow layer due to soil waterabsorption by the crop root system. The salt contents in the column with lower groundwater were mostly greater thanthose with high groundwater. The water contents in the soil columns increased from top to the bottom due to plant rootwater uptake. The changes in groundwater level had little influence on water content of the root zone in the soil columnswith crop planting. Comparison between the simulated and the determined values showed that model simulation resultswere ideal, so it is practicable to do numerical simulation of soil salt and water transport by the HYDRUS 1D model.Furthermore, if the actual movement of salt and water in fields is to be described in detail, much work needs to be done.The most important thing is to refine the parameters and select precise boundary conditions.

Multi-parameter numerical simulation of dynamic monitoring of rock deformation in deep mining

The level of deformation development of surrounding rocks is a vital predictor to evaluate impending coal mine disasters and it is important to establish accurate measurements of the deformed status to ensure coal mine safety. Traditional deformation monitoring methods are mostly based on single parameter, in this paper, multiple approaches are integrated: firstly, both electric and elastic models are established,from which electric field distribution and seismic wave recording are calculated and finally, the resistivity profiles and source position information are determined using inversion methods, from which then the deformation and failure of mine floor are evaluated. According to the inversion results of both electric and seismic field signals, multiple-parameter dynamic monitoring of surrounding rock deformation in deep mine can be performed. The methodology is validated using numerical simulation results which shows that the multi-parameter dynamic monitoring methods have better results for surrounding rock deformation in deep mine monitoring than single parameter methods.

Numerical simulation of dynamic process ofthe Tangshan earthquake

In this paper. :LDDA (Lagrangian Discontinuous Deformation Analysis) method is used in modeling thedynamic process of the M,=7.8 Tang shan earthquake on July 28, 1976 and obtain directly the dynamic and quasi static dislocations. shear stress drops, fracture velocities of the Tang shan earthquake fault. The simulation showsthai the slip history at each point of the fault is different. The displacement vectors at the concave side of the faultis greater than that of the convex side of the fault. The 'over shoot' of the fault slip is greatest at the middle part ofthe fault and attenuates to its ends. The rupture velocities of the fault from the epicenter towards south-west andtowards north-east are 3.08 m/s and 1. 18 m/s, respectively, the average one is 2.13 In/s. The maximum dynamic.m quasi-static dislocations are 7. 1 m and 6.2 m respectively. the average quasi-static one on the fault is 4.5 m.initial stress dynamic and quasi-static shear stress drops are 8.1 M Pa and 5.4 MP4 respectively, the averagequasi-static shear stress drop is 3.9 M Pa. We found that the rupture velocities and shear stress are related to theinitial stress states of the fault.

Model test and numerical simulation on the dynamic stability of the bedding rock slope under frequent microseisms

Shake table testing was performed to investigate the dynamic stability of a mid-dip bedding rock slope under frequent earthquakes. Then, numerical modelling was established to further study the slope dynamic stability under purely microseisms and the influence of five factors, including seismic amplitude, slope height, slope angle, strata inclination and strata thickness, were considered. The experimental results show that the natural frequency of the slope decreases and damping ratio increases as the earthquake loading times increase. The dynamic strength reduction method is adopted for the stability evaluation of the bedding rock slope in numerical simulation, and the slope stability decreases with the increase of seismic amplitude, increase of slope height, reduction of strata thickness and increase of slope angle. The failure mode of a mid-dip bedding rock slope in the shaking table test is integral slipping along the bedding surface with dipping tensile cracks at the slope rear edge going through the bedding surfaces. In the numerical simulation, the long-term stability of a mid-dip bedding slope is worst under frequent microseisms and the slope is at risk of integral sliding instability, whereas the slope rock mass is more broken than shown in the shaking table test. The research results are of practical significance to better understand the formation mechanism of reservoir landslides and prevent future landslide disasters.

Numerical simulation of dynamic fracture properties of rocks under different static stress conditions

When underground cavities are subjected to explosive stress waves,a uniquely damaged zone may appear due to the combined effect of dynamic loading and static pre-load stress.In this study,a rate-dependent two-dimensional rock dynamic constitutive model was established to investigate the dynamic fractures of rocks under different static stress conditions.The effects of the loading rate and peak amplitude of the blasting wave under different confining pressures and the vertical compressive coefficient(K_(0))were considered.The numerical simulated results reproduced the initiation and further propagation of primary radial crack fractures,which were in agreement with the experimental results.The dynamic loading rate,peak amplitude,static vertical compressive coefficient(K_(0))and confining pressure affected the evolution of fractures around the borehole.The heterogeneity parameter(m)plays an important role in the evolution of fractures around the borehole.The crack propagation path became more discontinuous and rougher in a smallerheterogeneity parameter case.

Three-dimensional transient numerical simulation for intake process in the engine intake port-valve-cylinder system

This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used in numerical study on internal combustion engine with vertical and inclined valves, and has higher calculation precision. A numerical simulation (on the intake process of a two-valve engine with a semi-sphere combustion chamber and a radial intake port) is provided for analysis of the velocity field and pressure field of different plane at different crank angles. The results revealed the formation of the tumble motion, the evolution of flow field parameters and the variation of tumble ratios as important information for the design of engine in-take system.

Review of numerical simulation on the dynamics of Qinghai-Xizang plateau

In recent twenty years, much numerical simulation work has been done on the evolution of Qinghai-Xizang (Tibetan) plateau. In this paper some principal numerical models and results are reviewed and analyzed. The earlier plane stress or plane strain model has much discrepancy with the actual deformation of Qinghai-Xizang plateau, such as the thickening of Tibetan crust and the lateral extrusion of Tibet along strike-slip faults. The thin viscous sheet model and the thin-plate model may simulate the change of the crustal thickness and the deformation pro-duced by gravitational force. It is suitable for studying the large-scale and long-time deformation. The influence of faults on the deformation of Tibetan plateau should be further studied.

Numerical Simulation of Dynamic Performance of the Molten Carbonate Fuel Cell

A three dimension of dynamic mathematical model of the molten carbonate fuel cell is established,in which the heat generation, mass transfer and electrochemical characteristics are described. The performance of the fuel cell including the distributions of the temperature and the velocity is predicted numerically. Then the experimental data including the output performance of the fuel cell generation system and the temperature distributions are compared. The numerical results are in agreement with the experiment results.

Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

Numerical Simulation of Dynamic Response of Towed Systems

In this paper, the dynamic response of undersea -+towed systems is numerically simulated. A two body towed system is especially considered in detail. The factors influencing the heave of towed-bodies, such as the weight of the towed-body (in sea water), the length and the weight (in sea water) per unit length of the cable between towed-bodies and towing ship, are investigated in detail. Calculations show that the two-body towed system can greatly increases the stability of the towed system.

Three-dimensional discrete element numerical simulation of Paleogene salt structures in the western Kuqa foreland thrust belt

Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.

Study on numerical simulation and dynamic mechanism of win-ter- time circulation in the eastern China seas

An MOM2 based 3-dimentional prognostic baroclinic Z-ordinate model was established to study the circulation in eastern China seas, considering the topography, inflow and outflow on the open boundary, wind stress, temperature and salinity exchange on the sea surface. The results were consistent with observation and showed that the Kuroshio intrudes in large scale into the East China Sea continental shelf East China, during which its water is exchanged ceaselessly with outer sea water along Ryukyu Island. The Tsushima Warm Current is derived from several sources, a branch of the Kuroshio, part of the Taiwan Warm Current, and Yellow Sea mixed water coming from the west of Cheju Island. The water from the west of Cheju Island contributes ap-proximately 13% of the Isushima Warm Current total transport through the Korea Strait. The circulation in the Bohai Sea and Yellow Sea is basically cyclonic circulation, and is comprised of coastal currents and the Yellow Sea Warm Current. Besides simulation of the real circulation, numerical experiments were conducted to study the dynamic mechanism. The numerical experiments indicated that wind directly drives the East China Sea and Yellow Sea Coastal Currents, and strengthens the Korea Coastal Current and Yellow Sea Warm Current. In the no wind case, the kinetic energy of the coastal current area and main YSWC area is only 1% of that of the wind case. Numerical experiments also showed that the Tsushima Warm Current is of great importance to the formation of the Korea Coastal Current and Yellow Sea Warm Current.

Dynamic Interference Photoelectron Spectra in Double Ionization:Numerical Simulation of 1D Helium

We report our numerical simulation on the dynamic interference photoelectron spectra for a one-dimensional （1D） He model exposed to intense ultrashort extreme ultraviolet （XUV） laser pulses. The results demonstrate an unambiguous interference feature in the photoelectron spectra, and the interference is unveiled to originate from the dynamic Stark effect. The interference photoelectron spectra are prompted for intense sub-femtosecond XUV laser pulses in double ionization. The stationary phase picture is corroborated qualitatively in the two-electron system. The ability of probing the dynamic Stark effect by the photoelectron spectra in a pragmatic experiment of single-photon double ionization of He may shed light on further investigation on multi-electron atoms and molecules.

Numerical simulation of the effects of the impact velocity on the particle deposition characteristics in cold gas dynamic spraying

In this study, the effects of the impact velocity on the particle deposition characteristics in cold gas dynamic spraying （CGDS） of 304 stainless steel （SS） on an interstitial free （IF） steel substrate are numerical simulated by means of a finite element analysis （FEA）. The results have illustrated that when the particle impact velocity exceeds a critical value at which adiabatic shear instability of the particle starts to occur. Meanwhile, the fatten ratio and impact crater depth （or the effective contacting area ） increase rapidly. The particle-substrate bonding and deposition mechanism can be attributed to such an adiabatic shear deformation induced by both the compressive force and the slide friction force of particle. The critical velocity can be predicted by numerical simulation, which is useful to optimize the CGDS processing parameters for various materials.