Determination method of mesh size for numerical simulation of blast load in near-ground detonation

In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.

Numerical Simulation of Combustion in 660MWTangentially Fired Pulverized Coal Boiler on Ultra-Low Load Operation

In this paper,the combustion conditions in the boiler furnace of a 660 MWtangential fired pulverized coal boiler are numerically simulated at 15%and 20%rated loads,to study the flexibility of coal-fired power units on ultra-low load operation.The numerical results show that the boiler can operate safely at 15%and 20%ultra-low loads,and the combustion condition in the furnace is better at 20%load,and the tangent circles formed by each characteristic section in the furnace are better,and when the boiler load is decreased to 15%,the tangent circles in the furnace begin to deteriorate.The average flue gas temperature of different areas of the furnace shows that when the boiler furnace operates under ultra-low load conditions,the average smoke temperature of the cold ash hopper at 20%load is basically the same as the average smoke temperature at 15%load;in the burner area,the average smoke temperature of the cold ash hopper at 20%load is about 50 K higher than that at 15%load;in the burned out area,the average smoke temperature of the cold ash hopper at 20%load is slightly higher than that at 15%load.The average temperature of flue gas in the furnace showed a tendency to increase rapidly with the height of the furnace,then slow down and fluctuate the temperature in the burner area,and finally increase slightly in the burnout area due to the further combustion of combustible components to release heat,and then began to decrease.

A Time-Domain Numerical Simulation for Free Motion Responses of Two Ships Advancing in Head Waves

The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.

Explosion resistance performance of reinforced concrete box girder coated with polyurea:Model test and numerical simulation

To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.

Numerical Simulation of Bed Load and Suspended Load Sediment Transport Using Well-Balanced Numerical Schemes

Sediment transport can be modelled using hydrodynamic models based on shallow water equations coupled with the sediment concentration conservation equation and the bed con-servation equation.The complete system of equations is made up of the energy balance law and the Exner equations.The numerical solution for this complete system is done in a seg-regated manner.First,the hyperbolic part of the system of balance laws is solved using a finite volume scheme.Three ways to compute the numerical flux have been considered,the Q-scheme of van Leer,the HLLCS approximate Riemann solver,and the last one takes into account the presence of non-conservative products in the model.The discretisation of the source terms is carried out according to the numerical flux chosen.In the second stage,the bed conservation equation is solved by using the approximation computed for the system of balance laws.The numerical schemes have been validated making comparisons between the obtained numerical results and the experimental data for some physical experiments.The numerical results show a good agreement with the experimental data.

State-of-the-art on the anchorage performance of rock bolts subjected to shear load

Rock bolts are extensively utilized in underground engineering as a means of offering support and stability to rock masses in tunnels,mines,and other underground structures.In environments of high ground stress,faults or weak zones can frequently arise in rock formations,presenting a significant challenge for engineering and potentially leading to underground engineering collapse.Rock bolts serve as a crucial structural element for the transmission of tensile stress and are capable of withstanding shear loads to prevent sliding of weak zones within rock mass.Therefore,a complete understanding of the behavior of rock bolts subjected to shear loads is essential.This paper presents a state-of-the-art review of the research progress of rock bolts subjected to shear load in three categories:experiment,numerical simulation,and analytical model.The review focuses on the research studies and developments in this area since the 1970s,providing a comprehensive overview of numerous factors that influence the anchorage performance of rock bolts.These factors include the diameter and angle of the rock bolt installation,rock strength,grouting material,bolt material,borehole diameter,rock bolt preload,normal stress,joint surface roughness and joint expansion angle.The paper reviews the improvement of mechanical parameter setting in numerical simulation of rock bolt shear.Furthermore,it delves into the optimization of the analytical model concerning rock bolt shear theory,approached from the perspectives of both Elastic foundation beam theory coupled with Elastoplasticity theory and Structural mechanic methods.The significance of this review lies in its ability to provide insights into the mechanical behavior of rock bolts.The paper also highlights the limitations of current research and guidelines for further research of rock bolts.

Numerical simulation on behavior of pile foundations under cyclic axial loads

