A generalized nonlinear three-dimensional Hoek‒Brown failure criterion

To understand the strengths of rocks under complex stress states,a generalized nonlinear threedimensional(3D)Hoek‒Brown failure(NGHB)criterion was proposed in this study.This criterion shares the same parameters with the generalized HB(GHB)criterion and inherits the parameter advantages of GHB.Two new parameters,b,and n,were introduced into the NGHB criterion that primarily controls the deviatoric plane shape of the NGHB criterion under triaxial tension and compression,respectively.The NGHB criterion can consider the influence of intermediate principal stress(IPS),where the deviatoric plane shape satisfies the smoothness requirements,while the HB criterion not.This criterion can degenerate into the two modified 3D HB criteria,the Priest criterion under triaxial compression condition and the HB criterion under triaxial compression and tension condition.This criterion was verified using true triaxial test data for different parameters,six types of rocks,and two kinds of in situ rock masses.For comparison,three existing 3D HB criteria were selected for performance comparison research.The result showed that the NGHB criterion gave better prediction performance than other criteria.The prediction errors of the strength of six types of rocks and two kinds of in situ rock masses were in the range of 2.0724%-3.5091%and 1.0144%-3.2321%,respectively.The proposed criterion lays a preliminary theoretical foundation for prediction of engineering rock mass strength under complex in situ stress conditions.

Estimating shear strength of high-level pillars supported with cemented backfilling using the HoekeBrown strength criterion

Deep metal mines are often mined using the high-level pillars with subsequent cementation backfilling(HLSCB)mining method.At the design stage,it is therefore important to have a reasonable method for determining the shear strength of the high-level pillars(i.e.cohesion and internal friction angle)when they are supported by cemented backfilling.In this study,a formula was derived for the upper limit of the confining pressure σ3max on a high-level pillar supported by cemented backfilling in a deep metal mine.A new method of estimating the shear strength of such pillars was then proposed based on the Hoek eBrown failure criterion.Our analysis indicates that the horizontal stress σhh acting on the cemented backfill pillar can be simplified by expressing it as a constant value.A reasonable and effective value for σ3max can then be determined.The value of s3max predicted using the proposed method is generally less than 3 MPa.Within this range,the shear strength of the high-level pillar is accurately calculated using the equivalent MohreCoulomb theory.The proposed method can effectively avoid the calculation of inaccurate shear strength values for the high-level pillars when the original HoekeBrown criterion is used in the presence of large confining pressures,i.e.the situation in which the cohesion value is too large and the friction angle is too small can effectively be avoided.The proposed method is applied to a deep metal mine in China that is being excavated using the HLSCB method.The shear strength parameters of the high-level pillars obtained using the proposed method were input in the numerical simulations.The numerical results show that the recommended level heights and sizes of the high-level pillars and rooms in the mine are rational.

Nonlinear empirical failure criterion for rocks under triaxial compression

Based on existing triaxial compression experimental data,a new empirical failure criterion with wide applicability was proposed considering hydrostatic pressure,second stress invariance,and maximum shear stress.Four fitting evaluation indicators were used to verify the consistency of the new failure criterion,and the differences with the other 6 failure criteria were discussed.The characteristics of the new failure criteria in the principal stress space were finally analyzed.The results indicate that(1)the new failure criterion exhibits strong predictive ability for triaxial experiments and has good applicability for both intact and jointed rocks;(2)the influence of hydrostatic pressure on the failure surface exhibits a non-linear trend,and different hydrostatic pressure also exhibits different distribution patterns on the deviatoric stress plane,with a distribution characteristic pattern of hexagonal snowflake-regular hexagon.The maximum shear stress has a torsional effect on the new criterion,in the three-dimensional failure surface.The parameters a and b of the rock have an impact on the failure surface morphology of the new criterion function on the offset surface.

Rockburst criterion and evaluation method for potential rockburst pit depth considering excavation damage effect

