Matrix Iteration. Giacomo Boffi.
|
|
- Shanna Hunter
- 5 years ago
- Views:
Transcription
1 Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 12, 2016
2 Outline Second -Ritz Method
3 Dynamic analysis of MDOF systems based on modal superposition is both simple and efficient simple: the modal response can be easily computed, analitically or numerically, with the techniques we have seen for SDOF systems, efficient: in most cases, only the modal responses of a few lower modes are required to accurately describe the structural response.
4 As the structural matrices are easily assembled using the FEM, our modal superposition procedure is ready to be applied to structures with tenth, thousands or millions of DOF s! except that we can compute the eigenpairs only when the analyzed structure has two, three or maybe four degrees of freedom...
5 As the structural matrices are easily assembled using the FEM, our modal superposition procedure is ready to be applied to structures with tenth, thousands or millions of DOF s! except that we can compute the eigenpairs only when the analyzed structure has two, three or maybe four degrees of freedom... We will discuss how it is possible to compute the eigenpairs of arbitrary dynamic systems using the so called s procedure and a number of variations derived from this fundamental idea.
6 Equilibrium First, we will see an iterative procedure whose outputs are the first, or fundamental, mode shape vector and the corresponding eigenvalue. When an undamped system freely vibrates with a harmonic time dependency of frequency ω i, the equation of motion, simplifying the time dependency, is K ψ i = ω 2 i M ψ i. In equilibrium terms, the elastic forces are equal to the inertial forces when the systems oscillates with frequency ω 2 i and mode shape ψ i
7 Proposal of an iterative procedure Our iterative procedure will be based on finding a new displacement vector x n+1 such that the elastic forces f S = K x i+1 are in equilibrium with the inertial forces due to the old displacement vector x n, f I = ω 2 i M x n, that is K x n+1 = ω 2 i M x n. Premultiplying by the inverse of K and introducing the Dynamic Matrix, D = K 1 M x n+1 = ω 2 i K 1 M x n = ω 2 i D x n.
8 Proposal of an iterative procedure Our iterative procedure will be based on finding a new displacement vector x n+1 such that the elastic forces f S = K x i+1 are in equilibrium with the inertial forces due to the old displacement vector x n, f I = ω 2 i M x n, that is K x n+1 = ω 2 i M x n. Premultiplying by the inverse of K and introducing the Dynamic Matrix, D = K 1 M x n+1 = ω 2 i K 1 M x n = ω 2 i D x n. In the generative equation above we miss a fundamental part, the square of the free vibration frequency ω 2 i.
9 The Procedure, 1 This problem is solved considering the x n as a sequence of normalized vectors and introducing the idea of an unnormalized new displacement vector, ˆx n+1, ˆx n+1 = D x n, note that we removed the explicit dependency on ω 2 i.
10 The Procedure, 1 The normalized vector is obtained applying to ˆx n+1 a normalizing factor, F n+1, x n+1 = ˆx n+1 F n+1, but x n+1 = ω 2 i D x n = ω 2 i ˆx n+1, 1 F = ω2 i
11 The Procedure, 1 The normalized vector is obtained applying to ˆx n+1 a normalizing factor, F n+1, x n+1 = ˆx n+1 F n+1, but x n+1 = ω 2 i D x n = ω 2 i ˆx n+1, 1 F = ω2 i If we agree that, near convergence, x n+1 x n, substituting in the previous equation we have x n+1 x n = ω 2 i ˆx n+1 ω 2 i x n ˆx n+1.
12 The Procedure, 1 The normalized vector is obtained applying to ˆx n+1 a normalizing factor, F n+1, x n+1 = ˆx n+1 F n+1, but x n+1 = ω 2 i D x n = ω 2 i ˆx n+1, 1 F = ω2 i If we agree that, near convergence, x n+1 x n, substituting in the previous equation we have x n+1 x n = ω 2 i ˆx n+1 ω 2 i x n ˆx n+1. Of course the division of two vectors is not an option, so we want to twist it into something useful.
13 Normalization First, consider x n = ψ i : in this case, for j = 1,..., N it is x n,j /ˆx n+1,j = ω 2 i. When x n ψ i it is possible to demonstrate that we can bound the eigenvalue { } { } xn,j min ω 2 i j=1,...,n ˆx max xn,j. n+1,j j=1,...,n ˆx n+1,j
14 Normalization First, consider x n = ψ i : in this case, for j = 1,..., N it is x n,j /ˆx n+1,j = ω 2 i. When x n ψ i it is possible to demonstrate that we can bound the eigenvalue { } { } xn,j min ω 2 i j=1,...,n ˆx max xn,j. n+1,j j=1,...,n ˆx n+1,j A more rational approach would make reference to a proper vector norm, so using our preferred vector norm we can write ω 2 i ˆxT n+1 M x n ˆx T n+1 M ˆx, n+1
15 Normalization First, consider x n = ψ i : in this case, for j = 1,..., N it is x n,j /ˆx n+1,j = ω 2 i. When x n ψ i it is possible to demonstrate that we can bound the eigenvalue { } { } xn,j min ω 2 i j=1,...,n ˆx max xn,j. n+1,j j=1,...,n ˆx n+1,j A more rational approach would make reference to a proper vector norm, so using our preferred vector norm we can write ω 2 i ˆxT n+1 M x n ˆx T n+1 M ˆx, n+1 (if memory helps, this is equivalent to the R 11 approximation, that we introduced studying quotient refinements).
16 Proof of Convergence, 1 Until now we postulated that the sequence x n converges to some, unspecified eigenvector ψ i, now we will demonstrate that the sequence converge to the first, or fundamental mode shape, lim x n = ψ 1. n 1. Expand x 0 in terms of eigenvectors an modal coordinates: x 0 = ψ 1 q 1,0 + ψ 2 q 2,0 + ψ 3 q 3, The inertial forces, assuming that the system is vibrating according to the fundamental frequency, are f I,n=0 = ω 2 1M (ψ 1 q 1,0 + ψ 2 q 2,0 + ψ 3 q 3,0 + ) ( = M ω 2 ω 2 1 1ψ 1 q 1,0 ω 2 + ω 2 ω 2 ) 1 2ψ 2 q 2,0 1 ω
17 Proof of Convergence, 2 3. The deflections due to these forces (no hat!, we have multiplied by ω 2 1) are ( ) x n=1 = K 1 M ω 2 ω 2 1 1ψ 1 q 1,0 + ω 2 ω 2 1 ω 2ψ 2 2 q 2,0 +, 1 ω 2 2 (note that every term has been multiplied and divided by the corresponding eigenvalue ω 2 i ). 4. With ω 2 j M ψ j = Kψ j, substituting and simplifying K 1 K = I, ( ( ) ω x n=1 = K ( ) Kψ 1 q 1 ω 2 1 ( ) 1,0 + Kψ ω 2 2 q 1 ω 2 1 ) 2,0 + Kψ 1 ω 2 3 q 1 3,0 + 2 ω 2 3 ω 2 1 ω 2 1 ω 2 1 = ψ 1 q 1,0 + ψ ω 2 2 q 2,0 + ψ 1 ω 2 3 q 3,0 +, 2 ω 2 3
18 Proof of Convergence, 3 5. applying again this procedure ( ( ) ω 2 2 ( ) x n=2 = ψ 1 q 1 ω 2 2 ( ) 1,0 + ψ ω 2 2 q 1 ω 2 2 ) 2,0 + ψ 1 ω 2 3 q 1 3,0 +, 2 ω applying the procedure n times ( ( ) ω 2 n ( ) x n = ψ 1 q 1 ω 2 n ( ) 1,0 + ψ ω 2 2 q 1 ω 2 n ) 2,0 + ψ 1 ω 2 3 q 1 3, ω 2 3
19 Proof of Convergence, 4 Going to the limit, because Consequently, lim n x n = ψ 1 q 1,0 lim n ( ω 2 1 ω 2 j ) n = δ 1j x n lim n ˆx n = ω2 1
20 Purified Vectors If we know ψ 1 and ω 2 1 from the matrix iteration procedure it is possible to compute the second eigenpair, following a slightly different procedure.
21 Purified Vectors If we know ψ 1 and ω 2 1 from the matrix iteration procedure it is possible to compute the second eigenpair, following a slightly different procedure. Express the initial iterate in terms of the (unknown) eigenvectors, x n=0 = Ψ q n=0 and premultiply by the (known) ψ T 1 M: solving for q 1,n=0 ψ T 1 M x n=0 = M 1 q 1,n=0 q1,n=0 = ψt 1 M x n=0 M 1.
22 Purified Vectors If we know ψ 1 and ω 2 1 from the matrix iteration procedure it is possible to compute the second eigenpair, following a slightly different procedure. Express the initial iterate in terms of the (unknown) eigenvectors, x n=0 = Ψ q n=0 and premultiply by the (known) ψ T 1 M: solving for q 1,n=0 ψ T 1 M x n=0 = M 1 q 1,n=0 q1,n=0 = ψt 1 M x n=0 M 1. Knowing the amplitude of the 1st modal contribution to x n=0 we can write a purified vector, y n=0 = x n=0 ψ 1 q 1,n=0.
23 Convergence (?) It is easy to demonstrate that using y n=0 as our starting vector lim y n = ψ 2 q 2,n=0, n y n lim n ŷ n = ω2 2. because the initial amplitude of the first mode is null.
24 Convergence (?) It is easy to demonstrate that using y n=0 as our starting vector lim y n = ψ 2 q 2,n=0, n y n lim n ŷ n = ω2 2. because the initial amplitude of the first mode is null. Due to numerical errors in the determination of fundamental mode and in the procedure itself, using a plain matrix iteration the procedure however converges to the 1st eigenvector, so to preserve convergence to the 2nd mode it is necessary that the iterated vector y n is purified at each step n.
25 Purification Procedure The purification procedure is simple, at each step the amplitude of the 1st mode is first computed, then removed from the iterated vector y n q 1,n = ψ T 1 My n /M 1, ( ŷ n+1 = D (y n ψ 1 q 1,n ) = D I 1 ) ψ 1 ψ T 1 M M 1 y n
26 Purification Procedure The purification procedure is simple, at each step the amplitude of the 1st mode is first computed, then removed from the iterated vector y n q 1,n = ψ T 1 My n /M 1, ( ŷ n+1 = D (y n ψ 1 q 1,n ) = D I 1 ) ψ 1 ψ T 1 M M 1 Introducing the sweeping matrix S 1 = I 1 M 1 ψ 1 ψ T 1 M and the modified dynamic matrix D 2 = DS 1, we can write y n ŷ n+1 = DS 1 y n = D 2 y n.
27 Purification Procedure The purification procedure is simple, at each step the amplitude of the 1st mode is first computed, then removed from the iterated vector y n q 1,n = ψ T 1 My n /M 1, ( ŷ n+1 = D (y n ψ 1 q 1,n ) = D I 1 ) ψ 1 ψ T 1 M M 1 Introducing the sweeping matrix S 1 = I 1 M 1 ψ 1 ψ T 1 M and the modified dynamic matrix D 2 = DS 1, we can write y n ŷ n+1 = DS 1 y n = D 2 y n. This is known as matrix iteration with sweeps.
28 Third Mode Using again the idea of purifying the iterated vector, starting with the knowledge of the first and the second eigenpair, with q n,1 as before and ŷ n+1 = D (y n ψ 1 q 1,n ψ 2 q 2,n ) q 2,n = ψ T 2 My n /M 2, substituting in the expression for the purified vector ŷ n+1 = D ( I 1 ψ 1 ψ T 1 M 1 ψ 2 ψ T 2 M ) M } 1 M {{} 2 S 1
29 Third Mode Using again the idea of purifying the iterated vector, starting with the knowledge of the first and the second eigenpair, with q n,1 as before and ŷ n+1 = D (y n ψ 1 q 1,n ψ 2 q 2,n ) q 2,n = ψ T 2 My n /M 2, substituting in the expression for the purified vector ŷ n+1 = D ( I 1 ψ 1 ψ T 1 M 1 ψ 2 ψ T 2 M ) M } 1 M {{} 2 S 1 The conclusion is that the sweeping matrix and the modified dynamic matrix to be used to compute the 3rd eigenvector are S 2 = S 1 1 M 2 ψ 2 ψ T 2 M, D 3 = D S 2.
30 Generalization to The results obtained for the third mode are easily generalised. It is easy to verify that the following procedure can be used to compute all the modes. Define S 0 = I, take i = 1, 1. compute the modified dynamic matrix to be used for mode i, D i = D S i i 2. compute ψ i using the modified dynamic matrix; 3. compute the modal mass M i = ψ T M ψ; 4. compute the sweeping matrix S i that sweeps the contributions of the first i modes from trial vectors, S i = S i 1 1 M i ψ i ψ T i M; 5. increment i, GOTO 1.
31 Generalization to The results obtained for the third mode are easily generalised. It is easy to verify that the following procedure can be used to compute all the modes. Define S 0 = I, take i = 1, 1. compute the modified dynamic matrix to be used for mode i, D i = D S i i 2. compute ψ i using the modified dynamic matrix; 3. compute the modal mass M i = ψ T M ψ; 4. compute the sweeping matrix S i that sweeps the contributions of the first i modes from trial vectors, S i = S i 1 1 M i ψ i ψ T i M; 5. increment i, GOTO 1. Well, we finally have a method that can be used to compute all the eigenpairs of our dynamic problems, full circle!
32 Discussion The method of matrix iteration with sweeping is not used in production because 1. D is a full matrix, even if M and K are banded matrices, and the matrix product that is the essential step in every iteration is computationally onerous, 2. the procedure is however affected by numerical errors, so, after having demonstrated that it is possible to compute all the eigenvectors of a large problem using an iterative procedure it is time to look for different, more efficient iterative procedures.
33 to Inverse iteration is based on the fact that the symmetric stiffness matrix has a banded structure, that is a relatively large triangular portion of the matrix is composed by zeroes. The banded structure is due to the FEM model: in every equation of equilibrium the only non zero elastic force coefficients are due to the degrees of freedom of the few FE s that contain the degree of freedom for which the equilibrium is written.
34 Definition of LU decomposition Every symmetric, banded matrix can be subjected to a so called LU decomposition, that is, for K we write K = L U where L and U are, respectively, a lower- and an upper-banded matrix. If we denote with b the bandwidth of K, we have L = [ { ] i < j l ij with l ij 0 for j < i b and U = [ u ij ] with u ij 0 for { i > j j > i + b
35 Twice the equations? In this case, with w n = M x n, the recursion can be written L U x n+1 = w n or as a system of equations, U x n+1 = z n+1 L z n+1 = w n Apparently, we have doubled the number of unknowns, but the z j s can be easily computed by the procedure of back substitution.
36 Back Substitution Temporarily dropping the n and n + 1 subscripts, we can write z 1 = (w 1 )/l 11 z 2 = (w 2 l 21 z 1 )/l 22 z 3 = (w 3 l 31 z 1 l 32 z 2 )/l 33 z j = (w j l jk z k )/l jj j 1 k=1
37 Back Substitution Temporarily dropping the n and n + 1 subscripts, we can write z 1 = (w 1 )/l 11 z 2 = (w 2 l 21 z 1 )/l 22 z 3 = (w 3 l 31 z 1 l 32 z 2 )/l 33 z j = (w j l jk z k )/l jj j 1 k=1 The x are then given by U x = z.
38 Back Substitution We have computed z by back substitution, we must solve U x = z but U is upper triangular, so we have x N = (z N )/u NN x N 1 = (z N 1 u N 1,N z N )/u N 1,N 1 x N 2 = (z N 2 u N 2,N z N u N 2,N 1 z N 1 )/u N 2,N 2 j 1 x N j = (z N j u N j,n k z N k )/u N j,n j, k=0
39 Back Substitution We have computed z by back substitution, we must solve U x = z but U is upper triangular, so we have x N = (z N )/u NN x N 1 = (z N 1 u N 1,N z N )/u N 1,N 1 x N 2 = (z N 2 u N 2,N z N u N 2,N 1 z N 1 )/u N 2,N 2 j 1 x N j = (z N j u N j,n k z N k )/u N j,n j, k=0 For moderately large systems, the reduction in operations count given by back substitution with respect to matrix multiplication is so large that the additional cost of the LU decomposition is negligible.
40 to Shifts Inverse iteration can be applied to each step of matrix iteration with sweeps, or to each step of a different procedure intended to compute all the eigenpairs, the matrix iteration with shifts.
41 , 1 If we write ω 2 i = µ + λ i, where µ is a shift and λ i is a shifted eigenvalue, the eigenvalue problem can be formulated as K ψ i = (µ + λ i )M ψ i or (K µm)ψ i = λ i M ψ i. If we introduce a modified stiffness matrix K = K µm, we recognize that we have a new problem, that has exactly the same eigenvectors and shifted eigenvalues, K φ i = λ i Mφ i, where φ i = ψ i, λ i = ω 2 i µ.
42 , 2 The shifted eigenproblem can be solved, e.g., by matrix iteration and the procedure will converge to the smallest absolute value shifted eigenvalue and to the associated eigenvector. After convergence is reached, ψ i = φ i, ω 2 i = λ i + µ.
43 , 2 The shifted eigenproblem can be solved, e.g., by matrix iteration and the procedure will converge to the smallest absolute value shifted eigenvalue and to the associated eigenvector. After convergence is reached, ψ i = φ i, ω 2 i = λ i + µ. The convergence of the method can be greatly enhanced if the shift µ is updated every few steps during the iterative procedure using the current best estimate of λ i, λ i,n+1 = ˆx n+1m x n, ˆx n+1 M ˆx n+1 to improve the modified stiffness matrix to be used in the following iterations, K = K λ i,n+1 M
44 , 2 The shifted eigenproblem can be solved, e.g., by matrix iteration and the procedure will converge to the smallest absolute value shifted eigenvalue and to the associated eigenvector. After convergence is reached, ψ i = φ i, ω 2 i = λ i + µ. The convergence of the method can be greatly enhanced if the shift µ is updated every few steps during the iterative procedure using the current best estimate of λ i, λ i,n+1 = ˆx n+1m x n, ˆx n+1 M ˆx n+1 to improve the modified stiffness matrix to be used in the following iterations, K = K λ i,n+1 M Much thought was spent on the problem of choosing the initial shifts, so that all the eigenvectors can be computed in sequence without missing any of them.
45 Quotient for Discrete Systems The matrix iteration procedures are usually used in conjunction with methods derived from the Quotient method. -Ritz Method
46 Quotient for Discrete Systems The matrix iteration procedures are usually used in conjunction with methods derived from the Quotient method. The Quotient method was introduced using distributed flexibilty systems and an assumed shape function, but we have seen also an example where the Quotient was computed for a discrete system using an assumed shape vector. -Ritz Method
47 Quotient for Discrete Systems The matrix iteration procedures are usually used in conjunction with methods derived from the Quotient method. The Quotient method was introduced using distributed flexibilty systems and an assumed shape function, but we have seen also an example where the Quotient was computed for a discrete system using an assumed shape vector. The procedure to be used for discrete systems can be summarized as x(t) = φz 0 sin ωt, ẋ(t) = ωφz 0 cos ωt, 2T max = ω 2 φ T M φ, 2V max = φ T K φ, -Ritz Method
48 Quotient for Discrete Systems The matrix iteration procedures are usually used in conjunction with methods derived from the Quotient method. The Quotient method was introduced using distributed flexibilty systems and an assumed shape function, but we have seen also an example where the Quotient was computed for a discrete system using an assumed shape vector. The procedure to be used for discrete systems can be summarized as x(t) = φz 0 sin ωt, equating the maxima, we have ẋ(t) = ωφz 0 cos ωt, 2T max = ω 2 φ T M φ, 2V max = φ T K φ, -Ritz Method ω 2 = φt K φ φ T M φ = k m, where φ is an assumed shape vector, not an eigenvector.
49 Ritz Coordinates For a N DOF system, an approximation to a displacement vector x can be written in terms of a set of M < N assumed shape, linearly independent vectors, φ i, i = 1,..., M < N and a set of Ritz coordinates z i, i 1,..., M < N: x = i φ i z i = Φ z. -Ritz Method We say approximation because a linear combination of M < N vectors cannot describe every point in a N-space.
50 Quotient in Ritz Coordinates We can write the quotient as a function of the Ritz coordinates, ω 2 (z) = zt Φ T K Φz z T φ T M φz = k(z) m(z), but this is not an explicit function for any modal frequency... -Ritz Method
51 Quotient in Ritz Coordinates We can write the quotient as a function of the Ritz coordinates, ω 2 (z) = zt Φ T K Φz z T φ T M φz = k(z) m(z), but this is not an explicit function for any modal frequency... On the other hand, we have seen that frequency estimates are always greater than true frequencies, so our best estimates are the the local minima of ω 2 (z), or the points where all the derivatives of ω 2 (z) with respect to z i are zero: -Ritz Method ω 2 m(z) k(z) k(z) m(z) (z) z = i z i = 0, for i = 1,..., M < N z j (m(z)) 2
52 Quotient in Ritz Coordinates Observing that k(z) = ω 2 (z)m(z) we can substitute into and simplify the preceding equation, k(z) z i ω 2 (z) m(z) z i = 0, for i = 1,..., M < N -Ritz Method
53 Quotient in Ritz Coordinates Observing that k(z) = ω 2 (z)m(z) we can substitute into and simplify the preceding equation, k(z) z i With the positions we have hence ω 2 (z) m(z) z i = 0, for i = 1,..., M < N K = Φ T K Φ, k(z) = z T Kz = i { k(z) z i } { = 2 j, M = Φ T M Φ k ij z j z i, j k ij z j } = 2Kz, -Ritz Method and, analogously { m(z) z i } = 2Mz.
54 Reduced Eigenproblem Substituting these results in k(z) z i ω 2 (z) m(z) z i = 0 we can write a new eigenvector problem, in the M DOF Ritz coordinates space, with reduced M M matrices: K z ω 2 M z = 0. -Ritz Method
55 Modal Superposition? After solving the reduced eigenproblem, we have a set of M eigenvalues ω 2 i and a corresponding set of M eigenvectors z i. What is the relation between these results and the eigenpairs of the original problem? The ω 2 i clearly are approximations from above to the real eigenvalues, and if we write ψ i = Φz i we see that, being ψ T i Mψ j = z T i ΦT }{{ MΦ} z j = M i δ ij, M -Ritz Method the approximated eigenvectors ψ i are orthogonal with respect to the structural matrices and can be used in ordinary modal superposition techniques.
56 A Last Question One last question: how many ω 2 i and ψ i are effective approximations to the true eigenpairs? Experience tells that an effective approximation is to be expected for the first M/2 eigenthings. -Ritz Method
Multi 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 informationGiacomo Boffi. April 7, 2016
http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 7, 2016 Outline Analysis The process of estimating the vibration
More informationStructural Matrices in MDOF Systems
in MDOF Systems http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 9, 2016 Outline Additional Static Condensation
More informationOutline. Structural Matrices. Giacomo Boffi. Introductory Remarks. Structural Matrices. Evaluation of Structural Matrices
Outline in MDOF Systems Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano May 8, 014 Additional Today we will study the properties of structural matrices, that is the operators that
More informationTruncation Errors Numerical Integration Multiple Support Excitation
Errors Numerical Integration Multiple Support Excitation http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 10,
More informationGeneralized Single Degree of Freedom Systems
Single Degree of Freedom http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March, 8 Outline Until now our were described
More informationReduction in number of dofs
Reduction in number of dofs Reduction in the number of dof to represent a structure reduces the size of matrices and, hence, computational cost. Because a subset of the original dof represent the whole
More informationFREE VIBRATION RESPONSE OF UNDAMPED SYSTEMS
Lecture Notes: STRUCTURAL DYNAMICS / FALL 2011 / Page: 1 FREE VIBRATION RESPONSE OF UNDAMPED SYSTEMS : : 0, 0 As demonstrated previously, the above Equation of Motion (free-vibration equation) has a solution
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationGeneralized Single Degree of Freedom Systems
Single Degree of Freedom http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April, 16 Outline Until now our were described
More information3. Mathematical Properties of MDOF Systems
3. Mathematical Properties of MDOF Systems 3.1 The Generalized Eigenvalue Problem Recall that the natural frequencies ω and modes a are found from [ - ω 2 M + K ] a = 0 or K a = ω 2 M a Where M and K are
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 informationOutline of parts 1 and 2
to Harmonic Loading http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March, 6 Outline of parts and of an Oscillator
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 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 informationABSTRACT Modal parameters obtained from modal testing (such as modal vectors, natural frequencies, and damping ratios) have been used extensively in s
ABSTRACT Modal parameters obtained from modal testing (such as modal vectors, natural frequencies, and damping ratios) have been used extensively in system identification, finite element model updating,
More informationMultiple Degree of Freedom Systems. The Millennium bridge required many degrees of freedom to model and design with.
Multiple Degree of Freedom Systems The Millennium bridge required many degrees of freedom to model and design with. The first step in analyzing multiple degrees of freedom (DOF) is to look at DOF DOF:
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 informationDynamic Response of Structures With Frequency Dependent Damping
Dynamic Response of Structures With Frequency Dependent Damping Blanca Pascual & S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
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 informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationStructural Dynamics A Graduate Course in Aerospace Engineering
Structural Dynamics A Graduate Course in Aerospace Engineering By: H. Ahmadian ahmadian@iust.ac.ir The Science and Art of Structural Dynamics What do all the followings have in common? > A sport-utility
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 informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationD && 9.0 DYNAMIC ANALYSIS
9.0 DYNAMIC ANALYSIS Introduction When a structure has a loading which varies with time, it is reasonable to assume its response will also vary with time. In such cases, a dynamic analysis may have to
More informationANALYSIS OF HIGHRISE BUILDING STRUCTURE WITH SETBACK SUBJECT TO EARTHQUAKE GROUND MOTIONS
ANALYSIS OF HIGHRISE BUILDING SRUCURE WIH SEBACK SUBJEC O EARHQUAKE GROUND MOIONS 157 Xiaojun ZHANG 1 And John L MEEK SUMMARY he earthquake response behaviour of unframed highrise buildings with setbacks
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi.
Periodic and Non-periodic http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 7, 218 Outline Undamped SDOF systems
More informationDesign of Structures for Earthquake Resistance
NATIONAL TECHNICAL UNIVERSITY OF ATHENS Design of Structures for Earthquake Resistance Basic principles Ioannis N. Psycharis Lecture 3 MDOF systems Equation of motion M u + C u + K u = M r x g(t) where:
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 informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 17 AA242B: MECHANICAL VIBRATIONS Solution Methods for the Generalized Eigenvalue Problem These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical
More informationME scope Application Note 28
App Note 8 www.vibetech.com 3/7/17 ME scope Application Note 8 Mathematics of a Mass-Spring-Damper System INTRODUCTION In this note, the capabilities of ME scope will be used to build a model of the mass-spring-damper
More informationA pragmatic approach to including complex natural modes of vibration in aeroelastic analysis
A pragmatic approach to including complex natural modes of vibration in aeroelastic analysis International Aerospace Symposium of South Africa 14 to 16 September, 215 Stellenbosch, South Africa Louw van
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 informationAA 242B / ME 242B: Mechanical Vibrations (Spring 2016)
AA 242B / ME 242B: Mechanical Vibrations (Spring 206) Solution of Homework #3 Control Tab Figure : Schematic for the control tab. Inadequacy of a static-test A static-test for measuring θ would ideally
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More informationPROJECT 1 DYNAMICS OF MACHINES 41514
PROJECT DYNAMICS OF MACHINES 454 Theoretical and Experimental Modal Analysis and Validation of Mathematical Models in Multibody Dynamics Ilmar Ferreira Santos, Professor Dr.-Ing., Dr.Techn., Livre-Docente
More informationIntroduction to structural dynamics
Introduction to structural dynamics p n m n u n p n-1 p 3... m n-1 m 3... u n-1 u 3 k 1 c 1 u 1 u 2 k 2 m p 1 1 c 2 m2 p 2 k n c n m n u n p n m 2 p 2 u 2 m 1 p 1 u 1 Static vs dynamic analysis Static
More informationDamping Matrix. Donkey2Ft
1 Damping Matrix DonkeyFt Damping in a single-degree-of-freedom (SDOF) system is well studied. Whether the system is under-damped, over-damped, or critically damped is well known. For an under-damped system,
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
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 informationRational Krylov methods for linear and nonlinear eigenvalue problems
Rational Krylov methods for linear and nonlinear eigenvalue problems Mele Giampaolo mele@mail.dm.unipi.it University of Pisa 7 March 2014 Outline Arnoldi (and its variants) for linear eigenproblems Rational
More informationTransient Response Analysis of Structural Systems
. 21 Transient Response Analysis of Structural Systems 21 1 Chapter 21: TRANSIENT RESPONSE ANALYSIS OF STRUCTURAL SYSTEMS 21 2 21.1 MODAL APPROACH TO TRANSIENT ANALYSIS Consider the following large-order
More informationSolution of Vibration and Transient Problems
. 9 Solution of Vibration and Transient Problems 9 Chapter 9: SOLUTION OF VIBRATION AND TRANSIENT PROBLEMS 9 9. MODAL APPROACH TO TRANSIENT ANALYSIS Consider the following large-order finite element model
More informationChapter 4 Analysis of a cantilever
Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3
ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3 ISSUED 24 FEBRUARY 2018 1 Gaussian elimination Let A be an (m n)-matrix Consider the following row operations on A (1) Swap the positions any
More informationIdentification of Damping Using Proper Orthogonal Decomposition
Identification of Damping Using Proper Orthogonal Decomposition M Khalil, S Adhikari and A Sarkar Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk
More informationThe quantum mechanics approach to uncertainty modeling in structural dynamics
p. 1/3 The quantum mechanics approach to uncertainty modeling in structural dynamics Andreas Kyprianou Department of Mechanical and Manufacturing Engineering, University of Cyprus Outline Introduction
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 12: Gauss for Linear Systems Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationDEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular
form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationGiacomo Boffi. Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano
http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 21, 2017 Outline of Structural Members Elastic-plastic Idealization
More informationDYNAMICS OF MACHINERY 41514
DYNAMICS OF MACHINERY 454 PROJECT : Theoretical and Experimental Modal Analysis and Validation of Mathematical Models in Multibody Dynamics Holistic Overview of the Project Steps & Their Conceptual Links
More informationStochastic Dynamics of SDOF Systems (cont.).
Outline of Stochastic Dynamics of SDOF Systems (cont.). Weakly Stationary Response Processes. Equivalent White Noise Approximations. Gaussian Response Processes as Conditional Normal Distributions. Stochastic
More informationSample Problems on Quantum Dynamics for PHYS301
MACQUARIE UNIVERSITY Department of Physics Division of ICS Sample Problems on Quantum Dynamics for PHYS30 The negative oxygen molecule ion O consists of a pair of oxygen atoms separated by a distance a
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationEigenvalue Problems CHAPTER 1 : PRELIMINARIES
Eigenvalue Problems CHAPTER 1 : PRELIMINARIES Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Preliminaries Eigenvalue problems 2012 1 / 14
More informationThe Eigenvalue Problem: Perturbation Theory
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 13 Notes These notes correspond to Sections 7.2 and 8.1 in the text. The Eigenvalue Problem: Perturbation Theory The Unsymmetric Eigenvalue Problem Just
More information( ) Chapter 3: Free Vibration of the Breakwater. 3.1 Introduction
Chapter : Free Vibration of the Breakwater. Introduction Before conducting forcing analyses on the breakwater, a free vibration study is necessary. This chapter will describe the analysis conducted to
More informationTWO-STAGE ISOLATION FOR HARMONIC BASE EXCITATION Revision A. By Tom Irvine February 25, 2008
TWO-STAGE ISOLATION FOR HARMONIC BASE EXCITATION Revision A By Tom Irvine Email: tomirvine@aol.com February 5, 008 Introduction Consider a base plate mass m and an avionics mass m modeled as two-degree-of-freedom.
More informationSHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine May 24, 2010
SHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine Email: tomirvine@aol.com May 4, 010 Introduction The primary purpose of this tutorial is to present the Modal Transient method
More informationLINEAR ALGEBRA SUMMARY SHEET.
LINEAR ALGEBRA SUMMARY SHEET RADON ROSBOROUGH https://intuitiveexplanationscom/linear-algebra-summary-sheet/ This document is a concise collection of many of the important theorems of linear algebra, organized
More informationOperating Deflection Shapes from Strain Measurement Data
Operating Deflection Shapes from Strain Measurement Data Timothy G. Hunter, Ph.D., P.E. President Wolf Star Technologies, LLC 3321 N. Newhall St., Milwaukee, WI 53211 Abstract Strain gauges are often more
More informationAppendix C. Modal Analysis of a Uniform Cantilever with a Tip Mass. C.1 Transverse Vibrations. Boundary-Value Problem
Appendix C Modal Analysis of a Uniform Cantilever with a Tip Mass C.1 Transverse Vibrations The following analytical modal analysis is given for the linear transverse vibrations of an undamped Euler Bernoulli
More informationLecture 4 Eigenvalue problems
Lecture 4 Eigenvalue problems Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationEffect of Mass Matrix Formulation Schemes on Dynamics of Structures
Effect of Mass Matrix Formulation Schemes on Dynamics of Structures Swapan Kumar Nandi Tata Consultancy Services GEDC, 185 LR, Chennai 600086, India Sudeep Bosu Tata Consultancy Services GEDC, 185 LR,
More informationAPPLICATIONS OF PURE AND COMBINED BUCKLING MODE CALCULATION OF THIN-WALLED MEMBERS USING THE FINITE ELEMENT METHOD
SDSS Rio 2010 STABILITY AND DUCTILITY OF STEEL STRUCTURES E. Batista, P. Vellasco, L. de Lima (Eds.) Rio de Janeiro, Brazil, September 8-10, 2010 APPLICATIONS OF PURE AND COMBINED BUCKLING MODE CALCULATION
More informationTheoretical Basis of Modal Analysis
American Journal of Mechanical Engineering, 03, Vol., No. 7, 73-79 Available online at http://pubs.sciepub.com/ajme//7/4 Science and Education Publishing DOI:0.69/ajme--7-4 heoretical Basis of Modal Analysis
More informationTuning TMDs to Fix Floors in MDOF Shear Buildings
Tuning TMDs to Fix Floors in MDOF Shear Buildings This is a paper I wrote in my first year of graduate school at Duke University. It applied the TMD tuning methodology I developed in my undergraduate research
More informationOutline. Introduction
Outline Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 5, 4 A periodic loading is characterized by the identity p(t) = p(t + T ) where T is the period
More informationDynamics of Structures
Dynamics of Structures Elements of structural dynamics Roberto Tomasi 11.05.2017 Roberto Tomasi Dynamics of Structures 11.05.2017 1 / 22 Overview 1 SDOF system SDOF system Equation of motion Response spectrum
More informationThe Finite Element Method
Page he Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. o. Dr. Michael Havbro bo Faber Dr. Nebojsa Mojsilovic Swiss Federal Institute of EH Zurich, Switzerland Contents of
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationAssignment 6. Using the result for the Lagrangian for a double pendulum in Problem 1.22, we get
Assignment 6 Goldstein 6.4 Obtain the normal modes of vibration for the double pendulum shown in Figure.4, assuming equal lengths, but not equal masses. Show that when the lower mass is small compared
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More information1792. On the conditions for synchronous harmonic free vibrations
1792. On the conditions for synchronous harmonic free vibrations César A. Morales Departamento de Mecánica, Universidad Simón Bolívar, Apdo. 89000, Caracas 1080A, Venezuela E-mail: cmorales@usb.ve (Received
More informationSpectral methods for fuzzy structural dynamics: modal vs direct approach
Spectral methods for fuzzy structural dynamics: modal vs direct approach S Adhikari Zienkiewicz Centre for Computational Engineering, College of Engineering, Swansea University, Wales, UK IUTAM Symposium
More informationMath 471 (Numerical methods) Chapter 4. Eigenvalues and Eigenvectors
Math 471 (Numerical methods) Chapter 4. Eigenvalues and Eigenvectors Overlap 4.1 4.3 of Bradie 4.0 Review on eigenvalues and eigenvectors Definition. A nonzero vector v is called an eigenvector of A if
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi. March 25, 2014
Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 25, 2014 Outline Introduction Fourier Series Representation Fourier Series of the Response Introduction
More informationNormal modes. where. and. On the other hand, all such systems, if started in just the right way, will move in a simple way.
Chapter 9. Dynamics in 1D 9.4. Coupled motions in 1D 491 only the forces from the outside; the interaction forces cancel because they come in equal and opposite (action and reaction) pairs. So we get:
More informationEmpirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems
Empirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems Antoni Ras Departament de Matemàtica Aplicada 4 Universitat Politècnica de Catalunya Lecture goals To review the basic
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationThe QR Decomposition
The QR Decomposition We have seen one major decomposition of a matrix which is A = LU (and its variants) or more generally PA = LU for a permutation matrix P. This was valid for a square matrix and aided
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationLab #2 - Two Degrees-of-Freedom Oscillator
Lab #2 - Two Degrees-of-Freedom Oscillator Last Updated: March 0, 2007 INTRODUCTION The system illustrated in Figure has two degrees-of-freedom. This means that two is the minimum number of coordinates
More information5.3 The Power Method Approximation of the Eigenvalue of Largest Module
192 5 Approximation of Eigenvalues and Eigenvectors 5.3 The Power Method The power method is very good at approximating the extremal eigenvalues of the matrix, that is, the eigenvalues having largest and
More informationMath 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Exam 2 will be held on Tuesday, April 8, 7-8pm in 117 MacMillan What will be covered The exam will cover material from the lectures
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationMathematical Properties of Stiffness Matrices
Mathematical Properties of Stiffness Matrices CEE 4L. Matrix Structural Analysis Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 0 These notes describe some of the
More informationChapter 23: Principles of Passive Vibration Control: Design of absorber
Chapter 23: Principles of Passive Vibration Control: Design of absorber INTRODUCTION The term 'vibration absorber' is used for passive devices attached to the vibrating structure. Such devices are made
More informationTHE subject of the analysis is system composed by
MECHANICAL VIBRATION ASSIGNEMENT 1 On 3 DOF system identification Diego Zenari, 182160, M.Sc Mechatronics engineering Abstract The present investigation carries out several analyses on a 3-DOF system.
More informationSTRUCTURAL DYNAMICS BASICS:
BASICS: STRUCTURAL DYNAMICS Real-life structures are subjected to loads which vary with time Except self weight of the structure, all other loads vary with time In many cases, this variation of the load
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
More informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
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 informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationIn this section again we shall assume that the matrix A is m m, real and symmetric.
84 3. The QR algorithm without shifts See Chapter 28 of the textbook In this section again we shall assume that the matrix A is m m, real and symmetric. 3.1. Simultaneous Iterations algorithm Suppose we
More informationEIGENVALUE PROBLEMS. Background on eigenvalues/ eigenvectors / decompositions. Perturbation analysis, condition numbers..
EIGENVALUE PROBLEMS Background on eigenvalues/ eigenvectors / decompositions Perturbation analysis, condition numbers.. Power method The QR algorithm Practical QR algorithms: use of Hessenberg form and
More information