THE STRESS SUBSPACE OF HYBRID STRESS ELEMENT AND THE DIAGONALIZATION METHOD FOR FLEXIBILITY MATRIX H

The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.

BCS Ground State and XXZ Antiferromagnetic Model as SU（2）,SU（1,1） Coherent States：AN Algebraic Diagonalization Method

An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byusing SU(2), SU(1,1) Lie algebraic method, respectively. Meanwhile, the eigenstates of the above two models are revealedto be SU(2), SU(1,1) coherent states, respectively. The relation between the usual Bogoliubov Valatin transformationand the algebraic method in a special case is also discussed.

Decoupling Control Research on Test System of Hydraulic Drive Unit of Quadruped Robot Based on Diagonal Matrix Method

The mathematical model of hydraulic drive unit of quadruped robot was built in this paper. According to the coupling characteristics between position control system and force control system, the decoupling control strategy was realized based on diagonal matrix method in AMESim?. The results of simulation show that using diagonal matrix method can achieve the decoupling control effectively and it can achieve the decoupling control more effectively with the method of not offset pole-zero in the S coordinate. This research can provide theoretical basis for the application of test system of hydraulic drive unit.

BCS Ground State and XXZ Antiferromagnetic Model as SU(2), SU(1,1) Coherent States: An Algebraic Diagonalization Method

An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byusing SU(2), SU(1,1) Lie algebraic method, respectively. Meanwhile, the eigenstates of the above two models are revealedto be SU(2), SU(1,1) coherent states, respectively. The relation between the usual Bogoliubov Valatin transformationand the algebraic method in a special case is also discussed.

DIAGONALIZATION METHOD FOR A SINGULARLY PERTURBED VECTOR ROBIN PROBLEM

In this paper the method and technique of the diagonalization are employed to transform a vector second-order nonlinear system into two first-order approximate diagonalized systems. The existence and the asymptotic behavior of the solutions are obtained for a vector second-order nonlinear Robin problem of singular perturbation type.

3D magnetotelluric inversions with unstructured finite-element and limited-memory quasi-Newton methods

Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.

FPGA-based Acceleration of Davidon-Fletcher-Powell Quasi-Newton Optimization Method

Quasi-Newton methods are the most widely used methods to find local maxima and minima of functions in various engineering practices. However, they involve a large amount of matrix and vector operations, which are computationally intensive and require a long processing time. Recently, with the increasing density and arithmetic cores, field programmable gate array(FPGA) has become an attractive alternative to the acceleration of scientific computation. This paper aims to accelerate Davidon-Fletcher-Powell quasi-Newton(DFP-QN) method by proposing a customized and pipelined hardware implementation on FPGAs. Experimental results demonstrate that compared with a software implementation, a speed-up of up to 17 times can be achieved by the proposed hardware implementation.

Completeness of the system of eigenvectors of off-diagonal operator matrices and its applications in elasticity theory

This paper deals with off-diagonal operator matrices and their applications in elasticity theory. Two kinds of completeness of the system of eigenvectors are proven, in terms of those of the compositions of two block operators in the off-diagonal operator matrices. Using these results, the double eigenfunction expansion method for solving upper triangular matrix differential systems is proposed. Moreover, we apply the method to the two-dimensional elasticity problem and the problem of bending of rectangular thin plates on elastic foundation.

Quasi-Newton Method for Optimal Blank Allowance Balancing

A balancing technique for casting or forging parts to be machined is presented in this paper.It allows an optimal part setup to make sure that no shortage of material(undercut)will occur during machining.Particularly in the heavy part in- dustry,where the resulting casting size and shape may deviate from expectations,the balancing process discovers whether or not the design model is totally enclosed in the actual part to be machined.The alignment is an iterative process involving nonlinear con- strained optimization,which forces data points to lie outside the nominal model under a specific order of priority.Newton methods for non-linear numerical minimization are rarely applied to this problem because of the high cost of computing.In this paper, Newton methods are applied to the balancing of blank part.The aforesaid algorithm is demonstrated in term of a marine propeller blade,and result shows that The Newton methods are more efficient and accurate than those implemented in past research and have distinct advantages compared to the registration methods widely used today.

OPTIMAL MOTION PLANNING FOR A RIGID SPACECRAFT WITH TWO MOMENTUM WHEELS USING QUASI-NEWTON METHOD

An optimal motion planning scheme based on the quasi-Newton method is proposed for a rigid spacecraft with two momentum wheels. A cost functional is introduced to incorporate the control energy, the final state errors and the constraints on states. The motion planning for determining control inputs to minimize the cost functional is formulated as a nonlinear optimal control problem. Using the control parametrization, one can transform the infinite dimensional optimal control problem to a finite dimensional one that is solved via the quasi-Newton methods for a feasible trajectory which satisfies the nonholonomic constraint. The optimal motion planning scheme was applied to a rigid spacecraft with two momentum wheels. The simulation results show the effectiveness of the proposed optimal motion planning scheme.

An Improved Quasi-Newton Method for Unconstrained Optimization