Excavation-induced disturbances in deep tunnels will lead to deterioration of rock properties and formation of excavation damaged zone(EDZ).This excavation damage effect may affect the potential rockburst pit depth.Taking two diversion tunnels of Jinping II hydropower station for example,the relationship between rockburst pit depth and excavation damage effect is first surveyed.The results indicate that the rockburst pit depth in tunnels with severe damage to rock masses is relatively large.Subsequently,the excavation-induced damage effect is characterized by disturbance factor D based on the Hoek-Brown criterion and wave velocity method.It is found that the EDZ could be further divided into a high-damage zone(HDZ)with D=1 and weak-damage zone(WDZ),and D decays from one to zero linearly.For this,a quantitative evaluation method for potential rockburst pit depth is established by presenting a three-element rockburst criterion considering rock strength,geostress and disturbance factor.The evaluation results obtained by this method match well with actual observations.In addition,the weakening of rock mass strength promotes the formation and expansion of potential rockburst pits.The potential rockburst pit depth is positively correlated with HDZ and WDZ depths,and the HDZ depth has a significant contribution to the potential rockburst pit depth.

Least-squares reverse time migration in visco-acoustic media based on symplectic stereo-modeling method

The authors proposed a symplectic stereo-modeling method(SSM)in the Birkhoffian dynam-ics and apply it to the visco-acoustic least-squares reverse time migration(LSRTM).The SSM adopts ste-reo-modeling operator in space and symplectic Runge-Kutta scheme in time,resulting in great ability in suppressing numerical dispersion and long-time computing.These advantages are further confirmed by numerical dispersion analysis,long-time computation test and computational efficiency comparison.After these theoretical analyses and experiments,acoustic and visco-acoustic LSRTM are tested and compared between SSM method and the conventional symplectic method(CSM)using the fault and marmousi models.Meanwhile,dynamic source encoding and exponential decay moving average gradients method are adopted to reduce the computation cost and improve the convergence rate.The imaging results show that LSRTM based on visco-acoustic wave equations effectively takes into account the influence of viscosity can therefore compensate for the amplitude attenuation.Besides,SSM method not only has high numerical accuracy and computational efficiency,but also performs effectively in LSRTM.

Enhanced Hoek-Brown(H-B)criterion for rocks exposed to chemical corrosion

Underground constructions often encounter water environments,where water–rock interaction can increase porosity,thereby weakening engineering rocks.Correspondingly,the failure criterion for chemically corroded rocks becomes essential in the stability analysis and design of such structures.This study enhances the applicability of the Hoek-Brown(H-B)criterion for engineering structures operating in chemically corrosive conditions by introducing a kinetic porosity-dependent instantaneous mi(KPIM).A multiscale experimental investigation,including nuclear magnetic resonance(NMR),X-ray diffraction(XRD),scanning electron microscopy(SEM),pH and ion chromatography analysis,and triaxial compression tests,is employed to quantify pore structural changes and their linkage with the strength responses of limestone under coupled chemical-mechanical(C-M)conditions.By employing ion chromatography and NMR analysis,along with incorporating the principles of free-face dissolution theory accounting for both congruent and incongruent dissolution,a kinetic chemical corrosion model is developed.This model aims to calculate the kinetic porosity alterations within rocks exposed to varying H+concentrations and durations.Subsequently,utilizing the generalized mixture rule(GMR),the kinetic porositydependent mi is formulated.Evaluation of the KPIM-enhanced H-B criterion using compression test data from 5 types of rocks demonstrated a high level of consistency between the criterion and the experimental results,with a coefficient of determination greater than 0.96,a mean absolute percentage error less than 4.84%,and a root-mean-square deviation less than 5.95 MPa.Finally,the physical significance of the porosity-dependent instantaneous mi is clarified:it serves as an indicator of a rock’s capacity to leverage the confining pressure effect.

Strength criterion for crystalline rocks considering grain size effect and tensile-compressive strength ratio

