Current situation and development trend of design methods for subgrade structure of high speed railways

Purpose–This method will become a new development trend in subgrade structure design for high speed railways.Design/methodology/approach–This paper summarizes the structural types and design methods of subgrade bed for high speed railways in China,Japan,France,Germany,the United States and other countries based on the study and analysis of existing literature and combined with the research results and practices of high speed railway subgrade engineering at home and abroad.Findings–It is found that in foreign countries,the layered reinforced structure is generally adopted for the subgrade bed of high speed railways,and the unified double-layer or multi-layer structure is adopted for the surface layer of subgrade bed,while the simple structure is adopted in China;in foreign countries,different inspection parameters are adopted to evaluate the compaction state of fillers according to their respective understanding and practice,while in China,compaction coefficient,subsoil coefficient and dynamic deformation modulus are adopted for such evaluation;in foreign countries,the subgrade top deformation control method,the subgrade bottom deformation control method,the subsurface fill strength control method are mainly adopted in subgrade bed structure design of high speed railways,while in China,dynamic deformation control of subgrade surface and dynamic strain control of subgrade bed bottom layer is adopted in the design.However,the cumulative deformation of subgrade caused by train cyclic vibration load is not considered in the existing design methods.Originality/value–This paper introduces a new subgrade structure design method based on whole-process dynamics analysis that meets subgrade functional requirements and is established on the basis of the existing research at home and abroad on prediction methods for cumulative deformation of subgrade soil.

Experimental and numerical study on dynamic mechanical behaviors of shale under true triaxial compression at high strain rate

High-energy gas fracturing of shale is a novel,high efficacy and eco-friendly mining technique,which is a typical dynamic perturbing behavior.To effectively extract shale gas,it is important to understand the dynamic mechanical properties of shale.Dynamic experiments on shale subjected to true triaxial compression at different strain rates are first conducted in this research.The dynamic stress-strain curves,peak strain,peak stress and failure modes of shale are investigated.The results of the study indicate that the intermediate principal stress and the minor principal stress have the significant influence on the dynamic mechanical behaviors,although this effect decreases as the strain rate increases.The characteristics of compression-shear failure primarily occur in shale subjected to triaxial compression at high strain rates,which distinguishes it from the fragmentation characteristics observed in shale under dynamic uniaxial compression.Additionally,a numerical three-dimensional Split Hopkinson Pressure Bar(3D-SHPB),which is established by coupling PFC3D and FLAC3D methods,is validated to replicate the laboratory characteristics of shale.The dynamic mechanical characteristics of shale subjected to different confining stresses are systematically investigated by the coupling PFC3D and FLAC3D method.The numerical results are in good agreement with the experimental data.

A Provable Positivity-Preserving Local Discontinuous Galerkin Method for the Viscous and Resistive MHD Equations

In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.

Research Progress of High Entropy Ceramic Materials

High-entropy materials(HEMs)have better mechanical,thermal,and electrical properties than traditional materials due to their special"high entropy effect".They can also adjust the performance of high entropy ceramics by adjusting the proportion of raw materials,and have broad application prospects in many fields.This article provides a review of the high entropy effect,preparation methods,and main applications of high entropy ceramic materials,especially exploring relevant research on high entropy perovskite ceramics.It is expected to provide reference for the promotion of scientific research and the development of further large-scale applications of high-entropy ceramic materials.

Mapping coalmine goaf using transient electromagnetic method and high density resistivity method in Ordos City, China

The research about subsurface characteristics by using transient electromagnetic method（TEM） and high density resistivity method（HDRM） were already conducted in Ordos. The objective of this research is to detect coalmine goaf areas based on rock resistivity. The data processing using wavelet transform, three point smoothing, RES2 DINV and Maxwell processing software to obtain 2D resistivity structure. The results showed that the layers with maximum resistivity values（30e33 U m on Line 1, 30e31 U m on Line 2, 32e40 U m on Line3） are founded at station 1e7, and 14e20 on Line 1,13e18 on Line 2, and 8e13 and 16e20 on Line 3 which is predicted as goaf layer, and the minimum resistivity values（20e26 U m of TEM, 45e75 U m of HDRM） at the other layers. This resistivity difference was caused by the geology and characteristics of the study area which is located close by the cleugh with rich coal, so the goaf area distinguishable with aquifer layer and coal seam. The results were also significant accidents and serious destruction of ecological environment.

Optimization of High-Protein Glutinous Rice Flour Production Using Response Surface Method