We present an improved method. If we assume that the objective function is twice continuously differentiable and uniformly convex, we discuss global and superlinear convergence of the improved quasi-Newton method.

The Prediction of Non-Life Claim Reserves under Inflation—An Analysis including Diagonal Effects

The extraction of various reserves is one of the most important measures that guarantee insurance companies’ solvency. Accurate assessment of non-life insurance claim reserves needs to consider the volatility risks of inflation. This paper presents a stochastic model of claim reserves including inflation factor and diagonal effects. By applying this model, we can predict the values of the claim reserves and evaluate predicting risks. Through analyzing actual data and using the bootstrap method, we can compare Bornhuetter-Ferguson method involving diagonal effects with chain ladder method. It is shown that the former is more efficient and robust than the latter.

A Study of BCI Signal Pattern Recognition by Using Quasi-Newton-SVM Method

The recognition of electroencephalogram (EEG) signals is the key of brain computer interface (BCI). Aimed at the problem that the recognition rate of EEG by using support vector machine (SVM) is low in BCI, based on the assumption that a well-defined physiological signal which also has a smooth form "hides" inside the noisy EEG signal, a Quasi-Newton-SVM recognition method based on Quasi-Newton method and SVM algorithm was presented. Firstly, the EEG signals were preprocessed by Quasi-Newton method and got the signals which were fit for SVM. Secondly, the preprocessed signals were classified by SVM method. The present simulation results indicated the Quasi-Newton-SVM approach improved the recognition rate compared with using SVM method; we also discussed the relationship between the artificial smooth signals and the classification errors.

BFGS quasi-Newton location algorithm using TDOAs and GROAs

With the emergence of location-based applications in various fields, the higher accuracy of positioning is demanded. By utilizing the time differences of arrival （TDOAs） and gain ratios of arrival （GROAs）, an efficient algorithm for estimating the position is proposed, which exploits the Broyden-Fletcher-Goldfarb-Shanno （BFGS） quasi-Newton method to solve nonlinear equations at the source location under the additive measurement error. Although the accuracy of two-step weighted-least-square （WLS） method based on TDOAs and GROAs is very high, this method has a high computational complexity. While the proposed approach can achieve the same accuracy and bias with the lower computational complexity when the signal-to-noise ratio （SNR） is high, especially it can achieve better accuracy and smaller bias at a lower SNR. The proposed algorithm can be applied to the actual environment due to its real-time property and good robust performance. Simulation results show that with a good initial guess to begin with, the proposed estimator converges to the true solution and achieves the Cramer-Rao lower bound （CRLB） accuracy for both near-field and far-field sources.

Structural Stress-Fatigue Life Curve Improvement of Spot Welding Based on Quasi‑Newton Method

ΔF-N curves are usually used to predict the fatigue life of spot welding in engineering,but they are time-consuming and laborious and not universal.For the purpose of predicting the fatigue life of spot welding accurately and efficiently,tensile-shear fatigue tests were conducted to obtain the fatigue life of spot-welded specimens with different sheet thicknesses combinations.These specimens were simulated by using the finite element method,and the structural stress was theoretically calculated.In the double logarithmic coordinate system,the structural stress-fatigue life(S-N)curve of spot welding was fitted by the least-squares method,based on the quasi-Newton method.The square of the correlation coefficient of the S-N curve was taken as the optimization objective,with the correction coefficients of force,bending moment,spot welding diameter,and sheet thickness as the variables.During the optimization process,three different ways were utilized to get three optimized spot welding S-N curves,which are suitable for different situations.The results show that the fitting effect of the S-N curve is improved,the data points are more compact,and the optimization effect is significant.These S-N curves can be used to predict the fatigue life,which provide the basis for practical engineering application.

Energy Spectrum for a Short-Range 1/r Singular Potential with a Non-Orbital Barrier Using the Asymptotic Iteration Method

Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.

A New Approximation to the Linear Matrix Equation AX = B by Modification of He’s Homotopy Perturbation Method

It is well known that the matrix equations play a significant role in engineering and applicable sciences. In this research article, a new modification of the homotopy perturbation method (HPM) will be proposed to obtain the approximated solution of the matrix equation in the form AX = B. Moreover, the conditions are deduced to check the convergence of the homotopy series. Numerical implementations are adapted to illustrate the properties of the modified method.

A parallel method for numerical solution of delay differential equations

A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the iteration parameters of the method are tuned in such a way that fast convergence to the value of corrector is achieved.

An Efficient Direct Method to Solve the Three Dimensional Poisson’s Equation

In this work, the three dimensional Poisson’s equation in Cartesian coordinates with the Dirichlet’s boundary conditions in a cube is solved directly, by extending the method of Hockney. The Poisson equation is approximated by 19-points and 27-points fourth order finite difference approximation schemes and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The efficiency of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results. It is shown that 19-point formula produces comparable results with 27-point formula, though computational efforts are more in 27-point formula.

Iterative convergence of boundary-volume integral equation method

The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.