Fast and Accurate Predictor-Corrector Methods Using Feedback-Accelerated Picard Iteration for Strongly Nonlinear Problems

Although predictor-corrector methods have been extensively applied,they might not meet the requirements of practical applications and engineering tasks,particularly when high accuracy and efficiency are necessary.A novel class of correctors based on feedback-accelerated Picard iteration(FAPI)is proposed to further enhance computational performance.With optimal feedback terms that do not require inversion of matrices,significantly faster convergence speed and higher numerical accuracy are achieved by these correctors compared with their counterparts;however,the computational complexities are comparably low.These advantages enable nonlinear engineering problems to be solved quickly and accurately,even with rough initial guesses from elementary predictors.The proposed method offers flexibility,enabling the use of the generated correctors for either bulk processing of collocation nodes in a domain or successive corrections of a single node in a finite difference approach.In our method,the functional formulas of FAPI are discretized into numerical forms using the collocation approach.These collocated iteration formulas can directly solve nonlinear problems,but they may require significant computational resources because of the manipulation of high-dimensionalmatrices.To address this,the collocated iteration formulas are further converted into finite difference forms,enabling the design of lightweight predictor-corrector algorithms for real-time computation.The generality of the proposed method is illustrated by deriving new correctors for three commonly employed finite-difference approaches:the modified Euler approach,the Adams-Bashforth-Moulton approach,and the implicit Runge-Kutta approach.Subsequently,the updated approaches are tested in solving strongly nonlinear problems,including the Matthieu equation,the Duffing equation,and the low-earth-orbit tracking problem.The numerical findings confirm the computational accuracy and efficiency of the derived predictor-corrector algorithms.

A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

Contact Stress Reliability Analysis Model for Cylindrical Gear with Circular Arc Tooth Trace Based on an Improved Metamodel

Although there is currently no unified standard theoretical formula for calculating the contact stress of cylindrical gears with a circular arc tooth trace(referred to as CATT gear),a mathematical model for determining the contact stress of CATT gear is essential for studying how parameters affect its contact stress and building the contact stress limit state equation for contact stress reliability analysis.In this study,a mathematical relationship between design parameters and contact stress is formulated using the KrigingMetamodel.To enhance the model’s accuracy,we propose a new hybrid algorithm that merges the genetic algorithm with the Quantum Particle Swarm optimization algorithm,leveraging the strengths of each.Additionally,the“parental inheritance+self-learning”optimization model is used to fine-tune the KrigingMetamodel’s parameters.Following this,amathematicalmodel for calculating the contact stress of Variable Hyperbolic Circular-Arc-Tooth-Trace(VH-CATT)gears using the optimized Kriging model was developed.We then examined how different gear parameters affect the VH-CATT gears’contact stress.Our simulation results show:(1)Improvements in R2,RMSE,and RMAE.R2 rose from0.9852 to 0.9974(a 1.22%increase),nearing 1,suggesting the optimized Kriging Metamodel’s global error is minimized.Meanwhile,RMSE dropped from3.9210 to 1.6492,a decline of 57.94%.The global error of the GA-IQPSO-Kriging algorithm was also reduced,with RMAE decreasing by 58.69%from 0.1823 to 0.0753,showing the algorithm’s enhanced precision.In a comparison of ten experimental groups selected randomly,the GA-IQPSO-Kriging and FEM-based contact analysis methods were used to measure contact stress.Results revealed a maximum error of 12.11667 MPA,which represents 2.85%of the real value.(2)Several factors,including the pressure angle,tooth width,modulus,and tooth line radius,are inversely related to contact stress.The descending order of their impact on the contact stress is:tooth line radius>modulus>pressure angle>tooth width.(3)Complex interactions are noted among various parameters.Specifically,when the tooth line radius interacts with parameters such as pressure angle,tooth width,and modulus,the resulting stress contour is nonlinear,showcasing amultifaceted contour plane.However,when tooth width,modulus,and pressure angle interact,the stress contour is nearly linear,and the contour plane is simpler,indicating a weaker coupling among these factors.

Stability Analysis of Nonlinear Models of Nose Landing Gear Shimmy

Shimmy can reduce the service life of the nose landing gear, affect ride comfort, and even cause fuselage damage leading to aircraft crashes. Taking a light aircraft as the research object, the torsional freedom of landing gear around strut axis and lateral deformation of tire are considered. Since the landing gear shimmy is a nonlinear system, a nonlinear mechanical model of the front landing gear shimmy is established. Sobol index method is proposed to analyze the influence of structural parameters on the stability region of the nose landing gear, and Routh-Huritz criterion is used to verify the reliability of the analysis results of Sobol index method. We analyse the effect of torsional stiffness of strut, caster length, rated initial tire inflation pressure, rake angle, and vertical force on the stability region of theront landing gear. And the research shows that the optimization of the torsional stiffness of the strut and the caster length of the nose landing gear should be emphasized, and the influence of vertical force on the stability region of the nose landing gear should be paid attention to.