A response surface method was employed to study the effect of α-amylase concentration, hydrolysis temperature and time on the production of high protein glutinous rice flour(HPGRF). The suspension of glutinous rice flour(15%) that contained 6.52% protein was gelatinized and subsequently hydrolyzed by thermostable α-amylase. The hydrolysis yielded 0.144–0.222 g/g HPGRF with 29.4%–45.4% protein content. Hydrolysis time exerted a significant effect, while enzyme concentration and hydrolysis temperature showed insignificant effect on the protein content and production yield of HPGRF. The result of response surface method showed that the optimum condition for the production of HPGRF that contained at least 36% protein was treating gelatinized 15% glutinous rice flour suspension with 0.90 Kilo Novo α-amylase Unit(KNU)/g α-amylase at 80 oC for 99 min. By carrying out the predicted hydrolysis condition, HPGRF with 35.9% protein and 61.8% carbohydrates was resulted. The process yielded 0.172 g/g HPGRF. HPGRF contained higher amount of essential amino acids compared to glutinous rice flour. HPGRF had higher solubility and lower swelling power, and also showed no pasting peak compared with glutinous rice flour.

Focused Wave Properties Based on A High Order Spectral Method with A Non-Periodic Boundary

In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.

On the interaction between bubbles and the free surface with high density ratio 3D lattice Boltzmann method

The bubbles rise up and burst at the free surface is a complex two-phase process.A free energy lattice Boltzmann method(LBM)model is adopted in this paper to study this phenomenon.The interface capturing technique[Zheng et al.,2006]is used to deal with the high density ratio problem.The Laplace law and the air-water interface capturing ability are validated for the multiphase model.The interaction between the single bubble or multiple bubbles and the free surface are studied by the multiphase model.The force acting on the bubble and the evolution of the free surface is studied.Meanwhile,effect of the initial distance between two adjacent bubbles on interaction effects of multiple bubbles is investigated as well.

High Efficient Methods of Content-based 3D Model Retrieval

Content-based 3D model retrieval is of great help to facilitate the reuse of existing designs and to inspire designers during conceptual design. However, there is still a gap to apply it in industry due to the low time efficiency. This paper presents two new methods with high efficiency to build a Content-based 3D model retrieval system. First, an improvement is made on the ＂Shape Distribution （D2）＂ algorithm, and a new algorithm named ＂Quick D2＂ is proposed. Four sample 3D mechanical models are used in an experiment to compare the time cost of the two algorithms. The result indicates that the time cost of Quick D2 is much lower than that of D2, while the descriptors extracted by the two algorithms are almost the same. Second, an expandable 3D model repository index method with high performance, namely, RBK index, is presented. On the basis of RBK index, the search space is pruned effectively during the search process, leading to a speed up of the whole system. The factors that influence the values of the key parameters of RBK index are discussed and an experimental method to find the optimal values of the key parameters is given. Finally, ＂3D Searcher＂, a content-based 3D model retrieval system is developed. By using the methods proposed, the time cost for the system to respond one query online is reduced by 75% on average. The system has been implemented in a manufacturing enterprise, and practical query examples during a case of the automobile rear axle design are also shown. The research method presented shows a new research perspective and can effectively improve the content-based 3D model retrieval efficiency.

High-precision solution to the moving load problem using an improved spectral element method

In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem.In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases.Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.

Hermite WENO-based limiters for high order discontinuous Galerkin method on unstructured grids

A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.

Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method

By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.

Re-study on Recurrence Period of Stokes Wave Train with High Order Spectral Method

Owing to the Benjamin-Feir instability, the Stokes wave train experiences a modulation-demodulation process, and presents a recurrence characteristics. Stiassnie and Shemer researched the unstable evolution process and provided a theoretical formulation for the recurrence period in 1985 on the basis of the nonlinear cubic Schrodinger equation （NLS）. However, NLS has limitations on the narrow band and the weak nonlinearity. The recurrence period is re-investigated in this paper by using a highly efficient High Order Spectral （HOS） method, which can be applied for the direct phase- resolved simulation of the nonlinear wave train evolution. It is found that the Stiassnie and Shemer＇s formula should be modified in the cases with most unstable initial conditions, which is important for such topics as the generation mechanisms of freak waves. A new recurrence period formula is presented and some new evolution characteristics of the Stokes wave train are also discussed in details.

Microstructure and Property of Hypereutectic High Chromium Cast Iron Prepared by Slope Cooling Body-Centrifugal Casting Method

In this paper, the ring-type ingot of hypereutectic high Cr cast iron was obtained by slope cooling bodycentrifugal casting method （SC-CCM）, and its microstructure and impact toughness were investigated, respectively. The results indicated that, first, the primary carbides in the microstructure are prominently finer than those in the hypereutectic high Cr cast iron prepared by conventional casting method. Second, in the ring-type ingot, the primary carbides near radial outer field are finer than those near radial inner field; furthermore, there is dividing field in the microstructure. Finally, the impact toughness values of the specimens impacted on the radial outer face and on the radial inner face are improved respectively about 36% and 138% more than that of the hypereutectic high Cr one prepared by conventional casting method.

Self-heating Probe Instrument and Method for Measuring High Temperature Melting Volume Change Rate of Material

