Numerical analysis on seismic performance of underground structures in liquefiable interlayer sites from centrifuge shaking table test

When an underground structure passes through a liquefiable soil layer,the soil liquefaction may pose a significant threat to the structure.A centrifuge shaking table test was performed to research the seismic response of underground structures in liquefiable interlayer sites,and a valid numerical model was obtained through simulation model test.Finally,the calibrated numerical model was used to perform further research on the influence of various distribution characteristics of liquefiable interlayers on the seismic reaction of underground structures.The key findings are as follows.The structure faces the most unfavorable condition once a liquefiable layer is located in the middle of the underground structure.When a liquefiable layer exists in the middle of the structure,the seismic reactions of both the underground structure and model site will increase with the rise of the thickness of the liquefiable interlayer.The inter-story drift of the structure in the non-liquefiable site is much smaller than that in the liquefiable interlayer site.The inter-story drift of the structure is not only associated with the site displacement and the soil-structure stiffness ratio but also closely associated with the slippage of the soil-structure contact interface under the condition of large deformation of the site.

Numerical investigation of the mechanical behavior of the backfill–rock composite structure under triaxial compression

To ensure safe and economical backfill mining,the mechanical response of the backfill–rock interaction system needs to be understood.The numerical investigation of the mechanical behavior of backfill–rock composite structure(BRCS)under triaxial compression,which includes deformation,failure patterns,strength characteristics,and acoustic emission(AE)evolution,was proposed.The models used in the tests have one rough interface,two cement–iron tailings ratios(CTRs),four interface angles(IAs),and three confining pressures(CPs).Results showed that the deformation,strength characteristics,and failure patterns of BRCS under triaxial compression depend on IA,CP,and CTR.The stress–strain curves of BRCS under triaxial compression could be divided into five stages,namely,compaction,elasticity,yield,strain softening,and residual stress.The relevant AE counts have corresponding relationships with different stages.The triaxial compressive strengths of composites increase linearly with the increase of the CP.Furthermore,the CP stress strengthening effect occurs.When the IAs are45°and 60°,the failure areas of composites appear in the interface and backfill.When the IAs are 75°and 90°,the failure areas of composites appear in the backfill,interface,and rock.Moreover,the corresponding failure modes yield the combined shear failure.The research results provide the basis for further understanding of the stability of the BRCS.

Design and mechanical properties analysis of a cellular Waterbomb origami structure

Cellular structures are commonly used to design energy-absorbing structures,and origami structures are becominga prevalent method of cellular structure design.This paper proposes a foldable cellular structure based on theWaterbomb origami pattern.The geometrical configuration of this structure is described.Quasi-static compressiontests of the origami tube cell of this cellular structure are conducted,and load-displacement relationship curvesare obtained.Numerical simulations are carried out to analyze the effects of aspect ratio,folding angle,thicknessand number of layers of origami tubes on initial peak force and specific energy absorption(SEA).Calculationformulas for initial peak force and SEA are obtained by the multiple linear regression method.The degree ofinfluence of each parameter on the mechanical properties of the single-layer tube cell is compared.The resultsshow that the cellular structure exhibits negative stiffness and periodic load-bearing capacity,as well as foldingangle has the most significant effect on the load-bearing and energy-absorbing capacity.By adjusting the designparameters,the stiffness,load-bearing capacity and energy absorption capacity of this cellular structure can beadjusted,which shows the programmable mechanical properties of this cellular structure.The foldability andthe smooth periodic load-bearing capacity give the structure potential for application as an energy-absorbing structure.

Structure Optimization of a Tesla Turbine Using an Organic Rankine Cycle Technology

The so-called ORC(Organic Rankine Cycle)heat recovery technology has attracted much attention with regard to medium and low temperature waste heat recovery.In the present study,it is applied to a Tesla turbine.At the same time,the effects of the disc speed,diameter and inter-disc gap on the internal flow field and output power of the turbine are also investigated by means of CFD(Computational Fluid Dynamics)numerical simulation,by which the pressure,velocity,and output efficiency of the internal flow field are obtained under different internal and external conditions.The highest efficiency(66.4%)is obtained for a number of nozzles equal to 4,a disk thickness of 1 mm,and a gap of 1 mm between the disks.The results of the study serve as a theoretical basis for the structural design and optimization of Tesla turbines.

