An overall vector radiative transfer theory was developed for numerical modeling, in both active and passive microwave remote sensing. The Theory and approaches are briefly summerized.To quantitatively understand scat...An overall vector radiative transfer theory was developed for numerical modeling, in both active and passive microwave remote sensing. The Theory and approaches are briefly summerized.To quantitatively understand scattering and thermal emission from targets in active and passive remote sensing, we have developed an overall vector radiative transfer theory for a set of theoretical models of discrete scatterer and continuous random media for the earth terrain (wet soil, vegetation, snow, sea-ice, etc.) and atmosphere, and numerical approaches for simulation, data analysis, and parameter sensitivity test. Our numerical results favorably agreed with experimental data in microwave re mote sensing of various earth surfaces. Main approaches are briefly summerized herewith.展开更多
With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ...With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.展开更多
The NATO agreement STANAG 4569 defines the protection levels for the occupants of logistic and light armored vehicle.The Allied Engineering Publication,AEP-55,Volume 2 document outlines the test conditions for underbe...The NATO agreement STANAG 4569 defines the protection levels for the occupants of logistic and light armored vehicle.The Allied Engineering Publication,AEP-55,Volume 2 document outlines the test conditions for underbelly improvised explosive device(IEDs),which must be buried in water-saturated sandy gravel.The use of sandy gravel has some drawbacks,for instance reproducibility,time consumption,and cost.This paper focuses on the investigation of four alternatives to sandy gravel,which could produce similar specific and cumulative impulses:a concrete pot filled with water,a concrete pot filled with quartz sand,a steel pot without filling and a concrete pot filled with glass spheres(diameter 200μm—300μm)and different water contents.The impulses are measured with a ring technology developed at the Fraunhofer EMI.A numerical soil model based on the work of Marrs,2014 and Fi serov a,2006 and considering the soil moisture was used to simulate the experiments with glass spheres at different water contents,showing much better agreement with the experiments than the classical Laine&Sandvik model,even for high saturation levels.These results can be used to create new test conditions at original scale that are more cost-effective,more reproducible and simpler to manage in comparison to the current tests carried out with STANAG sandy gravel.展开更多
Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical e...Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical experimental measurement and numerical simulation pose research challenges.This study focuses on the ice load of a cylinder structure breaking upward through the ice sheet form underneath in the Small Ice Model Basin of China Ship Scientific Research Center(CSSRC SIMB).A high-speed camera system was employed to observe the ice sheet failure during the tests,in which,with the loading position as center,local radial cracks and circumferential cracks were generated.A load sensor was used to measure the overall ice load during this process.Meanwhile,a numerical model was developed using LS-DYNA for validation and comparison.With this model,numerical simulation was conducted under various ice thicknesses and upgoing speeds to analyze the instantaneous curves of ice load.The calculation results were statistically analyzed under different working conditions to determine the influence of the factors on the ice load of the cylinder.The study explores the measurement method about ice load of objects vertically breaking through model ice sheet and is expected to provide some fundamental insights into the safety design of underwater structures operating in ice waters.展开更多
Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process...Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.展开更多
Drying is a complicated physical process which involves simultaneous heat and mass transfer in the removal of solvents inside propellants.Inappropriate drying techniques may result in the formation of a hard skin laye...Drying is a complicated physical process which involves simultaneous heat and mass transfer in the removal of solvents inside propellants.Inappropriate drying techniques may result in the formation of a hard skin layer near the surface to block the free access of most solvent through for long stick propellants with large web thickness,which lead to lower drying efficiency and worse drying quality.This study aims to gain a comprehensive understanding of drying process and clarify the mechanism of the blocked layer near the propellant surface.A new three-dimensional coupled heat and mass transfer(3D-CHMT)model was successfully developed under transient conditions.The drying experiment results show that the 3DCHMT model could be applied to describe the drying process well since the relative error of the content of solvent between simulation and experiment values is only 5.5%.The solvent behavior simulation demonstrates that the mass transfer process can be divided into super-fast(SF)and subsequent minorfast(MF)stages,and the SF stage is vital to the prevention of the blocked layer against the free access for solvent molecules inside propellant grains.The effective solvent diffusion coefficient(Deff)of the propellant surface initially increases from 3.4×10^(-6)to 5.3×10^(-6)m^(2)/s as the temperature increases,and then decreases to 4.1×10^(-8)m^(2)/s at 60-100 min.The value of Deffof surface between 0-1.4 mm has a unique trend of change compared with other regions,and it is much lower than that of the internal at100 min under simulation conditions.Meanwhile,the temperature of the propellant surface increases rapidly at the SF stage(0-100 min)and then very slowly thereafter.Both the evolution of Deffand temperature distribution demonstrate that the blocked layer near the propellant surface has been formed in the time period of approximately 0-100 min and its thickness is about 1.4 mm.To mitigate the formation of blocked layer and improve its drying quality of finial propellant products effectively,it should be initially dried at lower drying temperature(30-40℃)in 0-100 min and then dried at higher drying temperature(50-60℃)to reduce drying time for later drying process in double base gun propellants.The present results can provide theoretical guidance for drying process and optimization of drying parameters for long stick propellants with large web thickness.展开更多
A new complementary metal oxide semiconductor UV/blue-extended photodiode was presented for light detection in the UV/blue spectral range. Photoelectric characteristics of this presented photodiode were studied by num...A new complementary metal oxide semiconductor UV/blue-extended photodiode was presented for light detection in the UV/blue spectral range. Photoelectric characteristics of this presented photodiode were studied by numerical modeling and device simulation. Technology computer aided design simulation was done first to analyze its photoelectric characteristics. The structure characteristic and depletion situation of space between two adjacent P+ anodes were discussed. The reverse characteristic, spectral response characteristic and DC characteristic were discussed in detail. For the numerical modeling, dead layer effect is considered in the building of analytical mode. Dead layer is a space in which the boron doping profile decreases towards the surface due to high doping effects and boron redistribution, which affects the sensitivity of photodiode in the UV range seriously. Reverse characteristics and spectral response characteristics were modeled and analyzed typically. At last, silicon test results were given and compared with the simulated result, which shows reasonable match for each.展开更多
The structural performance of perforated steel columns(PSCs)is significantly more complex than the one of solid web I-shaped elements under the diversity of blast loading scenarios.The damage criterion of PSCs is not ...The structural performance of perforated steel columns(PSCs)is significantly more complex than the one of solid web I-shaped elements under the diversity of blast loading scenarios.The damage criterion of PSCs is not only related to initial deformation response during the blast but also the residual axial load capacity and it can be considered as a reliable index after the blast effects.Therefore,the PSCs damages will be studied in two stages;direct and post blast effects.In the present study,the dynamic response of PSCs was numerically evaluated under different levels of blast threats using LS-DYNA software.Extensive explicit finite element(FE)analyses are carried out to investigate the effect of various parameters,such as web opening shapes,boundary conditions and strengthening details on the damage index and toughness of the PSCs compared to the parent steel sections.The results of the comparative study show that the damage and toughness decrease when the support condition changes from pinned to fixed ends through the two stages of loadings.PSCs give high toughness compared to its parent sections during blast shock stage while,a remarkable decrease in toughness is observed during the application of axial gravity after blast shock.Furthermore,the web opening shapes have slight effects on the global dynamic behavior of PSCs,particularly in terms of residual capacity.On the contrary,the retrofitting strategy using both closed holes at end and vertical stiffeners have an effective enhancement to get higher toughness in case of the extreme blasts.展开更多
Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concret...Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.展开更多
Anchoring mechanism and failure characteristics of composite soft rock with weak interface usually exhibit remarkable difference from those in single rock mass.In order to fully understand the reinforcement mechanism ...Anchoring mechanism and failure characteristics of composite soft rock with weak interface usually exhibit remarkable difference from those in single rock mass.In order to fully understand the reinforcement mechanism of composite soft roof in western mining area of China,a mechanical model of composite soft rock with weak interface and rock bolt which considering the transverse shear sliding between different rock layers was established firstly.The anchoring effect was quantified by a factor defined as anchoring effect coefficient and its evolution equation was further deduced based on the deformation relationship and homogenized distribution assumption of stress acting on composite structure.Meanwhile,the numerical simulation model of composite soft rock with shear joint was prompted by finite element method.Then detailed analysis were carried out for the deformation features,stress distribution and failure behavior of rock mass and rock bolt near the joint under transverse load.The theoretical result indicates that the anchoring effect of rock-bolt through weak joint changes with the working status of rock mass and closely relates with the physical and geometric parameters of rock mass and rock bolt.From the numerical results,the bending deformation of rock bolt accurately characterized by Doseresp model is mainly concentrated between two plastic hinges near the shear joint.The maximum tensile and compression stresses distribute in the plastic hinge.However,the maximum shear stress appears at the positions of joint surface.The failure zones of composite rock are produced firstly at the joint surface due to the reaction of rock bolt.The above results laid a theoretical and computational foundation for further study of anchorage failure in composite soft rock.展开更多
Ammonium nitrate and fuel oil(ANFO)based explosive is a classic example of non-ideal high explosives.Its detonation is characterized by a strong dependence of detonation parameters on explosive charge diameter,presenc...Ammonium nitrate and fuel oil(ANFO)based explosive is a classic example of non-ideal high explosives.Its detonation is characterized by a strong dependence of detonation parameters on explosive charge diameter,presence and characteristics of confinement,as well as incomplete consumption of explosive at the sonic point.In this work we propose a detonation model based on the Wood-Kirkwood(WK)theory coupled with the thermochemical code EXPLO5 and supplemented with reaction rate models.Our objective is to analyze the validity of the model for highly non-ideal ANFO explosives,with emphasis on effect of reaction rate models.It was found that both single-step and two-step pressure-based models can be calibrated to reproduce experimental detonation velocity-charge radius data of ANFO at radii significantly above the failure radius(i.e.for D/D_(id)>~0.6).Single-step pressure-based model,with the pressure exponent equal to 1.4,proved to be the most accurate,even in the vicinity of the failure radius.The impact of the rate models is most evident on temporal(and spatial)distribution of flow parameters in detonation driving zone,especially when it comes to the conversion and width of detonation driving zone.展开更多
A 2-D numerical model was developed to predict the shape of weld pool in stationary GTA welding of commercial pure aluminium, without considering fluid flow in the weld pool. A Gaussian current density and heat input ...A 2-D numerical model was developed to predict the shape of weld pool in stationary GTA welding of commercial pure aluminium, without considering fluid flow in the weld pool. A Gaussian current density and heat input distribution on the surface of the workpiece were considered. The parameters of Gaussian distribution were modified by comparing calculated results with experimental ones. It was found that these distribution parameters are fimctions of applied current and arc length. Effects of arc length, applied current and welding time on the geometry of the weld pool were investigated. To check the validity of the model, a series of experiments were also conducted. In general, the agreement between calculated overall shape of the weld pool and the experimental one was acceptable, especially in low applied currents. Therefore, it can be concluded that in pure aluminium, the heat conduction is dominant mechanism of heat transfer in the weld pool.展开更多
The rapid development of high-speed transportation infrastructure such as highway and high-speed railway has resulted in the advancement of soft soil improvement techniques. Vacuum preloading combined with vertical dr...The rapid development of high-speed transportation infrastructure such as highway and high-speed railway has resulted in the advancement of soft soil improvement techniques. Vacuum preloading combined with vertical drains has been proved to be an effective method in the treatment of soft foundation. A three-dimensional numerical analysis of the coupled methods was presented, in which the smear zone and the well resistance were taken into account. The variations of the basic soil parameters including the permeability coefficient and the coefficient of volume compressibility were considered in the numerical model. The result of the numerical model was then compared to the measured value. The results indicate that the decrease of coefficient of volume compressibility accelerates the consolidation of the soil while the influence of hydraulic conductivity is insignificant. A cube drain presents the closest result to the real situation compared to the other equivalent methods of prefabricated vertical drain (PVD). The case study indicates that the numerical model with variation of soil parameters is closer to the measured value than the numerical model without variation of soil parameters.展开更多
Rapid dewatering and thickening of whole-tailings with ultrafine particles is one of the most important processes for the whole-tailings paste preparation. Deep-cone thickener, a kind of such process for the flocculat...Rapid dewatering and thickening of whole-tailings with ultrafine particles is one of the most important processes for the whole-tailings paste preparation. Deep-cone thickener, a kind of such process for the flocculation and settling of whole-tailings, is particularly necessary to study. However, there exist many problems in observing the flocculation and settling process of whole-tailings, as well as the particle size distribution(PSD) of whole-tailings floccules in deep-cone thickener. Population balance model(PBM) is applied to predict the PSD in deep-cone thickener, and LUO model and GHADIRI model are employed to study the aggregation and fragmentation mechanism of the whole-tailings particles, respectively. Through three-dimensional numerical simulation on the whole-tailings flocculation and settling in deep-cone thickener using computational fluid dynamics(CFD)-PBM, the distribution of density and turbulent kinetic energy in deep-cone thickener were obtained, at the same time the spatio-temporal changes of whole-tailings floccules particle size distribution are analyzed. Finally, the major flocculation position in deep-cone thickener is found and the flocculation settling rules of whole-tailings are achieved.展开更多
From the mathematical principles, the generalized potential theory can be employed to create constitutive model of geomaterial directly. The similar Cam-clay model, which is created based on the generalized potential ...From the mathematical principles, the generalized potential theory can be employed to create constitutive model of geomaterial directly. The similar Cam-clay model, which is created based on the generalized potential theory, has less assumptions,clearer mathematical basis, and better computational accuracy. Theoretically, it is more scientific than the traditional Cam-clay models. The particle flow code PFC3 D was used to make numerical tests to verify the rationality and practicality of the similar Cam-clay model. The verification process was as follows: 1) creating the soil sample for numerical test in PFC3 D, and then simulating the conventional triaxial compression test, isotropic compression test, and isotropic unloading test by PFC3D; 2)determining the parameters of the similar Cam-clay model from the results of above tests; 3) predicting the sample's behavior in triaxial tests under different stress paths by the similar Cam-clay model, and comparing the predicting results with predictions by the Cam-clay model and the modified Cam-clay model. The analysis results show that the similar Cam-clay model has relatively high prediction accuracy, as well as good practical value.展开更多
Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the c...Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the compaction processes are to obtain a compact with the geometrical requirements, without cracks, and with a uniform distribution of density. Design of such proc esses consist, essentially, in determine the sequence and relative displacements of die and punches in order to achieve such goals. A.B. Khoei presented a gener al framework for the finite element simulation of powder forming processes based on the following aspects; a large displacement formulation, centred on a total and updated Lagrangian formulation; an adaptive finite element strategy based on error estimates and automatic remeshing techniques; a cap model based on a hard ening rule in modelling of the highly non-linear behaviour of material; and the use of an efficient contact algorithm in the context of an interface element fo rmulation. In these references, the non-linear behaviour of powder was adequately desc ribed by the cap plasticity model. However, it suffers from a serious deficiency when the stress-point reaches a yield surface. In the flow theory of plasticit y, the transition from an elastic state to an elasto-plastic state appears more or less abruptly. For powder material it is very difficult to define the locati on of yield surface, because there is no distinct transition from elastic to ela stic-plastic behaviour. Results of experimental test on some hard met al powder show that the plastic effects were begun immediately upon loading. In such mater ials the domain of the yield surface would collapse to a point, so making the di rection of plastic increment indeterminate, because all directions are normal to a point. Thus, the classical plasticity theory cannot deal with such materials and an advanced constitutive theory is necessary. In the present paper, the constitutive equations of powder materials will be discussed via an endochronic theory of plasticity. This theory provides a unifi ed point of view to describe the elastic-plastic behaviour of material since it places no requirement for a yield surface and a ’loading function’ to disting uish between loading an unloading. Endochronic theory of plasticity has been app lied to a number of metallic materials, concrete and sand, but to the knowledge of authors, no numerical scheme of the model has been applied to powder material . In the present paper, a new approach is developed based on an endochronic rate independent, density-dependent plasticity model for describing the isothermal deformation behavior of metal powder at low homologous temperature. Although the concept of yield surface has not been explicitly assumed in endochronic theory, it is shown that the cone-cap plasticity yield surface (Fig.1), which is the m ost commonly used plasticity models for describing the behavior of powder materi al can be easily derived as a special case of the proposed endochronic theory. Fig.1 Trace of cone-cap yield function on the meridian pl ane for different relative density As large deformation is observed in powder compaction process, a hypoelastic-pl astic formulation is developed in the context of finite deformation plasticity. Constitutive equations are stated in unrotated frame of reference that greatly s implifies endochronic constitutive relation in finite plasticity. Constitutive e quations of the endochronic theory and their numerical integration are establish ed and procedures for determining material parameters of the model are demonstra ted. Finally, the numerical schemes are examined for efficiency in the model ling of a tip shaped component, as shown in Fig.2. Fig.2 A shaped tip component. a) Geometry, boundary conditio n and finite element mesh; b) density distribution at final stage of展开更多
A deep understanding of the internal ballistic process and the factors affecting it is of primary importance to efficiently design a gun system and ensure its safe management. One of the main goals of internal ballist...A deep understanding of the internal ballistic process and the factors affecting it is of primary importance to efficiently design a gun system and ensure its safe management. One of the main goals of internal ballistics is to estimate the gas pressure into the combustion chamber and the projectile muzzle velocity in order to use the propellant to its higher efficiency while avoiding over-pressure phenomena. Dealing with the internal ballistic problem is a complex undertaking since it requires handling the interaction between different constituents during a transient time lapse with very steep rise of pressure and temperature. Several approaches have been proposed in the literature, based on different assumptions and techniques. Generally, depending on the used mathematical framework, they can be classified into two categories: computational fluid dynamics-based models and lumped-parameter ones. By focusing on gun systems, this paper offers a review of the main contributions in the field by mentioning their advantages and drawbacks. An insight into the limitations of the currently available modelling strategies is provided,as well as some considerations on the choice of one model over another. Lumped-parameter models, for example, are a good candidate for performing parametric analysis and optimisation processes of gun systems, given their minimum requirements of computer resources. Conversely, CFD-based models have a better capacity to address more sophisticated phenomena like pressure waves and turbulent flow effects. The performed review also reveals that too little attention has been given to small calibre guns since the majority of currently available models are conceived for medium and large calibre gun systems.Similarly, aspects like wear phenomena, bore deformations or projectile-barrel interactions still need to be adequately addressed and our suggestion is to dedicate more effort on it.展开更多
Sudden earthquakes pose a threat to the running safety of trains on high-speed railway bridges,and the stiffness of piers is one of the factors affecting the dynamic response of train-track-bridge system.In this paper...Sudden earthquakes pose a threat to the running safety of trains on high-speed railway bridges,and the stiffness of piers is one of the factors affecting the dynamic response of train-track-bridge system.In this paper,a experiment of a train running on a high-speed railway bridge is performed based on a dynamic experiment system,and the corresponding numerical model is established.The reliability of the numerical model is verified by experiments.Then,the experiment and numerical data are analyzed to reveal the pier height effects on the running safety of trains on bridges.The results show that when the pier height changes,the frequency of the bridge below the 30 m pier height changes greater;the increase of pier height causes the transverse fundamental frequency of the bridge close to that of the train,and the shaking angle and lateral displacement of the train are the largest for bridge with 50 m pier,which increases the risk of derailment;with the pier height increases from 8 m to 50 m,the derailment coefficient obtained by numerical simulations increases by 75% on average,and the spectral intensity obtained by experiments increases by 120% on average,two indicators exhibit logarithmic variation.展开更多
A three-dimensional numerical torsion shear test is presented on hollow cylinder specimen which is performed on a spherical assemblage with fixed principal stress axes using the discrete element code PFC3D.Stack wall ...A three-dimensional numerical torsion shear test is presented on hollow cylinder specimen which is performed on a spherical assemblage with fixed principal stress axes using the discrete element code PFC3D.Stack wall technique boundary conditions are employed and optimized to reasonably capture the microstructure evolution.Parametric studies are conducted in terms of the ratio κ,normal and shear stiffness of particles,wall stiffness and friction coefficients.Afterwards,in comparison with physical test,numerical results for a fixed principal stress angle(α=45°) are presented.The results show that the numerical test could capture the macro-micro mechanical behavior of the spherical particle assembly.The evolution of the coordination number demonstrates that particles in shear banding undergo remarkable decrease.The effects of localization on specimens illustrate that global stress and strain recorded from a hollow cylinder apparatus could not represent the localized response.The shearing band initiation and evolution from porosity and shear rate are visualized by contour lines in different shear strains.展开更多
Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive stre...Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive strength (UCS) and cohesion) and acting forces (e.g., field stresses and mud pressure). Accurate determination of these parameters is time-consuming, expensive and sometimes even impossible. This work offers a systematic sensitivity analysis to quantify the amount of each parameter’s effect on the stability of a wellbore. Maximum wellbore wall displacement is used as a stability factor to study the stability of a wellbore. A 3D finite difference method with Mohr model is used for the numerical modeling. The numerical model is verified against an analytical solution. A dimensionless sensitivity factor is developed in order to compare the results of various parameters in the sensitivity analysis. The results show a different order of importance of parameters based on rock strength. The most sensitive properties for a weak rock are the maximum horizontal stress, internal friction angle and formation pressure, respectively, while for a strong rock, the most sensitive parameters are the maximum horizontal stress, mud pressure and pore pressure, respectively. The amount of error in wellbore stability analysis inflicted by the error in estimation of each parameter was also derived.展开更多
基金The Project supported by National National Science FoundationYing Tung Education Foundation
文摘An overall vector radiative transfer theory was developed for numerical modeling, in both active and passive microwave remote sensing. The Theory and approaches are briefly summerized.To quantitatively understand scattering and thermal emission from targets in active and passive remote sensing, we have developed an overall vector radiative transfer theory for a set of theoretical models of discrete scatterer and continuous random media for the earth terrain (wet soil, vegetation, snow, sea-ice, etc.) and atmosphere, and numerical approaches for simulation, data analysis, and parameter sensitivity test. Our numerical results favorably agreed with experimental data in microwave re mote sensing of various earth surfaces. Main approaches are briefly summerized herewith.
文摘With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.
基金Bundeswehr Technical Center for Weapons and Ammunition WTD-91 GF-440 in Meppen,Germany for funding this work。
文摘The NATO agreement STANAG 4569 defines the protection levels for the occupants of logistic and light armored vehicle.The Allied Engineering Publication,AEP-55,Volume 2 document outlines the test conditions for underbelly improvised explosive device(IEDs),which must be buried in water-saturated sandy gravel.The use of sandy gravel has some drawbacks,for instance reproducibility,time consumption,and cost.This paper focuses on the investigation of four alternatives to sandy gravel,which could produce similar specific and cumulative impulses:a concrete pot filled with water,a concrete pot filled with quartz sand,a steel pot without filling and a concrete pot filled with glass spheres(diameter 200μm—300μm)and different water contents.The impulses are measured with a ring technology developed at the Fraunhofer EMI.A numerical soil model based on the work of Marrs,2014 and Fi serov a,2006 and considering the soil moisture was used to simulate the experiments with glass spheres at different water contents,showing much better agreement with the experiments than the classical Laine&Sandvik model,even for high saturation levels.These results can be used to create new test conditions at original scale that are more cost-effective,more reproducible and simpler to manage in comparison to the current tests carried out with STANAG sandy gravel.
文摘Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical experimental measurement and numerical simulation pose research challenges.This study focuses on the ice load of a cylinder structure breaking upward through the ice sheet form underneath in the Small Ice Model Basin of China Ship Scientific Research Center(CSSRC SIMB).A high-speed camera system was employed to observe the ice sheet failure during the tests,in which,with the loading position as center,local radial cracks and circumferential cracks were generated.A load sensor was used to measure the overall ice load during this process.Meanwhile,a numerical model was developed using LS-DYNA for validation and comparison.With this model,numerical simulation was conducted under various ice thicknesses and upgoing speeds to analyze the instantaneous curves of ice load.The calculation results were statistically analyzed under different working conditions to determine the influence of the factors on the ice load of the cylinder.The study explores the measurement method about ice load of objects vertically breaking through model ice sheet and is expected to provide some fundamental insights into the safety design of underwater structures operating in ice waters.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.
基金supported by the National Natural Science Foundation of China(Grant No.22075146)。
文摘Drying is a complicated physical process which involves simultaneous heat and mass transfer in the removal of solvents inside propellants.Inappropriate drying techniques may result in the formation of a hard skin layer near the surface to block the free access of most solvent through for long stick propellants with large web thickness,which lead to lower drying efficiency and worse drying quality.This study aims to gain a comprehensive understanding of drying process and clarify the mechanism of the blocked layer near the propellant surface.A new three-dimensional coupled heat and mass transfer(3D-CHMT)model was successfully developed under transient conditions.The drying experiment results show that the 3DCHMT model could be applied to describe the drying process well since the relative error of the content of solvent between simulation and experiment values is only 5.5%.The solvent behavior simulation demonstrates that the mass transfer process can be divided into super-fast(SF)and subsequent minorfast(MF)stages,and the SF stage is vital to the prevention of the blocked layer against the free access for solvent molecules inside propellant grains.The effective solvent diffusion coefficient(Deff)of the propellant surface initially increases from 3.4×10^(-6)to 5.3×10^(-6)m^(2)/s as the temperature increases,and then decreases to 4.1×10^(-8)m^(2)/s at 60-100 min.The value of Deffof surface between 0-1.4 mm has a unique trend of change compared with other regions,and it is much lower than that of the internal at100 min under simulation conditions.Meanwhile,the temperature of the propellant surface increases rapidly at the SF stage(0-100 min)and then very slowly thereafter.Both the evolution of Deffand temperature distribution demonstrate that the blocked layer near the propellant surface has been formed in the time period of approximately 0-100 min and its thickness is about 1.4 mm.To mitigate the formation of blocked layer and improve its drying quality of finial propellant products effectively,it should be initially dried at lower drying temperature(30-40℃)in 0-100 min and then dried at higher drying temperature(50-60℃)to reduce drying time for later drying process in double base gun propellants.The present results can provide theoretical guidance for drying process and optimization of drying parameters for long stick propellants with large web thickness.
基金Projects(61233010,61274043)supported by the National Natural Science Foundation of ChinaProject(NCET-11-0975)supported by the Program for New Century Excellent Talents in University,China
文摘A new complementary metal oxide semiconductor UV/blue-extended photodiode was presented for light detection in the UV/blue spectral range. Photoelectric characteristics of this presented photodiode were studied by numerical modeling and device simulation. Technology computer aided design simulation was done first to analyze its photoelectric characteristics. The structure characteristic and depletion situation of space between two adjacent P+ anodes were discussed. The reverse characteristic, spectral response characteristic and DC characteristic were discussed in detail. For the numerical modeling, dead layer effect is considered in the building of analytical mode. Dead layer is a space in which the boron doping profile decreases towards the surface due to high doping effects and boron redistribution, which affects the sensitivity of photodiode in the UV range seriously. Reverse characteristics and spectral response characteristics were modeled and analyzed typically. At last, silicon test results were given and compared with the simulated result, which shows reasonable match for each.
文摘The structural performance of perforated steel columns(PSCs)is significantly more complex than the one of solid web I-shaped elements under the diversity of blast loading scenarios.The damage criterion of PSCs is not only related to initial deformation response during the blast but also the residual axial load capacity and it can be considered as a reliable index after the blast effects.Therefore,the PSCs damages will be studied in two stages;direct and post blast effects.In the present study,the dynamic response of PSCs was numerically evaluated under different levels of blast threats using LS-DYNA software.Extensive explicit finite element(FE)analyses are carried out to investigate the effect of various parameters,such as web opening shapes,boundary conditions and strengthening details on the damage index and toughness of the PSCs compared to the parent steel sections.The results of the comparative study show that the damage and toughness decrease when the support condition changes from pinned to fixed ends through the two stages of loadings.PSCs give high toughness compared to its parent sections during blast shock stage while,a remarkable decrease in toughness is observed during the application of axial gravity after blast shock.Furthermore,the web opening shapes have slight effects on the global dynamic behavior of PSCs,particularly in terms of residual capacity.On the contrary,the retrofitting strategy using both closed holes at end and vertical stiffeners have an effective enhancement to get higher toughness in case of the extreme blasts.
文摘Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.
基金Projects(51774196,41472280,51578327)supported by the National Natural Science Foundation of ChinaProject(2016M592221)supported by the China Postdoctoral Science FoundationProject(BJRC20160501)supported by the SDUST Young Teachers Teaching Talent Training Plan,China
文摘Anchoring mechanism and failure characteristics of composite soft rock with weak interface usually exhibit remarkable difference from those in single rock mass.In order to fully understand the reinforcement mechanism of composite soft roof in western mining area of China,a mechanical model of composite soft rock with weak interface and rock bolt which considering the transverse shear sliding between different rock layers was established firstly.The anchoring effect was quantified by a factor defined as anchoring effect coefficient and its evolution equation was further deduced based on the deformation relationship and homogenized distribution assumption of stress acting on composite structure.Meanwhile,the numerical simulation model of composite soft rock with shear joint was prompted by finite element method.Then detailed analysis were carried out for the deformation features,stress distribution and failure behavior of rock mass and rock bolt near the joint under transverse load.The theoretical result indicates that the anchoring effect of rock-bolt through weak joint changes with the working status of rock mass and closely relates with the physical and geometric parameters of rock mass and rock bolt.From the numerical results,the bending deformation of rock bolt accurately characterized by Doseresp model is mainly concentrated between two plastic hinges near the shear joint.The maximum tensile and compression stresses distribute in the plastic hinge.However,the maximum shear stress appears at the positions of joint surface.The failure zones of composite rock are produced firstly at the joint surface due to the reaction of rock bolt.The above results laid a theoretical and computational foundation for further study of anchorage failure in composite soft rock.
基金the Croatian Science Foundation(HRZZ),Croatia,under the projects IP-2019-04-1618 and I-2243-2017.
文摘Ammonium nitrate and fuel oil(ANFO)based explosive is a classic example of non-ideal high explosives.Its detonation is characterized by a strong dependence of detonation parameters on explosive charge diameter,presence and characteristics of confinement,as well as incomplete consumption of explosive at the sonic point.In this work we propose a detonation model based on the Wood-Kirkwood(WK)theory coupled with the thermochemical code EXPLO5 and supplemented with reaction rate models.Our objective is to analyze the validity of the model for highly non-ideal ANFO explosives,with emphasis on effect of reaction rate models.It was found that both single-step and two-step pressure-based models can be calibrated to reproduce experimental detonation velocity-charge radius data of ANFO at radii significantly above the failure radius(i.e.for D/D_(id)>~0.6).Single-step pressure-based model,with the pressure exponent equal to 1.4,proved to be the most accurate,even in the vicinity of the failure radius.The impact of the rate models is most evident on temporal(and spatial)distribution of flow parameters in detonation driving zone,especially when it comes to the conversion and width of detonation driving zone.
文摘A 2-D numerical model was developed to predict the shape of weld pool in stationary GTA welding of commercial pure aluminium, without considering fluid flow in the weld pool. A Gaussian current density and heat input distribution on the surface of the workpiece were considered. The parameters of Gaussian distribution were modified by comparing calculated results with experimental ones. It was found that these distribution parameters are fimctions of applied current and arc length. Effects of arc length, applied current and welding time on the geometry of the weld pool were investigated. To check the validity of the model, a series of experiments were also conducted. In general, the agreement between calculated overall shape of the weld pool and the experimental one was acceptable, especially in low applied currents. Therefore, it can be concluded that in pure aluminium, the heat conduction is dominant mechanism of heat transfer in the weld pool.
基金Project(2010THZ021)supported by Tsinghua University,ChinaProject(50978139)supported by the National Natural Science Foundation of ChinaProject(2012CB719804)supported by the National Basic Research Program of China
文摘The rapid development of high-speed transportation infrastructure such as highway and high-speed railway has resulted in the advancement of soft soil improvement techniques. Vacuum preloading combined with vertical drains has been proved to be an effective method in the treatment of soft foundation. A three-dimensional numerical analysis of the coupled methods was presented, in which the smear zone and the well resistance were taken into account. The variations of the basic soil parameters including the permeability coefficient and the coefficient of volume compressibility were considered in the numerical model. The result of the numerical model was then compared to the measured value. The results indicate that the decrease of coefficient of volume compressibility accelerates the consolidation of the soil while the influence of hydraulic conductivity is insignificant. A cube drain presents the closest result to the real situation compared to the other equivalent methods of prefabricated vertical drain (PVD). The case study indicates that the numerical model with variation of soil parameters is closer to the measured value than the numerical model without variation of soil parameters.
基金Project(51174032)supported by the National Natural Science Foundation of ChinaProject(NCET-10-0225)supported by the Program for New Century Excellent Talents in University,ChinaProject(FRF-TP-09-001A)supported by the Fundamental Research Funds for the Central Universities,China
文摘Rapid dewatering and thickening of whole-tailings with ultrafine particles is one of the most important processes for the whole-tailings paste preparation. Deep-cone thickener, a kind of such process for the flocculation and settling of whole-tailings, is particularly necessary to study. However, there exist many problems in observing the flocculation and settling process of whole-tailings, as well as the particle size distribution(PSD) of whole-tailings floccules in deep-cone thickener. Population balance model(PBM) is applied to predict the PSD in deep-cone thickener, and LUO model and GHADIRI model are employed to study the aggregation and fragmentation mechanism of the whole-tailings particles, respectively. Through three-dimensional numerical simulation on the whole-tailings flocculation and settling in deep-cone thickener using computational fluid dynamics(CFD)-PBM, the distribution of density and turbulent kinetic energy in deep-cone thickener were obtained, at the same time the spatio-temporal changes of whole-tailings floccules particle size distribution are analyzed. Finally, the major flocculation position in deep-cone thickener is found and the flocculation settling rules of whole-tailings are achieved.
基金Projects(51378131,51378403)supported by the National Natural Science Foundation of ChinaProject(2012210020203)supported by the Fundamental Research Funds for the Central Universities,China
文摘From the mathematical principles, the generalized potential theory can be employed to create constitutive model of geomaterial directly. The similar Cam-clay model, which is created based on the generalized potential theory, has less assumptions,clearer mathematical basis, and better computational accuracy. Theoretically, it is more scientific than the traditional Cam-clay models. The particle flow code PFC3 D was used to make numerical tests to verify the rationality and practicality of the similar Cam-clay model. The verification process was as follows: 1) creating the soil sample for numerical test in PFC3 D, and then simulating the conventional triaxial compression test, isotropic compression test, and isotropic unloading test by PFC3D; 2)determining the parameters of the similar Cam-clay model from the results of above tests; 3) predicting the sample's behavior in triaxial tests under different stress paths by the similar Cam-clay model, and comparing the predicting results with predictions by the Cam-clay model and the modified Cam-clay model. The analysis results show that the similar Cam-clay model has relatively high prediction accuracy, as well as good practical value.
文摘Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the compaction processes are to obtain a compact with the geometrical requirements, without cracks, and with a uniform distribution of density. Design of such proc esses consist, essentially, in determine the sequence and relative displacements of die and punches in order to achieve such goals. A.B. Khoei presented a gener al framework for the finite element simulation of powder forming processes based on the following aspects; a large displacement formulation, centred on a total and updated Lagrangian formulation; an adaptive finite element strategy based on error estimates and automatic remeshing techniques; a cap model based on a hard ening rule in modelling of the highly non-linear behaviour of material; and the use of an efficient contact algorithm in the context of an interface element fo rmulation. In these references, the non-linear behaviour of powder was adequately desc ribed by the cap plasticity model. However, it suffers from a serious deficiency when the stress-point reaches a yield surface. In the flow theory of plasticit y, the transition from an elastic state to an elasto-plastic state appears more or less abruptly. For powder material it is very difficult to define the locati on of yield surface, because there is no distinct transition from elastic to ela stic-plastic behaviour. Results of experimental test on some hard met al powder show that the plastic effects were begun immediately upon loading. In such mater ials the domain of the yield surface would collapse to a point, so making the di rection of plastic increment indeterminate, because all directions are normal to a point. Thus, the classical plasticity theory cannot deal with such materials and an advanced constitutive theory is necessary. In the present paper, the constitutive equations of powder materials will be discussed via an endochronic theory of plasticity. This theory provides a unifi ed point of view to describe the elastic-plastic behaviour of material since it places no requirement for a yield surface and a ’loading function’ to disting uish between loading an unloading. Endochronic theory of plasticity has been app lied to a number of metallic materials, concrete and sand, but to the knowledge of authors, no numerical scheme of the model has been applied to powder material . In the present paper, a new approach is developed based on an endochronic rate independent, density-dependent plasticity model for describing the isothermal deformation behavior of metal powder at low homologous temperature. Although the concept of yield surface has not been explicitly assumed in endochronic theory, it is shown that the cone-cap plasticity yield surface (Fig.1), which is the m ost commonly used plasticity models for describing the behavior of powder materi al can be easily derived as a special case of the proposed endochronic theory. Fig.1 Trace of cone-cap yield function on the meridian pl ane for different relative density As large deformation is observed in powder compaction process, a hypoelastic-pl astic formulation is developed in the context of finite deformation plasticity. Constitutive equations are stated in unrotated frame of reference that greatly s implifies endochronic constitutive relation in finite plasticity. Constitutive e quations of the endochronic theory and their numerical integration are establish ed and procedures for determining material parameters of the model are demonstra ted. Finally, the numerical schemes are examined for efficiency in the model ling of a tip shaped component, as shown in Fig.2. Fig.2 A shaped tip component. a) Geometry, boundary conditio n and finite element mesh; b) density distribution at final stage of
基金the support provided by the Royal Higher Institute for Defence (RHID) of the Belgian Defence, which has contributed to the progress of this ongoing research.
文摘A deep understanding of the internal ballistic process and the factors affecting it is of primary importance to efficiently design a gun system and ensure its safe management. One of the main goals of internal ballistics is to estimate the gas pressure into the combustion chamber and the projectile muzzle velocity in order to use the propellant to its higher efficiency while avoiding over-pressure phenomena. Dealing with the internal ballistic problem is a complex undertaking since it requires handling the interaction between different constituents during a transient time lapse with very steep rise of pressure and temperature. Several approaches have been proposed in the literature, based on different assumptions and techniques. Generally, depending on the used mathematical framework, they can be classified into two categories: computational fluid dynamics-based models and lumped-parameter ones. By focusing on gun systems, this paper offers a review of the main contributions in the field by mentioning their advantages and drawbacks. An insight into the limitations of the currently available modelling strategies is provided,as well as some considerations on the choice of one model over another. Lumped-parameter models, for example, are a good candidate for performing parametric analysis and optimisation processes of gun systems, given their minimum requirements of computer resources. Conversely, CFD-based models have a better capacity to address more sophisticated phenomena like pressure waves and turbulent flow effects. The performed review also reveals that too little attention has been given to small calibre guns since the majority of currently available models are conceived for medium and large calibre gun systems.Similarly, aspects like wear phenomena, bore deformations or projectile-barrel interactions still need to be adequately addressed and our suggestion is to dedicate more effort on it.
基金Projects(52022113,52278546)supported by the National Natural Science Foundation of ChinaProject(2020EEEVL0403)supported by the China Earthquake Administration。
文摘Sudden earthquakes pose a threat to the running safety of trains on high-speed railway bridges,and the stiffness of piers is one of the factors affecting the dynamic response of train-track-bridge system.In this paper,a experiment of a train running on a high-speed railway bridge is performed based on a dynamic experiment system,and the corresponding numerical model is established.The reliability of the numerical model is verified by experiments.Then,the experiment and numerical data are analyzed to reveal the pier height effects on the running safety of trains on bridges.The results show that when the pier height changes,the frequency of the bridge below the 30 m pier height changes greater;the increase of pier height causes the transverse fundamental frequency of the bridge close to that of the train,and the shaking angle and lateral displacement of the train are the largest for bridge with 50 m pier,which increases the risk of derailment;with the pier height increases from 8 m to 50 m,the derailment coefficient obtained by numerical simulations increases by 75% on average,and the spectral intensity obtained by experiments increases by 120% on average,two indicators exhibit logarithmic variation.
基金Project(41202186) supported by the National Natural Science Foundation of ChinaProject(LQ12E08007) supported by the Zhejiang Natural Science Foundation,ChinaProject(#11-KF-08) supported by the Partially Guangxi Key Laboratory of Geomechanics and Geotechnical Engineering,Guilin University of Technology,China
文摘A three-dimensional numerical torsion shear test is presented on hollow cylinder specimen which is performed on a spherical assemblage with fixed principal stress axes using the discrete element code PFC3D.Stack wall technique boundary conditions are employed and optimized to reasonably capture the microstructure evolution.Parametric studies are conducted in terms of the ratio κ,normal and shear stiffness of particles,wall stiffness and friction coefficients.Afterwards,in comparison with physical test,numerical results for a fixed principal stress angle(α=45°) are presented.The results show that the numerical test could capture the macro-micro mechanical behavior of the spherical particle assembly.The evolution of the coordination number demonstrates that particles in shear banding undergo remarkable decrease.The effects of localization on specimens illustrate that global stress and strain recorded from a hollow cylinder apparatus could not represent the localized response.The shearing band initiation and evolution from porosity and shear rate are visualized by contour lines in different shear strains.
文摘Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive strength (UCS) and cohesion) and acting forces (e.g., field stresses and mud pressure). Accurate determination of these parameters is time-consuming, expensive and sometimes even impossible. This work offers a systematic sensitivity analysis to quantify the amount of each parameter’s effect on the stability of a wellbore. Maximum wellbore wall displacement is used as a stability factor to study the stability of a wellbore. A 3D finite difference method with Mohr model is used for the numerical modeling. The numerical model is verified against an analytical solution. A dimensionless sensitivity factor is developed in order to compare the results of various parameters in the sensitivity analysis. The results show a different order of importance of parameters based on rock strength. The most sensitive properties for a weak rock are the maximum horizontal stress, internal friction angle and formation pressure, respectively, while for a strong rock, the most sensitive parameters are the maximum horizontal stress, mud pressure and pore pressure, respectively. The amount of error in wellbore stability analysis inflicted by the error in estimation of each parameter was also derived.