The castings defects are affected by the melting volume change rate of material. The change rate has an important effect on running safety of the high temperature thermal storage chamber, too. But the characteristics of existing measuring installations are complex structure, troublesome operation and low precision. In order to measure the melting volume change rate of material accurately and conveniently, a self-designed measuring instrument, self-heating probe instrument, and measuring method are described. Temperature in heating cavity is controlled by PID temperature controller; melting volume change rate υ and molten density are calculated based on the melt volume which is measured by the instrument. Positive and negative υ represent expansion and shrinkage of the sample volume after melting, respectively. Taking eutectic LiF＋CaF2 for example, its melting volume change rate and melting density at 1 123 K are -20.6% and 2 651 kg/m–3 measured by this instrument, which is only 0.71% smaller than literature value. Density and melting volume change rate of industry pure aluminum at 973 K and analysis pure NaCl at 1 123 K are detected by the instrument too. The measure results are agreed with report values. Measuring error sources are analyzed and several improving measures are proposed. In theory, the measuring errors of the change rate and molten density which are measured by the self-designed instrument is nearly 1/20-1/50 of that measured by the refitted mandril thermal expansion instrument. The self-designed instrument and method have the advantages of simple structure, being easy to operate, extensive applicability for material, relatively high accuracy, and most importantly, temperature and sample vapor pressure have little effect on the measurement accuracy. The presented instrument and method solve the problems of complicated structure and procedures, and large measuring errors for the samples with high vapor pressure by existing installations.

High-pressure phonon dispersion of copper by using the modified analytic embedded atom method

By using the Born-von Kfirmfin theory of lattice dynamics and the modified analytic embedded atom method, we reproduce the experimental results of the phonon dispersion in fcc metal Cu at zero pressure along three high symmetry directions and four oft-symmetry directions, and then simulate the phonon dispersion curves of Cu at high pressures of 50, 100, and 150 GPa. The results show that the shapes of dispersion curves at high pressures are very similar to that at zero pressure. All the vibration frequencies of Cu in all vibration branches at high pressures are larger than the results at zero pressure, and increase correspondingly as pressure reaches 50, 100, and 150 GPa sequentially. Moreover, on the basis of phonon dispersion, we calculate the values of specific heat of Cu at different pressures. The prediction of thermodynamic quantities lays a significant foundation for guiding and judging experiments of thermodynamic properties of solids under high pressures.

Synthesis of large diamond crystals containing high-nitrogen concentration at high pressure and high temperature using Ni-based solvent by temperature gradient method

This paper reprots that with Ni-based catalyst/solvent and with a dopant of NAN3, large green single crystal diamonds with perfect shape are successfully synthesized by temperature gradient method under high pressure and high temperature in a China-type cubic anvil high-pressure apparatus （SPD-6 × 1200）, and the highest nitrogen concentration reaches approximately 121-1257 ppm calculated by infrared absorption spectra. The synthesis conditions are about 5.5 CPa and 1240-1300 ℃. The growth behaviour of diamond with high-nitrogen concentration is investigated in detail. The results show that, with increasing the content of NaN3 added in synthesis system, the width of synthesis temperature region for growth high-quality diamonds becomes narrower, and the morphology of diamond crystal is changed from cube-octahedral to octahedral at same temperature and pressure, the crystal growth rate is slowed down, nevertheless, the nitrogen concentration doped in synthetic diamond increases.

HIGH PERFORMANCE SPARSE SOLVER FOR UNSYMMETRICAL LINEAR EQUATIONS WITH OUT-OF-CORE STRATEGIES AND ITS APPLICATION ON MESHLESS METHODS

A new direct method for solving unsymmetrical sparse linear systems（USLS） arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin （MLPG） method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.

THE STUDY ON MEASURING TECHNIQUE OF PARTIAL DISCHARGE IN GAS INSULATED SWITCHGEAR USING ULTRA HIGH FREQUENCY METHOD WITH EXTERNAL SENSORS

Objective In order to find early latent faults and prevent catastrophic failures, diagnosis of insulation condition by measuring technique of partial discharge(PD) in gas insulated switchgear (GIS) is applied in this paper, which is one of the most basic ways for diagnosis of insulation condition. Methods Ultra high frequency(UHF) PD detection method by using internal sensors has been proved efficient, because it may avoid the disturbance of corona, but the sensor installation of this method will be limited by the structure and operation condition of GIS. There are some of electromagnetic (E-M) waves leak from the place of insulation spacer, therefore, the external sensors UHF measuring PD technique is applied, which isn't limited by the operation condition of GIS. Results This paper analyzes propagated electromagnetic (E-M) waves of partial discharge pulse excited by using the finite-difference time-domain (FDTD) method. The signal collected at the outer point is more complex than that of the inner point, and the signals' amplitude of outer is about half of the inner, because it propagates through spacer and insulation slot. Set up UHF PD measuring system. The typical PD in 252kV GIS bus bar was measured using PD detection UHF technique with external sensors. Finally, compare the results of UHF measuring technique using external sensors with the results of FDTD method simulation and the traditional IEC60270 method detection. Conclusion The results of experiment shows that the UHF technique can realize the diagnosis of insulation condition, the results of FDTD method simulation and the result UHF method detection can demonstrate each other, which gives references to further researches and application for UHF PD measuring technique.