Mechanical responses of underground carriageway structures due to construction of metro tunnels beneath the existing structure:A case study

To understand the mechanical response pattern of the existing structure and ground due to the construction of metro tunnels underneath,the finite difference method is adopted to study the torsional deformation and stress variation of the existing structure and the effect of underground carriageway structures on the surface subsidence.The curves of the maximum differential subsidence,torsion angle,and distortion of the cross-section of the existing structure show two peaks in succession during traversing of two metro tunnels beneath it.The torsion angle of the existing structure changes when the two tunnels traverse beneath it in opposite directions.The first traversing of the shield tunnel mainly induces the magnitude variation in torsional deformation of the existing structure,but the second traversing of the subsurface tunnel may cause a dynamic change in the magnitude and form of torsional deformation in the existing structure.The shielding effect can reduce the surface subsidence caused by metro tunnel excavation to a certain extent,and the development trend of subsidence becomes slower as the excavation continues.

Optimization of Channel Structure of Alkaline Water Electrolyzer by Using an Expanded Mesh as a Bipolar Plate

Alkaline water electrolysis(AWE)is the most mature technology for hydrogen production by water electrolysis.Alkaline water electrolyzer consists of multiple electrolysis cells,and a single cell consists of a diaphragm,electrodes,bipolar plates and end plates,etc.The existing industrial bipolar plate channel is concave-convex structure,which is manufactured by complicated and high-cost mold punching.This structure still results in uneven electrolyte flow and low current density in the electrolytic cell,further increasing in energy consumption and cost of AWE.Thereby,in this article,the electrochemical and flow model is firstly constructed,based on the existing industrial concave and convex flow channel structure of bipolar plate,to study the current density,electrolyte flow and bubble distribution in the electrolysis cell.The reliability of the model was verified by comparison with experimental data in literature.Among which,the electrochemical current density affects the bubble yield,on the other hand,the generated bubbles cover the electrode surface,affecting the active specific surface area and ohmic resistance,which in turn affects the electrochemical reaction.The result indicates that the flow velocity near the bottom of the concave ball approaches zero,while the flow velocity on the convex ball surface is significantly higher.Additionally,vortices are observed within the flow channel structure,leading to an uneven distribution of electrolyte.Next,modelling is used to optimize the bipolar plate structure of AWE by simulating the electrochemistry and fluid flow performances of four kinds of structures,namely,concave and convex,rhombus,wedge and expanded mesh,in the bipolar plate of alkaline water electrolyzer.The results show that the expanded mesh channel structure has the largest current density of 3330 A/m^(2)and electrolyte flow velocity of 0.507 m/s in the electrolytic cell.Under the same current density,the electrolytic cell with the expanded mesh runner structure has the smallest potential and energy consumption.This work provides a useful guide for the comprehensive understanding and optimization of channel structures,and a theoretical basis for the design of large-scale electrolyzer.

Numerical Simulation Study of Vibration Characteristics of Cantilever Traffic Signal Support Structure under Wind Environment

Computational fluid dynamics(CFD)and the finite element method(FEM)are used to investigate the wind-driven dynamic response of cantilever traffic signal support structures as a whole.By building a finite element model with the same scale as the actual structure and performing modal analysis,a preliminary understanding of the dynamic properties of the structure is obtained.Based on the two-way fluid-structure coupling calculation method,the wind vibration response of the structure under different incoming flow conditions is calculated,and the vibration characteristics of the structure are analyzed through the displacement time course data of the structure in the crosswind direction and along-wind direction.The results show that the maximum response of the structure increases gradually with the increase of wind speed under 90°wind direction angle,showing a vibration dispersion state,and the vibration response characteristics are following the vibration phenomenon of galloping;under 270°wind direction angle,the maximum displacement response of the structure occurs at the lower wind speed of 5 and 6m/s,and the vibration generated by the structure is vortex vibration at this time;the displacement response of the structure in along-wind direction increaseswith the increase of wind speed.The along-wind displacement response of the structure will increase with increasing wind speed,and the effective wind area and shape characteristics of the structurewill also affect the vibration response of the structure.