The macroscopic mechanical properties of rocks are significantly influenced by their microstructure.As a material bonded by mineral grains,the grain morphology of crystalline rock is the primary factor influencing the strength.However,most strength criteria neglect the strength variations caused by different grain characteristics in rocks.Furthermore,the traditional linear criteria tend to overestimate tensile strength and exhibit apex singularity.To address these shortcomings,a piecewise strength criterion that considers the grain size effect has been proposed.A part of an ellipse was employed to construct the envelope of the tensive-shear region on the meridian plane,to accurately reproduce the low tensile-compressive strength ratio.Based on the analysis of experimental data,both linear and exponential modification functions that account for grain size effects were integrated into the proposed criterion.The corresponding finite element algorithm has been implemented.The accuracy and applicability of the proposed criterion were validated by comparing with the experimental data.

Maximum Correntropy Criterion-Based UKF for Loosely Coupling INS and UWB in Indoor Localization

Indoor positioning is a key technology in today’s intelligent environments,and it plays a crucial role in many application areas.This paper proposed an unscented Kalman filter(UKF)based on the maximum correntropy criterion(MCC)instead of the minimummean square error criterion(MMSE).This innovative approach is applied to the loose coupling of the Inertial Navigation System(INS)and Ultra-Wideband(UWB).By introducing the maximum correntropy criterion,the MCCUKF algorithm dynamically adjusts the covariance matrices of the system noise and the measurement noise,thus enhancing its adaptability to diverse environmental localization requirements.Particularly in the presence of non-Gaussian noise,especially heavy-tailed noise,the MCCUKF exhibits superior accuracy and robustness compared to the traditional UKF.The method initially generates an estimate of the predicted state and covariance matrix through the unscented transform(UT)and then recharacterizes the measurement information using a nonlinear regression method at the cost of theMCC.Subsequently,the state and covariance matrices of the filter are updated by employing the unscented transformation on the measurement equations.Moreover,to mitigate the influence of non-line-of-sight(NLOS)errors positioning accuracy,this paper proposes a k-medoid clustering algorithm based on bisection k-means(Bikmeans).This algorithm preprocesses the UWB distance measurements to yield a more precise position estimation.Simulation results demonstrate that MCCUKF is robust to the uncertainty of UWB and realizes stable integration of INS and UWB systems.

Three-dimensional limit variation analysis on the ultimate pullout capacity of the anchor cables based on the Hoek-Brown failure criterion

Only simplified two-dimensional model and a single failure mode are adopted to calculate the ultimate pullout capacity(UPC)of anchor cables in most previous research.This study focuses on a more comprehensive combination failure mode that consists of bond failure of an anchorage body and failure of an anchored rock mass.The three-dimensional ultimate pullout capacity of the anchor cables is calculated based on the Hoek-Brown failure criterion and variation analysis method.The numerical solution for the curvilinear function in fracture plane is obtained based on the finite difference theory,which more accurately reflects the failure state of the anchor cable,as opposed to that being assumed in advance.The results reveal that relying solely on a single failure mode for UPC calculations has limitations,as changes in parameter values not only directly impact the UPC value but also can alter the failure model and thus the calculation method.

Least-Squares及Galerkin谱元方法求解环形区域内的泊松方程

为研究基于Least-Squares变分及Galerkin变分两种形式的谱元方法的求解特性,推导了极坐标系中采用两种变分方法求解环形区域内Poisson方程时对应的弱解形式,采用Chebyshev多项式构造插值基函数进行空间离散,得到两种谱元方法对应的代数方程组,由此分析了系数矩阵结构的特点。数值计算结果显示:Least-Squares谱元方法为实现方程的降阶而引入新的求解变量,使得代数方程组形式更为复杂,但边界条件的处理比Galerkin谱元方法更为简单;两种谱元方法均能求解极坐标系中的Poisson方程且能获得高精度的数值解,二者绝对误差分布基本一致;固定单元内的插值阶数时,增加单元数可减小数值误差,且表现出代数精度的特点,误差降低速度较慢,而固定单元数时,在一定范围内数值误差随插值阶数的增加而减小的速度更快,表现出谱精度的特点;单元内插值阶数较高时,代数方程组系数矩阵的条件数急剧增多,方程组呈现病态,数值误差增大,这一特点限制了单元内插值阶数的取值。研究内容对深入了解两种谱元方法在极坐标系中求解Poisson方程时的特点、进一步采用相关分裂算法求解实际流动问题具有参考价值。

