Matrix Method for Vibration Analysis
Uploaded by
ScribdTranslationsMatrix Method for Vibration Analysis
Uploaded by
ScribdTranslationsMATRIX ITERATION
1 CONCEPT.-It is an iterative procedure that leads to the main modes of vibration.
of a system and its natural frequencies.
The movements of the masses are calculated approximately and based on this.
displacements the matrix equation of the system is written. The influence coefficients of the
systems are replaced in the matrix equation, which is then developed. The
Normalization of the displacement and the development of the matrix are repeated. This process continues.
until the first mode repeats, with the degree of accuracy desired.
For the higher natural modes and frequencies, the principle of orthogonality is used with
the end of obtaining a new matrix equation, which is free from any lower order.
the iterative procedure is repeated.
In this work 'Matrix method to calculate natural frequencies and normal modes of vibration
of deformable systems" we have developed a method that allows calculating the modes
normal modes of vibration using the dynamic matrix [D] and by explaining it through a case
typical, we have estimated the errors of the method, which turn out to be significant (and not satisfactory) to
as the order of the mode increases.
An influence coefficient, denoted by a12is defined as the static deflection of the system e
position 1, due to a unit force applied at position 2 when this unit force
it is the only force that acts. Therefore, the method of influence coefficients is
advantageous for calculating all the deflections produced by the applied forces and for
establish the differential equations of motion of the system.
It can be proven that the following expression is true:
a ij= deflection in the position due to a unit force applied in the proportion j.
ahideflection in the position due to a unit force applied in proportion.
2 STATEMENT OF THE EXERCISE
3 MECHANICAL IMPEDANCE METHOD
FK= K (X + γcos) X1X2X3
∑K γcos = ∑W = 0
∑F = m X1°°
m X1°°equal to -4KX1–3K(X1–X2)
m X1°°+ 4KX1+ 3K (X1-X2) = 0
m X1°°+ 7KX13KX2= 0 I
∑F = m X2°°
m X2°°= 3K(X1–X2)–2K(X2–X3)
m X2°°- 3Kx1+ 3KX2+ 2KX22KX3= 0
m X2°°+ 5KX2-3KX1-2KX3= 0 II
∑F = m X3°°
m X3°°= 2K(X2–X3)
m X3°°+ 2KX3-2KX2= 0 III
MIM: XJ= XJeiWt; X°J= iWXJeiWt ; X°°J = - W2XJeiWt
-m W2X1eiWt+ 7KX1eiWt-3 KX2eiWt= 0 ÷ eiWt
(7K–m W2)X1-3KX2= 0 I
- M W2X2eiWt+5 KX2eiWt-3 KX1eiWt-2 KX3eiWt= 0 ÷ eiWt
-3KX1+ (5K–m W2)X2–2KX3= 0 II
-M W2X3+ 2 KX3eiWt-2KX3eiWt= 0 ÷ eiWt
-2KX2+ (2K–m W2)X3= 0 III
=
Δdet.(7K–m W2) (5K–m W2) (2K–m W2) + 0 + 0–4K2(7K–m W2) -9K2(2K–m W2)
Δdet.(35 K2-7K m W2-5 m W2+ m2W4) (2K–m W2)–4K2(7K–m W2)–9K2(2K–m W2)
Δdet.70K3-24K2m W2 + 2 km2W4-35K2m W2+ 12K m2W4–m3W6-4K2(7K–m W2)–
9K2……...... …….(2K –m W2)
Δdet.= - m3W6+ 14K m2W4-59K2m W2+ 70 K3-28K3+ 4K2m W2-18K3+ 9K2m W2
Δdet.-m3W6+ 14K m2W4-46K2m W2+ 24K3
Δdet.= 0
-m3W6+ 14K m2W4–46 K2m W2+ 24K3= 0
-m3(W2)3+ 14K m2(W2)2-46K2m (W2)1+ 24K3= 0
W2n1= 9.36 K/m W1= 3.06 √K/m
W2n2= 4 K/m W2= 2 √K/m
W2n3= 0.64 K/m W3= 0.8 √K/m
= W2m/K
X1= 0 ; X2= 2 ; X3= 3
4 MATRIX ITERATION METHOD
Mass analysis 1
Sum of Forces in the X direction1=0 Sum of Forces in the X Direction2=0
4KX11-3K(X11X21) +F = 0 3K(X11–X21)–2K(X21–X31) = 0
-7KX11+ 3K X21= - F 3KX11-3KX21-2KX21+ 2KX31= 0
7KX113K X21 = 1 I 3KX11+ 5KX21+ 2KX31= 0 II
Sum of Forces3=0
2K(X21–X31)
2KX21-2KX31= 0 III
=
=1
X11= ¼
X21= ¼
X31= ¼
Mass analysis 2 m 2; J = 2 x2x1>x3 Note: 1.1 2.2 3.3 = Positive
Sum of Forces1=0 Sum of FX2=0
-4KX123K(X22–X12) = 0 -3K(X22–X12)–2K(X22–X32) + F = 0
7KX12-3KX22= 0 -5X22+ 3KX12+ 2KX32-1
5KX22-3KX12-2KX32= 1
∑FX3=0
2K(X22–X32) = 0
2KX22-2KX32= 0
= =
X12= ¼
X22= 7/12
X32= 7/12
Mass Analysis 3 m 1J = 3 x3>x2x1
Sum of Forces1=0 sum of forces in the X direction2=0
-4KX13+3K(X23-X13) = 0 -3K(X23-X13+ 2K(X33–X23) = 0
-7KX13+ 3KX23= 0 -5KX23+ 3KX13+ 2KX33= 0
Sum of FX3=0
-2K(X33–X23) + F = 0
-2KX33+ 2KX23= -F
-2K23+ 2KX33= F
X13= ¼
X23= 7/12
X33= 1.0833