Furnace structure analysis for copper flash continuous smelting based on numerical simulation

According to the innate characteristic of four types of furnace, the copper flash continuous smelting （CFCS） furnace can be considered a synthetic reactor of two relatively independent processes： flash matte smelting process （FMSP） and copper continuous converting process （CCCP）. Then, the CFCS thermodynamic model was proposed by establishing the multi-phase equilibrium model of FMSP and the local-equilibrium model of CCCP, respectively, and by combining them through the smelting intermediates. Subsequently, the influences of the furnace structures were investigated using the model on the formation of blister copper, the Fe3O4 behavior, the copper loss in slag and the copper recovery rate. The results show that the type D furnace, with double flues and a slag partition wall, is an ideal CFCS reactor compared with the other three types furnaces. For CFCS, it is effective to design a partition wall in the furnace to make FMSP and CCCP perform in two relatively independent zones, respectively, and to make smelting gas and converting gas discharge from respective flues.

NUMERICAL MODEL OF WAVE FORCES ON CONTINUOUS CYLINDER STRUCTURES

A numerical model of wave force upon continuous cylinder structures with a large diameter using the boundary element method (BEM) is presented. A numerical model of reflecting wave upon continuous cylinders was established on the basis of linear wave theory.The fundamental solution to the Helmholtz equation within an infinite strip area that explicitly satisfies two infinite parallel boundaries is used together with Radiation condition rather than the solution of an infinite area.According to the proposed theory and method,the computer programs have been composed in Visual C ++ Development Studio.Several examples show that the technique and its program are feasible and efficient.And the wave forces upon continuous cylinders can be decreased by as much as 14%～24% under a ratio of D/L= 0.09～0.19 compared with the square caissons.

VIBRATION NUMERICAL ANALYSIS OF COUNTER-ROTATING TURBINE WITH WAKE-FLOW USING FLUID-STRUCTURE INTERACTION METHOD

The design of counter-rotating turbine is one of new techniques to improve the thrust-weight ratio of jet propulsion engines.Numerical analysis of a low pressure（LP）counter-rotating turbine rotor blade is presented by using ANSYS/CFX software.Interaction of aerodynamics and solid mechanics coupling in the computation is applied.In some rating of turbine,stress distribution and vibration characteristics of low pressure turbine（LPT）blade are computed.The wake aerodynamic forces and LPT blade vibration are transformed in frequency domain using fast Fourier transform（FFT）method.The results show that under wake aerodynamic force excitation,the first order modal vibration is more easily aroused and the higher order response cannot be ignored.Moreover,with different temperature fields,the vibration responses of blade are also different.

Prediction of atomization characteristics of pressure swirl nozzle with different structures

The structure of the pressure swirl nozzle is an important factor affecting its spray performance.This work aims to study pressure swirl nozzles with different structures by experiment and simulation.In the experiment,10 nozzles with different structures are designed to comprehensively cover various geometric factors.In terms of simulation,steady-state simulation with less computational complexity is used to study the flow inside the nozzle.The results show that the diameter of the inlet and outlet,the direction of the inlet,the diameter of the swirl chamber,and the height of the swirl chamber all affect the atomization performance,and the diameter of the inlet and outlet has a greater impact.It is found that under the same flow rate and pressure,the geometric differences do have a significant impact on the atomization characteristics,such as spray angle and SMD(Sauter mean diameter).Specific nozzle structures can be customized according to the actual needs.Data analysis shows that the spray angle is related to the swirl number,and the SMD is related to turbulent kinetic energy.Through data fitting,the equations for predicting the spray angle and the SMD are obtained.The error range of the fitting equation for the prediction of spray angle and SMD is within 15% and 10% respectively.The prediction is expected to be used in engineering to estimate the spray performance at the beginning of a real project.

