Numerical simulations of full-wave fi elds and analysis of channel wave characteristics in 3-D coal mine roadway models

Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.

3-D Numerical Simulation on the Chip Machining Process of a Metal Block

In this paper, the cutting process of a metal block is numerically simulated by the dynamic explicit FE code ABAQUS. Taking thermo-mechanical coupling effect into consideration, the simulation presents the variation of temperature, stress and strain distribution in the workpiece and chip. The effective plastic strain failure criterion is applied to modeling the chip separation and plastic formation. And the phenomenon of the contact and friction between the workpiece and the cutting tool are described in the paper.

SIMULATION OF 3-D DEFORMATION AND MATERIAL FLOW DURING ROLL FORGING PROCESS USING SYSTEM OF OVAL-ROUND GROOVE

Basing on the analysis of the traits of the roll forging process, a system-model of computer simulation has been established. Three-dimensional rigid-plastic FEM has been used for the simulation of the deformation process in the oval and round pass rolling, including the entering, rolling, and separating stages. The analysis was conducted using the Deform-3D ver.5.0 code. The important information concerned with the deformation area characteristic, material flow, and velocity field has been presented. Otherwise, the location of the neutral plane in the deformation area was shown clearly.

3-D fracture network dynamic simulation based on error analysis in rock mass of dam foundation

Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.

HEXAHEDRAL ELEMENT REFINEMENT FOR THE PREDICTION- CORRECTION ALE FEM SIMULATION OF 3D BULKINGFORMING PROCESS

Based on the characteristics of 3D bulk forming process, the arbitrary Lagrangian-Eulerian (ALE) formulation-based FEM is studied, and a prediction-correction ALE-based FEM is proposed which integrates the advantages of precisely predicting the boundary configuration of the deformed material, and of efficiently avoiding hexahedron remeshing processes. The key idea of the prediction-correction ALE FEM is elaborated in detail. Accordingly, the strategy of mesh quality control, one of the key enabling techniques for the 3D bulk forming process numerical simulation by the prediction-correction ALE FEM is carefully investigated, and the algorithm for hexahedral element refinement is formulated based on the mesh distortion energy.

Simulation of subcritical flow pattern in 180° uniform and convergent open-channel bends using SSIIM 3-D model

In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180v bend, 0.6 m to 0.45 m wide, was simulated using the SSI1M 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widtbwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel＇s centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.

3-D Simulation Study on Seismic Response of Bridge Piles in Landslide

The anti-slide support structure is widely used in the anti-seismic reinforcement of bridge foundations,but related experimental research was processing slowly. Based on the prototype of the Jiuzhaigou bridge at the Chengdu-Lanzhou Railway,a 3-D simulation model was established on the basis of the shaking table model test,and the rationality of the dynamic analysis model was verified by indicators such as the bending moment of the bridge piles,peak soil pressure,and PGA amplification factors. The results show that the inertia force of the bridge pier has an important influence on the deformation of the pile foundation. The bending moment and shearing force are larger in lateral bridge piles,and the maximum value is near the pile top. The PGA amplification factor is stronger in the back of the rear anti-slide piles and so is it in front of the bridge pier,and the soil is prone to slip and damage. The bedrock is rigid and the dynamic response is maintained at a low level. The anti-slide piles in the rear row play a major role in the anti-seismic reinforcement design,and the anti-slide piles in the front row can be used as an auxiliary support structure.

3-D Finite-Element Numerical Simulation of Centrifugal Induction Electrosalg Casting Solidification Process

A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.

Filling Simulation of Three-dimension Braided Composite in Resin Transfer Molding

This paper measured permeability of three-dimension braided preform by radial technology. The results show that principal permeability tensor coincided with their braiding axial direction. The software of one dimensional flow filling mold was designed using Visual C++ language. Filling time is predicted and validated. The result showed that the filling time of the mold centerline agrees with the prediction value. The filling time of the mould edge is shorter than that of the prediction. An actual plate of 3D braided preform/ modified polyarylacetylene composite is produced according to prediction value and validation analysis.

THE SYSTEM SIMULATION OF THREE-DIMENSIONAL RADAR

To provide a test platform for Electronic Warfare (EW) system, it is needed to simulate the radar received Intermediate Frequency (IF) signals and radar system functions.This letter gives a description of a radar system simulation software developed for frequencyphase scanning three-dimensional (3-D) radar. Experimental results prove that the software could be used for system evaluation and for training purposes as an attractive alternative to real EW system.

Magnetic field of mathematical modeling and simulation of 3D magnetic pole array spherical actuator

In this paper a new spherical actuator is designed and its advantages are compared to an existing spherical actuator, which function is limited by several design bottlenecks. First the output torque is too small. Second, the attitude is difficult to be accurately detected. The new three-dimen- sional magnetic pole array can solve these major problems. The new actuator features an outer rotor with multiple permanent magnet （PM） poles. Using an analytical solution and the finite element so- lution simulation, the feasibility of the approach is verified. A prototype was developed, tested, and experiments were conducted to obtain the practical value of the magnetic flux density.

Three-dimensional numerical simulation of geothermal fi eld in space-wavenumber domain

Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.

Three-dimensional Building Damage Simulation Based on GIS

GIS technology has been applied to building damage analysis around the world. However, most previous studies focused on the application of 2-D GIS technology, and the results from traditional earthquake damage prediction are displayed in 2-D figures and charts, which is incapable of demonstrating the 3-D spatial characteristics of buildings. Taking brick-concrete building as an example, we study the characteristics of building damage, and effectively combine the information of building textures and earthquake damage. Then, we apply Google SketchUp techniques to create building models and display them with seismic damage texture in the ArcGIS Engine software development environment. In this paper we propose a solid idea for 3-D simulation of earthquake damage, which is helpful in earthquake damage prediction, virtual emergency rescue practice and earthquake knowledge education.

