The proposed robust reversible watermarking algorithm addresses the compatibility challenges between robustness and reversibility in existing video watermarking techniques by leveraging scene smoothness for frame grou...The proposed robust reversible watermarking algorithm addresses the compatibility challenges between robustness and reversibility in existing video watermarking techniques by leveraging scene smoothness for frame grouping videos.Grounded in the H.264 video coding standard,the algorithm first employs traditional robust watermark stitching technology to embed watermark information in the low-frequency coefficient domain of the U channel.Subsequently,it utilizes histogram migration techniques in the high-frequency coefficient domain of the U channel to embed auxiliary information,enabling successful watermark extraction and lossless recovery of the original video content.Experimental results demonstrate the algorithm’s strong imperceptibility,with each embedded frame in the experimental videos achieving a mean peak signal-to-noise ratio of 49.3830 dB and a mean structural similarity of 0.9996.Compared with the three comparison algorithms,the performance of the two experimental indexes is improved by 7.59%and 0.4%on average.At the same time,the proposed algorithm has strong robustness to both offline and online attacks:In the face of offline attacks,the average normalized correlation coefficient between the extracted watermark and the original watermark is 0.9989,and the average bit error rate is 0.0089.In the face of online attacks,the normalized correlation coefficient between the extracted watermark and the original watermark is 0.8840,and the mean bit error rate is 0.2269.Compared with the three comparison algorithms,the performance of the two experimental indexes is improved by 1.27%and 18.16%on average,highlighting the algorithm’s robustness.Furthermore,the algorithm exhibits low computational complexity,with the mean encoding and the mean decoding time differentials during experimental video processing being 3.934 and 2.273 s,respectively,underscoring its practical utility.展开更多
We introduce the notion of K-very smoothness which is a generalization of very smoothness in Banach spaces. A necessary and sufficient condition for a Banach space to be K-very smooth is obtained. We also consider som...We introduce the notion of K-very smoothness which is a generalization of very smoothness in Banach spaces. A necessary and sufficient condition for a Banach space to be K-very smooth is obtained. We also consider some relations between K-very smoothness and other geometrical notions.展开更多
Fault degradation prognostic, which estimates the time before a failure occurs and process breakdowns, has been recognized as a key component in maintenance strategies nowadays. Fault degradation processes are, in gen...Fault degradation prognostic, which estimates the time before a failure occurs and process breakdowns, has been recognized as a key component in maintenance strategies nowadays. Fault degradation processes are, in general,slowly varying and can be modeled by autoregressive models. However, industrial processes always show typical nonstationary nature, which may bring two challenges: how to capture fault degradation information and how to model nonstationary processes. To address the critical issues, a novel fault degradation modeling and online fault prognostic strategy is developed in this paper. First, a fault degradation-oriented slow feature analysis(FDSFA) algorithm is proposed to extract fault degradation directions along which candidate fault degradation features are extracted. The trend ability assessment is then applied to select major fault degradation features. Second, a key fault degradation factor(KFDF) is calculated to characterize the fault degradation tendency by combining major fault degradation features and their stability weighting factors. After that, a time-varying regression model with temporal smoothness regularization is established considering nonstationary characteristics. On the basis of updating strategy, an online fault prognostic model is further developed by analyzing and modeling the prediction errors. The performance of the proposed method is illustrated with a real industrial process.展开更多
Smoothness prior approach for spectral smoothing is investigated using Fourier frequency filter analysis.We show that the regularization parameter in penalized least squares could continuously control the bandwidth of...Smoothness prior approach for spectral smoothing is investigated using Fourier frequency filter analysis.We show that the regularization parameter in penalized least squares could continuously control the bandwidth of low-pass filter.Besides,due to its property of interpolating the missing values automatically and smoothly,a spectral baseline correction algorithm based on the approach is proposed.This algorithm generally comprises spectral peak detection and baseline estimation.First,the spectral peak regions are detected and identified according to the second derivatives.Then,generalized smoothness prior approach combining identification information could estimate the baseline in peak regions.Results with both the simulated and real spectra show accurate baseline-corrected signals with this method.展开更多
Let h, be a measurable function defined on R^+ ×R^+. Let Ω ∈ L(log L^+)^υq (S^n1-1 × S^n2-1) (1≤ υq ≤ 2) be homogeneous of degree zero and satisfy certain cancellation conditions. We show that...Let h, be a measurable function defined on R^+ ×R^+. Let Ω ∈ L(log L^+)^υq (S^n1-1 × S^n2-1) (1≤ υq ≤ 2) be homogeneous of degree zero and satisfy certain cancellation conditions. We show that the singular integral Tf(x1,x2)=p.v.∫∫R^n1+n2 Ω(y′1,y′2)h(|y1|,|y2|)/|y1|^n1|y2|^n2 f(x1-y1,x2-y2)dy1dy2maps from Sp,q^α1,α2F(R^n1×R^n2)boundedly to itself for 1 〈 p, q 〈 ∞, α1, α2 ∈R.展开更多
Based on the traditional fifth-order weighted essentially non-oscillatory(WENO)scheme,a smoothness indicator is introduced to improve the capability of WENO schemes for resolving short waves.In the construction of the...Based on the traditional fifth-order weighted essentially non-oscillatory(WENO)scheme,a smoothness indicator is introduced to improve the capability of WENO schemes for resolving short waves.In the construction of the new smoothness indicator,the proportion of the first-order term in the original smoothness indicator is reduced by replacing the square of the first-order term with the product of the first-order and the third-order terms.To preserve the fifth-order of convergence rate,the smoothness indicator is combined with the method of Borges,et al.The numerical results show that the proposed schemes are more suitable for simulating turbulent flows or aeroacoustics problems than the previous fifth-order WENO schemes,thanks to its improved resolution on short waves.展开更多
Extending the results of [4] in the univariate case, in this paper we prove that the bivariate interpolation polynomials of Hermite-Fejér based on the Chebyshev nodes of the first kind, those of Lagrange based o...Extending the results of [4] in the univariate case, in this paper we prove that the bivariate interpolation polynomials of Hermite-Fejér based on the Chebyshev nodes of the first kind, those of Lagrange based on the Chebyshev nodes of second kind and ±1, and those of bivariate Shepard operators, have the property of partial preservation of global smoothness, with respect to various bivariate moduli of continuity.展开更多
Nonlinear m-term approximation plays an important role in machine learning, signal processing and statistical estimating. In this paper by means of a nondecreasing dominated function, a greedy adaptive compression num...Nonlinear m-term approximation plays an important role in machine learning, signal processing and statistical estimating. In this paper by means of a nondecreasing dominated function, a greedy adaptive compression numerical algorithm in the best m -term approximation with regard to tensor product wavelet-type basis is pro-posed. The algorithm provides the asymptotically optimal approximation for the class of periodic functions with mixed Besov smoothness in the L q norm. Moreover, it depends only on the expansion of function f by tensor pro-duct wavelet-type basis, but neither on q nor on any special features of f.展开更多
Automatically assessing fabric smoothness grade is very important in the evaluation of fabric appearance.A system for objectively evaluating the fabric smoothness grade based on a grating projection unit and double co...Automatically assessing fabric smoothness grade is very important in the evaluation of fabric appearance.A system for objectively evaluating the fabric smoothness grade based on a grating projection unit and double colored CCD(short form of charge coupled device) was constructed in this paper.Two images captured by different CCD compensated each other which reduced the influence of noises.The application of the four-step phase-shifting method enabled the calculation of the exact phase in a point easy and quick.A large amount of 3D points with three coordinates X,Y and Z were obtained precisely making the definition and calculation of fabric smoothness characters easy.Then four parameters which intuitively denoted the fabric smoothness degree were obtained.Finally,a proper neural network was built,which successfully performed the fabric smoothness classification.The experimental results show that the system is applicable for all the fabric whatever pattern or color.The experimental grades provided by this grating projection system are also highly consistent with the subjective results.展开更多
For garment or fabric appearance, the cloth smoothness grade is one of the most important performance factors in textile and garment community. In this paper, on the base of Rough Set Theory,a new objective method for...For garment or fabric appearance, the cloth smoothness grade is one of the most important performance factors in textile and garment community. In this paper, on the base of Rough Set Theory,a new objective method for fabric smoothness grade evaluation was constructed. The objective smoothness grading model took the parameters of 120 AATCC replicas' point-sampled models as the conditional attributes and formed the smoothness grading decision table. Then, NS discretization method and genetic algorithm reduction method were used in the attributes discretization and feature reduction. Finally, the grading model was expressed as simple and intuitional classification rules. The simulation results show the validity of the fabric smoothness grading system which is built on the use of rough sets.展开更多
A method of real-time data smoothness which is applied in a hardware-in-the-loop (HIL) simulation platform for a plug-in hybrid electric vehicle synthetical power device is described. The input signal of the platfor...A method of real-time data smoothness which is applied in a hardware-in-the-loop (HIL) simulation platform for a plug-in hybrid electric vehicle synthetical power device is described. The input signal of the platform comes from a AC/DC switch power with output containing noises. A linear slide average arithmetic is applied to smooth the noises. To average a certain number input sample signals, this method can decrease the noises voltage level, which meet the requirement of the simulation platform. The efficiency and signal delay time are presented to describe the result of this method, and a statistical index is used to judge the arithmetic' s efficiency. The tests results show that the arithmetic fit the requirement of the HIL simulation platform.展开更多
For the step-weight function , we prove that the Holder spaces ∧a,p on the interval [-1,1], defined in terms of moduli of smoothness with the step-weight function ,are linearly isomorphic to some sequence spaces, an...For the step-weight function , we prove that the Holder spaces ∧a,p on the interval [-1,1], defined in terms of moduli of smoothness with the step-weight function ,are linearly isomorphic to some sequence spaces, and the isomorphism is given by the cofficients of function with respect to a system of orthonormal splines with knots uniformly distributed according to the measure with density . In case ∧a,p is contained in the space of continuous functions, we give a discrete characterization of this space, using only values of function at the appropriate knots. Application of these results to characterize the order of polynomial approximation is presented.展开更多
The criterion for k-smooth points of the Orlicz sequence space endowed with the Orlicz norm is proved. The necessary and sufficient conditions of k-smoothness of l M and l (M ) are obtained, respectively. Finally, w...The criterion for k-smooth points of the Orlicz sequence space endowed with the Orlicz norm is proved. The necessary and sufficient conditions of k-smoothness of l M and l (M ) are obtained, respectively. Finally, we give the counterexamples which show that previous results are not true.展开更多
In this note, we establish a new formulation of smoothness conditions for piecewise polynomial (: =pp) functions in terms of the B-net representation in the general n-dimensional setting. It plays an important role fo...In this note, we establish a new formulation of smoothness conditions for piecewise polynomial (: =pp) functions in terms of the B-net representation in the general n-dimensional setting. It plays an important role for 2-dimensional setting in the constructive proof of the fact that the spaces of polynomial splines with smoothness rand total degree k≥3r+2 over arbitrary triangulations achieve the optimal approximation order with the approximation constant depending only on k and the smallest angle of the partition in [5].展开更多
In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SP...In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.展开更多
Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strat...Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.展开更多
基金supported in part by the National Natural Science Foundation of China under Grants 62202496,62272478the Basic Frontier Innovation Project of Engineering university of People Armed Police under Grants WJY202314,WJY202221.
文摘The proposed robust reversible watermarking algorithm addresses the compatibility challenges between robustness and reversibility in existing video watermarking techniques by leveraging scene smoothness for frame grouping videos.Grounded in the H.264 video coding standard,the algorithm first employs traditional robust watermark stitching technology to embed watermark information in the low-frequency coefficient domain of the U channel.Subsequently,it utilizes histogram migration techniques in the high-frequency coefficient domain of the U channel to embed auxiliary information,enabling successful watermark extraction and lossless recovery of the original video content.Experimental results demonstrate the algorithm’s strong imperceptibility,with each embedded frame in the experimental videos achieving a mean peak signal-to-noise ratio of 49.3830 dB and a mean structural similarity of 0.9996.Compared with the three comparison algorithms,the performance of the two experimental indexes is improved by 7.59%and 0.4%on average.At the same time,the proposed algorithm has strong robustness to both offline and online attacks:In the face of offline attacks,the average normalized correlation coefficient between the extracted watermark and the original watermark is 0.9989,and the average bit error rate is 0.0089.In the face of online attacks,the normalized correlation coefficient between the extracted watermark and the original watermark is 0.8840,and the mean bit error rate is 0.2269.Compared with the three comparison algorithms,the performance of the two experimental indexes is improved by 1.27%and 18.16%on average,highlighting the algorithm’s robustness.Furthermore,the algorithm exhibits low computational complexity,with the mean encoding and the mean decoding time differentials during experimental video processing being 3.934 and 2.273 s,respectively,underscoring its practical utility.
文摘We introduce the notion of K-very smoothness which is a generalization of very smoothness in Banach spaces. A necessary and sufficient condition for a Banach space to be K-very smooth is obtained. We also consider some relations between K-very smoothness and other geometrical notions.
基金Project(U1709211) supported by NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization,ChinaProject(ICT2021A15) supported by the State Key Laboratory of Industrial Control Technology,Zhejiang University,ChinaProject(TPL2019C03) supported by Open Fund of Science and Technology on Thermal Energy and Power Laboratory,China。
文摘Fault degradation prognostic, which estimates the time before a failure occurs and process breakdowns, has been recognized as a key component in maintenance strategies nowadays. Fault degradation processes are, in general,slowly varying and can be modeled by autoregressive models. However, industrial processes always show typical nonstationary nature, which may bring two challenges: how to capture fault degradation information and how to model nonstationary processes. To address the critical issues, a novel fault degradation modeling and online fault prognostic strategy is developed in this paper. First, a fault degradation-oriented slow feature analysis(FDSFA) algorithm is proposed to extract fault degradation directions along which candidate fault degradation features are extracted. The trend ability assessment is then applied to select major fault degradation features. Second, a key fault degradation factor(KFDF) is calculated to characterize the fault degradation tendency by combining major fault degradation features and their stability weighting factors. After that, a time-varying regression model with temporal smoothness regularization is established considering nonstationary characteristics. On the basis of updating strategy, an online fault prognostic model is further developed by analyzing and modeling the prediction errors. The performance of the proposed method is illustrated with a real industrial process.
基金Supported by the National Basic Research Program of China(61178072)
文摘Smoothness prior approach for spectral smoothing is investigated using Fourier frequency filter analysis.We show that the regularization parameter in penalized least squares could continuously control the bandwidth of low-pass filter.Besides,due to its property of interpolating the missing values automatically and smoothly,a spectral baseline correction algorithm based on the approach is proposed.This algorithm generally comprises spectral peak detection and baseline estimation.First,the spectral peak regions are detected and identified according to the second derivatives.Then,generalized smoothness prior approach combining identification information could estimate the baseline in peak regions.Results with both the simulated and real spectra show accurate baseline-corrected signals with this method.
文摘Let h, be a measurable function defined on R^+ ×R^+. Let Ω ∈ L(log L^+)^υq (S^n1-1 × S^n2-1) (1≤ υq ≤ 2) be homogeneous of degree zero and satisfy certain cancellation conditions. We show that the singular integral Tf(x1,x2)=p.v.∫∫R^n1+n2 Ω(y′1,y′2)h(|y1|,|y2|)/|y1|^n1|y2|^n2 f(x1-y1,x2-y2)dy1dy2maps from Sp,q^α1,α2F(R^n1×R^n2)boundedly to itself for 1 〈 p, q 〈 ∞, α1, α2 ∈R.
基金Supported by the National Natural Science Foundation of China(50830201,11102179)the Nanjing University of Aeronautics and Astronautics Research Funding(NP 2011033)
文摘Based on the traditional fifth-order weighted essentially non-oscillatory(WENO)scheme,a smoothness indicator is introduced to improve the capability of WENO schemes for resolving short waves.In the construction of the new smoothness indicator,the proportion of the first-order term in the original smoothness indicator is reduced by replacing the square of the first-order term with the product of the first-order and the third-order terms.To preserve the fifth-order of convergence rate,the smoothness indicator is combined with the method of Borges,et al.The numerical results show that the proposed schemes are more suitable for simulating turbulent flows or aeroacoustics problems than the previous fifth-order WENO schemes,thanks to its improved resolution on short waves.
文摘Extending the results of [4] in the univariate case, in this paper we prove that the bivariate interpolation polynomials of Hermite-Fejér based on the Chebyshev nodes of the first kind, those of Lagrange based on the Chebyshev nodes of second kind and ±1, and those of bivariate Shepard operators, have the property of partial preservation of global smoothness, with respect to various bivariate moduli of continuity.
基金Supported by National Natural Science Foundation of China (No. 60872161, 10501026, 60675010 and 10626029)Natural Science Foundation of Tianjin (No. 08JCYBJC09600)China Postdoctoral Science Foundation ( No. 20070420708).
文摘Nonlinear m-term approximation plays an important role in machine learning, signal processing and statistical estimating. In this paper by means of a nondecreasing dominated function, a greedy adaptive compression numerical algorithm in the best m -term approximation with regard to tensor product wavelet-type basis is pro-posed. The algorithm provides the asymptotically optimal approximation for the class of periodic functions with mixed Besov smoothness in the L q norm. Moreover, it depends only on the expansion of function f by tensor pro-duct wavelet-type basis, but neither on q nor on any special features of f.
文摘Automatically assessing fabric smoothness grade is very important in the evaluation of fabric appearance.A system for objectively evaluating the fabric smoothness grade based on a grating projection unit and double colored CCD(short form of charge coupled device) was constructed in this paper.Two images captured by different CCD compensated each other which reduced the influence of noises.The application of the four-step phase-shifting method enabled the calculation of the exact phase in a point easy and quick.A large amount of 3D points with three coordinates X,Y and Z were obtained precisely making the definition and calculation of fabric smoothness characters easy.Then four parameters which intuitively denoted the fabric smoothness degree were obtained.Finally,a proper neural network was built,which successfully performed the fabric smoothness classification.The experimental results show that the system is applicable for all the fabric whatever pattern or color.The experimental grades provided by this grating projection system are also highly consistent with the subjective results.
文摘For garment or fabric appearance, the cloth smoothness grade is one of the most important performance factors in textile and garment community. In this paper, on the base of Rough Set Theory,a new objective method for fabric smoothness grade evaluation was constructed. The objective smoothness grading model took the parameters of 120 AATCC replicas' point-sampled models as the conditional attributes and formed the smoothness grading decision table. Then, NS discretization method and genetic algorithm reduction method were used in the attributes discretization and feature reduction. Finally, the grading model was expressed as simple and intuitional classification rules. The simulation results show the validity of the fabric smoothness grading system which is built on the use of rough sets.
文摘A method of real-time data smoothness which is applied in a hardware-in-the-loop (HIL) simulation platform for a plug-in hybrid electric vehicle synthetical power device is described. The input signal of the platform comes from a AC/DC switch power with output containing noises. A linear slide average arithmetic is applied to smooth the noises. To average a certain number input sample signals, this method can decrease the noises voltage level, which meet the requirement of the simulation platform. The efficiency and signal delay time are presented to describe the result of this method, and a statistical index is used to judge the arithmetic' s efficiency. The tests results show that the arithmetic fit the requirement of the HIL simulation platform.
文摘For the step-weight function , we prove that the Holder spaces ∧a,p on the interval [-1,1], defined in terms of moduli of smoothness with the step-weight function ,are linearly isomorphic to some sequence spaces, and the isomorphism is given by the cofficients of function with respect to a system of orthonormal splines with knots uniformly distributed according to the measure with density . In case ∧a,p is contained in the space of continuous functions, we give a discrete characterization of this space, using only values of function at the appropriate knots. Application of these results to characterize the order of polynomial approximation is presented.
基金Project supported by the National Natural Science Foundation of China (Grant No.10971129)
文摘The criterion for k-smooth points of the Orlicz sequence space endowed with the Orlicz norm is proved. The necessary and sufficient conditions of k-smoothness of l M and l (M ) are obtained, respectively. Finally, we give the counterexamples which show that previous results are not true.
文摘In this note, we establish a new formulation of smoothness conditions for piecewise polynomial (: =pp) functions in terms of the B-net representation in the general n-dimensional setting. It plays an important role for 2-dimensional setting in the constructive proof of the fact that the spaces of polynomial splines with smoothness rand total degree k≥3r+2 over arbitrary triangulations achieve the optimal approximation order with the approximation constant depending only on k and the smallest angle of the partition in [5].
基金supported by the National Natural Science Foundation of China(Grant No.52201323).
文摘In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.
基金funded by the National Key Research and Development Program of China(Grant No.2023YFC3008300,Grant No.2019YFC1509702)the National Natural Science Foundation of China(Grant No.42172296).
文摘Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.