Search results
1 – 10 of over 1000Mehdi Dehghan and Masoud Hajarian
Solving the non‐linear equation f(x)=0 has nice applications in various branches of physics and engineering. Sometimes the applications of the numerical methods to solve…
Abstract
Purpose
Solving the non‐linear equation f(x)=0 has nice applications in various branches of physics and engineering. Sometimes the applications of the numerical methods to solve non‐linear equations depending on the second derivatives are restricted in physics and engineering. The purpose of this paper is to propose two new modified Newton's method for solving non‐linear equations. Convergence results show that the order of convergence of the proposed iterative methods for a simple root is four. The iterative methods are free from second derivative and can be used for solving non‐linear equations without computing the second derivative. Finally, several numerical examples are given to illustrate that proposed iterative algorithms are effective.
Design/methodology/approach
In this paper, first the authors introduce two new approximations for the definite integral arising from Newton's theorem. Then by considering these approximations, two new iterative methods are provided with fourth‐order convergence which can be used for solving non‐linear equations without computing second derivatives.
Findings
In this paper, the authors propose two new iterative methods without second derivatives for solving the non‐linear equation f(x)=0. From numerical results, it is observed that the new methods are comparable with various iterative methods. Also numerical results corroborate the theoretical analysis.
Originality/value
The best property of these schemes is that they are second derivative free. Also from numerical results, it is observed that the new methods are comparable with various iterative methods. The numerical results corroborate the theoretical analysis.
Details
Keywords
The initial stiffness method has been extensively adopted for elasto‐plastic finite element analysis. The main problem associated with the initial stiffness method, however, is…
Abstract
Purpose
The initial stiffness method has been extensively adopted for elasto‐plastic finite element analysis. The main problem associated with the initial stiffness method, however, is its slow convergence, even when it is used in conjunction with acceleration techniques. The Newton‐Raphson method has a rapid convergence rate, but its implementation resorts to non‐symmetric linear solvers, and hence the memory requirement may be high. The purpose of this paper is to develop more advanced solution techniques which may overcome the above problems associated with the initial stiffness method and the Newton‐Raphson method.
Design/methodology/approach
In this work, the accelerated symmetric stiffness matrix methods, which cover the accelerated initial stiffness methods as special cases, are proposed for non‐associated plasticity. Within the computational framework for the accelerated symmetric stiffness matrix techniques, some symmetric stiffness matrix candidates are investigated and evaluated.
Findings
Numerical results indicate that for the accelerated symmetric stiffness methods, the elasto‐plastic constitutive matrix, which is constructed by mapping the yield surface of the equivalent material to the plastic potential surface, appears to be appealing. Even when combined with the Krylov iterative solver using a loose convergence criterion, they may still provide good nonlinear convergence rates.
Originality/value
Compared to the work by Sloan et al., the novelty of this study is that a symmetric stiffness matrix is proposed to be used in conjunction with acceleration schemes and it is shown to be more appealing; it is assembled from the elasto‐plastic constitutive matrix by mapping the yield surface of the equivalent material to the plastic potential surface. The advantage of combining the proposed accelerated symmetric stiffness techniques with the Krylov subspace iterative methods for large‐scale applications is also emphasized.
Details
Keywords
The purpose of this paper is to propose three iterative finite element methods for equations of thermally coupled incompressible magneto-hydrodynamics (MHD) on 2D/3D bounded…
Abstract
Purpose
The purpose of this paper is to propose three iterative finite element methods for equations of thermally coupled incompressible magneto-hydrodynamics (MHD) on 2D/3D bounded domain. The detailed theoretical analysis and some numerical results are presented. The main results show that the Stokes iterative method has the strictest restrictions on the physical parameters, and the Newton’s iterative method has the higher accuracy and the Oseen iterative method is stable unconditionally.
Design/methodology/approach
Three iterative finite element methods have been designed for the thermally coupled incompressible MHD flow on 2D/3D bounded domain. The Oseen iterative scheme includes solving a linearized steady MHD and Oseen equations; unconditional stability and optimal error estimates of numerical approximations at each iterative step are established under the uniqueness condition. Stability and convergence of numerical solutions in Newton and Stokes’ iterative schemes are also analyzed under some strong uniqueness conditions.
Findings
This work was supported by the NSF of China (No. 11971152).
Originality/value
This paper presents the best choice for solving the steady thermally coupled MHD equations with different physical parameters.
Details
Keywords
The linear regression technique is widely used to determine empirical parameters of fatigue life profile while the results may not continuously depend on experimental data. Thus…
Abstract
Purpose
The linear regression technique is widely used to determine empirical parameters of fatigue life profile while the results may not continuously depend on experimental data. Thus Tikhonov-Morozov method is utilized here to regularize the linear regression results and consequently reduces the influence of measurement noise without notably distorting the fatigue life distribution. The paper aims to discuss these issues.
Design/methodology/approach
Tikhonov-Morozov regularization method would be shown to effectively reduce the influences of measurement noise without distorting the fatigue life distribution. Moreover since iterative regularization methods are known to be an attractive alternative to Tikhonov regularization, four gradient iterative methods called as simple iteration, minimum error, steepest descent and conjugate gradient methods are examined with an appropriate initial guess of regularized coefficients.
Findings
It has been shown that in case of sparse fatigue life measurements, linear regression results may not have continuous dependence on experimental data and measurement error could lead to misinterpretations of the solution. Therefore from engineering safety point of view, utilizing regularization method could successfully reduce the influence of measurement noise without significantly distorting the fatigue life distribution.
Originality/value
An excellent initial guess for mixed iterative-direct algorithm is introduced and it has been shown that the combination of Newton iterative approach and Morozov discrepancy principle is one of the interesting strategies for determination of regularization parameter having an excellent rate of convergence. Moreover since iterative methods are known to be an attractive alternative to Tikhonov regularization, four gradient descend methods are examined here for regularization of the linear regression problem. It has been found that all of gradient decent methods with an appropriate initial guess of regularized coefficients have an excellent convergence to Tikhonov-Morozov regularization results.
Details
Keywords
In a model resulting from Maxwell's equations with a constitutive law using Preisach operators for incorporating magnetization hysteresis, this paper aims at identifying the…
Abstract
Purpose
In a model resulting from Maxwell's equations with a constitutive law using Preisach operators for incorporating magnetization hysteresis, this paper aims at identifying the hysteresis operator, i.e. the Preisach weight function, from indirect measurements.
Design/methodology/approach
Dealing with a nonlinear inverse problem, one has to apply iterative methods for its numerical solution. For this purpose several approaches are proposed based on fixed point or Newton type ideas. In the latter case, one has to take into account nondifferentiability of the hysteresis operator. This is done by using differentiable substitutes or quasi‐Newton methods.
Findings
Numerical tests with synthetic data show that fixed point methods based on fitting after a full forward sweep (alternating iteration) and Newton type iterations using the hysteresis centerline or commutation curve exhibit a satisfactory convergence behavior, while fixed point iterations based on subdividing the time interval (Kaczmarz) suffer from instability problems and quasi Newton iterations (Broyden) are too slow in some cases.
Research limitations/implications
Application of the proposed methods to measured data will be the subject of future research work.
Practical implications
The proposed methodologies allow to determine material parameters in hysteresis models from indirect measurements.
Originality/value
Taking into account the full PDE model, one can expect to get accurate and reliable results in this model identification problem. Especially the use of Newton type methods – taking into account nondifferentiability – is new in this context.
Details
Keywords
Wenhua Ling and Henryk K. Stolarski
Some frictional contact problems are characterized by significant variations in the location and size of the contact area occurring in the process of deformation. When this…
Abstract
Some frictional contact problems are characterized by significant variations in the location and size of the contact area occurring in the process of deformation. When this feature is combined with strongly non‐linear, path‐dependent material behaviour, difficulties with convergence of the typically used iterative processes can be encountered. Demonstrates this by analysis of press‐fit connection, a typical problem in which both of those characteristics can be present. Offers an explanation as to the possible source of those difficulties. Suggests in support of this explanation, two simple modifications of the usual iterative schemes. In spite of their simplicity, they are found to be more robust than those usual schemes which are normally used in numerical analysis of similar problems.
Details
Keywords
Manpreet Kaur, Sanjeev Kumar and Munish Kansal
The purpose of the article is to construct a new class of higher-order iterative techniques for solving scalar nonlinear problems.
Abstract
Purpose
The purpose of the article is to construct a new class of higher-order iterative techniques for solving scalar nonlinear problems.
Design/methodology/approach
The scheme is generalized by using the power-mean notion. By applying Neville's interpolating technique, the methods are formulated into the derivative-free approaches. Further, to enhance the computational efficiency, the developed iterative methods have been extended to the methods with memory, with the aid of the self-accelerating parameter.
Findings
It is found that the presented family is optimal in terms of Kung and Traub conjecture as it evaluates only five functions in each iteration and attains convergence order sixteen. The proposed family is examined on some practical problems by modeling into nonlinear equations, such as chemical equilibrium problems, beam positioning problems, eigenvalue problems and fractional conversion in a chemical reactor. The obtained results confirm that the developed scheme works more adequately as compared to the existing methods from the literature. Furthermore, the basins of attraction of the different methods have been included to check the convergence in the complex plane.
Originality/value
The presented experiments show that the developed schemes are of great benefit to implement on real-life problems.
Details
Keywords
M.A. Gutierrez, Y.M. Ojanguren and J.J. Anza
The numerical simulation of metal forming processes approximated by means of finite element techniques, require large computational effort, which contradicts the need of…
Abstract
The numerical simulation of metal forming processes approximated by means of finite element techniques, require large computational effort, which contradicts the need of interactivity for industrial applications. This work analyses the computational efficiency of algorithms combining elastoplasticity with finite deformation and contact mechanics, and in particular, the optimum solution of the linear systems to be solved through the incremental‐iterative schemes associated with non linear implicit analysis. A method based on domain decomposition techniques especially adapted to contact problems is presented, as well as the improved performance obtained in the application to hot rolling simulation, as a consequence of bandwidth reduction and the differentiated treatment of subdomains along the non linear analysis.
Details
Keywords
Liang Zhang, Changzhu Wei, Yin Diao and Naigang Cui
This paper aims to investigate the problem of on-line orbit planning and guidance for an advanced upper stage.
Abstract
Purpose
This paper aims to investigate the problem of on-line orbit planning and guidance for an advanced upper stage.
Design/methodology/approach
The double impulse optimal transfer orbit is planned by the Lambert algorithm and the improved particle swarm optimization (IPSO) method, which can reduce the total velocity increment of the transfer orbit. More specially, a simplified formula is developed to obtain the working time of the main engine for two phases of flight based on the theorem of impulse. Subsequently, the true anomalies of the start position and the end position for both two phases are planned by the Newton iterative algorithm and the Kepler equation. Finally, the first phase of flight is guided by a novel iterative guidance (NIG) law based on the true anomaly update with respect to the geometrical relationship. Also, a completely analytical powered explicit guidance (APEG) law is presented to realize orbital injection for the second phase of flight.
Findings
Simulations including Monte Carlo and three typical orbit transfer missions are carried out to demonstrate the efficiency of the proposed scheme.
Originality/value
A novel on-line orbit planning algorithm is developed based on the Lambert problem, IPSO optimization method and Newton iterative algorithm. The NIG and APEG are presented to realize the designed transfer orbit for the first and second phases of flight. Both two guidance laws achieve higher orbit injection accuracies than traditional guidance laws.
Details
Keywords
Fabio De Angelis and Robert L. Taylor
The purpose of this paper is to present an efficient return mapping algorithm for elastoplastic constitutive problems of ductile metals with an exact closed form solution of the…
Abstract
Purpose
The purpose of this paper is to present an efficient return mapping algorithm for elastoplastic constitutive problems of ductile metals with an exact closed form solution of the local constitutive problem in the small strain regime. A Newton Raphson iterative method is adopted for the solution of the boundary value problem.
Design/methodology/approach
An efficient return mapping algorithm is illustrated which is based on an elastic predictor and a plastic corrector scheme resulting in an implicit and accurate numerical integration method. Nonlinear kinematic hardening rules and linear isotropic hardening rules are used to describe the components of the hardening variables. In the adopted algorithmic approach the solution of the local constitutive equations reduces to only one straightforward nonlinear scalar equation.
Findings
The presented algorithmic scheme naturally leads to a particularly simple form of the nonlinear scalar equation which ultimately scales down to an algebraic (polynomial) equation with a single variable. The straightforwardness of the present approach allows to find the analytical solution of the algebraic equation in a closed form. Further, the consistent tangent operator is derived as associated with the proposed algorithmic scheme and it is shown that the proposed computational procedure ensures a quadratic rate of asymptotic convergence when used with a Newton Raphson iterative method for the global solution procedure.
Originality/value
In the present approach the solution of the algebraic nonlinear equation is found in a closed form and accordingly no iterative method is required to solve the problem of the local constitutive equations. The computational procedure ensures a quadratic rate of asymptotic convergence for the global solution procedure typical of computationally efficient solution schemes. In the paper it is shown that the proposed algorithmic scheme provides an efficient and robust computational solution procedure for elastoplasticity boundary value problems. Numerical examples and computational results are reported which illustrate the effectiveness and robustness of the adopted integration algorithm for the finite element analysis of elastoplastic structures also under elaborate loading conditions.
Details