Time-Varying Mesh Stiffness Calculation and Dynamic Modeling of Spiral Bevel Gear with Spalling Defects

Time-varying mesh stiffness(TVMS)is a vital internal excitation source for the spiral bevel gear(SBG)transmission system.Spalling defect often causes decrease in gear mesh stiffness and changes the dynamic characteristics of the gear system,which further increases noise and vibration.This paper aims to calculate the TVMS and establish dynamic model of SBG with spalling defect.In this study,a novel analytical model based on slice method is proposed to calculate the TVMS of SBG considering spalling defect.Subsequently,the influence of spalling defect on the TVMS is studied through a numerical simulation,and the proposed analytical model is verified by a finite element model.Besides,an 8-degrees-of-freedom dynamic model is established for SBG transmission system.Incorporating the spalling defect into TVMS,the dynamic responses of spalled SBG are analyzed.The numerical results indicate that spalling defect would cause periodic impact in time domain.Finally,an experiment is designed to verify the proposed dynamic model.The experimental results show that the spalling defect makes the response characterized by periodic impact with the rotating frequency of spalled pinion.

Statistical Modification Analysis of Helical Planetary Gears based on Response Surface Method and Monte Carlo Simulation

Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random errors on gear modification effects. In order to investigate the uncertainties of tooth modification amount variations on system＇s dynamic behaviors of a helical planetary gears, an analytical dynamic model including tooth modification parameters is proposed to carry out a deterministic analysis on the dynamics of a helical planetary gear. The dynamic meshing forces as well as the dynamic transmission errors of the sun-planet 1 gear pair with and without tooth modifications are computed and compared to show the effectiveness of tooth modifications on gear dynamics enhancement. By using response surface method, a fitted regression model for the dynamic transmission error（DTE） fluctuations is established to quantify the relationship between modification amounts and DTE fluctuations. By shifting the inevitable random errors arousing from manufacturing and installing process to tooth modification amount variations, a statistical tooth modification model is developed and a methodology combining Monte Carlo simulation and response surface method is presented for uncertainty analysis of tooth modifications. The uncertainly analysis reveals that the system＇s dynamic behaviors do not obey the normal distribution rule even though the design variables are normally distributed. In addition, a deterministic modification amount will not definitely achieve an optimal result for both static and dynamic transmission error fluctuation reduction simultaneously.

Modification Design Method for an Enveloping Hourglass Worm Gear with Consideration of Machining and Misalignment Errors

The influences of machining and misalignment errors play a very critical role in the performance of the anti-backlash double-roller enveloping hourglass worm gear(ADEHWG).However,a corresponding efficient method for eliminating or reducing these errors on the tooth profile of the ADEHWG is seldom reported.The gear engagement equation and tooth profile equation for considering six different errors that could arise from the machining and gear misalignment are derived from the theories of differential geometry and gear meshing.Also,the tooth contact analysis(TCA) is used to systematically investigate the influence of the machining and misalignment errors on the contact curves and the tooth profile by means of numerical analysis and three-dimensional solid modeling.The research results show that vertical angular misalignment of the worm wheel(Δβ) has the strongest influences while the tooth angle error(Δα) has the weakest influences on the contact curves and the tooth profile.A novel efficient approach is proposed and used to minimize the effect of the errors in manufacturing by changing the radius of the grinding wheel and the approaching point of contact.The results from the TCA and the experiment demonstrate that this tooth profile design modification method can indeed reduce the machining and misalignment errors.This modification design method is helpful in understanding the manufacturing technology of the ADEHWG.

Tooth surface geometry optimization of spiral bevel and hypoid gears generated by duplex helical method with circular profile blade

In order to effectively improve meshing performance of spiral bevel and hypoid gears generated by the duplex helical method, the effects of straight lined and circular cutting edges profile on meshing and contact of spiral bevel and hypoid gears were investigated analytically. Firstly, a mathematical model of spiral bevel and hypoid gears with circular blade profile was established according to the cutting characteristics of the duplex helical method. Based on a hypoid gear drive, the tooth bearings and the functions of transmission errors of four design cases were analyzed respectively by the use of the tooth contact analysis(TCA), and the contact stresses of the four design cases were analyzed and compared using simulation software. Finally, the curvature radius of the circular profile blade was optimized. The results show that the contact stresses are availably reduced, and the areas of edge contact and severe contact stresses can be avoided by selecting appropriate circular blade profile. In addition, the convex and concave sides are separately modified by the use of different curvature radii of inside and outside blades, which can increase the flexibility of the duplex helical method.