NUMERICAL SIMULATION ON THE WIND FIELD STRUCTURE OF A MOUNTAINOUS AREA BESIDE SOUTH CHINA SEA DURING THE LANDFALL OF TYPHOON MOLAVE

Leveraging the commercial CFD software FLUENT,the fine-scale three-dimensional wind structure over the Paiya Mountains on the Dapeng Peninsula near Shenzhen,a city on the seashore of South China Sea,during the landfall of Typhoon Molave has been simulated and analyzed.Through the study,a conceptual wind structure model for mountainous areas under strong wind condition is established and the following conclusions are obtained as follows:(1)FLUENT can reasonably simulate a three-dimensional wind structure over mountainous areas under strong wind conditions;(2)the kinetic effect of a mountain can intensify wind speed in the windward side of the mountain and the area over the mountain peak;and(3)in the leeward side of the mountain,wind speed is relatively lower with relatively stronger wind shear and turbulence.

Mesh Size Effect in Numerical Simulation of Blast Wave Propagation and Interaction with Structures

Numerical method is popular in analysing the blast wave propagation and interaction with structures.However,because of the extremely short duration of blast wave and energy trans-mission between different grids,the numerical results are sensitive to the finite element mesh size.Previous numerical simulations show that a mesh size acceptable to one blast scenario might not be proper for another case,even though the difference between the two scenarios is very small,indicating a simple numerical mesh size convergence test might not be enough to guarantee accu-rate numerical results.Therefore,both coarse mesh and fine mesh were used in different blast scenarios to investigate the mesh size effect on numerical results of blast wave propagation and interaction with structures.Based on the numerical results and their comparison with field test re-sults and the design charts in TM5-1300,a numerical modification method was proposed to correct the influence of the mesh size on the simulated results.It can be easily used to improve the accu-racy of the numerical results of blast wave propagation and blast loads on structures.

Numerical simulation and control of welding distortion for double floor structure of high speed train

The welding heat source models and the plastic tension zone sizes of a typical weld joint involved in the double floor structure of high speed train under different welding parameters were calculated by a thermal-elastic-plastic FEM analysis based on SYSWELD code.Then,the welding distortion of floor structure was predicted using a linear elastic FEM and shrinkage method based on Weld Planner software.The effects of welding sequence,clamping configuration and reverse deformation on welding distortion of floor structure were examined numerically.The results indicate that the established elastic FEM model for floor structure is reliable for predicting the distribution of welding distortion in view of the good agreement between the calculated results and the measured distortion for real double floor structure.Compared with the welding sequence,the clamping configuration and the reverse deformation have a significant influence on the welding distortion of floor structure.In the case of30 mm reverse deformation,the maximum deformation can be reduced about 70%in comparison to an actual welding process.

Basin Structure and Numerical Simulation for the Mechanisms of Seismic Disasters

In the present study seismic wave propagation in heterogeneous media is numerically simulated by using the pseudospectral method with the staggered grid RFFT differentiation in order to clarify the cause for the complicated distribution characteristics of strong ground motion in regions with basin structure. The results show that the maximum amplitudes of simulated ground acceleration waveforms are closely related to the basin structure. Interference of seismic waves in the basin strongly affects the distribution of maximum seismic waveforms, which may result in peak disasters during earthquakes. Peak disasters might be away from basin boundaries or earthquake faults. Seismic energy transmitted into the basin from the bedrock can hardly penetrate the bottom of the basin and then travel back into the bedrock region. The seismic energy is absorbed by basin media, and transferred into the kinematical energy of seismic waves with great amplitude in the basin. Seismic waves between basins may result in serious damage to buildings over the basin. This is significant for aseismatic research. Geological surveys in and around urban areas would benefit aseismatic research and mitigation of seismic disasters of a city. Such geological surveys should involve seismic velocity structure in the media above the bedrock besides such subjects as active faults and geological structure.

