In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least s...In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.展开更多
Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
RNAs have important biological functions and the functions of RNAs are generally coupled to their structures, especiallytheir secondary structures. In this work, we have made a comprehensive evaluation of the performa...RNAs have important biological functions and the functions of RNAs are generally coupled to their structures, especiallytheir secondary structures. In this work, we have made a comprehensive evaluation of the performances of existingtop RNA secondary structure prediction methods, including five deep-learning (DL) based methods and five minimum freeenergy (MFE) based methods. First, we made a brief overview of these RNA secondary structure prediction methods.Afterwards, we built two rigorous test datasets consisting of RNAs with non-redundant sequences and comprehensivelyexamined the performances of the RNA secondary structure prediction methods through classifying the RNAs into differentlength ranges and different types. Our examination shows that the DL-based methods generally perform better thanthe MFE-based methods for RNAs with long lengths and complex structures, while the MFE-based methods can achievegood performance for small RNAs and some specialized MFE-based methods can achieve good prediction accuracy forpseudoknots. Finally, we provided some insights and perspectives in modeling RNA secondary structures.展开更多
The optimization of the waverider is constrained by the reversely designed leading edge and the constant shock angle distribution. This paper proposes a design method called the variable Leading-Edge Cone (vLEC) metho...The optimization of the waverider is constrained by the reversely designed leading edge and the constant shock angle distribution. This paper proposes a design method called the variable Leading-Edge Cone (vLEC) method to address these limitations. In the vLEC method, the waverider is directly designed from the preassigned leading edge and the variable shock angle distribution based on the Leading-Edge Cone (LEC) concept. Since the vLEC method is an approximate method, two test waveriders are designed and evaluated using numerical simulations to validate the shock design accuracy and the effectiveness of the vLEC method. The results show that the shocks of the test waveriders coincide well with the preassigned positions. Furthermore, four specifically designed application cases are conducted to analyze the performance benefits of the vLEC waveriders. The results of these cases indicate that, due to their variable shock angle distributions, the vLEC waveriders exhibit higher lift-to-drag ratios and better longitudinal static stability than conventional waveriders. Additionally, the vLEC waveriders demonstrate superior volumetric capacities near the symmetry plane, albeit with a minor decrease in volumetric efficiency.展开更多
The visual noise of each light intensity area is different when the image is drawn by Monte Carlo method.However,the existing denoising algorithms have limited denoising performance under complex lighting conditions a...The visual noise of each light intensity area is different when the image is drawn by Monte Carlo method.However,the existing denoising algorithms have limited denoising performance under complex lighting conditions and are easy to lose detailed information.So we propose a rendered image denoising method with filtering guided by lighting information.First,we design an image segmentation algorithm based on lighting information to segment the image into different illumination areas.Then,we establish the parameter prediction model guided by lighting information for filtering(PGLF)to predict the filtering parameters of different illumination areas.For different illumination areas,we use these filtering parameters to construct area filters,and the filters are guided by the lighting information to perform sub-area filtering.Finally,the filtering results are fused with auxiliary features to output denoised images for improving the overall denoising effect of the image.Under the physically based rendering tool(PBRT)scene and Tungsten dataset,the experimental results show that compared with other guided filtering denoising methods,our method improves the peak signal-to-noise ratio(PSNR)metrics by 4.2164 dB on average and the structural similarity index(SSIM)metrics by 7.8%on average.This shows that our method can better reduce the noise in complex lighting scenesand improvethe imagequality.展开更多
The separation-of-variable(SOV)methods,such as the improved SOV method,the variational SOV method,and the extended SOV method,have been proposed by the present authors and coworkers to obtain the closed-form analytica...The separation-of-variable(SOV)methods,such as the improved SOV method,the variational SOV method,and the extended SOV method,have been proposed by the present authors and coworkers to obtain the closed-form analytical solutions for free vibration and eigenbuckling of rectangular plates and circular cylindrical shells.By taking the free vibration of rectangular thin plates as an example,this work presents the theoretical framework of the SOV methods in an instructive way,and the bisection–based solution procedures for a group of nonlinear eigenvalue equations.Besides,the explicit equations of nodal lines of the SOV methods are presented,and the relations of nodal line patterns and frequency orders are investigated.It is concluded that the highly accurate SOV methods have the same accuracy for all frequencies,the mode shapes about repeated frequencies can also be precisely captured,and the SOV methods do not have the problem of missing roots as well.展开更多
To ensure the dimensional accuracy of the final blade profile,it is necessary for precision Electrochemical Machining(ECM)of blade profile to come into an equilibrium state.However,after Electrochemical Trepanning(ECT...To ensure the dimensional accuracy of the final blade profile,it is necessary for precision Electrochemical Machining(ECM)of blade profile to come into an equilibrium state.However,after Electrochemical Trepanning(ECTr),the cascade channel of the blisk is narrow,and the blank allowance distribution is uneven,making it difficult for the precision ECM to become balanced.In blisk production,the two-step method cannot make precision ECM enter equilibrium for some blisk types.A three-step processing method is proposed to overcome this problem.The threestep method adds Electrochemical Homogenizing Machining(ECHM)between the ECTr and precision ECM steps so that the blank allowance can be homogenized quickly without unduly affecting the minimum allowance.Comparative machining experiments of the two-and three-step methods were performed to verify the improvement to blade machining accuracy.The processing results show that the contour parameters of the blade after the three-step method implementation are much better.The allowance difference of the concave(convex)side decreased by 70.5%(65%).In addition,the current in the three-step method is stable at 110 A at the end of precision ECM,verifying successfully entering the equilibrium state.展开更多
Purpose–For the commonly used concrete mix for railway tunnel linings,concrete model specimens were made,and springback and core drilling tests were conducted at different ages.The springback strength was measured to...Purpose–For the commonly used concrete mix for railway tunnel linings,concrete model specimens were made,and springback and core drilling tests were conducted at different ages.The springback strength was measured to the compressive strength of the core sample with a diameter of 100mm and a height-to-diameter ratio of 1:1.By comparing the measured strength values,the relationship between the measured values under different strength measurement methods was analyzed.Design/methodology/approach–A comparative test of the core drilling method and the rebound method was conducted on the side walls of tunnel linings in some under-construction railways to study the feasibility of the rebound method in engineering quality supervision and inspection.Findings–Tests showed that the rebound strength was positively correlated with the core drill strength.The core drill test strength was significantly higher than the rebound test strength,and the strength still increased after 56 days of age.The rebound method is suitable for the general survey of concrete strength during the construction process and is not suitable for direct supervision and inspection.Originality/value–By studying the correlation of test strength of tunnel lining concrete using two methods,the differences in test results of different methods are proposed to provide a reference for the test and evaluation of tunnel lining strength in railway engineering.展开更多
Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using in...The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.展开更多
Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and ...Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.展开更多
Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of un...Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.展开更多
The complex vibration directly affects the dynamic safety of drill string in ultra-deep wells and extra-deep wells.It is important to understand the dynamic characteristics of drill string to ensure the safety of dril...The complex vibration directly affects the dynamic safety of drill string in ultra-deep wells and extra-deep wells.It is important to understand the dynamic characteristics of drill string to ensure the safety of drill string.Due to the super slenderness ratio of drill string,strong nonlinearity implied in dynamic analysis and the complex load environment,dynamic simulation of drill string faces great challenges.At present,many simulation methods have been developed to analyze drill string dynamics,and node iteration method is one of them.The node iteration method has a unique advantage in dealing with the contact characteristics between drill string and borehole wall,but its drawback is that the calculation consumes a considerable amount of time.This paper presents a dynamic simulation method of drilling string in extra-deep well based on successive over-relaxation node iterative method(SOR node iteration method).Through theoretical analysis and numerical examples,the correctness and validity of this method were verified,and the dynamics characteristics of drill string in extra-deep wells were calculated and analyzed.The results demonstrate that,in contrast to the conventional node iteration method,the SOR node iteration method can increase the computational efficiency by 48.2%while achieving comparable results.And the whirl trajectory of the extra-deep well drill string is extremely complicated,the maximum rotational speed downhole is approximately twice the rotational speed on the ground.The dynamic torque increases rapidly at the position of the bottom stabilizer,and the lateral vibration in the middle and lower parts of drill string is relatively intense.展开更多
In this article,we develop the Laplace transform(LT)based Chebyshev spectral collocation method(CSCM)to approximate the time fractional advection-diffusion equation,incorporating the Atangana-Baleanu Caputo(ABC)deriva...In this article,we develop the Laplace transform(LT)based Chebyshev spectral collocation method(CSCM)to approximate the time fractional advection-diffusion equation,incorporating the Atangana-Baleanu Caputo(ABC)derivative.The advection-diffusion equation,which governs the transport of mass,heat,or energy through combined advection and diffusion processes,is central to modeling physical systems with nonlocal behavior.Our numerical scheme employs the LT to transform the time-dependent time-fractional PDEs into a time-independent PDE in LT domain,eliminating the need for classical time-stepping methods that often suffer from stability constraints.For spatial discretization,we employ the CSCM,where the solution is approximated using Lagrange interpolation polynomial based on the Chebyshev collocation nodes,achieving exponential convergence that outperforms the algebraic convergence rates of finite difference and finite element methods.Finally,the solution is reverted to the time domain using contour integration technique.We also establish the existence and uniqueness of the solution for the proposed problem.The performance,efficiency,and accuracy of the proposed method are validated through various fractional advection-diffusion problems.The computed results demonstrate that the proposed method has less computational cost and is highly accurate.展开更多
Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an...Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.展开更多
Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of thes...Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.展开更多
Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extens...Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extensive attention due to their advantages of small size,high multiplexing efficiency,convenient mass production,and low cost.An RGB beam combiner based on directional couplers is designed,with a core-cladding relative refractive index difference of 0.75%.The RGB beam combiner is optimized from the perspective of parameter optimization.Using the beam propagation method(BPM),the relationship between the performance of the RGB beam combiner and individual parameters is studied,achieving preliminary optimization of the device’s performance.The key parameters of the RGB beam combiner are optimized using the entropy weight-technique for order preference by similarity to an ideal solution TOPSIS method,establishing the optimal parameter scheme and further improving the device’s performance indicators.The results show that after optimization,the multiplexing efficiencies for red,green,and blue lights,as well as the average multiplexing efficiency,reached 99.17%,99.76%,96.63%and 98.52%,respectively.The size of the RGB beam combiner is 4.768 mm×0.062 mm.展开更多
In the aerospace field, residual stress directly affects the strength, fatigue life and dimensional stability of thin-walled structural components, and is a key factor to ensure flight safety and reliability. At prese...In the aerospace field, residual stress directly affects the strength, fatigue life and dimensional stability of thin-walled structural components, and is a key factor to ensure flight safety and reliability. At present, research on residual stress at home and abroad mainly focuses on the optimization of traditional detection technology, stress control of manufacturing process and service performance evaluation, among which research on residual stress detection methods mainly focuses on the improvement of the accuracy, sensitivity, reliability and other performance of existing detection methods, but it still faces many challenges such as extremely small detection range, low efficiency, large error and limited application range.展开更多
NLLoc is a nonlinear search positioning method.In this study,we use simulated arrival time data to quantitatively evaluate the NLLoc method from three aspects:arrival time picking accuracy,station distribution,and vel...NLLoc is a nonlinear search positioning method.In this study,we use simulated arrival time data to quantitatively evaluate the NLLoc method from three aspects:arrival time picking accuracy,station distribution,and velocity model.The results show that the NLLoc method exhibits high positioning accuracy and stability in terms of arrival time picking accuracy and station distribution;however,it is sensitive to the velocity model.The positioning accuracy is higher when the velocity model is smaller than the true velocity.We combined absolute and relative positioning methods.First,we use the NLLoc method for absolute positioning of seismic data and then the double difference positioning method for relative positioning to obtain a more accurate relocation result.Furthermore,we used the combined method to locate the earthquake sequence after collecting dense seismic array data on the Luanzhou M_(S)4.3 earthquake that occurred on April 16,2021,in Hebei Province.By fitting the fault plane with the relocated earthquake sequences,the results show that the strike and dip angles of the seismogenic fault of the Luanzhou M_(S)4.3 earthquake are 208.5°and 85.6°,respectively.This indicates a high-dip angle fault with North-North-East strike and North-West dip directions.Furthermore,we infer that the seismogenic fault of the Luanzhou M_(S)4.3 earthquake is the Lulong fault.展开更多
基金supported by the Graduate Student Scientific Research Innovation Project through Research Innovation Fund for Graduate Students in Shanxi Province(Project No.2024KY648).
文摘In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
基金supported by grants from the National Science Foundation of China(Grant Nos.12375038 and 12075171 to ZJT,and 12205223 to YLT).
文摘RNAs have important biological functions and the functions of RNAs are generally coupled to their structures, especiallytheir secondary structures. In this work, we have made a comprehensive evaluation of the performances of existingtop RNA secondary structure prediction methods, including five deep-learning (DL) based methods and five minimum freeenergy (MFE) based methods. First, we made a brief overview of these RNA secondary structure prediction methods.Afterwards, we built two rigorous test datasets consisting of RNAs with non-redundant sequences and comprehensivelyexamined the performances of the RNA secondary structure prediction methods through classifying the RNAs into differentlength ranges and different types. Our examination shows that the DL-based methods generally perform better thanthe MFE-based methods for RNAs with long lengths and complex structures, while the MFE-based methods can achievegood performance for small RNAs and some specialized MFE-based methods can achieve good prediction accuracy forpseudoknots. Finally, we provided some insights and perspectives in modeling RNA secondary structures.
基金supported by grants from the National Natural Science Foundation of China(No.U20B2006)the Guangdong Basic and Applied Basic Research Foundation(No.2022A1515110145)Young Elite Scientists Sponsorship Program by China Association for Science and Technology(No.2022QNRC001).
文摘The optimization of the waverider is constrained by the reversely designed leading edge and the constant shock angle distribution. This paper proposes a design method called the variable Leading-Edge Cone (vLEC) method to address these limitations. In the vLEC method, the waverider is directly designed from the preassigned leading edge and the variable shock angle distribution based on the Leading-Edge Cone (LEC) concept. Since the vLEC method is an approximate method, two test waveriders are designed and evaluated using numerical simulations to validate the shock design accuracy and the effectiveness of the vLEC method. The results show that the shocks of the test waveriders coincide well with the preassigned positions. Furthermore, four specifically designed application cases are conducted to analyze the performance benefits of the vLEC waveriders. The results of these cases indicate that, due to their variable shock angle distributions, the vLEC waveriders exhibit higher lift-to-drag ratios and better longitudinal static stability than conventional waveriders. Additionally, the vLEC waveriders demonstrate superior volumetric capacities near the symmetry plane, albeit with a minor decrease in volumetric efficiency.
基金supported by the National Natural Science(No.U19A2063)the Jilin Provincial Development Program of Science and Technology (No.20230201080GX)the Jilin Province Education Department Scientific Research Project (No.JJKH20230851KJ)。
文摘The visual noise of each light intensity area is different when the image is drawn by Monte Carlo method.However,the existing denoising algorithms have limited denoising performance under complex lighting conditions and are easy to lose detailed information.So we propose a rendered image denoising method with filtering guided by lighting information.First,we design an image segmentation algorithm based on lighting information to segment the image into different illumination areas.Then,we establish the parameter prediction model guided by lighting information for filtering(PGLF)to predict the filtering parameters of different illumination areas.For different illumination areas,we use these filtering parameters to construct area filters,and the filters are guided by the lighting information to perform sub-area filtering.Finally,the filtering results are fused with auxiliary features to output denoised images for improving the overall denoising effect of the image.Under the physically based rendering tool(PBRT)scene and Tungsten dataset,the experimental results show that compared with other guided filtering denoising methods,our method improves the peak signal-to-noise ratio(PSNR)metrics by 4.2164 dB on average and the structural similarity index(SSIM)metrics by 7.8%on average.This shows that our method can better reduce the noise in complex lighting scenesand improvethe imagequality.
基金supported by the National Natural Science Foundation of China(12172023).
文摘The separation-of-variable(SOV)methods,such as the improved SOV method,the variational SOV method,and the extended SOV method,have been proposed by the present authors and coworkers to obtain the closed-form analytical solutions for free vibration and eigenbuckling of rectangular plates and circular cylindrical shells.By taking the free vibration of rectangular thin plates as an example,this work presents the theoretical framework of the SOV methods in an instructive way,and the bisection–based solution procedures for a group of nonlinear eigenvalue equations.Besides,the explicit equations of nodal lines of the SOV methods are presented,and the relations of nodal line patterns and frequency orders are investigated.It is concluded that the highly accurate SOV methods have the same accuracy for all frequencies,the mode shapes about repeated frequencies can also be precisely captured,and the SOV methods do not have the problem of missing roots as well.
基金supported by the National Natural Science Foundation of China(No.52075253)the Innovation Research Team of the National Natural Science Foundation of China(No.51921003)。
文摘To ensure the dimensional accuracy of the final blade profile,it is necessary for precision Electrochemical Machining(ECM)of blade profile to come into an equilibrium state.However,after Electrochemical Trepanning(ECTr),the cascade channel of the blisk is narrow,and the blank allowance distribution is uneven,making it difficult for the precision ECM to become balanced.In blisk production,the two-step method cannot make precision ECM enter equilibrium for some blisk types.A three-step processing method is proposed to overcome this problem.The threestep method adds Electrochemical Homogenizing Machining(ECHM)between the ECTr and precision ECM steps so that the blank allowance can be homogenized quickly without unduly affecting the minimum allowance.Comparative machining experiments of the two-and three-step methods were performed to verify the improvement to blade machining accuracy.The processing results show that the contour parameters of the blade after the three-step method implementation are much better.The allowance difference of the concave(convex)side decreased by 70.5%(65%).In addition,the current in the three-step method is stable at 110 A at the end of precision ECM,verifying successfully entering the equilibrium state.
文摘Purpose–For the commonly used concrete mix for railway tunnel linings,concrete model specimens were made,and springback and core drilling tests were conducted at different ages.The springback strength was measured to the compressive strength of the core sample with a diameter of 100mm and a height-to-diameter ratio of 1:1.By comparing the measured strength values,the relationship between the measured values under different strength measurement methods was analyzed.Design/methodology/approach–A comparative test of the core drilling method and the rebound method was conducted on the side walls of tunnel linings in some under-construction railways to study the feasibility of the rebound method in engineering quality supervision and inspection.Findings–Tests showed that the rebound strength was positively correlated with the core drill strength.The core drill test strength was significantly higher than the rebound test strength,and the strength still increased after 56 days of age.The rebound method is suitable for the general survey of concrete strength during the construction process and is not suitable for direct supervision and inspection.Originality/value–By studying the correlation of test strength of tunnel lining concrete using two methods,the differences in test results of different methods are proposed to provide a reference for the test and evaluation of tunnel lining strength in railway engineering.
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金funded by Vice Chancellor of Research at Shiraz University(grant 3GFU2M1820).
文摘The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.
基金supported by the National Natural Science Foundation of China(No.42372287 and No.U24A20178)the Fundamental Research Funds for the Central Universities CHD(No.2024SHEEAR002)+3 种基金the Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shaanxi Province(No.2020024)the China Postdoctoral Science Foundation(GZC20232955,2024M753472,and 2024MD763937)the Science-Technology Foundation for Young Scientists of Gansu Province,China(No.24JRRA097)the Study of biodiversity survey and limiting factor analysis of Yinkentala(2023ZL01).
文摘Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.
基金supported by the National Natural Science Foundation of China(No.92252201)the Fundamental Research Funds for the Central Universitiesthe Academic Excellence Foundation of Beihang University(BUAA)for PhD Students。
文摘Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.
基金supported by the National Natural Science Foundation of China(52174003,52374008).
文摘The complex vibration directly affects the dynamic safety of drill string in ultra-deep wells and extra-deep wells.It is important to understand the dynamic characteristics of drill string to ensure the safety of drill string.Due to the super slenderness ratio of drill string,strong nonlinearity implied in dynamic analysis and the complex load environment,dynamic simulation of drill string faces great challenges.At present,many simulation methods have been developed to analyze drill string dynamics,and node iteration method is one of them.The node iteration method has a unique advantage in dealing with the contact characteristics between drill string and borehole wall,but its drawback is that the calculation consumes a considerable amount of time.This paper presents a dynamic simulation method of drilling string in extra-deep well based on successive over-relaxation node iterative method(SOR node iteration method).Through theoretical analysis and numerical examples,the correctness and validity of this method were verified,and the dynamics characteristics of drill string in extra-deep wells were calculated and analyzed.The results demonstrate that,in contrast to the conventional node iteration method,the SOR node iteration method can increase the computational efficiency by 48.2%while achieving comparable results.And the whirl trajectory of the extra-deep well drill string is extremely complicated,the maximum rotational speed downhole is approximately twice the rotational speed on the ground.The dynamic torque increases rapidly at the position of the bottom stabilizer,and the lateral vibration in the middle and lower parts of drill string is relatively intense.
基金extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through Large Research Project under grant number RGP2/174/46.
文摘In this article,we develop the Laplace transform(LT)based Chebyshev spectral collocation method(CSCM)to approximate the time fractional advection-diffusion equation,incorporating the Atangana-Baleanu Caputo(ABC)derivative.The advection-diffusion equation,which governs the transport of mass,heat,or energy through combined advection and diffusion processes,is central to modeling physical systems with nonlocal behavior.Our numerical scheme employs the LT to transform the time-dependent time-fractional PDEs into a time-independent PDE in LT domain,eliminating the need for classical time-stepping methods that often suffer from stability constraints.For spatial discretization,we employ the CSCM,where the solution is approximated using Lagrange interpolation polynomial based on the Chebyshev collocation nodes,achieving exponential convergence that outperforms the algebraic convergence rates of finite difference and finite element methods.Finally,the solution is reverted to the time domain using contour integration technique.We also establish the existence and uniqueness of the solution for the proposed problem.The performance,efficiency,and accuracy of the proposed method are validated through various fractional advection-diffusion problems.The computed results demonstrate that the proposed method has less computational cost and is highly accurate.
基金Supported by the National Natural Science Foundation of China (Grant Nos.52088102 and 51879287)National Key Research and Development Program of China (Grant No.2022YFB2602301)。
文摘Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.
基金supported by the Open Funds for Hubei Key Laboratory of Marine Geological Resources,China University of Geosciences(No.MGR202308)the Natural Science Foundation of Shandong Province(No.ZR2020MD085)+3 种基金the National Natural Science Foundation of China(No.41821004)the Taishan Scholar Program(No.tstp2022114)the Shandong Provincial Natural Science Foundation(No.DKXZZ202206)the National Key Research and Development Program of China(No.2016YFC1402404).
文摘Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.
基金Project(52175445)supported by the National Natural Science Foundation of ChinaProject(2022JJ30743)supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(2023GK2024)supported by the Key Research and Development Program of Hunan Province,ChinaProject(2023ZZTS0391)supported by the Fundamental Research Funds for the Central Universities of China。
文摘Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extensive attention due to their advantages of small size,high multiplexing efficiency,convenient mass production,and low cost.An RGB beam combiner based on directional couplers is designed,with a core-cladding relative refractive index difference of 0.75%.The RGB beam combiner is optimized from the perspective of parameter optimization.Using the beam propagation method(BPM),the relationship between the performance of the RGB beam combiner and individual parameters is studied,achieving preliminary optimization of the device’s performance.The key parameters of the RGB beam combiner are optimized using the entropy weight-technique for order preference by similarity to an ideal solution TOPSIS method,establishing the optimal parameter scheme and further improving the device’s performance indicators.The results show that after optimization,the multiplexing efficiencies for red,green,and blue lights,as well as the average multiplexing efficiency,reached 99.17%,99.76%,96.63%and 98.52%,respectively.The size of the RGB beam combiner is 4.768 mm×0.062 mm.
文摘In the aerospace field, residual stress directly affects the strength, fatigue life and dimensional stability of thin-walled structural components, and is a key factor to ensure flight safety and reliability. At present, research on residual stress at home and abroad mainly focuses on the optimization of traditional detection technology, stress control of manufacturing process and service performance evaluation, among which research on residual stress detection methods mainly focuses on the improvement of the accuracy, sensitivity, reliability and other performance of existing detection methods, but it still faces many challenges such as extremely small detection range, low efficiency, large error and limited application range.
基金Supported by the Foundation:This research project is jointly supported by Hebei Provincial Science and Technology Program(No.22375406D)The Earthquake Science and Technology Program of Hebei Province(No.DZ2023120500009,DZ2024120500001).
文摘NLLoc is a nonlinear search positioning method.In this study,we use simulated arrival time data to quantitatively evaluate the NLLoc method from three aspects:arrival time picking accuracy,station distribution,and velocity model.The results show that the NLLoc method exhibits high positioning accuracy and stability in terms of arrival time picking accuracy and station distribution;however,it is sensitive to the velocity model.The positioning accuracy is higher when the velocity model is smaller than the true velocity.We combined absolute and relative positioning methods.First,we use the NLLoc method for absolute positioning of seismic data and then the double difference positioning method for relative positioning to obtain a more accurate relocation result.Furthermore,we used the combined method to locate the earthquake sequence after collecting dense seismic array data on the Luanzhou M_(S)4.3 earthquake that occurred on April 16,2021,in Hebei Province.By fitting the fault plane with the relocated earthquake sequences,the results show that the strike and dip angles of the seismogenic fault of the Luanzhou M_(S)4.3 earthquake are 208.5°and 85.6°,respectively.This indicates a high-dip angle fault with North-North-East strike and North-West dip directions.Furthermore,we infer that the seismogenic fault of the Luanzhou M_(S)4.3 earthquake is the Lulong fault.