On the basis of the two dimensional finite element analysis model, the pile foundations' mechanical effect of the rigid pile composite foundation under the dynamic load was researched. Through the research, the development law and deformation property of axial force of pile body, shaft resistance of pile, and cumulative settlement of pile head under vertical cyclic dynamic loads were concluded. Through the comparison and analysis of the test results of dynamic models, the test results of Poulos(1989) and cumulative settlement model of the single pile under cyclic loads were confirmed. Based on the above research, Fortran language was adopted to introduce the soil attenuation factor, the secondary development of relevant modules of ABAQUS was carried out, and the effect of soil attenuation factor on dynamic property of pile-soil was discussed further.

Novel Models and Algorithms of Load Balancing for Variable-structured Collaborative Simulation under HLA/RTI

High level architecture(HLA) is the open standard in the collaborative simulation field. Scholars have been paying close attention to theoretical research on and engineering applications of collaborative simulation based on HLA/RTI, which extends HLA in various aspects like functionality and efficiency. However, related study on the load balancing problem of HLA collaborative simulation is insufficient. Without load balancing, collaborative simulation under HLA/RTI may encounter performance reduction or even fatal errors. In this paper, load balancing is further divided into static problems and dynamic problems. A multi-objective model is established and the randomness of model parameters is taken into consideration for static load balancing, which makes the model more credible. The Monte Carlo based optimization algorithm(MCOA) is excogitated to gain static load balance. For dynamic load balancing, a new type of dynamic load balancing problem is put forward with regards to the variable-structured collaborative simulation under HLA/RTI. In order to minimize the influence against the running collaborative simulation, the ordinal optimization based algorithm(OOA) is devised to shorten the optimization time. Furthermore, the two algorithms are adopted in simulation experiments of different scenarios, which demonstrate their effectiveness and efficiency. An engineering experiment about collaborative simulation under HLA/RTI of high speed electricity multiple units(EMU) is also conducted to indentify credibility of the proposed models and supportive utility of MCOA and OOA to practical engineering systems. The proposed research ensures compatibility of traditional HLA, enhances the ability for assigning simulation loads onto computing units both statically and dynamically, improves the performance of collaborative simulation system and makes full use of the hardware resources.

Numerical simulation of sand load applied on high-speed train in sand environment

High-speed train running in the sand environment is different from the general environment. In the former situation, there will be sand load applied on high-speed train(SLAHT) caused by sand particles hitting train surface. This will have a great impact on the train stability, running drag and surface corrosion. Numerical simulation method of SLAHT in sand environment is studied. The velocity and mass flow rate models of saltation and suspension sand particles and the calculation model of SLAHT caused by sand particles hitting train surface are established. The discrete phase method is adopted for numerical simulating the process of saltation and suspension sand particles moving to train surface and generating sand load. By comparison with the field tests, the numerical simulation reliability is analysed. The theoretical formula of SLAHT changing with cross-wind and train speed is proposed. SLAHT changing law is analyzed. Research results indicate that SLAHT changing with cross-wind and train speed is a quadratic relationship. When train speed is constant, SLAHT increases quadratically with cross-wind speed improvement. When cross-wind speed is constant, SLAHT increases quadratically with train speed improvement.

Effects of loading waveforms on rock damage using particle simulation method

The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at the mesoscopic scale, such as loading frequency, stress amplitude, mean stress, confining pressure and loading sequence, are also investigated with sinusoidal waveform in detail. The related numerical results have demonstrated that： 1） the loading waveform has a certain effect on rock failure processes. The square waveform has the most damage within these waveforms, while the triangle waveform has less damage than sinusoidal waveform. In each cycle, the number of microscopic cracks increases in the loading stage, while it keeps nearly constant in the unloading stage. 2） The loading frequency, stress amplitude, mean stress, confining pressure and loading sequence have considerable effects on rock damage subjected to cyclic loading. The higher the loading frequency, stress amplitude and mean stress, the greater the damage the rock accumulated; in contrast, the lower the confining pressure, the greater the damage the rock has accumulated. 3） There is a threshold value of mean stress and stress amplitude, below which no further damage accumulated after the first few cycle loadings. 4） The high-to-low loading sequence has more damage than the low-to-high loading sequence, suggesting that the rock damage is loading-path dependent.

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.

Experiment and simulation on the pyrite removal from the recirculating load of pulverizer with a dilute phase gas-solid fluidized bed

