Seminar Assignments Structural Dynamics Assignments
Seminar Assignments Structural Dynamics Assignments
Seminar Assignments Structural Dynamics Assignments
Assignment 1 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on February 15, 2015).
Derive the steady-state response of SDOF (Single Degree of Freedom System) whose
dynamic behavior is governed by the following equations:
d 2x dx
m 2
+ c =
+ kx f p d ( t nt0 ) t > 0, t0 > 0
dt dt n =0
x ( 0 ) = 0,
dx = v0
dt t =0
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_1_2015_YourName.pdf. For example:
Assignment_1_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
(1)
The solution consists of a homogeneous and a particular part, no initial conditions are applied yet1.
(2)
(3)
With
, , (4)
However, in this assignment, only the steady-state response is asked. When taking the limit of time
going to infinity of any particular solution, it changes to the unique steady-state response.
To rewrite the particular solution the sum rule in integration is used, i.e. the integral is split into a
sum of three integrals.
(5)
1
See also slides of lecture 3, slide 9
(6)
(7)
Because of the principle of superposition, the integral of a sum of Dirac delta pulses is the same as a
sum of integrals of Dirac delta pulses, so the first part of the integral can be rewritten to.
(8)
(10)
Assignment 2 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on February 18, 2015).
Consider a mass-spring-dashpot system shown in Figure 1. The left ends of the spring
and of the dashpot are attached to a movable rigid wall, whose displacement in the
horizontal direction is given by u ( t ) .
u(t)
k x(t)
x
Figure 1.
t
2 3
Figure 2.
Questions:
1. Express the time dependence of u ( t ) on time analytically in accordance with Figure 2.
2. Derive an equation of motion that governs small horizontal displacement of the mass
caused by the motion of the wall.
3. Derive an analytical expression for the displacement of the mass x ( t ) assuming that
at t = 0 both the deflection and velocity of the mass are zero (an expression in the
form of an integral would suffice).
4. Plot the normalized displacement of the mass x = x u0 as function of time using the
following figures for the system n
parameters: = =
k m 10 rad s ,
= c 2= (
=
km 0.3, )
1 n . To plot the dependence, MAPLE or any other
suitable software may be used.
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_2_2015_YourName.pdf. For example:
Assignment_2_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
3. First the equation of motion is rewritten to have less parameters. Using the natural frequency
and the damping factor , the following is obtained:
The homogeneous solution to the rewritten equation of motion describing a free motion is1:
with
Assuming zero initial conditions, the homogeneous solution becomes zero and the total solution is
equal to the particular solution. Splitting up in three domains and substituting in the force
with u(t') for different domains results in:
For
For
3 0)) sin 1
1
In correspondence with Slides lecture 3, slide 9 where n = n
For
>
>
>
>
>
>
Assignment 3 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on February 21, 2015).
Derive analytically the steady-state solution of the equation below. It is mandatory to use
the complex-valued form of the solution along the lines explained on Slide 18 of Lecture
2&3 and on page 36-37 of the lecture notes to CT 4145.
d 3x d 2x dx
a3 3
+ a2 2
F sin (t ) ,
+ a1 + a0 x =
dt dt dt
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_3_2015_YourName.pdf. For example:
Assignment_3_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
(3)
(4)
(5)
(6)
Thus, an expression for the complex amplitude is found which can be filled in xpart:
(11)
(12)
To get the imaginary part from the denominator to the numerator one can multiply with the complex
conjugate of the denominator after rearranging into the form
(13)
(14)
(15)
Now the equation is in such a form that the imaginary part can be taken out (i.e. )
(16)
(17)
With:
(18)
(19)
Since the steady-state solution becomes xpart in the limit of t to infinity, xsteady-state = xpart.
Assignment 4 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on February 25, 2015).
Assume that a structure may perform coupled vertical-angular vibrations in wind and the
equations of motion governing these vibrations are given as
1 2
mz( t ) + kz ( t ) = v BLC3 ( t )
2
1
I ( t ) + k ( t ) = B 3 LC7 z ( t )
2
= =
m 10 kg, I 5 kg=
m 2 , k 1000 N =
m , k 2000 N m
=B 5=
m, L 10 m
= 1.29
= kg m3 , C3 1
Question.
Find the range of magnitudes of the coefficient C7 , which would assure that the structure
is stable at wind velocities smaller than 0.3m s .
Use of MAPLE or any other symbolic software is not permitted.
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_4_2015_YourName.pdf. For example:
Assignment_4_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
(6)
(7)
This is an algebraic equation which can be solved for with the abc-formula.
with , a = 50, b = , c = 2 000 000 (12)
with (13)
Much depends on the value of the discriminant, which can split in three domains. For D = 0 the roots of p can
be found with the abc-formula again, after first rewriting D.
(15)
(16)
(17)
which gives p1 = 1.731 and p2 = 0.192
Now the stability can be analyzed. The system is stable if the real part of is negative or zero.1 For
the various values of the discriminant the roots p can be filled in the characteristic equation and the
real part of can be determined.
1
Actually, the real part of sn.
From the graph the value of D for different values of p can be seen.
(18)
(19)
(20)
So for stability p must be smaller or equal to 0.192. Given that v is smaller than 0.3 m/s, the
coefficient C7 can be found:
(21)
(22)
(23)
Assignment 5 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on 28th February 2015).
Consider a square-shaped block of mass M as depicted in the figure below. Assume that
it can move only in the plane of the figure, i.e. it can rotate around the x3 -axis and move
along the x1 and x2 axes.
a
k2 a 4
a x2
M
x1 x3
X 0 cos ( t )
k1 c k1
The block is supported by three springs and one dashpot as shown in the figure. A point
mass M is attached to the lower left corner of the block, as shown in the figure. The
mass of the body is the same as that of the block.
Questions.
1. Derive equations of motion that govern small vibrations of the block+mass
system with respect to the centre of inertia of this system.
2. Assuming the following numerical values for the system parameters:
= =
M 20000 kg, k1 50000 = N m , k2 30000 = N m , a 2 m,
find the natural frequencies of the undamped system, its eigenmatrix and the
modal mass and stiffness matrices. (MAPLE may be used to find these).
3. Derive the amplitude-frequency response matrix assuming that c = 400 Ns m .
4. Plot the amplitude of the horizontal, x2 ( t ) , steady-state displacement of the block
as a function of the frequency of the external moment.
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_5_2015_YourName.pdf. For example:
Assignment_5_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
(1)
(2)
(3)
(6)
In matrix form
Taking a coordinate system in the centre of inertia, for the polar moment of inertia J one finds:
(7)
>
>
>
>
>
>
>
>
3. >
C* is not diagonal, so the modal analysis won't work. The damping can still be assumed to be
uncoupled if it is uniformly distributed and not too large, but that could be a farfetched assumption.
Another option is to consider proportional damping, but then parameters a0 and a1 must be known.
The amplitude-frequency response matrix HxF cannot be found with given information.
4. No plot available.
Assignment 6 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on March 10, 2015).
Consider a simply supported beam, whose motion is governed by the following set of
equations:
2w 4w
A 2 + EI 4 f ( x, t ) ,
= 0 < x < L, 0<t <
t x
x
f sin
f ( x, t ) = , 0 < t <t
L
0, t >t
w ( 0, x=
) 0, w t ( x, t )t ==
0
0
w (t, 0) =
=
2 2
x 0=
w x
2 2
x L
( x, t ) w (t, L ) =
= w x ( x, t ) =
0
Question.
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_6_2015_YourName.pdf. For example:
Assignment_6_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
(2)
Initial conditions:
(3)
(4)
So there is zero initial displacement and initial velocity.
Boundary conditions:
(5)
(6)
So the beam has zero deflection and zero moments at both ends, and thus is pinned at both ends.
First the free vibration must be analyzed, to obtain the natural frequencies and modes of the system.
The (free vibration form of the) equation of motion can be rewritten to:
with (7)
(23a)
(24a)
(23b)
(24b)
Summing (23b) and (24b) results in , which means C3 = 0.
Now the external force can be taken into account, split in two time domains.
0<t<
It is attempted to find the solution in the form of the normal modes
(27)
(33)
Everything which does not depend on time can be put in front of the Duhamel integral
(34)
(35)
To find and equation (35) is substituted in equation (27) and then the initial conditions
(3) and (4) are applied.
(36)
(37)
(40)
(41)
(42)
Now equation (27) can be written, with equation (35) substituted in it with and known:
(43)
(44)
The integral is only non-zero and equal to L/2 for , so when n = 1, the solution is:
(45)
with and
t>
The moment t is larger than the forcing is zero and the equation of motion for free vibrations (7)
holds:
The initial conditions are found at the end of the first time domain, for t = .
(47)
(48)
(51)
(52)
with and
Assignment 7 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on March 13, 2015).
Consider a string of the length L as shown in the figure below. The left end of the string
is fixed, whereas its right end is attached to a mass-spring-dashpot element, as shown in
the figure. The string is subject to a harmonic force which is uniformly distributed over
the interval 0 x L 2 . At the later interval the string is supported by a uniform elastic
foundation as shown in the figure below.
p sin ( t )
, A, T
kd k c0
x
0 L 2 L
1. Formulate governing equations to the problem (an equation of motion for the string,
the boundary and interface conditions).
2. Derive the steady-state response of the string to the applied force (solution in the form
of a matrix equation would suffice).
3. Using the following numerical values for the system parameters
= 8000 kg=
m3 , A 5 = =
mm 2 , T 500 N, L 1m
= =
m 0.2 =
kg, k 100 N =
m, c0 600 kg s , kd 200 N m 2
p = 40 N m
plot the amplitude of the string displacement at the point x = L / 3 as a function of the
frequency of the external force. Use of MAPLE is advisable here.
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_7_2015_YourName.pdf. For example:
Assignment_7_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
with , (1)
Left of L/2 will be called w1 and to the right of L/2 will be called w2 for clarity.
Interface conditions at x=L/2:
(4)
(5)
(7)
(8)
(9)
(10)
(11)
with (14)
(16)
Substituting the general solution in the boundary conditions and initial conditions gives:
Boundary condition 1:
(17)
Boundary condition 2:
(18)
Interface condition 1:
(21)
Interface condition 2:
(22)
Assignment 8 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on March 19, 2015).
Consider a structure shown in the figure below. It consists of an axially tensioned beam
and a string. The beam occupies the interval 0 x L 2 , whereas the string occupies the
interval L 2 x L . At x = L 2 the beam and the string are fixed to each other. The left
end of the beam and the right end of the string are fixed.
M cos ( t ) P sin ( t )
m string
, Ab , E , I , T , As , T
beam c0
x
0 L 4 L 2 3L 4 L
The interface between the beam and the string ( x = L 2 ) is supported by a dashpot and a
point-mass is attached to the system at this point. The system is loaded by an external
moment applied at x = L 4 and an external force applied at the interface.
Questions.
1. Formulate a system of equations that govern small vibration of the system. These
equations should include an equation of motion of the beam, that of the string,
boundary conditions at x = 0 and x = L , and interface conditions at x = L 4 and
x = L 2.
2. State a form in which the steady-state response of the system can be found and, by
substituting this form into the governing equations, rewrite those such that no
time is present in the formulation.
3. Write an expression that gives the amplitude of the steady-state displacement of
the beam at the point x = 3L / 4 .
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_8_2013_YourName.pdf. For example:
Assignment_8_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
(1)
1 and 2: (2)
3: (3)
Boundary conditions:
x=0 (4)
(5)
x=L (6)
Interface conditions:
x = L/4 (7)
(8)
(9)
(10)
x = L/2 (11)
(12)
(13)
3: (16)
(18)
x=L (19)
(21)
(22)
(23)
x = L/2 (24)
(25)
(26)
Note:
(28)
(29)
with
These 3 equations have 10 unknowns, and there are also 10 equations (3 boundary conditions, 4
interface conditions and 3 interface conditions), creating a inhomogeneous (nonzero right hand side)
system of equations which can be solved.
Assignment 9 (bonus 0.1; to get the bonus, please deliver your solution not later than
at 23.55 on March 25, 2013).
Consider a semi-infinite string at the interval < x 0 . The end of the string is
(1)
attached to a dashpot with the damping coefficient cdp . Additionally, at x = 3L , the
( 2)
string is supported by another dashpot with the damping coefficient cdp .
Assume that an incident pulse propagates towards the end of the string. The string
deflection in this pulse is given as
0, x ct < 6 L
( x=
, t )t < L c w0 , 6 L < x ct < 4 L
0, 4 L < x ct
where c is the wave speed in the string. Thus, the situation at t = 0 is as shown in the
figure below
w
w0
c
( )
2 ()
1
cdp cdp
x
7L 6L 5L 4L 3L 2L L 0
Question.
Derive an expression for and plot the string deflection w ( x, t ) at t = 12 L c assuming that
(1) ( 2)
=cdp T=
c , cdp 3T c where T is the tension in the string.
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_9_2013_YourName.pdf. For example:
Assignment_9_2013_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
The damping coefficients are: cdp,1 = T/c which corresponds to non-reflective and cdp,2 = 3T/c.
Initial conditions:
(2)
(3)
Boundary conditions:
x=0
(4)
Interface conditions:
x = -3L
(5)
(6)
The values of the dashpots absorption are: cdp,1 = T/c and cdp,2 = 3T/c. Their influence on the wave
must be investigated.
Assuming a solution in the form of the modified D' Alembert solution form:
(7)
In the first interface condition (5) the left and right part are, using :
(incident and reflected) (8)
(transmitted) (9)
Since
(11)
From the first interface conditions one can write the reflected wave in terms of the incident and
transmitted wave :
(12)
Assuming 1
(13)
(14)
So 2/5 of the incident wave will be transmitted and 3/5 will be reflected.
Now looking at the boundary condition at x = 0, one sees this is standard viscous boundary.
Formulating and
This means the wave is not reflected at all, at the boundary at x = 0.
1
valid for
L
i
W
>
<
-3
UK
-
IA/.
L
<
-3'
\
tl
-3 L
No
c
t
L (1
ff
67 L
"TT'
i
i
4
.-
4- c
>
-3
l
i--'
>
>
Hh 7
0
-1' t L -f
*>
>
Assignment 10 (bonus 0.1; to get the bonus, please deliver your solution not later
than at 23.55 on March 31st , 2015).
Consider an infinite beam that performs bending vibrations. A point mass m with mass
moment of inertia J and a rotational spring (this spring resists to the rotational beam
motion only and does not restrain the vertical displacement) are attached to the beam at
x = 0 as shown in the figure below.
m, J
w0 cos ( t x )
kr
x
0
Questions.
1. Formulate governing equations to the problem (an equation of motion for the beam
and interface conditions at x = 0 ).
2. Derive the steady-state displacement of the beam caused by the incident wave.
3. Using the following numerical values for the system parameters
= =
m 10 3
kg, kr 106 N m, J=104 N m 2
3
= =
w0 10 m, L 10 m.
plot the amplitude of the beam displacement at the point x = 3L 2 as a function of the
frequency of the incident wave. Use of MAPLE is advisable to answer this
question.
Your solution should be submitted electronically in the pdf format and the file name
should be of the following form: Assignment_10_2015_YourName.pdf. For example:
Assignment_10_2015_AVMetrikine.pdf
The solution may be handwritten but then scanned (in a good quality), converted to the
pdf format and sent as ONE file.
Solutions which will be delivered differently wont be considered.
Interface conditions:
x=0
(2)
(3)
(4)
(5)
2. The incident wave will cause a reflection and a transmission at the interface. The reflected and
transmitted wave should have the same frequency as the incident wave. The total displacement is
the superposition of the reflected and transmitted waves:
(6)
(8)
(9)
From equation (10) the wavenumber of the wave that is excited in the beam must be found.
with (11)
and or (12)
or and or (13)
Because the incident wave is of a cosine form the real part of the obtained result must be taken.