A novel method to predict static transmission error for spur gear pair based on accuracy grade

This paper proposes a novel method to predict the spur gear pair’s static transmission error based on the accuracy grade,in which manufacturing errors(MEs),assembly errors(AEs),tooth deflections(TDs)and profile modifications(PMs)are considered.For the prediction,a discrete gear model for generating the error tooth profile based on the ISO accuracy grade is presented.Then,the gear model and a tooth deflection model for calculating the tooth compliance on gear meshing are coupled with the transmission error model to make the prediction by checking the interference status between gear and pinion.The prediction method is validated by comparison with the experimental results from the literature,and a set of cases are simulated to study the effects of MEs,AEs,TDs and PMs on the static transmission error.In addition,the time-varying backlash caused by both MEs and AEs,and the contact ratio under load conditions are also investigated.The results show that the novel method can effectively predict the range of the static transmission error under different accuracy grades.The prediction results can provide references for the selection of gear design parameters and the optimization of transmission performance in the design stage of gear systems.

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.

Tooth surface error correction of hypoid gears machined by duplex helical method

In this work,synchronous cutting of concave and convex surfaces was achieved using the duplex helical method for the hypoid gear,and the problem of tooth surface error correction was studied.First,the mathematical model of the hypoid gears machined by the duplex helical method was established.Second,the coordinates of discrete points on the tooth surface were obtained by measurement center,and the normal errors of the discrete points were calculated.Third,a tooth surface error correction model is established,and the tooth surface error was corrected using the Levenberg-Marquard algorithm with trust region strategy and least square method.Finally,grinding experiments were carried out on the machining parameters obtained by Levenberg-Marquard algorithm with trust region strategy,which had a better effect on tooth surface error correction than the least square method.After the tooth surface error is corrected,the maximum absolute error is reduced from 30.9μm before correction to 6.8μm,the root mean square of the concave error is reduced from 15.1 to 2.1μm,the root mean square of the convex error is reduced from 10.8 to 1.8μm,and the sum of squared errors of the concave and convex surfaces was reduced from 15471 to 358μm^(2).It is verified that the Levenberg-Marquard algorithm with trust region strategy has a good accuracy for the tooth surface error correction of hypoid gear machined by duplex helical method.

Quasi Ellipsoid Gear Surface Reconstruction Based on Meshless Local Petrov-Galerkin Method and Transmission Characteristic

Special transmission 3D model simulation must be based on surface discretization and reconstruction, but special transmission usually has complicated tooth shape and movement, so present software can＇t provide technical support for special transmission 3D model simulation. Currently, theoretical calculation and experimental method are difficult to exactly solve special transmission contact analysis problem. How to reduce calculation and computer memories consume and meet calculation precision is key to resolve special transmission contact analysis problem. According to 3D model simulation and surface reconstruction of quasi ellipsoid gear is difficulty, this paper employes meshless local Petrov-Galerkin （MLPG） method. In order to reduce calculation and computer memories consume, we disperse tooth mesh into finite points--sparseness points cloud or grid mesh, and then we do interpolation reconstruction in some necessary place of the 3D surface model during analysis. Moving least square method （MLSM） is employed for tooth mesh interpolation reconstruction, there are some advantages to do interpolation by means of MLSM, such as high precision, good flexibility and no require of tooth mesh discretization into units. We input the quasi ellipsoid gear reconstruction model into simulation software, we complete tooth meshing simulation. Simulation transmission ratio during meshing period was obtained, compared with theoretical transmission ratio, the result inosculate preferably. The method using curve reconstruction realizes surface reconstruction, reduce simulation calculation enormously, so special gears simulation can be realized by minitype computer. The method provides a novel solution for special transmission 3D model simulation analysis and contact analysis.

New Predictor-Corrector Methods Based on Exponential Time Differencing forSystems of Nonlinear Differential Equations

We present the new predictor-corrector methods for systems of nonlinear differential equations, based on the method of exponential time differencing. We compare the present schemes with the explicit multistep exponential time differencing and Adams–Bashforth–Moulton method. The numerical results show that the schemes are more accurate and more efficient than Adams predictor-corrector method. The exponential time differencing method has been developed and perfected by the present studies.

Calculation method of sliding ratios for conjugate-curve gear pair and its application