Numerical modeling of deep-seated landslides interacting with man-made structures

This paper describes the interaction between deep-seated landslides and man-made structures such as dams, penstocks, viaducts, and tunnels. Selected case studies are reported first with the intent to gain insights into the complexities associated with the interaction of these structures with deep-seated landslides(generally referred to as deep-seated gravity slope deformations, DSGSDs). The main features, which characterize these landslides, are mentioned together with the interaction problems encountered in each case. Given the main objective of this paper, the numerical modeling methods adopted are outlined as means for increase in the understanding of the interaction problems being investigated. With the above in mind, the attention moves to an important and unique case history dealing with the interaction of a large-size twin-tunnel excavated with an earth pressure balance(EPB)tunnel boring machine(TBM) and a deep-seated landslide, which was reactivated due to the stress changes induced by tunnel excavation in landslide shear zone. The geological and geotechnical conditions are described together with the available monitoring data on the landslide movements, based on the advanced and conventional monitoring tools used. Numerical modeling is illustrated as an aid to back-analyze the monitored surface and subsurface deformations and to assist in finding the appropriate engineering solution for putting the tunnel into service and as a follow-up means for future understanding and control of the interaction problems. The simulation is based on a novel time-dependent model representing the landslide behavior.

Numerical Simulation of Dynamic Response and Collapse for Steel Frame Structures Subjected to Blast Load

The progressive collapse of steel frame structures under the blast load was investigated using LS-DYNA. The multi-material Eulerian and Lagrangian coupling algorithm was adopted. A flu-id-structure coupling finite element model was established which consists of Lagrange element for simulating steel frame structures and concrete ground, multiple ALE element for simulating air and TNT explosive material. Numerical simulations of the blast pressure wave propagation, struc-tural dynamic responses and deformation, and progressive collapse of a five-story steel frame structure in the event of an explosion near above ground were performed. The numerical analysis showed that the Lagrangian and Eulerian coupling algorithm gave good simulations of the shock wave propagation in the mediums and blast load effects on the structure. The columns subjected to blast load may collapse by shear yielding rather than by flexural deformation. The columns and joints of steel beam to column in the front steel frame structure generated enormous plastic defor-mation subjected to intensive blast waves, and columns lost carrying capacity, subsequently lead-ing to the collapse of the whole structure. The approach coupling influence between struc-tural deformation and fluid load well simulated the progressive collapse process of structures, and provided an effective tool for analyzing the collapse mechanism of the steel frame structure under blast load.

Investigation on the Mechanical Properties and Shape Memory Effect of Landing Buffer Structure Based on NiTi Alloy Printing

With the deepening of human research on deep space exploration,our research on the soft landing methods of landers has gradually deepened.Adding a buffer and energy-absorbing structure to the leg structure of the lander has become an effective design solution.Based on the energy-absorbing structure of the leg of the interstellar lander,this paper studies the appearance characteristics of the predatory feet of the Odontodactylus scyllarus.The predatory feet of the Odontodactylus scyllarus can not only hit the prey highly when preying,but also can easily withstand the huge counter-impact force.The predatory feet structure of the Odontodactylus scyllarus,like a symmetrical cone,shows excellent rigidity and energy absorption capacity.Inspired by this discovery,we used SLM technology to design and manufacture two nickel-titanium samples,which respectively show high elasticity,shape memory,and get better energy absorption capacity.This research provides an effective way to design and manufacture high-mechanical energy-absorbing buffer structures using bionic 3D printing technology and nickel-titanium alloys.

Numerical Simulation of Large Diameter Cylindrical Structure Slamming

The water entry of large diameter cylindrical structure is studied by applying numerical simulation method. The processes of different diameter cyhndrical structures impacting water with various constant velocities are calculated numerically. Thereafter, analyzed are the distribution of slamming pressure on structure during slamming course and the influence of slamming velocity and cylindrical diameter on slamming process. Furthermore, presented herein is an equation being used to forecast the peak slamming force on a large diameter cylindrical structure.

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.