Regularized least-squares migration of simultaneous-source seismic data with adaptive singular spectrum analysis

Simultaneous-source acquisition has been recog- nized as an economic and efficient acquisition method, but the direct imaging of the simultaneous-source data produces migration artifacts because of the interference of adjacent sources. To overcome this problem, we propose the regularized least-squares reverse time migration method （RLSRTM） using the singular spectrum analysis technique that imposes sparseness constraints on the inverted model. Additionally, the difference spectrum theory of singular values is presented so that RLSRTM can be implemented adaptively to eliminate the migration artifacts. With numerical tests on a fiat layer model and a Marmousi model, we validate the superior imaging quality, efficiency and convergence of RLSRTM compared with LSRTM when dealing with simultaneoussource data, incomplete data and noisy data.

Characterizing and estimating rice brown spot disease severity using stepwise regression,principal component regression and partial least-square regression

Detecting plant health conditions plays a key role in farm pest management and crop protection. In this study, measurement of hyperspectral leaf reflectance in rice crop (Oryzasativa L.) was conducted on groups of healthy and infected leaves by the fungus Bipolaris oryzae (Helminthosporium oryzae Breda. de Hann) through the wavelength range from 350 to 2 500 nm. The percentage of leaf surface lesions was estimated and defined as the disease severity. Statistical methods like multiple stepwise regression, principal component analysis and partial least-square regression were utilized to calculate and estimate the disease severity of rice brown spot at the leaf level. Our results revealed that multiple stepwise linear regressions could efficiently estimate disease severity with three wavebands in seven steps. The root mean square errors (RMSEs) for training (n=210) and testing (n=53) dataset were 6.5% and 5.8%, respectively. Principal component analysis showed that the first principal component could explain approximately 80% of the variance of the original hyperspectral reflectance. The regression model with the first two principal components predicted a disease severity with RMSEs of 16.3% and 13.9% for the training and testing dataset, respec-tively. Partial least-square regression with seven extracted factors could most effectively predict disease severity compared with other statistical methods with RMSEs of 4.1% and 2.0% for the training and testing dataset, respectively. Our research demon-strates that it is feasible to estimate the disease severity of rice brown spot using hyperspectral reflectance data at the leaf level.

Single Point Positioning with Sequential Least-Squares Filter and Estimated Real-Time Stochastic Model

To obtain higher accurate position estimates, the stochastic model is estimated by using residual of observations, hence, the stochastic model describes the noise and bias in measurements more realistically. By using GPS data and broadcast ephemeris, the numerical results indicating the accurate position estimates at sub-meter level are obtainable.

AN ASSESSMENT OF THE MESHLESS WEIGHTED LEAST-SQUARE METHOD

The meshless weighted least-square (MWLS) method was developed based on the weighted least-square method. The method possesses several advantages, such as high accuracy, high stability and high e?ciency. Moreover, the coe?cient matrix obtained is symmetric and semi- positive de?nite. In this paper, the method is further examined critically. The e?ects of several parameters on the results of MWLS are investigated systematically by using a cantilever beam and an in?nite plate with a central circular hole. The numerical results are compared with those obtained by using the collocation-based meshless method (CBMM) and Galerkin-based meshless method (GBMM). The investigated parameters include the type of approximations, the type of weight functions, the number of neighbors of an evaluation point, as well as the manner in which the neighbors of an evaluation point are determined. This study shows that the displacement accuracy and convergence rate obtained by MWLS is comparable to that of the GBMM while the stress accuracy and convergence rate yielded by MWLS is even higher than that of GBMM. Furthermore, MWLS is much more e?cient than GBMM. This study also shows that the instability of CBMM is mainly due to the neglect of the equi- librium residuals at boundary nodes. In MWLS, the residuals of all the governing equations are minimized in a weighted least-square sense.