In order to reduce the energy consumption and subsequent air pollution of coal-fired power station, based on the analysis to size and density distribution of particles from the recirculating load of the classifier of pulverizer, the separation experiment on sampling material from power plant with a dilute phase fluidized bed to remove pyrite and other minerals and numerical simulation on the separation process were done. The results show that the minimum fluidization velocity is 1.62 cm/s. Pyrite and other minerals in the material are separated. Ash of the upper and bottom layer material account for 33.34% and 73.42% respectively and sulfur content occupy 1.12% and 8.96% respectively. Scanning electron microscopy and spectroscopy tests show that sulfur in the bottom material exist in the form of pyrite. Numerical simulation on the flow field form of the dilute phase separation bed with gas-solid two phase and particle motion verifies the experimental results.

Finite element simulation of aluminum alloy cross valve forming by multi-way loading

Deformation behavior,temperature evolution and coupled effects have a significant influence on forming process and quality of component formed,which are very complex in forming process of aluminum alloy 7075 cross valve under multi-way loading due to the complexity of loading path and the multiplicity of associated processing parameters.A model of the process was developed under DFEORM-3D environment based on the coupled thermo-mechanical finite element method.The comparison between two process models,the conventional isothermal process model and the non-isothermal process model developed in this study,was carried out,and the results indicate that the thermal events play an important role in the aluminum alloy forming process under multi-way loading.The distributions and evolutions of the temperature field and strain filed are obtained by non-isothermal process simulation.The plastic zone and its extension in forming process of cross valve were analyzed.The results may provide guidelines for the determination of multi-way loading forming scheme and loading conditions of the forming cross valve components.

Numerical Simulation for Prediction of Wind Loads on Building Surfaces

The reliability of the numerical K ε model for determining wind pressure on building surfaces is evaluated. The solution algorithm is based on a body fitted non orthogonal curvilinear coordinate system and a staggered grid arrangement. The covariant velocity components are chosen as dependent variables. Convective fluxes are described by the Power Law Scheme. The grids are generated with an elliptic grid generator using control functions. The results compare favorably with those by Oxford wind tunnel measurements.

Investigation on automated loading of dynamic 3D heat source model for welding simulation

Since programing complex and dynamic heat source model for welding simulation is a complex job,the parametric methods are studied in this paper.Firstly,an overall flow to achieve automatically modeling welding was introduced.Secondly,an expert module rule for selecting welding heat source model was founded,which is based on simulation knowledge and experiences.Thirdly,a modularity routine method was investigated using writing with C++programing,which automatically creates subroutines of 3D dynamic heat source model for user.To realize the dynamic weld path,the local weld path coordinate system was moved in the global coordinate system and it is used to model the direction of weld gun,welding path and welding pose.The weld path data file was prepared by the automatic tool for the welding heat source subroutines.All above functions were integrated in the user interface and the connection with architecture was introduced.At last,a laser beam welding heat source modeling was automatically modeled and the weld pool geometry was compared with the reported literature.It demonstrated that the automated tool is valid for welding simulation.Since modeling became convenient for welding simulation using the tool proposed,it could be easy and useful for welding engineers to acquire the needed information.

Macro Meso Response and Stress Wave Propagation Characteristics of MCT High-Voltage Switch Under Shock load

In order to study the dynamic and electrical coupling response characteristics of Metal Oxide Semiconductor Controlled Thyristor(MCT)high-voltage switch under the synergic action of mechanical load and high voltage,the separated Hopkinson pressure bar(SHPB)test system was used to simulate different impact load environments,and combined with the multi-layer high-voltage ceramic capacitor charging and discharging system,the instantaneous electrical signals of MCT high-voltage switch were collected.Combined with numerical simulation and theoretical analysis,the failure mode and stress wave propagation characteristics of MCT high voltage switch were determined.The mechanical and electrical coupling response characteristics and failure mechanism of MCT high voltage switch under dynamic load were revealed from macroscopic and microscopic levels.The results show that the damage modes of MCT high-voltage switches can be divided into non-functional damage,recoverable functional damage,non-recoverable damage and structural damage.Due to the gap between the metal gate and the oxide layer,the insulating oxide layer was charged.After placing for a period of time,the elastic deformation of the metal gate partially recovered and the accumulated charge disappeared,which induced the recoverable functional damage failure of the device.In addition,obvious cracks appeared on both sides of the monocrystalline silicon inside the MCT high-voltage switch,leading to unrecoverable damage of the device.

Micro defects formation and dynamic response analysis of steel plate of quasi-cracking area subjected to explosive load

