To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection techniq...To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linea...How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).展开更多
A composite beam is symmetric if both the material property and support are symmetric with respect to the middle point. In order to study the free vibration performance of the symmetric composite beams with different ...A composite beam is symmetric if both the material property and support are symmetric with respect to the middle point. In order to study the free vibration performance of the symmetric composite beams with different complex nonsmooth/discontinuous interfaces, we develop an R(x)-orthonormal theory, where R(x) is an integrable flexural rigidity function. The R(x)-orthonormal bases in the linear space of boundary functions are constructed, of which the second-order derivatives of the boundary functions are asked to be orthonormal with respect to the weight function R(x). When the vibration modes of the symmetric composite beam are expressed in terms of the R(x)-orthonormal bases we can derive an eigenvalue problem endowed with a special structure of the coefficient matrix A :=[aij ],aij= 0 if i + j is odd. Based on the special structure we can prove two new theorems, which indicate that the characteristic equation of A can be decomposed into the product of the characteristic equations of two sub-matrices with dimensions half lower. Hence, we can sequentially solve the natural frequencies in closed-form owing to the specialty of A. We use this powerful new theory to analyze the free vibration performance and the vibration modes of symmetric composite beams with three different interfaces.展开更多
A double optimal solution of an n-dimensional system of linear equations Ax=b has been derived in an affine m-dimensional Krylov subspace with m <<n.We further develop a double optimal iterative algorithm(DOIA),...A double optimal solution of an n-dimensional system of linear equations Ax=b has been derived in an affine m-dimensional Krylov subspace with m <<n.We further develop a double optimal iterative algorithm(DOIA),with the descent direction z being solved from the residual equation Az=r0 by using its double optimal solution,to solve ill-posed linear problem under large noise.The DOIA is proven to be absolutely convergent step-by-step with the square residual error ||r||^2=||b-Ax||^2 being reduced by a positive quantity ||Azk||^2 at each iteration step,which is found to be better than those algorithms based on the minimization of the square residual error in an m-dimensional Krylov subspace.In order to tackle the ill-posed linear problem under a large noise,we also propose a novel double optimal regularization algorithm(DORA)to solve it,which is an improvement of the Tikhonov regularization method.Some numerical tests reveal the high performance of DOIA and DORA against large noise.These methods are of use in the ill-posed problems of structural health-monitoring.展开更多
In order to recover unknown space-dependent function G(x)or unknown time-dependent function H(t)in the wave source F(x;t)=G(x)H(t),we develop a technique of homogenized function and differencing equations,which can si...In order to recover unknown space-dependent function G(x)or unknown time-dependent function H(t)in the wave source F(x;t)=G(x)H(t),we develop a technique of homogenized function and differencing equations,which can significantly reduce the difficulty in the inverse wave source recovery problem,only needing to solve a few equations in the problem domain,since the initial condition/boundary conditions and a supplementary final time condition are satisfied automatically.As a consequence,the eigenfunctions are used to expand the trial solutions,and then a small scale linear system is solved to determine the expansion coefficients from the differencing equations.Because the ill-posedness of the inverse wave source problem is greatly reduced,the present method is accurate and stable against a large noise up to 50%,of which the numerical tests confirm the observation.展开更多
Earthquake is a violent and irregular ground motion that can severely damage structures. In this paper we subject a single-degree-of-freedom system, consisting of spring and damper, to an earthquake excitation, and me...Earthquake is a violent and irregular ground motion that can severely damage structures. In this paper we subject a single-degree-of-freedom system, consisting of spring and damper, to an earthquake excitation, and meanwhile investigate the response behavior from a novel theory about the dynamical system, by viewing the time-varying signum function of It can reflect the characteristic property of each earthquake through and the second component of f, where is a time-sampling record of the acceleration of a ground motion. The barcode is formed by plotting with respect to time. We analyze the complex jumping behavior in a barcode and an essential property of a high percentage occupation of the first set of dis-connectivity in the barcode from four strong earthquake records: 1940 El Centro earthquake, 1989 Loma earthquake, and two records of 1999 Chi-Chi earthquake. Through the comparisons of four earthquakes, we can observe that strong earthquake leads to large percentage of the first set of dis-connectivity.展开更多
基金supported by the the National Science and Technology Council(Grant Number:NSTC 112-2221-E239-022).
文摘To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
基金support provided by the Ministry of Science and Technology,Taiwan,ROC under Contract No.MOST 110-2221-E-019-044.
文摘How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).
文摘A composite beam is symmetric if both the material property and support are symmetric with respect to the middle point. In order to study the free vibration performance of the symmetric composite beams with different complex nonsmooth/discontinuous interfaces, we develop an R(x)-orthonormal theory, where R(x) is an integrable flexural rigidity function. The R(x)-orthonormal bases in the linear space of boundary functions are constructed, of which the second-order derivatives of the boundary functions are asked to be orthonormal with respect to the weight function R(x). When the vibration modes of the symmetric composite beam are expressed in terms of the R(x)-orthonormal bases we can derive an eigenvalue problem endowed with a special structure of the coefficient matrix A :=[aij ],aij= 0 if i + j is odd. Based on the special structure we can prove two new theorems, which indicate that the characteristic equation of A can be decomposed into the product of the characteristic equations of two sub-matrices with dimensions half lower. Hence, we can sequentially solve the natural frequencies in closed-form owing to the specialty of A. We use this powerful new theory to analyze the free vibration performance and the vibration modes of symmetric composite beams with three different interfaces.
文摘A double optimal solution of an n-dimensional system of linear equations Ax=b has been derived in an affine m-dimensional Krylov subspace with m <<n.We further develop a double optimal iterative algorithm(DOIA),with the descent direction z being solved from the residual equation Az=r0 by using its double optimal solution,to solve ill-posed linear problem under large noise.The DOIA is proven to be absolutely convergent step-by-step with the square residual error ||r||^2=||b-Ax||^2 being reduced by a positive quantity ||Azk||^2 at each iteration step,which is found to be better than those algorithms based on the minimization of the square residual error in an m-dimensional Krylov subspace.In order to tackle the ill-posed linear problem under a large noise,we also propose a novel double optimal regularization algorithm(DORA)to solve it,which is an improvement of the Tikhonov regularization method.Some numerical tests reveal the high performance of DOIA and DORA against large noise.These methods are of use in the ill-posed problems of structural health-monitoring.
文摘In order to recover unknown space-dependent function G(x)or unknown time-dependent function H(t)in the wave source F(x;t)=G(x)H(t),we develop a technique of homogenized function and differencing equations,which can significantly reduce the difficulty in the inverse wave source recovery problem,only needing to solve a few equations in the problem domain,since the initial condition/boundary conditions and a supplementary final time condition are satisfied automatically.As a consequence,the eigenfunctions are used to expand the trial solutions,and then a small scale linear system is solved to determine the expansion coefficients from the differencing equations.Because the ill-posedness of the inverse wave source problem is greatly reduced,the present method is accurate and stable against a large noise up to 50%,of which the numerical tests confirm the observation.
文摘Earthquake is a violent and irregular ground motion that can severely damage structures. In this paper we subject a single-degree-of-freedom system, consisting of spring and damper, to an earthquake excitation, and meanwhile investigate the response behavior from a novel theory about the dynamical system, by viewing the time-varying signum function of It can reflect the characteristic property of each earthquake through and the second component of f, where is a time-sampling record of the acceleration of a ground motion. The barcode is formed by plotting with respect to time. We analyze the complex jumping behavior in a barcode and an essential property of a high percentage occupation of the first set of dis-connectivity in the barcode from four strong earthquake records: 1940 El Centro earthquake, 1989 Loma earthquake, and two records of 1999 Chi-Chi earthquake. Through the comparisons of four earthquakes, we can observe that strong earthquake leads to large percentage of the first set of dis-connectivity.