Mixed-Weights Least-Squares Stable Predictive Control Algorithm with Soft and Hard Constraints

Mixed-weight least-squares (MWLS) predictive control algorithm, compared with quadratic programming (QP) method, has the advantages of reducing the computer burden, quick calculation speed and dealing with the case in which the optimization is infeasible. But it can only deal with soft constraints. In order to deal with hard constraints and guarantee feasibility, an improved algorithm is proposed by recalculating the setpoint according to the hard constraints before calculating the manipulated variable and MWLS algorithm is used to satisfy the requirement of soft constraints for the system with the input constraints and output constraints. The algorithm can not only guarantee stability of the system and zero steady state error, but also satisfy the hard constraints of input and output variables. The simulation results show the improved algorithm is feasible and effective.

Multisource least-squares reverse-time migration with structure-oriented filtering

The technology of simultaneous-source acquisition of seismic data excited by several sources can significantly improve the data collection efficiency. However, direct imaging of simultaneous-source data or blended data may introduce crosstalk noise and affect the imaging quality. To address this problem, we introduce a structure-oriented filtering operator as preconditioner into the multisource least-squares reverse-time migration （LSRTM）. The structure-oriented filtering operator is a nonstationary filter along structural trends that suppresses crosstalk noise while maintaining structural information. The proposed method uses the conjugate-gradient method to minimize the mismatch between predicted and observed data, while effectively attenuating the interference noise caused by exciting several sources simultaneously. Numerical experiments using synthetic data suggest that the proposed method can suppress the crosstalk noise and produce highly accurate images.

Quantitative energy-dispersive X-ray fluorescence analysis for unknown samples using full-spectrum least-squares regression

The full-spectrum least-squares(FSLS) method is introduced to perform quantitative energy-dispersive X-ray fluorescence analysis for unknown solid samples.Based on the conventional least-squares principle, this spectrum evaluation method is able to obtain the background-corrected and interference-free net peaks, which is significant for quantization analyses. A variety of analytical parameters and functions to describe the features of the fluorescence spectra of pure elements are used and established, such as the mass absorption coefficient, the Gi factor, and fundamental fluorescence formulas. The FSLS iterative program was compiled in the C language. The content of each component should reach the convergence criterion at the end of the calculations. After a basic theory analysis and experimental preparation, 13 national standard soil samples were detected using a spectrometer to test the feasibility of using the algorithm. The results show that the calculated contents of Ti, Fe, Ni, Cu, and Zn have the same changing tendency as the corresponding standard content in the 13 reference samples. Accuracies of 0.35% and 14.03% are obtained, respectively, for Fe and Ti, whose standard concentrations are 8.82% and 0.578%, respectively. However, the calculated results of trace elements (only tens of lg/g) deviate from the standard values. This may be because of measurement accuracy and mutual effects between the elements.

LEAST-SQUARES METHOD-BASED FEATURE FITTING AND EXTRACTION IN REVERSE ENGINEERING

The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.

Solution of shallow-water equations using least-squares finite-element method

A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.

Adaptive Control of a Two-Link Robot Using Batch Least-Square Identifier

We design a regulation-triggered adaptive controller for robot manipulators to efficiently estimate unknown parameters and to achieve asymptotic stability in the presence of coupled uncertainties.Robot manipulators are widely used in telemanipulation systems where they are subject to model and environmental uncertainties.Using conventional control algorithms on such systems can cause not only poor control performance,but also expensive computational costs and catastrophic instabilities.Therefore,system uncertainties need to be estimated through designing a computationally efficient adaptive control law.We focus on robot manipulators as an example of a highly nonlinear system.As a case study,a 2-DOF manipulator subject to four parametric uncertainties is investigated.First,the dynamic equations of the manipulator are derived,and the corresponding regressor matrix is constructed for the unknown parameters.For a general nonlinear system,a theorem is presented to guarantee the asymptotic stability of the system and the convergence of parameters'estimations.Finally,simulation results are discussed for a two-link manipulator,and the performance of the proposed scheme is thoroughly evaluated.