As the protective component,steel plate had attracted extensive attention because of frequently threats of explosive loads.In this paper,the evolution of microstructure and the mechanism of damage in the quasi-cracking area of steel plate subjected to explosive load were discussed and the relationships between micro defects and dynamic mechanical response were revealed.After the explosion experiment,five observation points were selected equidistant from the quasi-cracking area of the section of the steel plate along the thickness direction,and the characteristics of micro defects at the observation points were analyzed by optical microscope(OM),scanning electron microscope(SEM) and electron backscattered diffraction(EBSD).The observation result shows that many slip bands(SBs) appeared,and the grain orientation changed obviously in the steel plate,the two were the main damage types of micro defects.In addition,cracks,peeling pits,grooves and other lager micro defects were appeared in the lower area of the plate.The stress parameters of the observation points were obtained through an effective numerical model.The mechanism of damage generation and crack propagation in the quasicracking area were clarified by comparing the specific impulse of each observation point with the corresponding micro defects.The result shows that the generation and expansion of micro defects are related to the stress area(i.e.the upper compression area,the neutral plane area,and the lower tension area).The micro defects gather and expand at the grain boundary,and will become macroscopic damage under the continuous action of tensile stress.Besides,the micro defects at the midpoint of the section of the steel plate in the direction away from the explosion center(i.e.the horizontal direction) were also studied.It was found that the specific impulse at these positions were much smaller than that in the thickness direction,the micro defects were only SBs and a few micro cracks,and the those decreased with the increase of the distance from the explosion center.

Random vibration analysis of FGM plates subjected to moving load using a refined stochastic finite element method

The article introduces a finite element procedure using the bilinear quadrilateral element or four-node rectangular element(namely Q4 element) based on a refined first-order shear deformation theory(rFSDT) and Monte Carlo simulation(MCS), so-called refined stochastic finite element method to investigate the random vibration of functionally graded material(FGM) plates subjected to the moving load.The advantage of the proposed method is to use r-FSDT to improve the accuracy of classical FSDT, satisfy the stress-free condition at the plate boundaries, and combine with MCS to analyze the vibration of the FGM plate when the parameter inputs are random quantities following a normal distribution. The obtained results show that the distribution characteristics of the vibration response of the FGM plate depend on the standard deviation of the input parameters and the velocity of the moving load.Furthermore, the numerical results in this study are expected to contribute to improving the understanding of FGM plates subjected to moving loads with uncertain input parameters.

Numerical and theoretical study of load transfer behavior during cascading pillar failure

To further study the load transfer mechanism of roofemulti-pillarefloor system during cascading pillar failure(CPF),numerical simulation and theoretical analysis were carried out to study the three CPF modes according to the previous experimental study on treble-pillar specimens,e.g.successive failure mode(SFM),domino failure mode(DFM)and compound failure mode(CFM).Based on the finite element code rock failure process analysis(RFPA^(2D)),numerical models of treble-pillar specimen with different mechanical properties were established to reproduce and verify the experimental results of the three CPF modes.Numerical results show that the elastic rebound of roofefloor system induced by pillar instability causes dynamic disturbance to adjacent pillars,resulting in sudden load increases and sudden jump displacement of adjacent pillars.The phenomena of load transfer in the roofemulti-pillarefloor system,as well as the induced accelerated damage behavior in adjacent pillars,were discovered and studied.In addition,based on the catastrophe theory and the proposed mechanical model of treble-pillar specimen edisc spring group system,a potential function that characterizes the evolution characteristics of roof emulti-pillarefloor system was established.The analytical expressions of sudden jump and energy release of treble-pillar specimenedisc spring group system of the three CPF modes were derived according to the potential function.The numerical and theoretical results show good agreement with the experimental results.This study further reveals the physical essence of load transfer during CPF of roof emulti-pillarefloor system,which provides references for mine design,construction and disaster prevention.

基于Simulink对机械臂可视化上料研究

自动上料设备可以通过现代化技术实现数据采集和分析,以优化制造流程,提高生产效率和产品质量,从而更好地实现智能制造。为了实现链板等零件的自动上料过程,文章介绍了机械臂的三维可视化人机交互平台模型的构建过程,以Gluon机械臂为分析对象,进行正运动学分析和逆运动学验证,并通过Simulink搭建了控制系统模型,实现了对Gluon机械臂的联合仿真。仿真过程中,对各关节的旋转角度、角速度以及末端执行器的齐次变换矩阵实时显示。通过验证,为后期搭建实物平台提供了可行性支撑,为进一步优化和完善系统功能奠定了基础。