The current theory in NF EN 1995-1-1/NA of Eurocode 5, which is based on maximum deflection, has been investigated on softwoods. Therefore, this theory is not adapted for slender glulam beam columns made of tropical h...The current theory in NF EN 1995-1-1/NA of Eurocode 5, which is based on maximum deflection, has been investigated on softwoods. Therefore, this theory is not adapted for slender glulam beam columns made of tropical hardwood species from the Congo Basin. This maximum deflection is caused by a set of loads applied to the structure. However, Eurocode 5 doesn’t provide how to predict this deflection in case of long-term load for such structures. This can be done by studying load-displacement (P-Δ) behaviour of these structures while taking into account second order effects. To reach this goal, a nonlinear analysis has been performed on a three-dimensional beam column embedded on both ends. Since conducting experimental investigations on large span structural products is time-consuming and expensive especially in developing countries, a numerical model has been implemented using the Newton-Raphson method to predict load-displacement (P-Δ) curve on a slender glulam beam column made of tropical hardwood species. On one hand, the beam has been analyzed without wood connection. On the other hand, the beam has been analyzed with a bolted wood connection and a slotted-in steel plate. The load cases considered include self-weight and a uniformly applied long-term load. Combinations of serviceability limit states (SLS) and ultimate limit states (ULS) have also been considered, among other factors. A finite-element software RFEM 5 has been used to implement the model. The results showed that the use of steel can reduce displacement by 20.96%. Additionally, compared to the maximum deflection provided by Eurocode 5 for softwoods, hardwoods can exhibit an increasing rate of 85.63%. By harnessing the plastic resistance of steel, the bending resistance of wood can be increased by 32.94%.展开更多
An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small phys...An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.展开更多
Based on the nonlinear constitutive equation,a piezoelectric semiconductor(PSC)fiber under axial loads and Ohmic contact boundary conditions is investigated.The analytical solutions of electromechanical fields are der...Based on the nonlinear constitutive equation,a piezoelectric semiconductor(PSC)fiber under axial loads and Ohmic contact boundary conditions is investigated.The analytical solutions of electromechanical fields are derived by the homotopy analysis method(HAM),indicating that the HAM is efficient for the nonlinear analysis of PSC fibers,along with a rapid rate of convergence.Furthermore,the nonlinear characteristics of electromechanical fields are discussed through numerical results.It is shown that the asymmetrical distribution of electromechanical fields is obvious under a symmetrical load,and the piezoelectric effect is weakened by an applied electric field.With the increase in the initial carrier concentration,the electric potential decreases,and owing to the screen-ing effect of electrons,the distribution of electromechanical fields tends to be symmetrical.展开更多
The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elas...The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.展开更多
The offshore reinforced concrete structures are always subject to cyclic load, such as wave load.In this paper a new finite element analysis model is developed to analyze the stress and strain state of reinforced conc...The offshore reinforced concrete structures are always subject to cyclic load, such as wave load.In this paper a new finite element analysis model is developed to analyze the stress and strain state of reinforced concrete structures including offshore concrete structures, subject to any number of the cyclic load. On the basis of the anal ysis of the experimental data,this model simplifies the number of cycles-total cyclic strain curve of concrete as three straight line segments,and it is assumed that the stress-strain curves of different cycles in each segment are the same, thus the elastoplastic analysis is only needed for the first cycle of each segment, and the stress or strain corresponding to any number of cycles can be obtained by superposition of stress or strain obtained by the above e lastoplastic analysis based on the cyclic numbers in each segment.This model spends less computer time,and can obtain the stress and strain states of the structures after any number of cycles.The endochronic-damage and ideal offshore concrete platform subject to cyclic loading are experimented and analyzed by the finite element method based on the model proposed in this paper. The results between the experiment and the finite element analysis are in good agreement,which demonstrates the validity and accuracy of the proposed model.展开更多
The nonlinear analysis of reinforced concrete rectangular slabs undermonotonic transverse loads is performed by finite element method.The layered rectangu-lar element with 4 nodes and 20 degrees of freedom is develope...The nonlinear analysis of reinforced concrete rectangular slabs undermonotonic transverse loads is performed by finite element method.The layered rectangu-lar element with 4 nodes and 20 degrees of freedom is developed,in whichbending-stretching coupling effect is taken into account.An orthotropic equivalentuniaxial stress-strain constitutive model of concrete is used.A program is worked out andused to calculate two reinforced concrete slabs.The results of calculation are in goodconformity with the corresponding test results.In addition,the influence of tension stif-fening effect of cracked concrete on the results of calculation is discussed.展开更多
This paper first analyzes the features of two classes of numerical methods for global analysis of nonlinear dynamical systems, which regard state space respectively as continuous and discrete ones. On basis of this un...This paper first analyzes the features of two classes of numerical methods for global analysis of nonlinear dynamical systems, which regard state space respectively as continuous and discrete ones. On basis of this understanding it then points out that the previously proposed method of point mapping under cell reference (PMUCR), has laid a frame work for the development of a two scaled numerical method suitable for the global analysis of high dimensional nonlinear systems, which may take the advantages of both classes of single scaled methods but will release the difficulties induced by the disadvantages of them. The basic ideas and main steps of implementation of the two scaled method, namely extended PMUCR, are elaborated. Finally, two examples are presented to demonstrate the capabilities of the proposed method.展开更多
A high order single step β algorithm, a new direct integration algorithm is proposed for solution of equations of motion. Whenβ=0.5, the accuracy of displacement, velocity and acceleration is of forth order (a trunc...A high order single step β algorithm, a new direct integration algorithm is proposed for solution of equations of motion. Whenβ=0.5, the accuracy of displacement, velocity and acceleration is of forth order (a truncation error of Δ t 5), and the algorithm is unconditionally stable and has no arithmetic damping and no overshooting. When >0.5, and an arithmetic damping is adopted, the algorithm is again unconditionally stable with a third order accuracy (a truncation error of Δ t 4). The analyses run with typical examples show that the algorithm proposed has higher speed, higher precision and better properties than other direct integration methods, such as Wilson θ method and Newmark β method in analysing linear elastic responses and nonlinear earthquake responses.展开更多
The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quad...The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.展开更多
In this paper, the homotopy analysis method is applied to deduce the periodic solutions of a conservative nonlinear oscillator for which the elastic force term is proportional to<em> u</em><sup>1/3&l...In this paper, the homotopy analysis method is applied to deduce the periodic solutions of a conservative nonlinear oscillator for which the elastic force term is proportional to<em> u</em><sup>1/3</sup>. By introducing the auxiliary linear operator and the initial guess of solution, the homotopy analysis solving is set up. By choosing the suitable convergence-control parameter, the accurate high-order approximations of solution and frequency for the whole range of initial amplitudes can easily be obtained. Comparison of the results obtained using this method with those obtained by different methods reveals that the former is more accurate, effective and convenient for these types of nonlinear oscillators.展开更多
In this study, the homotopy analysis method (HAM) is used to solve the generalized Duffing equation. Both the frequencies and periodic solutions of the nonlinear Duffing equation can be explicitly and analytically for...In this study, the homotopy analysis method (HAM) is used to solve the generalized Duffing equation. Both the frequencies and periodic solutions of the nonlinear Duffing equation can be explicitly and analytically formulated. Accuracy and validity of the proposed techniques are then verified by comparing the numerical results obtained based on the HAM and numerical integration method. Numerical simulations are extended for even very strong nonlinearities and very good correlations which achieved between the results. Besides, the optimal HAM approach is introduced to accelerate the convergence of solutions.展开更多
Based on the upper bound limit analysis theorem and the shear strength reduction technique, the equation for expressing critical limit-equilibrium state was employed to define the safety factor of a given slope and it...Based on the upper bound limit analysis theorem and the shear strength reduction technique, the equation for expressing critical limit-equilibrium state was employed to define the safety factor of a given slope and its corresponding critical failure mechanism by means of the kinematical approach of limit analysis theory. The nonlinear shear strength parameters were treated as variable parameters and a kinematically admissible failure mechanism was considered for calculation schemes. The iterative optimization method was adopted to obtain the safety factors. Case study and comparative analysis show that solutions presented here agree with available predictions when nonlinear criterion reduces to linear criterion, and the validity of present method could be illuminated. From the numerical results, it can also be seen that nonlinear parameter rn, slope foot gradient ,β, height of slope H, slope top gradient a and soil bulk density γ have significant effects on the safety factor of the slope.展开更多
We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on t...We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on the Laplace trans- form with the homotopy analysis method （HAM）. This method is a powerful tool for solving a large amount of problems. This technique provides a series of functions which may converge to the exact solution of the problem. A good agreement between the obtained solution and some well-known results is obtained.展开更多
This paper introduces an orthogonal expansion method for general stochastic processes. In the method, a normalized orthogonal function of time variable t is first introduced to carry out the decomposition of a stochas...This paper introduces an orthogonal expansion method for general stochastic processes. In the method, a normalized orthogonal function of time variable t is first introduced to carry out the decomposition of a stochastic process and then a correlated matrix decomposition technique, which transforms a correlated random vector into a vector of standard uncorrelated random variables, is used to complete a double orthogonal decomposition of the stochastic processes. Considering the relationship between the Hartley transform and Fourier transform of a real-valued function, it is suggested that the first orthogonal expansion in the above process is carried out using the Hartley basis function instead of the trigonometric basis function in practical applications. The seismic ground motion is investigated using the above method. In order to capture the main probabilistic characteristics of the seismic ground motion, it is proposed to directly carry out the orthogonal expansion of the seismic displacements. The case study shows that the proposed method is feasible to represent the seismic ground motion with only a few random variables. In the second part of the paper, the probability density evolution method （PDEM） is employed to study the stochastic response of nonlinear structures subjected to earthquake excitations. In the PDEM, a completely uncoupled one-dimensional partial differential equation, the generalized density evolution equation, plays a central role in governing the stochastic seismic responses of the nonlinear structure. The solution to this equation will yield the instantaneous probability density function of the responses. Computational algorithms to solve the probability density evolution equation are described. An example, which deals with a nonlinear frame structure subjected to stochastic ground motions, is illustrated to validate the above approach.展开更多
Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for...Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for high tailings dam was developed according to the results of engineering geological exploration, static and dynamic test and stability analysis on Baizhishan tailing dam 113.5 m high. The law of generation, diffusion and dissipation of seismic pore water pressure during and after earthquake was investigated, and the results of tailings dam’s acceleration, seismic dynamic stress and pore water pressure were obtained. The results show that the seismic stability and liquefaction resistance of high tailings dam are strengthened remarkably, and the scope and depth of liquefaction area at the top of dam are reduced greatly. The interior stress is compressive stress, the stress level of every element is less than 1.0 and the safety coefficient of every element is greater than 1.0. The safety coefficient against liquefaction of every element of tailing dam is greater than 1.5 according to the effective stress analysis of seismic response by finite element method. The calculated results prove that liquefaction is the main reason of seismic failure of high tailing dams, and the effect of seismic inertia forces on high tailing dams’ stability during earthquake is secondary reason.展开更多
A homotopy analysis method（HAM）is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existenc...A homotopy analysis method（HAM）is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existence of small parameters in the considered equation.The HAM provides a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter.Two examples are presented to show that the HAM solutions agree well with the results of the modified Linstedt-Poincar＇e method and the incremental harmonic balance method.展开更多
Soft nonlinear support is a major engineering project,but there are few relevant studies.In this paper,a dynamic pipeline model with soft nonlinear supports at both ends is established.By considering the influence of ...Soft nonlinear support is a major engineering project,but there are few relevant studies.In this paper,a dynamic pipeline model with soft nonlinear supports at both ends is established.By considering the influence of the Coriolis force and centrifugal force,the dynamical coupling equation of fluid-structure interaction is derived with extended Hamilton’s principle.Then,the approximate analytical solutions are sought via the harmonic balance method.The amplitude-frequency response curves show that different effects can be determined by approximate analysis.It is demonstrated that the increase in the fluid velocity can increase the amplitude of the pipeline system.The frequency range of unstable response increases when the fluid pressure raises.The combination of the soft nonlinear clamp and the large geometrical deformation of the pipeline affects the nonlinear vibration characteristic of the system,and the external excitation force and damping have significant effects on the stability.展开更多
Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic s...Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.展开更多
Under the hypotheses that the second-order and third-order derivatives of a function are bounded, an estimate of the radius of the convergence ball of Ostrowski-Traub＇s method is obtained. An error analysis is given ...Under the hypotheses that the second-order and third-order derivatives of a function are bounded, an estimate of the radius of the convergence ball of Ostrowski-Traub＇s method is obtained. An error analysis is given which matches the convergence order of the method. Finally, two examples are provided to show applications of our theorem.展开更多
Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ...Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ordinary solution techniques lead to instability near the limit points and also have problems in case of snap-through and snap-back. Thus they fail to predict the complete load-displacement response. The arc-length method serves the purpose well in principle, received wide acceptance in finite element analysis, and has been used extensively. However modifications to the basic idea are vital to meet the particular needs of the analysis. This paper reviews some of the recent developments of the method in the last two decades, with particular emphasis on nonlinear finite element analysis of reinforced concrete structures.展开更多
文摘The current theory in NF EN 1995-1-1/NA of Eurocode 5, which is based on maximum deflection, has been investigated on softwoods. Therefore, this theory is not adapted for slender glulam beam columns made of tropical hardwood species from the Congo Basin. This maximum deflection is caused by a set of loads applied to the structure. However, Eurocode 5 doesn’t provide how to predict this deflection in case of long-term load for such structures. This can be done by studying load-displacement (P-Δ) behaviour of these structures while taking into account second order effects. To reach this goal, a nonlinear analysis has been performed on a three-dimensional beam column embedded on both ends. Since conducting experimental investigations on large span structural products is time-consuming and expensive especially in developing countries, a numerical model has been implemented using the Newton-Raphson method to predict load-displacement (P-Δ) curve on a slender glulam beam column made of tropical hardwood species. On one hand, the beam has been analyzed without wood connection. On the other hand, the beam has been analyzed with a bolted wood connection and a slotted-in steel plate. The load cases considered include self-weight and a uniformly applied long-term load. Combinations of serviceability limit states (SLS) and ultimate limit states (ULS) have also been considered, among other factors. A finite-element software RFEM 5 has been used to implement the model. The results showed that the use of steel can reduce displacement by 20.96%. Additionally, compared to the maximum deflection provided by Eurocode 5 for softwoods, hardwoods can exhibit an increasing rate of 85.63%. By harnessing the plastic resistance of steel, the bending resistance of wood can be increased by 32.94%.
文摘An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.
基金supported by the National Natural Science Foundation of China(Nos.11702251,12002316)。
文摘Based on the nonlinear constitutive equation,a piezoelectric semiconductor(PSC)fiber under axial loads and Ohmic contact boundary conditions is investigated.The analytical solutions of electromechanical fields are derived by the homotopy analysis method(HAM),indicating that the HAM is efficient for the nonlinear analysis of PSC fibers,along with a rapid rate of convergence.Furthermore,the nonlinear characteristics of electromechanical fields are discussed through numerical results.It is shown that the asymmetrical distribution of electromechanical fields is obvious under a symmetrical load,and the piezoelectric effect is weakened by an applied electric field.With the increase in the initial carrier concentration,the electric potential decreases,and owing to the screen-ing effect of electrons,the distribution of electromechanical fields tends to be symmetrical.
基金supported by the Iranian Nanotechnology Development Committee and the University of Kashan(No.363452/10)
文摘The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.
文摘The offshore reinforced concrete structures are always subject to cyclic load, such as wave load.In this paper a new finite element analysis model is developed to analyze the stress and strain state of reinforced concrete structures including offshore concrete structures, subject to any number of the cyclic load. On the basis of the anal ysis of the experimental data,this model simplifies the number of cycles-total cyclic strain curve of concrete as three straight line segments,and it is assumed that the stress-strain curves of different cycles in each segment are the same, thus the elastoplastic analysis is only needed for the first cycle of each segment, and the stress or strain corresponding to any number of cycles can be obtained by superposition of stress or strain obtained by the above e lastoplastic analysis based on the cyclic numbers in each segment.This model spends less computer time,and can obtain the stress and strain states of the structures after any number of cycles.The endochronic-damage and ideal offshore concrete platform subject to cyclic loading are experimented and analyzed by the finite element method based on the model proposed in this paper. The results between the experiment and the finite element analysis are in good agreement,which demonstrates the validity and accuracy of the proposed model.
文摘The nonlinear analysis of reinforced concrete rectangular slabs undermonotonic transverse loads is performed by finite element method.The layered rectangu-lar element with 4 nodes and 20 degrees of freedom is developed,in whichbending-stretching coupling effect is taken into account.An orthotropic equivalentuniaxial stress-strain constitutive model of concrete is used.A program is worked out andused to calculate two reinforced concrete slabs.The results of calculation are in goodconformity with the corresponding test results.In addition,the influence of tension stif-fening effect of cracked concrete on the results of calculation is discussed.
基金supported by the National Natural Science Foundation of China (NSFC) (10872155)
文摘This paper first analyzes the features of two classes of numerical methods for global analysis of nonlinear dynamical systems, which regard state space respectively as continuous and discrete ones. On basis of this understanding it then points out that the previously proposed method of point mapping under cell reference (PMUCR), has laid a frame work for the development of a two scaled numerical method suitable for the global analysis of high dimensional nonlinear systems, which may take the advantages of both classes of single scaled methods but will release the difficulties induced by the disadvantages of them. The basic ideas and main steps of implementation of the two scaled method, namely extended PMUCR, are elaborated. Finally, two examples are presented to demonstrate the capabilities of the proposed method.
文摘A high order single step β algorithm, a new direct integration algorithm is proposed for solution of equations of motion. Whenβ=0.5, the accuracy of displacement, velocity and acceleration is of forth order (a truncation error of Δ t 5), and the algorithm is unconditionally stable and has no arithmetic damping and no overshooting. When >0.5, and an arithmetic damping is adopted, the algorithm is again unconditionally stable with a third order accuracy (a truncation error of Δ t 4). The analyses run with typical examples show that the algorithm proposed has higher speed, higher precision and better properties than other direct integration methods, such as Wilson θ method and Newmark β method in analysing linear elastic responses and nonlinear earthquake responses.
文摘The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.
文摘In this paper, the homotopy analysis method is applied to deduce the periodic solutions of a conservative nonlinear oscillator for which the elastic force term is proportional to<em> u</em><sup>1/3</sup>. By introducing the auxiliary linear operator and the initial guess of solution, the homotopy analysis solving is set up. By choosing the suitable convergence-control parameter, the accurate high-order approximations of solution and frequency for the whole range of initial amplitudes can easily be obtained. Comparison of the results obtained using this method with those obtained by different methods reveals that the former is more accurate, effective and convenient for these types of nonlinear oscillators.
文摘In this study, the homotopy analysis method (HAM) is used to solve the generalized Duffing equation. Both the frequencies and periodic solutions of the nonlinear Duffing equation can be explicitly and analytically formulated. Accuracy and validity of the proposed techniques are then verified by comparing the numerical results obtained based on the HAM and numerical integration method. Numerical simulations are extended for even very strong nonlinearities and very good correlations which achieved between the results. Besides, the optimal HAM approach is introduced to accelerate the convergence of solutions.
基金Project(2006318802111) supported by West Traffic Construction Science and Technology of ChinaProject(2008yb004) supported by Excellent Doctorate Dissertations of Central South University, China Project(2008G032-3) supported by Key Item of Science and Technology Research of Railway Ministry of China
文摘Based on the upper bound limit analysis theorem and the shear strength reduction technique, the equation for expressing critical limit-equilibrium state was employed to define the safety factor of a given slope and its corresponding critical failure mechanism by means of the kinematical approach of limit analysis theory. The nonlinear shear strength parameters were treated as variable parameters and a kinematically admissible failure mechanism was considered for calculation schemes. The iterative optimization method was adopted to obtain the safety factors. Case study and comparative analysis show that solutions presented here agree with available predictions when nonlinear criterion reduces to linear criterion, and the validity of present method could be illuminated. From the numerical results, it can also be seen that nonlinear parameter rn, slope foot gradient ,β, height of slope H, slope top gradient a and soil bulk density γ have significant effects on the safety factor of the slope.
文摘We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on the Laplace trans- form with the homotopy analysis method （HAM）. This method is a powerful tool for solving a large amount of problems. This technique provides a series of functions which may converge to the exact solution of the problem. A good agreement between the obtained solution and some well-known results is obtained.
基金National Natural Science Foundation of China for Innovative Research Groups Under Grant No.50321803 & 50621062National Natural Science Foundation of China Under Grant No.50808113 & 10872148
文摘This paper introduces an orthogonal expansion method for general stochastic processes. In the method, a normalized orthogonal function of time variable t is first introduced to carry out the decomposition of a stochastic process and then a correlated matrix decomposition technique, which transforms a correlated random vector into a vector of standard uncorrelated random variables, is used to complete a double orthogonal decomposition of the stochastic processes. Considering the relationship between the Hartley transform and Fourier transform of a real-valued function, it is suggested that the first orthogonal expansion in the above process is carried out using the Hartley basis function instead of the trigonometric basis function in practical applications. The seismic ground motion is investigated using the above method. In order to capture the main probabilistic characteristics of the seismic ground motion, it is proposed to directly carry out the orthogonal expansion of the seismic displacements. The case study shows that the proposed method is feasible to represent the seismic ground motion with only a few random variables. In the second part of the paper, the probability density evolution method （PDEM） is employed to study the stochastic response of nonlinear structures subjected to earthquake excitations. In the PDEM, a completely uncoupled one-dimensional partial differential equation, the generalized density evolution equation, plays a central role in governing the stochastic seismic responses of the nonlinear structure. The solution to this equation will yield the instantaneous probability density function of the responses. Computational algorithms to solve the probability density evolution equation are described. An example, which deals with a nonlinear frame structure subjected to stochastic ground motions, is illustrated to validate the above approach.
基金Projects(03JJY3078, 04JJ40032) supported by the Natural Science Foundation of Hunan Province, China project(03A006) supported by Scientific Research Fund of Hunan Provincial Education Department, China
文摘Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for high tailings dam was developed according to the results of engineering geological exploration, static and dynamic test and stability analysis on Baizhishan tailing dam 113.5 m high. The law of generation, diffusion and dissipation of seismic pore water pressure during and after earthquake was investigated, and the results of tailings dam’s acceleration, seismic dynamic stress and pore water pressure were obtained. The results show that the seismic stability and liquefaction resistance of high tailings dam are strengthened remarkably, and the scope and depth of liquefaction area at the top of dam are reduced greatly. The interior stress is compressive stress, the stress level of every element is less than 1.0 and the safety coefficient of every element is greater than 1.0. The safety coefficient against liquefaction of every element of tailing dam is greater than 1.5 according to the effective stress analysis of seismic response by finite element method. The calculated results prove that liquefaction is the main reason of seismic failure of high tailing dams, and the effect of seismic inertia forces on high tailing dams’ stability during earthquake is secondary reason.
基金supported by the Fundamental Research Funds for the Central Universities(No.N090405009)
文摘A homotopy analysis method（HAM）is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existence of small parameters in the considered equation.The HAM provides a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter.Two examples are presented to show that the HAM solutions agree well with the results of the modified Linstedt-Poincar＇e method and the incremental harmonic balance method.
基金supported by the National Natural Science Foundation of China(No.11972112)the Fundamental Research Funds for the Central Universities of China(Nos.N2103024 and N2003014)the National Science and Technology Major Project of China(No.J2019-I-0008-0008)。
文摘Soft nonlinear support is a major engineering project,but there are few relevant studies.In this paper,a dynamic pipeline model with soft nonlinear supports at both ends is established.By considering the influence of the Coriolis force and centrifugal force,the dynamical coupling equation of fluid-structure interaction is derived with extended Hamilton’s principle.Then,the approximate analytical solutions are sought via the harmonic balance method.The amplitude-frequency response curves show that different effects can be determined by approximate analysis.It is demonstrated that the increase in the fluid velocity can increase the amplitude of the pipeline system.The frequency range of unstable response increases when the fluid pressure raises.The combination of the soft nonlinear clamp and the large geometrical deformation of the pipeline affects the nonlinear vibration characteristic of the system,and the external excitation force and damping have significant effects on the stability.
文摘Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.
基金Supported by the National Natural Science Foundation of China (10871178)the Natural Science Foundation(Y606154)the Foundation of the Eduction Department of Zhejiang Province of China (Y200804008)
文摘Under the hypotheses that the second-order and third-order derivatives of a function are bounded, an estimate of the radius of the convergence ball of Ostrowski-Traub＇s method is obtained. An error analysis is given which matches the convergence order of the method. Finally, two examples are provided to show applications of our theorem.
文摘Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ordinary solution techniques lead to instability near the limit points and also have problems in case of snap-through and snap-back. Thus they fail to predict the complete load-displacement response. The arc-length method serves the purpose well in principle, received wide acceptance in finite element analysis, and has been used extensively. However modifications to the basic idea are vital to meet the particular needs of the analysis. This paper reviews some of the recent developments of the method in the last two decades, with particular emphasis on nonlinear finite element analysis of reinforced concrete structures.