The calculation method of sliding ratios for conjugate-curve gear pair, generated based on the theory of conjugate curves,is proposed. The theoretical model of conjugate-curve gear drive is briefly introduced. The general calculation formulas of sliding ratios are developed according to the conjugate curves. The applications to the circular arc gears based on conjugate curves and the novel involute-helix gears are studied. A comparison on the sliding coefficient with the conventional corresponding gear drive is also carried out. The influences of gear parameters such as spiral parameter, gear ratio and modulus on the sliding ratios of gear drive are discussed. Brief description of manufacturing method for conjugate-curve gear pair is given. The research results show that the sliding ratios of gear pair become smaller with the increase of spiral parameter and gear ratio, respectively. And it will be greater with the increase of modulus for the tooth profiles. The meshing characteristics of conjugate-curve gears are further reflected and the optimization design of tooth profiles with high performance may be obtained.

Finite Volume Element Predictor-corrector Method for a Class of Nonlinear Parabolic Systems

A finite volume element predictor-corrector method for a class of nonlinear parabolic system of equations is presented and analyzed. Suboptimal L^2 error estimate for the finite volume element predictor-corrector method is derived. A numerical experiment shows that the numerical results are consistent with theoretical analysis.

A Method for Gear Fault Diagnosis Based on the Empirical Mode Decomposition

According to the characteristics of gear fault vibration signals, a methodfor gear fault diagnosis based upon the empirical mode decomposition (EMD) is proposed in thispaper. By using EMD, any complicated signal can be decomposed into a finite and often small numberof intrinsic mode functions (IMFs) , which are based upon the local characteristic time scale of thesignal. Thus, EMD is perfectly suitable for non-stationary signal processing and faultcharacteristics extracting. It is well known that a gear vibration signal consists of a number offrequency family components, each of which is a modulated signal. Thus, we can use EMD to decomposea gear fault vibration signal into a number of IMF components, some of which correspond to thefrequency families, and the others are noises. Therefore, the frequency families can be separatedand the noise can be decreased at the same time. The proposed method has been applied to gear faultdiagnosis. The results show that both the sensitivity and the reliability of this method aresatisfactory.

An Exact Analytical Method for Magnetic Field Computation and Electromagnetic Torque in a Concentric Magnetic Gear

准确计算磁力齿轮电磁转矩是设计、分析磁力齿轮的关键,采用二维全局解析法计算同心式磁力齿轮气隙磁场。求解场域划分为内外转子永磁体、内外两层气隙和调磁定子的槽形区域,3类子区域的拉普拉斯方程和泊松方程通过边界连续条件建立联系。得到内外两层气隙区域的矢量磁位磁通密度解析表达式,有利于方便、快速、精确地计算任意转子位置的电磁转矩。计算了内外两层气隙磁场和内外转子电磁转矩,将气隙磁场波形和内外转子电磁转矩波形分别与二维有限元法计算波形作比较,结果吻合,证明了方法的正确性和有效性。

Perturbation Theory Combined with Boundary Element Method for Analysis of Gear Contact Problems

This paper combines the perturbation theory with the boundary element methodfor contact problems of three-dimensional elasticity mechanism to analyse the effect oferrors on the shape of the contact area and pressure distribution in gear drive through theperturbation of a cubic order geometry,there by greatly bringing down both computationwork volume and cost and providing a powerful tool for engineering study on the effectof errors on structural strength.

Calculating method of contact stress for non-circular gears

According to Hertz theory, the difference of contact stress for non-circular gears and equivalent gears is compared in the paper, a calculating method of contact stress for non-circular gears by using equivalent gears is researched, and computing formulas of power and rotation speed for equivalent gears are deduced. A numerical simulation of contact stress for non-circular gears has also been conducted based on the finite element method. By the comparison of fitting curves, the feasibility of using equivalent gears instead of non-circular gears to calculate the contact stress is testified.

Chaos control of gear system with elastomeric web based on multi-parameter multi-step method

This paper employs a multi-parameter multi-step chaos control method, which is built up on the OGY method, to stabilize desirable UPOs of a gear system with elastomeric web as a high-dimensional and non-hyperbolic chaotic system, and the analyses are carried out. Three types of relations between components of a certain control parameter combination are defined in a certain control process. Special emphasis is put on the comparison of control efficiencies of the multi-parameter multi-step method and single-parameter multi-step method. The numerical experiments show the ability to switch between different orbits and the method can be a good chaos control alternative since it provides a more effective UPOs stabilization of high-dimensional and non-hyperbolic chaotic systems than the single-parameter chaos control, and according to the relation between components of each parameter combination, the best combination for chaos control in a certain UPO stabilization process are obtained.