A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcom...The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.展开更多
An improved precise integration method (IPIM) for solving the differential Riccati equation (DRE) is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise in...An improved precise integration method (IPIM) for solving the differential Riccati equation (DRE) is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method (PIM) for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the展开更多
The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary diffe...The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.展开更多
In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the met...In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations(ODEs) are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.展开更多
General purpose graphic processing unit (GPU) calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has ...General purpose graphic processing unit (GPU) calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps. In this study, we introduce a GPU-based parallel calculation method of a precise integration method (PIM) for seismic forward modeling. Compared with CPU single-core calculation, GPU parallel calculating perfectly keeps the features of PIM, which has small bandwidth, high accuracy and capability of modeling complex substructures, and GPU calculation brings high computational efficiency, which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.展开更多
This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matr...This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.展开更多
A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in t...A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.展开更多
In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DN...In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DNN)method.The NPIM is employed to modify the exponential matrix and loading vector based on the DNN integral method.This involves incorporating the basic assumption of the Newmark-βmethod into the dynamic equation and eliminating the acceleration term from the dynamic equilibrium equation.As a result,the equation is reduced to a first-order linear equation system.Subsequently,the PIM is applied to integrate the system step by step within the NPIM.The DNN method is adopted to solve the inhomogeneous term through fitting the integrand and the original function with a pair of neural networks,and the integral term is solved using the Newton–Leibniz formula.Numerical examples demonstrate that the proposed methodology significantly improves computing efficiency and provides sufficient precision compared to the DNN method.This is particularly evident when analyzing large-scale structures under blast loading conditions.展开更多
This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analys...This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analysis so that the parameter N, which has been takenas a constant and an independent pa- rameter without consideration ofthe problems in the original method, can be generated automaticallyby the algorithm itself.展开更多
Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix d...Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.展开更多
This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order...This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a "major compo- nent" and a "high order small quantity" and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.展开更多
The vortex-induced vibrations (VIV) of barrel-shaped mast are calculated by three numerical metods, i.e. New- mark- β, HPD-L (High Precision Direct integration scheme-Linear form), and HPD-S (High Precision Dire...The vortex-induced vibrations (VIV) of barrel-shaped mast are calculated by three numerical metods, i.e. New- mark- β, HPD-L (High Precision Direct integration scheme-Linear form), and HPD-S (High Precision Direct integration scheme-Sinusoidal form). According to the measured value, the accuracy curves are given to show the advantages of HPD-S method over others. Based on the comparison above, HPD-S method is used to calculate the influence of ice covering on the mast to VIV responses. It has been proved that the vortex-induced responses of barrel-shaped mast are changed along with ice thicknesses and types.展开更多
Based on the governing equation of vibration of a kind of cylindrical shells written in a matrix differential equation of the first order, a new matrix method is presented for steady-state vibration analysis of a nonc...Based on the governing equation of vibration of a kind of cylindrical shells written in a matrix differential equation of the first order, a new matrix method is presented for steady-state vibration analysis of a noncircular cylindrical shell simply sup- ported at two ends and circumferentially stiffened by rings under harmonic pressure. Its difference from the existing works by Yamada and Irie is that the matrix differential equation is solved by using the extended homogeneous capacity precision integration' approach other than the Runge-Kutta-Gill integration method. The transfer matrix can easily be determined by a high precision integration scheme. In addition, besides the normal interacting forces, which were commonly adopted by researchers earlier, the tangential interacting forces between the cylindrical shell and the rings are considered at the same time by means of the Dirac-δ function. The effects of the exciting frequencies on displacements and stresses responses have been investigated. Numerical results show that the proposed method is more efficient than the aforementioned method.展开更多
The pseudo-viscous frictional energy dissipator(PVFED) is a new energy dissipator. This dissipator can be widely used in engineering for not only the friction is in direct ratio to velocity, but also the problem of ...The pseudo-viscous frictional energy dissipator(PVFED) is a new energy dissipator. This dissipator can be widely used in engineering for not only the friction is in direct ratio to velocity, but also the problem of viscous energy dissipator mucilage easily leaked has been overcome. The problem of how to get response of the PVFED sys- tem need to be solved before this dissipator can be used widely in engineering. The response calculation methods of the PVFED system on sina load was researched. Wilson-θ,Newmark-β and a precise integration algorithm was used separately to solve the system response and the calculation result in a different time step was compared. It was found from comparison that three calculation results were almost equivalent in a small time step. Calculation precision of Newmark-β and Wilson-θ was reduced and high calculation precision of a precise integration algorithm was kept in a large time step. The results show that it is an effective way to solve the response of a PVFED system by a precise integration method.展开更多
Based on the nonlocal continuum theory, the nonlinear vibration of an embedded single-walled carbon nanotube (SWCNT) subjected to a harmonic load is in- vestigated. In the present study, the SWCNT is assumed to be a...Based on the nonlocal continuum theory, the nonlinear vibration of an embedded single-walled carbon nanotube (SWCNT) subjected to a harmonic load is in- vestigated. In the present study, the SWCNT is assumed to be a curved beam, which is unlike previous similar work. Firstly, the governing equations of motion are derived by the Hamilton principle, meanwhile, the Galerkin approach is carried out to convert the nonlinear integral-differential equation into a second-order nonlinear ordinary differ- ential equation. Then, the precise integration method based on the local linearzation is appropriately designed for solving the above dynamic equations. Besides, the numerical example is presented, the effects of the nonlocal parameters, the elastic medium constants, the waviness ratios, and the material lengths on the dynamic response are analyzed. The results show that the above mentioned effects have influences on the dynamic behavior of the SWCNT.展开更多
The Non-uniform rational B-spline (NURBS) enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this pap...The Non-uniform rational B-spline (NURBS) enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.展开更多
The governing equations of elasticity theory for natural vibration and buck- ling of anisotropic plate are derived from Hellinger-Reissner's variational principle with nonlinear strain-displacement relations. Simply ...The governing equations of elasticity theory for natural vibration and buck- ling of anisotropic plate are derived from Hellinger-Reissner's variational principle with nonlinear strain-displacement relations. Simply supported rectangular hybrid plates are studied with a precise integration method. This method, in contrast to the traditional finite difference approximation, gives highly precise numerical results that approach the full computer precision. So the results for natural vibration and stability of hybrid plates presented in the paper can be riewed as approximate analytical solutions. Furthermore, several types of coupling effects such as coupling between bending and twisting, and coupling between extension and bending, when the layer stacking sequence is asymmetric, are considered by only one set of governing equations.展开更多
This paper presents a new type of base isolation system, i. e. , slide-limited friction (S-LF) base isolation system . Based on this system, the harmonic and subharmonic periodic response of S-LF subjected to harmonic...This paper presents a new type of base isolation system, i. e. , slide-limited friction (S-LF) base isolation system . Based on this system, the harmonic and subharmonic periodic response of S-LF subjected to harmonic motions is investigated by using Fourier-Galerkin-Newton (FGN) method with Flo-quet theory. The dynamic response of S-LF subjected to earthquake ground motions is calculated with a high order precision direct integration method, and the numerical results are presented in maximum acceleration response spectra of superstructure and maximum sliding displacement response spectrum form. The comparison of isolating effects of S-LF, pure-friction base isolation system (P-F) and resilient-friction base isolation system (R-FBI) shows that the isolating property of S-LF is superior to those of P-F and R-FBI. Finally, by analyzing an engineering example, it is observed that the distribution of the maximum shear between floors and absolute acceleration of S-LF to earthquake ground motion is very different from that of traditional structures.展开更多
Based on the classical theory of thin plate and Biot theory, a precise model of the transverse vibrations of a thin rectangular porous plate is proposed. The first order differential equations of the porous plate are ...Based on the classical theory of thin plate and Biot theory, a precise model of the transverse vibrations of a thin rectangular porous plate is proposed. The first order differential equations of the porous plate are derived in the frequency domain. By considering the coupling effect between the solid phase and the fluid phase and without any hypothesis for the fluid displacement, the model presented here is rigorous and close to the real materials. Owing to the use of extended homogeneous capacity precision integration method and precise element method, the model can be applied in higher frequency range than pure numerical methods. This model also easily adapts to various boundary conditions. Numerical results are given for two different porous plates under different excitations and boundary conditions.展开更多
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
文摘A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
文摘The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method (IPIM) for solving the differential Riccati equation (DRE) is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method (PIM) for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
文摘The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.
基金supported by the National Natural Science Foundation of China (11132004 and 51078145)the Natural Science Foundation of Guangdong Province (9251064101000016)
文摘In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations(ODEs) are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.
基金supported by the National Natural Science Foundation of China (Nos 40974066 and 40821062)National Basic Research Program of China (No 2007CB209602)
文摘General purpose graphic processing unit (GPU) calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps. In this study, we introduce a GPU-based parallel calculation method of a precise integration method (PIM) for seismic forward modeling. Compared with CPU single-core calculation, GPU parallel calculating perfectly keeps the features of PIM, which has small bandwidth, high accuracy and capability of modeling complex substructures, and GPU calculation brings high computational efficiency, which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.
基金Project supported by the National Natural Science Foundation of China(No.10672194)the China-Russia Cooperative Project(the National Natural Science Foundation of China and the Russian Foundation for Basic Research)(No.10811120012)
文摘This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.
基金the National Natural Science Foundation of China(No.11672064)。
文摘A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.
基金supported by the National Natural Science Foundation of China(Grant Nos.12072288,U2241274,and 12272319).
文摘In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DNN)method.The NPIM is employed to modify the exponential matrix and loading vector based on the DNN integral method.This involves incorporating the basic assumption of the Newmark-βmethod into the dynamic equation and eliminating the acceleration term from the dynamic equilibrium equation.As a result,the equation is reduced to a first-order linear equation system.Subsequently,the PIM is applied to integrate the system step by step within the NPIM.The DNN method is adopted to solve the inhomogeneous term through fitting the integrand and the original function with a pair of neural networks,and the integral term is solved using the Newton–Leibniz formula.Numerical examples demonstrate that the proposed methodology significantly improves computing efficiency and provides sufficient precision compared to the DNN method.This is particularly evident when analyzing large-scale structures under blast loading conditions.
基金the National Natural Science Foundation of China (No.19872016,19872017)the National Key Basic Research Special Foundation (G1999032805)the Foundation for University Key Teachers by the Ministry of Education of China
文摘This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analysis so that the parameter N, which has been takenas a constant and an independent pa- rameter without consideration ofthe problems in the original method, can be generated automaticallyby the algorithm itself.
基金Project supported by the National Natural Science Foundation of China(Nos.60273048and60174023).
文摘Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.
基金supported by the National Natural Science Foun-dation of China (11172334)
文摘This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a "major compo- nent" and a "high order small quantity" and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.
文摘The vortex-induced vibrations (VIV) of barrel-shaped mast are calculated by three numerical metods, i.e. New- mark- β, HPD-L (High Precision Direct integration scheme-Linear form), and HPD-S (High Precision Direct integration scheme-Sinusoidal form). According to the measured value, the accuracy curves are given to show the advantages of HPD-S method over others. Based on the comparison above, HPD-S method is used to calculate the influence of ice covering on the mast to VIV responses. It has been proved that the vortex-induced responses of barrel-shaped mast are changed along with ice thicknesses and types.
基金Project supported by the Doctoral Foundation of the National Education Ministry of China(No.20040487013)
文摘Based on the governing equation of vibration of a kind of cylindrical shells written in a matrix differential equation of the first order, a new matrix method is presented for steady-state vibration analysis of a noncircular cylindrical shell simply sup- ported at two ends and circumferentially stiffened by rings under harmonic pressure. Its difference from the existing works by Yamada and Irie is that the matrix differential equation is solved by using the extended homogeneous capacity precision integration' approach other than the Runge-Kutta-Gill integration method. The transfer matrix can easily be determined by a high precision integration scheme. In addition, besides the normal interacting forces, which were commonly adopted by researchers earlier, the tangential interacting forces between the cylindrical shell and the rings are considered at the same time by means of the Dirac-δ function. The effects of the exciting frequencies on displacements and stresses responses have been investigated. Numerical results show that the proposed method is more efficient than the aforementioned method.
文摘The pseudo-viscous frictional energy dissipator(PVFED) is a new energy dissipator. This dissipator can be widely used in engineering for not only the friction is in direct ratio to velocity, but also the problem of viscous energy dissipator mucilage easily leaked has been overcome. The problem of how to get response of the PVFED sys- tem need to be solved before this dissipator can be used widely in engineering. The response calculation methods of the PVFED system on sina load was researched. Wilson-θ,Newmark-β and a precise integration algorithm was used separately to solve the system response and the calculation result in a different time step was compared. It was found from comparison that three calculation results were almost equivalent in a small time step. Calculation precision of Newmark-β and Wilson-θ was reduced and high calculation precision of a precise integration algorithm was kept in a large time step. The results show that it is an effective way to solve the response of a PVFED system by a precise integration method.
基金Project supported by the National Basic Research Program of China (No. 2011CB610300)the National Natural Science Foundation of China (Nos. 10972182, 11172239, and 10902089)+3 种基金the 111 Project of China (No. B07050)the Ph. D. Programs Foundation of Ministry of Education of China (No. 20106102110019)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment (No. GZ0802)the Doctorate Foundation of Northwestern Polytechnical University (No. CX201224)
文摘Based on the nonlocal continuum theory, the nonlinear vibration of an embedded single-walled carbon nanotube (SWCNT) subjected to a harmonic load is in- vestigated. In the present study, the SWCNT is assumed to be a curved beam, which is unlike previous similar work. Firstly, the governing equations of motion are derived by the Hamilton principle, meanwhile, the Galerkin approach is carried out to convert the nonlinear integral-differential equation into a second-order nonlinear ordinary differ- ential equation. Then, the precise integration method based on the local linearzation is appropriately designed for solving the above dynamic equations. Besides, the numerical example is presented, the effects of the nonlocal parameters, the elastic medium constants, the waviness ratios, and the material lengths on the dynamic response are analyzed. The results show that the above mentioned effects have influences on the dynamic behavior of the SWCNT.
基金support by the National Natural Science Foundation of China(grant No.51779033,51409038)the National Key Research and Development Plan(grant No.2016YFB0201001)the National Natural Science Foundation of China(grant No.51421064)
文摘The Non-uniform rational B-spline (NURBS) enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper. The scaled boundary finite element method is a semi-analytical technique, which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction. In this method, only the boundary is discretized in the finite element sense leading to a re- duction of the spatial dimension by one with no fundamental solution required. Neverthe- less, in case of the complex geometry, a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often un- avoidable in the conventional finite element approach, which leads to huge computational efforts and loss of accuracy. NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape. In the proposed methodology, the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions, while the straight part of the boundary is discretized by the conventional Lagrange shape functions. Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analy- sis and the solution is obtained using the modified precise integration method. The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion. Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method. The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.
文摘The governing equations of elasticity theory for natural vibration and buck- ling of anisotropic plate are derived from Hellinger-Reissner's variational principle with nonlinear strain-displacement relations. Simply supported rectangular hybrid plates are studied with a precise integration method. This method, in contrast to the traditional finite difference approximation, gives highly precise numerical results that approach the full computer precision. So the results for natural vibration and stability of hybrid plates presented in the paper can be riewed as approximate analytical solutions. Furthermore, several types of coupling effects such as coupling between bending and twisting, and coupling between extension and bending, when the layer stacking sequence is asymmetric, are considered by only one set of governing equations.
文摘This paper presents a new type of base isolation system, i. e. , slide-limited friction (S-LF) base isolation system . Based on this system, the harmonic and subharmonic periodic response of S-LF subjected to harmonic motions is investigated by using Fourier-Galerkin-Newton (FGN) method with Flo-quet theory. The dynamic response of S-LF subjected to earthquake ground motions is calculated with a high order precision direct integration method, and the numerical results are presented in maximum acceleration response spectra of superstructure and maximum sliding displacement response spectrum form. The comparison of isolating effects of S-LF, pure-friction base isolation system (P-F) and resilient-friction base isolation system (R-FBI) shows that the isolating property of S-LF is superior to those of P-F and R-FBI. Finally, by analyzing an engineering example, it is observed that the distribution of the maximum shear between floors and absolute acceleration of S-LF to earthquake ground motion is very different from that of traditional structures.
基金Project supported by the National Natural Science Foundation of China(nos.11162001,11502056 and 51665006)
文摘Based on the classical theory of thin plate and Biot theory, a precise model of the transverse vibrations of a thin rectangular porous plate is proposed. The first order differential equations of the porous plate are derived in the frequency domain. By considering the coupling effect between the solid phase and the fluid phase and without any hypothesis for the fluid displacement, the model presented here is rigorous and close to the real materials. Owing to the use of extended homogeneous capacity precision integration method and precise element method, the model can be applied in higher frequency range than pure numerical methods. This model also easily adapts to various boundary conditions. Numerical results are given for two different porous plates under different excitations and boundary conditions.