A new simple method of implicit time integration for dynamic problems of engineering structures
|
|
- Audra Page
- 6 years ago
- Views:
Transcription
1 Acta Mech Sin :91 99 DOI /s RESEARCH PAPER A new simple method of implicit time integration for dynamic problems of engineering structures Jun Zhou Youhe Zhou Received: 17 April 2006 / Revised: 10 July 2006 / Accepted: 6 September 2006 / Published online: 16 January 2007 Springer-Verlag 2007 Abstract This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems. Keywords Initial-value problems Time integration Implicit method Higher accuracy Time step and stability The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China G , 2003CB716707, the National Science Fund for Distinguished Young Scholars , the National Natural Science Foundation of China Key Program , the Research Fund for Oversea Chinese The English text was polished by Keren Wang. J. Zhou Y. Zhou B Department of Mechanics, Lanzhou University, Lanzhou , China zhouyh@lzu.edu.cn 1 Introduction The simulation of a complex dynamic system requires a high efficient algorithm of time integration, with high accuracy and limited amount of computation. This requirement has attracted many researchers 1,2]. Many approaches have been established to solve different dynamic problems, such as structure, fluid flowing, and multi-rigid-body systems 3 ]. They are mainly classified into two categories: explicit and implicit algorithms. The advantage of an explicit approach is that it avoids solving algebraic equations at each time step and leads to less computation. Recently, some comparative studies of explicit time integration methods were conducted 13, 1]. These methods were improved for solving linear or nonlinear dynamic problems 15, 16], and some of them were applied to large-scale dynamic systems such as railway ballast vibration systems 16, 17]. However, the inherent deficiency of an explicit approach is that it is conditionally stable and may require a tradeoff between the accuracy, stability and the computation time in simulating dynamics of large-scale structures. The trade-off sometimes depends on the researchers experience. In order to overcome this deficiency, some implicit approaches are proposed to achieve unconditional stability for large time-step sizes, e.g., Wilson-θ method 18], Newmark-β method 19] and Runge Kutta method. Of course, implicit methods may also need a balance between the accuracy and the computation time in practice when the stability is to be assured. In order to get higher order accuracy most studies were based on the Newmark method for different applications in dynamic systems 20 26]. To solve nonlinear dynamic systems with energy conservation, for example, the basic Newmark-β method was modified by keeping
2 92 J. Zhou, Y. Zhou the total energy conservation at the cost of reducing accuracy 2]. Thus the accuracy of the modified algorithm is degenerated to order one for some circumstances and a small time step is required to ensure the desired accuracy. In order to get an accuracy of arbitrary order, a complex-time-step Newmark method was established by operating complex matrices at each time step, with a number of parameters chosen 20]. It is obvious that this method may render emerge computations comparable with other ones when it is applied to a dynamic system with large number of degrees of freedom. In Ref. 25], a multi-time-step technique was developed on the basis of the basic formulae of the Newmark method when the formulae are re-written in the matrix form, and the accuracy and stability remain unchanged. In this paper, a new simple two-parameter implicit time integration method is developed. The new implicit scheme possesses satisfactory integration properties and may serve as a possible basis for further modifications to achieve better performances on accuracy, stability and numerical dissipation and dispersion. It has a second-order accuracy in general, and reaches a thirdorder accuracy as its control parameters are set to their optimal values. Its stability is analysized and its amplitude decay rate AD and the relative period error TD are lower than those of the Newmark-β method and the Zhai-method 16]. In the following section, the time integration algorithm is briefly introduced, while the analysis of accuracy and stability of the algorithm are carried out in Sect. 3. After that, the evaluations of amplitude decay rate AD and relative error of period TD insimulating a harmonic system are conducted in Sect., and comparisons with other efficient methods are also presented. Finally, some concluding remarks are made in Sect Essence of the new implicit time integration method The dynamic systems described by second order ordinary differential equations can be written in a compact form of matrix Mẍ + Cẋ + Kx = F, 1 or MA + CV + KX = F, 2 where M, C, and K stand for the mass, damping, and stiffness matrices, respectively; x, V= ẋ, and A= ẍ are the vectors of displacement, velocity, and acceleration, respectively; and F represents the external force. The symbol dot over a quantity means differentiation of the quantity with respect to time variable t. For the initial-value problem of a dynamic system, the initial conditions are formulated as X t=0 = X 0 and V t=0 = V 0, 3 in which X 0 and V 0 are specific vectors. The goal of a numerical integration of Eqs. 1 or2 associated with the initial conditions of Eq. 3 is to find an efficient algorithm when the matrices in the equations are very large. Applying the Taylor expansion of Vt n + t up to order 2 and Xt n + t up to the order 3 at the time step t n = n t, we have Vt n + t Vt n + At n t + 1 2!Ȧt n t 2, Xt n + t Xt n + Vt n t + 1 2! At n t !Ȧt n t 3, 5 where A n, V n, and X n are, respectively, the acceleration, velocity, and displacement at time step t n n = 0, 1, 2,..., i.e., A n = A t=tn, V n = V t=tn, and X n = X t=tn. Further, the central difference is employed to the calculation of Ȧt n, i.e., Ȧt n = A n+1 A n t Substituting Eq. 6 into Eqs. and 5, we obtain the functions of V n+1 and X n+1 in terms of variables of A n 1, A n, and A n+1 except for V n and X n. Once we introduce six weighted factors to the former three variables in the resulting two formulae, we have V n+1 = V n + ξ 1 A n t + 1 δ 1A n+1 t 1 ζ 1A n 1 t, 7 X n+1 = X n + V n t ξ 2A n t δ 2A n+1 t 2 1 ζ 2A n 1 t 2. 8 For A n 1 and A n+1 in the above equations, we take their Taylor expansions at A n up to order 3 and have A n 1 A n Ȧ n t Än t 1 A n t 3, 6 9 A n+1 A n + Ȧ n t Än t + 1 A n t
3 A new simple method of implicit time integration 93 Then, the substitution of Eqs. 9 and 10 into Eqs. 7 and 8 leads to V n+1 = V n + ξ 1 + δ 1 ζ 1 A n t + δ 1 + ζ 1 Ȧ n t 2 + δ 1 ζ 1 Ä n t δ 1 + ζ 1 A n t + O t 5, 11 2 ξ2 X n+1 = X n + V n t δ 2 ζ 2 A n t 2 + δ 2 + ζ 2 Ȧ n t 3 v δ 2 ζ 2 Ä n t 2 + δ 2 + ζ 2 A n t 5 + O t Comparing Eqs. 11 and witheqs. and 5, respectively, the accuracy of order 2 of the scheme in fact, the velocity has an accuracy of order 2 and the displacement has a accuracy of order 3 may be ensured by using the following relations ξ 1 = ζ 1, δ 1 = 2 ζ 1, ξ 2 = ζ 2, δ 2 = 2 ζ Thus two parameters in the scheme are independent of each other. Substituting the relations in Eq. 13 into Eqs. 11 and, we get the implicit algorithm of time integration with two control parameters in the form 1 V n+1 = V n ζ ζ 1 2 X n+1 = X n + V n t ζ 2 A n t A n+1 t 1 ζ 1A n 1 t, ζ 2 A n t 2 6 A n+1 t 2 1 ζ 2A n 1 t Substituting Eqs. 1 and 15 into Eq. 2 at time step t n+1,or MA n+1 + CV n+1 + KX n+1 = F n+1, 16 we obtain A n+1 = M 1 F n+1, 17 here, M = M + ζ C t + ζ K t 2, 18 6 F n+1 = F n+1 KX n C + K tv n + ζ1 C t + ζ ζ1 C t + ζ 2 K t ] K t 2 A n 3 ] A n From Eqs. 1 19, one sees that the calculation of time integration for Eqs. 1 or2 may start once the initial conditions of Eq. 2 are known and the values of acceleration A at the initial instant and at one time step before that time are specified. Without loss of generality, the initial acceleration may be obtained by substituting the initial conditions of Eq. 2 into the differential Eq. 2, and assuming A 1 = 0. FromEq.18, we find that the revised mass matrix is a constant matrix once the parameters ζ 1 and ζ 2 are specified and the equal time step t is taken, so the calculations of inversion matrix M 1 can be avoided at all time steps except for the first time step. Thus a large amount of CPU time may be saved when it is applied to complex dynamic systems with emerge degrees of freedom. Suppose X R m and a fixed time step size is used in the computation, then the computation of inverse matrix M 1 needs to be done only once. And subsequently, one needs about O2m 2 operations to compute vector F n+1, and another Om 2 operations to calculate vector A n+1 in each time step, as can be seen in Eqs. 17 and 19. Thus, the total number of operations is about O3m 2 in each time step, which is about the same as the Newmark method and the Zhai s method. Of course, the new method needs a bit more computation for vectors V n+1 and X n+1 in each time step, as compared with the Newmark method and the Zhai s method. 3 Analysis of accuracy and stability As shown in the previous section, the implicit method of time integration with Eqs has an accuracy of order 2 when the relations in Eq. 13 are satisfied. In this section, we present some conditions for higher accuracy of the scheme and stability of the calculations.
4 9 J. Zhou, Y. Zhou 3.1 Higher order accuracy The substitution of the expansions of A n 1 and A n+1 expressed in Eqs. 9 and 10 into Eqs. 1 and 15 leads to the following formulae V n+1 = V n + A n t + 1 2Ȧn t ζ 1Ä n t A n t +O t 5, 20 X n+1 = X n + V n t A n t Ȧn t ζ 2Ä n t + 1 A n t 5 +O t It is obvious that the accuracy of the scheme depends on the two independent parameters ζ 1 and ζ 2. Subtracting the Taylor expansions of V n+1 and X n+1 of one order higher than the formulae of Eqs. and 5 from Eqs. 20 and 21, respectively, we get the local truncation errors of the form EV n+1 = ζ 1 Ä n t 3 + O t, 22 EX n+1 = ζ 2 Ä n t + O t When ζ 1 = 1/3, 2 we have EV n+1 = O t and EX n+1 = O t, which tells us that the scheme has an accuracy of order 3 higher than that of simple methods of the Newmark method and the Wilson-θ method. 3.2 Stability analysis For simplicity, here, we present a stability analysis of the proposed scheme of implicit time integration with two parameters for the linear dynamic systems. In this case, the system of linear dynamic equations formulated by Eq. 2 may be always decoupled into a set of independent single-degree-of-freedom ordinary differential equations by using the modal decomposition technique. Thus, the stability analysis is usually performed on the basis of a dynamic system of a single-degree-of-freedom without damping in the dimensionless form, i.e., ẍ + ω 2 x = 0, 25 or equivalently in the difference form, a n+1 + ω 2 x n+1 = 0, 26 where ω= k/m is a natural frequency of the system; a = ẍ, and m and k are the mass and the stiffness, respectively. By means of Eqs. 1 and 15 for the case of the single-degree-of-freedom system, and substituting the formulas into Eq. 26, then eliminating variables V n and V n+1, we get a difference equation of Eq. 2 as follows ζ ] p ζ 2 ζ 1 2 x n+2 + ] p ζ 2 ζ ] 1 p x n+1 ζ2 x n + ζ ] 1 p x n 1 =0, 27 here, p = ω t 2 and n 1. Let x n+1 = λx n, then one obtains an eigenvalue equation of Eq. 27 in terms of λ as follows ζ 2p + 1 ] ζ 2 1 ] ζ 1 p λ 2 ] p λ+ p = 0. λ ζ ζ 1 1 ζ 2 1 ζ 1 28 Denote the roots of Eq. 28byλ i i = 1, 2, 3. Following the theory of difference 27], we know that the general solution of Eq. 27 can be expressed as x n = 3 c i λ n i, n = 1, 2, 3,..., 29 i=1 in which the coefficients c i i = 1, 2, 3 are determined by the initial conditions of Eq. 3. The stability of the scheme requires that the general solution x n should be bounded as n, which means that the roots should all have a magnitude less than or equal to 1, i.e., λ i 1, i ={1, 2, 3}, 30 which determine the regions of parameters. In order to get the formulae of stability conditions, we introduce the complex variable transformation of the form λ = 1 + z/1 z. 31 From the theory of complex variables, we know that the conditions in Eq. 30 are equivalent to Rez <0. Substituting Eq. 31 into28, we get ã 3 z 3 + ã 2 z 2 + ã 1 + ã 0 = 0, 32 where ã 0 = ã 1 = 3p, ã 2 = 3ζ 1 + 1p, and ã 3 = 2ζ 2 3ζ 1 + 1p. Thus, the stability conditions become Rez i <0, i = 1, 2, 3, 33
5 A new simple method of implicit time integration 95 in which z i i = 1, 2, 3 are the root of Eq. 32. According to the Routh Hurwitz criterion 28] of stability for Eq. 32, the stability conditions correspond to the case of no sign change of the first column in the Routh criterion table as follows λ 3 λ 2 ã 3 ã 1 ã 2 ã 0 λ 1, 3 b 2 0 λ 0 c 2 0 where b 2 = ã 3 ã 0 ã 2 ã 1 /ã 2 and c 2 = ã 0.Duetothe fact that p > 0or c 2 < 0, the stability conditions of the scheme reduce to the form 2ζ 2 3ζ 1 + 1p < 0, 3ζ 1 + 1p < 0, 35 and ζ 2 3ζ 1 p 2 > 0. From the inequalities in Eq. 35, we see that the unconditional stability of the scheme is obtained when the parameters satisfy the following inequalities 2ζ 2 3ζ , 3ζ , 36 and ζ 2 3ζ 1 > 0, or explicitly ζ 1 1/3, and 3ζ 1 <ζ 2 3ζ 1 1/2. 37 If the parameters are in the regions of 2ζ 2 3ζ > 0 and 3ζ > 0, the conditions of conditional stability of the scheme are obtained, ζ 1 > 1/3, ζ 2 > 3ζ 1, 38 and p < /2ζ 2 3ζ It is evident that the optimal values of ζ 1 expressed in Eq. 2 do not satisfy the first two conditions of Eq. 37 but 38, which implies that the stability of the new scheme with accuracy of order 3 is conditional. The time step should satisfy p = ω t 2 < 6/ζ 2 or t < 2.9/ω ζ 2 ζ 2 > 1. Usually, the time step is not greater than 0.2T, orf t 0.2, where T = 2π/ω and f = 1/T represent a period of the dynamic system and corresponding frequency, respectively. Corresponding to the criterion value of f t max = 0.2, we get q 0 = 625/π , where q = 2ζ 2 3ζ 1 +1 = q 0. For the case of q 625/π 2, the following inequality is obtained by the last inequality of Eq. 38 2ζ 2 3ζ /π Fig. 1 The conditional and unconditional stability regions of parameters ζ 1 and ζ 2 in the implicit scheme of time integration for the dynamic systems. Here, q = 2ζ 2 3ζ 1 + 1andq 0 = Then, selecting different value of q q 0, we can draw the lines of 2ζ 2 3ζ = q in the plane of ζ 1 ζ 2. According to the conditions of Eqs. 37 and 38, the regions of unconditional and conditional stability of the scheme are plotted in Fig. 1 for the parameters ζ 1 and ζ 2 of the scheme in the plane of ζ 1 ζ 2. When the parameters are in these two regions, the implicit method of time integration proposed in this paper for a dynamic system has an accuracy of at least order 2, and it has a third order accuracy when the parameters are chosen on the line of ζ 1 = 1/3 in the conditional stability region. When the parameters are on the line of 2ζ 2 3ζ = q, when q q 0, in the conditional stability region, the third inequality of Eq. 38 may be rewritten as f t /q /q 0, which gives a restriction of the time step with respect to the frequency of the dynamic system concerned. Figure 2 shows the curves of t versus frequency for different parameter q q 0.Asq increases, the curve drops down, and the time step should be in the region under and/or on the line marked by the value of q. For the parameters in the unconditional stability region and those under line of 2ζ 2 3ζ = q 0 in the conditional stability region, the time step should be selected in the region under the curve marked by q = q 0 in Fig. 2. Evaluation of dissipation and dispersion In order to further show the efficiency of the scheme, we need to evaluate the dissipation and dispersion of the new scheme. That is, we should evaluate the amplitude and period for a periodical dynamics. From the algebraic
6 96 J. Zhou, Y. Zhou Fig. 2 The diagram of log log curves of time step t versus frequency f Fig. 3 Curves of relative error of period TD versus the dimensionless time step t/t for the case of AD = 0 theory, we know that Eq. 28 has at least one real root and may have two conjugate complex roots. Denote the complex roots by λ 1 and λ 2, and the real root by λ 3.Let the conjugate roots be λ 1,2 = α ± iβ, where α and β are real and i = 1. When ζ 2 = 3ζ 1, we get λ 3 = 0. In the stability regions of parameters near the line of ζ 2 = 3ζ 1, we know that the conjugate complex roots are principal, or λ 3 0. In this case, the general solution of Eq. 27 may be rewritten as x n = eᾱt n c 1 cos ω t n + c 2 sin ωt n ], 0 in which ᾱ = lnα2 + β 2 2 t, ω = tan 1 β/α. 1 t After that, the measures of dissipation and dispersion are, respectively, defined by the amplitude decay function AD and the relative error of period TD of the form AD := 1 e 2πᾱ/ω, 2 TD := T T = 1, 3 T ω ω where T = 2π/ω and T = 2π/ ω. For the case of ζ 2 = 3ζ 1, it is shown that λ 1 λ 2 = 1, or α 2 + β 2 = 1, and { 2 λ 1,2 = ] 2 ζ 1 p / ±i } ] 2 3ζ 1p, here, = p p 2 /3 p 2 ζ 1. The complex roots of λ 1 and λ 2 satisfy the condition of >0 which means the second inequality of Eq. 35. In addition, the result α 2 + β 2 = 1 leads to AD = 0 from the first formula of Eq. 1 and Eq. 2. These results show that the line of ζ 2 = 3ζ 1 in the plane of parameters shown in Fig. 1 is in the stability regions, and the harmonic response may be well obtained by the new algorithm proposed in this paper for a harmonic dynamic system. Comparing this method with others in literature, e.g., the Newmarkβ method and the Zhai s method 16], one sees that AD = 0 is held only for a set of parameters in them, for example, γ = 1/ and β = 1/2 in the Newmark-β method, and ψ = φ = 1/2 in the Zhai s method. Our method, however, gives the parameters ζ 1 and ζ 2 in the line ζ 2 = 3ζ 1 with ζ 1 > 1/3. It is evident that the selected parameters for AD = 0 in the method of this paper are much more convenient than those in other methods. Considering the relative error of period TD for the case of AD = 0, we find that the relative error depends on both the independent parameter ζ 1 and the time step t. Figure 3 plots the curves of the relative error of period versus the ratio of dimensional time step for different parameter ζ 1 in the region of 1/8, 1/2]. From this figure, one sees that all relative errors for the case in the regions mentioned above are within 6%, and they are always less than those of the Newmark-β method and the Zhai s method, where the maximum relative errors are about 16 and 10%, respectively. When the new scheme has accuracy of order 3 and AD = 0, i.e., ζ 1 = 1/3 and ζ 2 = 1, the curve indicates that the TD is only within 1.3%. Figure shows the curves of amplitude dissipation ratio versus the dimensional time step for different parameters ζ 1 and ζ 2 in the vicinity of the stability regions near the line of ζ 2 = 3ζ 1. When ζ 1 is in the region of 1/8, 1/2], for example, the value of ζ 2 is chosen 20% higher than 3ζ 1,orζ 2 = 3.6ζ 1. This departure of parameters from those when AD = 0 in this paper is greater than those selected in other methods. For the curves of the Newmark method and the Zhai s method, there are 5%
7 A new simple method of implicit time integration 97 Fig. Characteristic curves of amplitude dissipation ratio AD versus dimensionless time step for different parameters ζ 1 and ζ 2 in the stability region Fig. 5 Comparison of the simulation results with analytical solution for the displacement response of a harmonic motion. The curves marked by values of ζ 1 and ζ 2 are obtained by the new scheme of time integration proposed in this paper departures from their perfect parameters. From Fig., we find that the amplitude dissipation ratio depends on the selection of the parameters and the time step, and most of the relative errors of the new method are less than those of other methods. When t/t < 0.15, the amplitude dissipation ratio of this method is always less than those of others, and the ratio generated from this new method is insensitive to the ratio of time step to period when the dimensionless time step is not greater than 0.1, but the dissipation ratio generated from other methods shown in this figure is strongly sensitive to the dimensionless time step even if the dimensionless time step is small, or close to zero. 5 Case study In this section, two numerical examples of applying the new scheme of time integration are displayed to show the efficiency of the scheme. The first one is the harmonic motion of Eq. 25 with the initial conditions of x0 = 0.1 and ẋ0 = 0, and with frequency ω = 10 rad/s. To this example, there is no difficulty to obtain its analytical solution xt = 0.1 cos10t. In the numerical simulation of this dynamic problem, we choose the time step t = T/. Figure 5 shows the simulation curves of displacement response of the problem during 5 10 periods. By comparing them with those obtained by the Newmark method and the Zhai s approach, we also plot their simulation results when their parameters in their schemes are taken as the perfect ones for the case of AD = 0. From this figure, one sees that the simulation results of this new scheme with accuracy of order 3, when either AD = 0see the curve marked by ζ 1 = 1/3 and ζ 2 = 3ζ 1 or AD = 0 see the curve marked by ζ 1 = 1/3 and ζ 2 = 1.1 are consistent with the analytical solution, while the curve with parameters ζ 1 = 1/6 and ζ 2 = 3ζ 1 when the scheme has accuracy of order 2 with AD = 0 is better than the curves of other methods with AD = 0 as compared to the analytic solution. The other example is a multi-degree-of-freedom dynamic system, a similar system of discrete axial vibration of one-dimensional straight bar, without damping and excitation, M 1 ẍ 1 + k 1 + k 2 x 1 k 2 x 2 = 0, M i ẍ i k i x i 1 + k i + k i+1 x i k i+1 x i+1 = 0, M n ẍ n k n x n 1 + k n x n = 0, 2 i n 1, 5 where M i and k i are the mass and spring stiffness, respectively; and the integer n 1 denotes the number of degrees of freedom. Here, x i represent the coordinates of mass M i relative to their equilibrium position. In this case, the parameters are taken as M i = kg and k i = N/m corresponding to a practical bar. In order to get a response consisting of all order vibrations, here, the initial conditions of displacement disturbance are chosen in the form x i 0 = 1 10 m, ẋ i 0 = 0, i = 1, 2,..., n. 6 From the modal analysis, it is not difficult for us to get the natural frequencies or periods of the dynamic system. Denote the maximum and minimum periods by T max and T min, respectively. In the simulation of this dynamic system, we take ζ 1 = 1/3 and ζ 2 = 1 in the new
8 98 J. Zhou, Y. Zhou 6 Conclusions A new simple algorithm of implicit time integration of dynamic systems with multi-degrees-of-freedom is proposed in this paper. It is found that this new algorithm with two control parameters has accuracy of at least order 2 when the parameters are in the stability regions, and of order 3 when the parameters are optimally taken as ζ 1 = 1/3 and ζ 2 ζ 1. The stability regions of the parameters and the time step are obtained by analysis of stability behavior of the scheme calculations and amplitude dissipation and dispersion. The numerical examples show that the proposed new algorithm can achieve satisfactory accuracy. Fig. 6 Numerical solutions of time response of the 11th mass in the duration of 30T max, which are obtained by the Newmark method, Zhai s method and the new scheme, respectively Fig. 7 Higher order components of the numerical solutions of the 11th mass scheme with third order accuracy and t = T min /6, and the simulation of response is conducted in the duration of 30T max. To make comparison, the explicit Runge Kutta method of order 29,30] is employed as the third-party method to yield a solution of the system with sufficient accuracy by setting a small time-step size, then numerical solutions are obtained by the Newmark method, Zhai-approach and the new method, respectively. Figure 6 shows the displacement response of the 11th mass. As can be seen in Fig. 6, the solution obtained by the new scheme is in good agreement with the ones obtained by the Newmark-method, Zhai s method and the Runge Kutta method. In addition, the new scheme can also achieve satisfactory accuracy for higher order components of the solution, which is shown in Fig. 7 by zooming out the figure of the displacement response. References 1. Belytschko, T., Hughes, T.J.R.: Computational Methods for Transient Analysis Computational Methods in Mechanics vol. 1. NorthHooland Elsevier, Amsterdam Wang, X.C., Shao, M.: The Fundamentals of the Finite Element Method and Numerical Methods, 2nd edn. Tsinghua University Press, Beijing Hilber, H.M., Hughes, T.J., Taylor, R.L.: Improved numerical dissipation for time integration algorithms in structural dynamics. Earthq. Eng. Struct. Dyn. 5, Stewart, D.E., Trinkle, J.C.: An implicit time-stepping scheme for rigid body dynamics with inelastic collisions and coulomb friction. Int. J. Numer. Methods Eng. 39, Villasenor, R.: A fully coupled, implicit, numerical scheme for laminar and turbulent parabolic flows. Int. J. Numer. Methods Eng. 0, Chen, S., Hansen, J.M., Tortorelli, D.A.: Unconditionally energy stable implicit time integration: application to multibody system analysis and design. Int. J. Numer. Methods Eng. 8, Modak, S., Sotelino, E.D.: The iterative group implicit algorithm for parallel transient finite element analysis. Int. J. Numer. Methods Eng. 7, Mugan, A., Hulbert, G.M.: Frequency-domain analysis of time-integration methods for semidiscrete finite element equations part I: Parabolic problems. Int. J. Numer. Methods Eng. 51, Mugan, A., Hulbert, G.M.: Frequency-domain analysis of time-integration methods for semidiscrete finite element equations part II: Hyperbolic and parabolic-hyperbolic problems. Int. J. Numer. Methods Eng. 51, Greyvenstein, G.P.: An implicit method for the analysis of transient flows in pipe networks. Int. J. Numer. Methods Eng. 53, Schreyer, H.L.: On time integration of viscoplastic constitutive models suitable for creep. Int. J. Numer. Methods Eng. 53, Mugan, A.: Discrete equivalent time integration methods for transient analysis. Int. J. Numer. Methods Eng. 57, Rio, G., Soive, A., Grolleau, V.: Comparative study of numerical explicit time integration algorithms. Adv. Eng. Softw. 36,
9 A new simple method of implicit time integration Rougier, E., Munjiza, A., John, N.W.M.: Numerical comparison of some explicit time integration schemes used in DEM, FEM/DEM and molecular dynamics. Int. J. Numer. Methods Eng. 61, Chung, J., Lee, J.M.: A new family of explicit time integration methods for linear and non-linear structural dynamics. Int. J. Numer. Methods Eng. 37, Zhai, W.M.: Two simple fast integration methods for largescale dynamic problems in engineering. Int. J. Numer. Methods Eng. 39, Zhai, W.M., Wang, K.Y., Lin, J.H.: Modelling and experiment of railway ballast vibration. J. Sound Vib. 270, Wilson, E.L., Farhoomand, I., Bathe, K.J.: Nonlinear dynamic analysis of complex structure. Earthq. Eng. Struct. Dyn. 1, Newmark, N.A.: A method of computation for structural dynamics. J. Eng. Mech.-ASCE, 85EM3, Fung, T.C.: Complex-time-step Newmark methods with controllable numerical dissipation. Int. J. Numer. Methods Eng. 1, Fish, J., Chen, W.: On accuracy, stability and efficiency of the Newmark method with incomplete solution by multilevel methods. Int. J. Numer. Methods Eng. 6, Kane, C., Marsden, J.E., Ortiz, M., West, M.: Variational integrators and the Newmark algorithm for conservative and dissipative mechanical systems. Int. J. Numer. Methods Eng. 9, Gravouil, A., Combescure, A.: Multi-time-step explicitimplicit method for non-linear structural dynamics. Int. J. Numer. Methods Eng. 50, Bui, Q.V.: Modified Newmark family for non-linear dynamic analysis. Int. J. Numer. Methods Eng. 61, Prakash, A., Hjelmstad, K.D.: A FETI-based multi-time-step coupling method for Newmark schemes in structural dynamics. Int. J. Numer. Methods Eng. 61, Krysl, P., Endres, L.: Explicit Newmark/Verlet algorithm for time integration of the rotational dynamics of rigid bodies. Int. J. Numer. Methods Eng. 62, Lick, W.J.: Difference Equations from Differential Equations. Springer, Berlin Dorf, R.C., Bishop, R.H.: Modern Control Systems. Addison Wesley Longman, Inc., New York Hairer, E., Nùrest, S.P., Wanner, G.: Solving Ordinary Differential Equations, vol. 1. Springer, Berlin Hairer, E., Wanner, G.: Solving Ordinary Differential Equations, vol. 2. Springer, Berlin 1991
New implicit method for analysis of problems in nonlinear structural dynamics
Applied and Computational Mechanics 5 (2011) 15 20 New implicit method for analysis of problems in nonlinear structural dynamics A. A. Gholampour a,, M. Ghassemieh a a School of Civil Engineering, University
More informationDAMPING MODELLING AND IDENTIFICATION USING GENERALIZED PROPORTIONAL DAMPING
DAMPING MODELLING AND IDENTIFICATION USING GENERALIZED PROPORTIONAL DAMPING S. Adhikari Department of Aerospace Engineering, University of Bristol, Queens Building, University Walk, Bristol BS8 1TR (U.K.)
More informationON THE ACCURACY AND COMPUTATIONAL COST OF TIME HISTORY ANALYSIS OF RESIDENTIAL BUILDINGS BY THE QUASI- WILSON-THETA METHOD
11 th International Conference on Vibration Problems Z. Dimitrovová et al. (eds.) Lisbon, Portugal, 9-12 September 2013 ON THE ACCURACY AND COMPUTATIONAL COST OF TIME HISTORY ANALYSIS OF RESIDENTIAL BUILDINGS
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 5/Part A - 23 November,
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationFree Vibration Analysis of Uniform Beams with Arbitrary Number of Cracks by using Adomian Decomposition Method
World Applied Sciences Journal 19 (1): 171-17, 1 ISSN 1818? 95 IDOSI Publications, 1 DOI: 1.589/idosi.was.1.19.1.581 Free Vibration Analysis of Uniform Beams with Arbitrary Number of Cracks by using Adomian
More informationCALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD
Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationActive control for a flexible beam with nonlinear hysteresis and time delay
THEORETICAL & APPLIED MECHANICS LETTERS 3, 635 (23) Active control for a flexible beam with nonlinear hysteresis and time delay Kun Liu, Longxiang Chen,, 2, a) and Guoping Cai ) Department of Engineering
More informationS12 PHY321: Practice Final
S12 PHY321: Practice Final Contextual information Damped harmonic oscillator equation: ẍ + 2βẋ + ω0x 2 = 0 ( ) ( General solution: x(t) = e [A βt 1 exp β2 ω0t 2 + A 2 exp )] β 2 ω0t 2 Driven harmonic oscillator
More informationNumerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
More informationImproving convergence of incremental harmonic balance method using homotopy analysis method
Acta Mech Sin (2009) 25:707 712 DOI 10.1007/s10409-009-0256-4 RESEARCH PAPER Improving convergence of incremental harmonic balance method using homotopy analysis method Yanmao Chen Jike Liu Received: 10
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationDynamics of Structures: Theory and Analysis
1. Free vibrations 2. Forced vibrations 3. Transient response 4. Damping mechanisms Dynamics of Structures: Theory and Analysis Steen Krenk Technical University of Denmark 5. Modal analysis I: Basic idea
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationSimulating Two-Dimensional Stick-Slip Motion of a Rigid Body using a New Friction Model
Proceedings of the 2 nd World Congress on Mechanical, Chemical, and Material Engineering (MCM'16) Budapest, Hungary August 22 23, 2016 Paper No. ICMIE 116 DOI: 10.11159/icmie16.116 Simulating Two-Dimensional
More information1 Introduction. Minho Lee 1 Jihoon Lee 1 Gunhee Jang 1
DOI 10.1007/s005-015-5-5 TECHNICAL PAPER Stability analysis of a whirling rigid rotor supported by stationary grooved FDBs considering the five degrees of freedom of a general rotor bearing system Minho
More informationRen-He s method for solving dropping shock response of nonlinear packaging system
Chen Advances in Difference Equations 216 216:279 DOI 1.1186/s1662-16-17-z R E S E A R C H Open Access Ren-He s method for solving dropping shock response of nonlinear packaging system An-Jun Chen * *
More informationSome Aspects of Structural Dynamics
Appendix B Some Aspects of Structural Dynamics This Appendix deals with some aspects of the dynamic behavior of SDOF and MDOF. It starts with the formulation of the equation of motion of SDOF systems.
More informationReduced-dimension Models in Nonlinear Finite Element Dynamics of Continuous Media
Reduced-dimension Models in Nonlinear Finite Element Dynamics of Continuous Media Petr Krysl, Sanjay Lall, and Jerrold E. Marsden, California Institute of Technology, Pasadena, CA 91125. pkrysl@cs.caltech.edu,
More informationThe student will experimentally determine the parameters to represent the behavior of a damped oscillatory system of one degree of freedom.
Practice 3 NAME STUDENT ID LAB GROUP PROFESSOR INSTRUCTOR Vibrations of systems of one degree of freedom with damping QUIZ 10% PARTICIPATION & PRESENTATION 5% INVESTIGATION 10% DESIGN PROBLEM 15% CALCULATIONS
More informationConserving energy and momentum in nonlinear dynamics: A simple implicit time integration scheme
Computers and Structures 5 (7) 37 5 www.elsevier.com/locate/compstruc Conserving energy and momentum in nonlinear dynamics: A simple implicit integration scheme Klaus-Jürgen Bathe * Massachusetts Institute
More informationSection 3.7: Mechanical and Electrical Vibrations
Section 3.7: Mechanical and Electrical Vibrations Second order linear equations with constant coefficients serve as mathematical models for mechanical and electrical oscillations. For example, the motion
More informationEXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro
EXTENDED ABSTRACT Dynamic analysis of elastic solids by the finite element method Vítor Hugo Amaral Carreiro Supervisor: Professor Fernando Manuel Fernandes Simões June 2009 Summary The finite element
More informationTransactions on Modelling and Simulation vol 20, 1998 WIT Press, ISSN X
A BEM-formulation applied to the solution of nonlinear dynamic problems H.J. Holl, A.K. Belyaev & H. Irschik Division of Technical Mechanics, University Linz, Austria Email: holl@mechatronik. uni-linz.
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationFigure XP3.1 (a) Mass in equilibrium, (b) Freebody diagram, (c) Kinematic constraint relation Example Problem 3.1 Figure XP3.1 illustrates a mass m
LECTURE 7. MORE VIBRATIONS ` Figure XP3.1 (a) Mass in equilibrium, (b) Freebody diagram, (c) Kinematic constraint relation Example Problem 3.1 Figure XP3.1 illustrates a mass m that is in equilibrium and
More informationA NEW METHOD FOR VIBRATION MODE ANALYSIS
Proceedings of IDETC/CIE 25 25 ASME 25 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference Long Beach, California, USA, September 24-28, 25 DETC25-85138
More informationThe Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry Friction
Send Orders for Reprints to reprints@benthamscience.ae 308 The Open Mechanical Engineering Journal, 2014, 8, 308-313 Open Access The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry
More informationCOMPOSITE MATERIAL WITH NEGATIVE STIFFNESS INCLUSION FOR VIBRATION DAMPING: THE EFFECT OF A NONLINEAR BISTABLE ELEMENT
11 th International Conference on Vibration Problems Z. Dimitrovová et.al. (eds.) Lisbon, Portugal, 9 12 September 2013 COMPOSITE MATERIAL WITH NEGATIVE STIFFNESS INCLUSION FOR VIBRATION DAMPING: THE EFFECT
More informationDEVELOPMENT OF A REAL-TIME HYBRID EXPERIMENTAL SYSTEM USING A SHAKING TABLE
DEVELOPMENT OF A REAL-TIME HYBRID EXPERIMENTAL SYSTEM USING A SHAKING TABLE Toshihiko HORIUCHI, Masahiko INOUE And Takao KONNO 3 SUMMARY A hybrid experimental method, in which an actuator-excited vibration
More information3 Mathematical modeling of the torsional dynamics of a drill string
3 Mathematical modeling of the torsional dynamics of a drill string 3.1 Introduction Many works about torsional vibrations on drilling systems [1, 12, 18, 24, 41] have been published using different numerical
More informationON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS
8th International Congress on Computational Mechanics, Volos, 1-15 July 015 ON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS E.G. Kakouris, V.K. Koumousis Institute of Structural
More informationIDENTIFICATION OF SEMI-RIGID JOINTS IN FRAME STRUCTURES
6th ECCOMAS Conference on Smart Structures and Materials SMART213 Politecnico di Torino, 24-26 June 213 E. Carrera, F. Miglioretti and M. Petrolo (Editors) www.smart213.com IDENTIFICATION OF SEMI-RIGID
More informationNonlinear Stability and Bifurcation of Multi-D.O.F. Chatter System in Grinding Process
Key Engineering Materials Vols. -5 (6) pp. -5 online at http://www.scientific.net (6) Trans Tech Publications Switzerland Online available since 6//5 Nonlinear Stability and Bifurcation of Multi-D.O.F.
More informationAseismic design of structure equipment systems using variable frequency pendulum isolator
Aseismic design of structure equipment systems using variable frequency pendulum isolator Pranesh Murnal a, Ravi Sinha b, a Department of Applied Mechanics, Government College of Engineering, Karad 415124,
More informationAn improved brake squeal source model in the presence of kinematic and friction nonlinearities
An improved brake squeal source model in the presence of kinematic and friction nonlinearities Osman Taha Sen, Jason T. Dreyer, and Rajendra Singh 3 Department of Mechanical Engineering, Istanbul Technical
More information1. Multiple Degree-of-Freedom (MDOF) Systems: Introduction
1. Multiple Degree-of-Freedom (MDOF) Systems: Introduction Lesson Objectives: 1) List examples of MDOF structural systems and state assumptions of the idealizations. 2) Formulate the equation of motion
More informationDr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationImaginary. Axis. Real. Axis
Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 25 3:3-6:3 p.m. ffl Print your name and sign the honor code statement ffl You may use your course notes,
More informationFinal Exam April 30, 2013
Final Exam Instructions: You have 120 minutes to complete this exam. This is a closed-book, closed-notes exam. You are allowed to use a calculator during the exam. Usage of mobile phones and other electronic
More informationAn explicit time-domain finite-element method for room acoustics simulation
An explicit time-domain finite-element method for room acoustics simulation Takeshi OKUZONO 1 ; Toru OTSURU 2 ; Kimihiro SAKAGAMI 3 1 Kobe University, JAPAN 2 Oita University, JAPAN 3 Kobe University,
More informationImaginary. Axis. Real. Axis
Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 2005 3:30-6:30 p.m. ffl Print your name and sign the honor code statement ffl You may use your course
More informationInteraction of Soil Static and Dynamic Stiffness and Buried Pipeline Under Harmonic Vibration
American J. of Engineering and Applied Sciences 3 (): 44-448, ISSN 94-7 Science Publications Interaction of Soil Static and Dynamic Stiffness and Buried Pipeline Under Harmonic Vibration Alireza Mirza
More informationAdditive resonances of a controlled van der Pol-Duffing oscillator
Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University
More informationWEEKS 8-9 Dynamics of Machinery
WEEKS 8-9 Dynamics of Machinery References Theory of Machines and Mechanisms, J.J.Uicker, G.R.Pennock ve J.E. Shigley, 2011 Mechanical Vibrations, Singiresu S. Rao, 2010 Mechanical Vibrations: Theory and
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationMODIFIED RUNGE_KUTTA METHOD FOR SOLVING NONLINEAR VIBRATION OF AXIALLY TRAVELLING STRING SYSTEM
The 21 st International Congress on Sound and Vibration 13-17 July, 214, Beijing/China MODIFIED RUNGE_KUTTA METHOD FOR SOLVING NONLINEAR VIBRATION OF AXIALLY TRAVELLING STRING SYSTEM Qun Wu, Enwei Chen*,
More informationSolutions of Nonlinear Oscillators by Iteration Perturbation Method
Inf. Sci. Lett. 3, No. 3, 91-95 2014 91 Information Sciences Letters An International Journal http://dx.doi.org/10.12785/isl/030301 Solutions of Nonlinear Oscillators by Iteration Perturbation Method A.
More informationThursday, August 4, 2011
Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.
More informationDynamic Loads CE 543. Examples. Harmonic Loads
CE 543 Structural Dynamics Introduction Dynamic Loads Dynamic loads are time-varying loads. (But time-varying loads may not require dynamic analysis.) Dynamics loads can be grouped in one of the following
More information1. Introduction
15. Pounding of adjacent buildings considering pile-soil-structure interaction Lihua Zou, Kai Huang, Liyuan Wang, Laiqing Fang 15. POUNDING OF ADJACENT BUILDINGS CONSIDERING PILE-SOIL-STRUCTURE INTERACTION.
More informationFinite Element Analysis Lecture 1. Dr./ Ahmed Nagib
Finite Element Analysis Lecture 1 Dr./ Ahmed Nagib April 30, 2016 Research and Development Mathematical Model Mathematical Model Mathematical Model Finite Element Analysis The linear equation of motion
More informationStructural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake
Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme
More informationBSc/MSci EXAMINATION. Vibrations and Waves. Date: 4 th May, Time: 14:30-17:00
BSc/MSci EXAMINATION PHY-217 Vibrations and Waves Time Allowed: 2 hours 30 minutes Date: 4 th May, 2011 Time: 14:30-17:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from
More informationMobility and Impedance Methods. Professor Mike Brennan
Mobility and Impedance Methods Professor Mike Brennan ibration control ibration Problem Understand problem Modelling (Mobility and Impedance Methods) Solve Problem Measurement Mobility and Impedance The
More informationLyapunov exponent calculation of a two-degreeof-freedom vibro-impact system with symmetrical rigid stops
Chin. Phys. B Vol. 20 No. 4 (2011) 040505 Lyapunov exponent calculation of a two-degreeof-freedom vibro-impact system with symmetrical rigid stops Li Qun-Hong( ) and Tan Jie-Yan( ) College of Mathematics
More informationDETERMINATION OF THE FREQUENCY-AMPLITUDE RELATION FOR NONLINEAR OSCILLATORS WITH FRACTIONAL POTENTIAL USING HE S ENERGY BALANCE METHOD
Progress In Electromagnetics Research C, Vol. 5, 21 33, 2008 DETERMINATION OF THE FREQUENCY-AMPLITUDE RELATION FOR NONLINEAR OSCILLATORS WITH FRACTIONAL POTENTIAL USING HE S ENERGY BALANCE METHOD S. S.
More informationResearch on the iterative method for model updating based on the frequency response function
Acta Mech. Sin. 2012) 282):450 457 DOI 10.1007/s10409-012-0063-1 RESEARCH PAPER Research on the iterative method for model updating based on the frequency response function Wei-Ming Li Jia-Zhen Hong Received:
More informationANALYSIS OF NONUNIFORM BEAMS ON ELASTIC FOUNDATIONS USING RECURSIVE DIFFERENTATION METHOD
Engineering MECHANICS, Vol. 22, 2015, No. 2, p. 83 94 83 ANALYSIS OF NONUNIFORM BEAMS ON ELASTIC FOUNDATIONS USING RECURSIVE DIFFERENTATION METHOD Mohamed Taha Hassan*, Samir Abo Hadima* Analytical solutions
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationGeneralized projective synchronization between two chaotic gyros with nonlinear damping
Generalized projective synchronization between two chaotic gyros with nonlinear damping Min Fu-Hong( ) Department of Electrical and Automation Engineering, Nanjing Normal University, Nanjing 210042, China
More information1-DOF Vibration Characteristics. MCE371: Vibrations. Prof. Richter. Department of Mechanical Engineering. Handout 7 Fall 2017
MCE371: Vibrations Prof. Richter Department of Mechanical Engineering Handout 7 Fall 2017 Free Undamped Vibration Follow Palm, Sect. 3.2, 3.3 (pp 120-138), 3.5 (pp 144-151), 3.8 (pp. 167-169) The equation
More informationVibrations Qualifying Exam Study Material
Vibrations Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering vibrations topics. These topics are listed below for clarification. Not all instructors
More informationA Second Order Extension of the Generalized α Method for Constrained Systems in Mechanics
A Second Order Extension of the Generalized α Method for Constrained Systems in Mechanics Laurent O. Jay 1 and Dan Negrut 2 1 Department of Mathematics, The University of Iowa, 14 MacLean Hall, Iowa City,
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationVibration analysis of concrete bridges during a train pass-by using various models
Journal of Physics: Conference Series PAPER OPEN ACCESS Vibration analysis of concrete bridges during a train pass-by using various models To cite this article: Qi Li et al 2016 J. Phys.: Conf. Ser. 744
More informationVibrations and Waves Physics Year 1. Handout 1: Course Details
Vibrations and Waves Jan-Feb 2011 Handout 1: Course Details Office Hours Vibrations and Waves Physics Year 1 Handout 1: Course Details Dr Carl Paterson (Blackett 621, carl.paterson@imperial.ac.uk Office
More informationIntroduction to Vibration. Professor Mike Brennan
Introduction to Vibration Professor Mie Brennan Introduction to Vibration Nature of vibration of mechanical systems Free and forced vibrations Frequency response functions Fundamentals For free vibration
More informationNonlinear Dynamic Systems Homework 1
Nonlinear Dynamic Systems Homework 1 1. A particle of mass m is constrained to travel along the path shown in Figure 1, which is described by the following function yx = 5x + 1x 4, 1 where x is defined
More informationHopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators
The Open Acoustics Journal 8 9-3 9 Open Access Hopf ifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators Jianping Cai *a and Jianhe Shen b a Department of
More informationDIFFERENTIAL TRANSFORMATION METHOD TO DETERMINE TEMPERATURE DISTRIBUTION OF HEAT RADIATING FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY
Published by Global Research Publications, New Delhi, India DIFFERENTIAL TRANSFORMATION METHOD TO DETERMINE TEMPERATURE DISTRIBUTION OF HEAT RADIATING FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY
More informationDynamical behaviour of a controlled vibro-impact system
Vol 17 No 7, July 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(07)/2446-05 Chinese Physics B and IOP Publishing Ltd Dynamical behaviour of a controlled vibro-impact system Wang Liang( ), Xu Wei( ), and
More informationSection 4.9; Section 5.6. June 30, Free Mechanical Vibrations/Couple Mass-Spring System
Section 4.9; Section 5.6 Free Mechanical Vibrations/Couple Mass-Spring System June 30, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This Session: (1) Free
More information1 Review of simple harmonic oscillator
MATHEMATICS 7302 (Analytical Dynamics YEAR 2017 2018, TERM 2 HANDOUT #8: COUPLED OSCILLATIONS AND NORMAL MODES 1 Review of simple harmonic oscillator In MATH 1301/1302 you studied the simple harmonic oscillator:
More information18.12 FORCED-DAMPED VIBRATIONS
8. ORCED-DAMPED VIBRATIONS Vibrations A mass m is attached to a helical spring and is suspended from a fixed support as before. Damping is also provided in the system ith a dashpot (ig. 8.). Before the
More informationPhysics 351, Spring 2017, Homework #2. Due at start of class, Friday, January 27, 2017
Physics 351, Spring 2017, Homework #2. Due at start of class, Friday, January 27, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at
More informationK. Pyragas* Semiconductor Physics Institute, LT-2600 Vilnius, Lithuania Received 19 March 1998
PHYSICAL REVIEW E VOLUME 58, NUMBER 3 SEPTEMBER 998 Synchronization of coupled time-delay systems: Analytical estimations K. Pyragas* Semiconductor Physics Institute, LT-26 Vilnius, Lithuania Received
More informationLab 1: Damped, Driven Harmonic Oscillator
1 Introduction Lab 1: Damped, Driven Harmonic Oscillator The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationNew Homoclinic and Heteroclinic Solutions for Zakharov System
Commun. Theor. Phys. 58 (2012) 749 753 Vol. 58, No. 5, November 15, 2012 New Homoclinic and Heteroclinic Solutions for Zakharov System WANG Chuan-Jian ( ), 1 DAI Zheng-De (à ), 2, and MU Gui (½ ) 3 1 Department
More informationIn-Plane and Out-of-Plane Dynamic Responses of Elastic Cables under External and Parametric Excitations
Applied Mathematics 5, 5(6): -4 DOI:.59/j.am.556. In-Plane and Out-of-Plane Dynamic Responses of Elastic Cables under External and Parametric Excitations Usama H. Hegazy Department of Mathematics, Faculty
More informationME 563 HOMEWORK # 7 SOLUTIONS Fall 2010
ME 563 HOMEWORK # 7 SOLUTIONS Fall 2010 PROBLEM 1: Given the mass matrix and two undamped natural frequencies for a general two degree-of-freedom system with a symmetric stiffness matrix, find the stiffness
More informationKEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY OSCILLATIONS AND WAVES PRACTICE EXAM
KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY-10012 OSCILLATIONS AND WAVES PRACTICE EXAM Candidates should attempt ALL of PARTS A and B, and TWO questions from PART C. PARTS A and B should be answered
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More informationChaos suppression of uncertain gyros in a given finite time
Chin. Phys. B Vol. 1, No. 11 1 1155 Chaos suppression of uncertain gyros in a given finite time Mohammad Pourmahmood Aghababa a and Hasan Pourmahmood Aghababa bc a Electrical Engineering Department, Urmia
More informationLab 1: damped, driven harmonic oscillator
Lab 1: damped, driven harmonic oscillator 1 Introduction The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationAn Estimation of Error-Free Frequency Response Function from Impact Hammer Testing
85 An Estimation of Error-Free Frequency Response Function from Impact Hammer Testing Se Jin AHN, Weui Bong JEONG and Wan Suk YOO The spectrum of impulse response signal from the impact hammer testing
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 13 http://acousticalsociety.org/ ICA 13 Montreal Montreal, Canada - 7 June 13 Signal Processing in Acoustics Session 1pSPc: Miscellaneous Topics in Signal
More informationReal-Time Hybrid Simulation of Single and Multiple Tuned Liquid Column Dampers for Controlling Seismic-Induced Response
6 th International Conference on Advances in Experimental Structural Engineering 11 th International Workshop on Advanced Smart Materials and Smart Structures Technology August 1-2, 215, University of
More informationHigh accuracy numerical and signal processing approaches to extract flutter derivatives
High accuracy numerical and signal processing approaches to extract flutter derivatives *NakHyun Chun 1) and Hak-eun Lee 2) 1), 2) School of Civil, Environmental and Architectural Engineering, Korea University,
More informationEffects of viscosity and varying gravity on liquid sloshing in a carrier subjected to external excitations
Int. J. Dynam. Control (204) 2:52 532 DOI 0.007/s40435-04-0072-y Effects of viscosity and varying gravity on liquid sloshing in a carrier subjected to external excitations Liming Dai Xiaojie Wang Received:
More informationTheory of Vibrations in Stewart Platforms
Theory of Vibrations in Stewart Platforms J.M. Selig and X. Ding School of Computing, Info. Sys. & Maths. South Bank University London SE1 0AA, U.K. (seligjm@sbu.ac.uk) Abstract This article develops a
More informationKink, singular soliton and periodic solutions to class of nonlinear equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10 Issue 1 (June 015 pp. 1 - Applications and Applied Mathematics: An International Journal (AAM Kink singular soliton and periodic
More informationEffects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device
Effects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device A. Krishnamoorthy Professor, Department of Civil Engineering Manipal
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More information