大型高架桥无支架盖梁施工过程模板受力特性的3D-FEM模拟分析

以上海S32高速公路闵浦大桥西侧引桥段的盖梁施工为例,通过合理的边界条件处理,基于ABAQUS有限元分析软件建立了复杂结构的3D-FEM模型,并以混凝土弹性模量和泊松比的变化模拟实际施工,进而分析无支架盖梁的受力状态.计算结果与现场实测的结果基本吻合,为无支架盖梁施工技术的研究分析提供了参考方法.研究表明:混凝土浇注完成约10 h左右,即具有90%的设计强度,基本能承受自身重力荷载,盖梁的竖向及横向荷载逐渐向混凝土转移;同时,该时间的确定为现场确定拆模时间、提高施工进度提供了依据.

Reconstruction of 3-D digital cores using a hybrid method

A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method which combines sedimentation simulation and simulated annealing （SA） method to generate 3-D digital cores based on 2-D images of rocks. The method starts with the sedimentation simulation to build a 3-D digital core, which is the initial configuration for the SA method. We update the initial digital core using the SA method to match the auto-correlation function of the 2-D rock image and eventually build the final 3-D digital core. Compared with the typical SA method, the hybrid method has significantly reduced the computation time. Local porosity theory is applied to quantitatively compare the reconstructed 3-D digital cores with the X-ray micro CT 3-D images. The results indicate that the 3-D digital cores reconstructed by the hybrid method have homogeneity and geometric connectivity similar to those of the X-ray micro CT image. The formation factors and permeabilities of the reconstructed 3-D digital cores are estimated using the finite element method （FEM） and lattice Boltzmann method （LBM）, respectively. The simulated results are in good agreement with the experimental measurements. Comparison of the simulation results suggests that the digital cores reconstructed by the hybrid method more closely reflect the true transport properties than the typical SA method alone.

Numerical simulation and analysis of solid-liquid two-phase threedimensional unsteady flow in centrifugal slurry pump

Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.

Simulation of Polymer Melt Flow Fields in Intermeshing Co-Rotating Three-Screw Extruders

Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package POLYFLOW. Based on the velocity fields calculated, the particle trajectories in both machines are visualized using particle tracking technique. The numerical results indicate that the flow patterns in three-screw extruders are similar to those in twin-screw extruders. The triangularly-arranged three-screw extruder has the largest pumping capacity and also the highest extrusion stability in terms of flowrate fluctuation with screw rotation. The instantaneous mixing and cumulative residence time distribution （RTD） characteristics are also analyzed and compared with traditional intermeshing co-rotating twin-screw extruders. It is shown that the start section of the cumulative RTD curve for the triangularly-arranged machine has a small shoulder, which is attributed to the faster flow in the central region of this type of extruder.

ANALYSIS OF SLAB EDGING BY A 3-D RIGID VISCO-PLASTIC FINITE ELEMENT METHOD

3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Cholesky decomposition of the stiffness matrix with thesolution of the equations for velocity increment by the conjugate gradient method is combined. Thistechnique, termed the shifted ICCG method, is then employed to solve the slab edging problem. Theperformance of this algorithm in terms of the number of iterations, friction variation, shiftedparameter psi and the results of simulation for processing parameters are analysed. Numerical testsand application of this technique verify the efficiency and stability of the shifted ICCG method inthe analysis of slab edging.

Comparison of Fully 3-D,Simplified 3-D and 2-D Numerical Simulations on a Circulating Fluidized Bed Boiler

Based on the fully three-dimensional(3-D)and two-dimensional(2-D)comprehensive CFD(Computational Fluid Dynamics)combustion models for a circulating fluidized bed boiler,a simplified 3-D computational domain considering the corrections of furnace side wall openings is proposed.It aims to compensate for the deficiencies of the large amount of computation in the fully 3-D model and improve the air and gas flow treatments at the openings in the simplified 2-D model.Three different computational domains,named as the fully 3-D model,simplified 3-D model and 2-D model,were implemented to perform a comparative CFD analysis in an ultra-supercritical circulating fluidized bed boiler including the hydrodynamics,penetration depth of secondary air,temperature and species distribution.The simulation results computed by the simplified3-D model yield better agreement with the fully 3-D simulation results than those of the 2-D model.The simplified 3-D model is recommended as an alternative computational domain for the conventional 2-D model in the numerical simulation of large-scale circulating fluidized bed boiler.

Aseismic ridge subduction and flat subduction:Insights from three-dimensional numerical models

Flat subduction can significantly influence the distribution of volcanism,stress state,and surface topography of the overriding plate.However,the mechanisms for inducing flat subduction remain controversial.Previous two-dimensional(2-D)numerical models and laboratory analogue models suggested that a buoyant impactor(aseismic ridge,oceanic plateau,or the like)may induce flat subduction.However,three-dimensional(3-D)systematic studies on the relationship between flat subduction and buoyant blocks are still lacking.Here,we use a 3-D numerical model to investigate the influence of the aseismic ridge,especially its width(which is difficult to consider in 2-D numerical models),on the formation of flat subduction.Our model results suggest that the aseismic ridge needs to be wide and thick enough to induce flat subduction,a condition that is difficult to satisfy on the Earth.We also find that the subduction of an aseismic ridge parallel to the trench or a double aseismic ridge normal to the trench has a similar effect on super-wide aseismic ridge subduction in terms of causing flat subduction,which can explain the flat subduction observed beneath regions such as Chile and Peru.