Engineering Mechanics
Force Vectores and Resultants
Resultant Force
Resulatant Moments
Resultant of a distribution Force
Normally force are ditributed along the structure, but can be simplified into a point
To do this we use Integrals
𝑭𝑹 = Integral(𝒘(𝒙)) 𝒅𝒙
M=Fd
The force has to perpendicular to the line of action
Right hand rules
F=ma
These can be broken down into X and Y components. Then use trig to find the resultant
If the force is parraller to the line of action then the moments is 0
M=Fd
M=F*0
M=0
Equilibrium
EQUILIBRIUM OF A PARTICLE
This is when the resultant force of a particle is zero. This means that the formula is Fr=0
Meaning that the Fx=0 and Fy=0
The resultant force being 0 means that the particle is not experiencing any acceleration. This could mean that the object is not moving or that he object is moving at a constant velocity.
A free body diagram is used to see cleary where are all the forces coming from and how they are interacting with the particle.
useful applications
Opposed to this, larger objects (beams etc) can not be represented as particles. This kind of objects are called rigid bodies
These analysis often requie us to look at tension on string and forces on spring. So it good to also know Hook's law, F=kx
Kinematics
Kinematics of Particles
When a particle is only moving in one direction (let’s say along s-axis), it is said to be in rectilinear motion
With kinematics, there are four main variables that we will be working with. Velocity, displacement, time and acceleration.
The relations between these quantities are described as time derivatives
If our motion is not in just one direction but two directions (let’s say x and y), then our motion is said to be curvilinear
In these cases, both our velocity and acceleration become Pythagorean products
of their x- and y-components, so their magnitudes are:
Direction angles can be calculated from x- and y-components using inverse
tangent:
The most common application for this would be for projectiles but the (x,y)-component method works well for cases, where the curvilinear nature of motion is a product of things happening in both directions – i.e. the motion must not be curvilinear by definition
If our situation is different in such a way that the curvilinear nature is forced because of our particle is moving along a path that is known beforehand, then it is usually simpler to perform our calculations in a plane which consists of normal and tangential components (n and t)
Here at represents the acceleration caused by change in the magnitude of the velocity and an represents the acceleration caused by change in the direction of the velocity
Tangential and normal components of the acceleration can be easily calculated
using the following formulas
Tangential acceleration is the time derivative of velocity
Normal acceleration is velocity squared
and divided by radius of curvature ρ
If our path is circular, radius of curvature is constant (= radius of circle)
If our path can be expressed as y = f(x), the radius of curvature is
Third possible coordinate system that we can use is the system of cylindrical coordinates (or, in plane, polar coordinates r and θ )
This is useful in applications where our particle moves at certain angular velocity, The radius of the path of our particle changes respect to time
Very common in machine design:
Slider-connecting rod –mechanisms
Robot arm movements
In polar coordinates, our velocity equation takes the form
Note:we often use dots to imply time (t) derivatives
• Likewise, our acceleration equation takes the form
NOTE! If radius r is expressed as a function of angle 𝜃 instead of time
(i.e. r = r(θ)), we must remember to apply chain rule in derivation!
𝑟' = 𝑟′(𝜃) ∙ 𝜃'
Force and acceleration
Acceleration is linked to the resultant force acting on the particle as well as the mass
of the particle via Newton’s 2nd law
Mass is also present in Newton’s 4th law – the law of gravitation. If we consider a particle with mass m that is falling near surface of Earth (mass 𝑚2),
we can write that the particle and Earth act upon each other via gravitation, so
If an object has mass, it also has inertia, which means that the object will resist changes in its state of movement. Accelerations and decelerations are often measured in g’s
All calculations must be performed in inertial coordinate system. i.e. a coordinate system which has an origin that does not rotate or have any acceleration. Unfortunately, a coordinate system like this has not (yet) been discovered
In “traditional” engineering tasks we can use an inertial coordinate system which
is fixed on the surface of Earth
Just like in kinematics, also kinetics problems can be handled in different
coordinate systems
Friction
For better accuracy, it is possible to separate static (𝜇𝑠) and dynamic (𝜇𝑘) COF. Becaurful with N as it doesn't always equal to the mass of the object.
x, y coordinates
In X,Y-coordinates, we can divide Newton’s 2nd law to x- and y-components and hence get two equations. Sometimes it is clever to set your coordinate system “tilted”
N, T coordinates
If our particle is moving along a known curved path, it is often wiser to write the
equations in N,T-coordinates
Cylindrical coordinates
Work and Energy
If a particle is subjected to force F and the particle undergoes a finite displacement along its path from s1 to s2 , the work done by force F is the product of force component in the direction of the path and displacement
we must remember, that our forces usually act on tangential and normal directions – which are not the same as θ and r
Hence, we will introduce a new angle ψ which means the angle between the
extension of the radial coordinate r and tangent of the path curve
If our force component in the direction of movement is constant, we can just multiply: U = F∆s
If the displacement is in the same direction as the force, work U is positive, and if it's opposite then U is negative.
The force of a spring is kx, and by performing the integration we get:
But becuase the force that is exerted on the particle is in the opposite direction, the final formula is negative
If we consider a particle moving along its path in an inertial coordinate system, the direction of movement is naturally tangential direction we can use the kinematic equation and intigrate to find the work done. Which results in:
This principle can be applied not only to a single particle but also a system of particles – with following limitations:
All particles must be a part of a rigid body which is either at rest or in translation (so, no rotation)
OR particles must be connected to each other with
ropes that have zero stretch
Note! Normal forces don’t do work, so when examining forces in normal direction, we have to use the equation of motion
Power is the time-rate of doing work – so, time derivative of work:. 𝑷 = 𝑭v
The concept of power is often needed when we need to select a suitable motor for
certain applications
s it is essential to differentiate between power input and power
output. Between these two we can form the relationship, This relationship is called mechanical efficiency of the machine
Conservation of energy
Forces are said to be conservative if the work done by them is independent of the
path followed by the particle
Conservative forces are for example:
- Weight
- sping
On the other hand, for example friction is a non-conservative force.
If our system only consists of forces which are conservative, we can write the
principle of work and energy into the following form:
T = Kinetic energy, V = potential energy (gravetational and spring)
click to edit
Impulse and Momentum
Linear Impulse and Momentum: If we consider the already familiar Newton’s 2nd law and remember that
acceleration is the time derivative of velocity, we get
F=ma=mv dv/dt
If we multiply this by dt and put in bounds for integrals, we get
This equation comes in handy when we have problems which are dealing with force, velocity and time. The equation can be applied separately in x- and y-directions, also be applied to a system of particles.
When thinking about dynamics processes, we can sometimes simplify our calculation by neglecting forces which have small impulses
Forces which produce a very small impulse are called
non-impulsive forces.
Forces which are large (in comparison to other forces) but act for a small time period are called impulsive forces
CENTRAL DIRECT IMPACT: The collision of two objects during a very short time is called an impact. The tangential plane of the contact surfaces is called plane of contact. The normal drawn to this plane at the contact point is called line of impact
If the mass centers of both objects lie on the line of impact, the impact in
question is a central impact
If the velocity vectors of contact points of both objects are parallel to the line of
impact, the impact in question is direct impact
In the beginning of the collision, objects exert deformation impulses integral P dt to
each other.
For a very short instant of maximum deformation, the objects will move with a
common velocity v
After this instant, the objects exert restitution impulses integral R dt to each other. As a result the velocity of the two objects changes
Because both deformation and restitution impulses are internal and cancel each
other out, the momentum of our system of two particles is conserved, so
If initial velocities are known, this leaves us with two unknowns. Therefore we need another equation in order to solve for final velocities. The second equation we can get by considering both particles separately at both
phases
The ratio of impulses is called coefficient of restitution:
By performing the same analysis for particle B, we get:
Now, the common velocity v (at the instant of maximum deformation) is only valid for a very short time and it’s not of our interest, so we eliminate that from equations for e
Value of e is always something between 0 and 1. In experimental tests it has been observed that e has correlation with initial
velocities, material, shape, mass of particles etc. If we want to calculate accurate results, this should be defined separately
Elastic impact: e = 1 – the particles “bounce” off each other
Plastic impact: e = 0 – the particles stick to each other (no restitution)
If e is not given, we can get reasonable results by simplifying our impact as
elastic or plastic
CENTRAL OBLIQUE IMPACT: If the velocities of colliding particles are not aligned with the line of impact, the
impact is called oblique
If we fix y-axis along contact plane and x-axis along line of impact, we can see that the impulses only affect in x-direction. The y-components of velocities remain unchanged.
Applying the coefficient of restitution and conservation of momentum for the x-components of velocities.
After the x-components of final velocities have been solved, final velocities and their direction angles can be solved using the Pythagorean theorem & inverse tangent
ANGULAR IMPULSE AND MOMENTUM: Angular momentum of a particle about point O is the “moment of momentum”
Here d is the perpendicular distance from O to the line
of action of mv
It can be mathematically derived that the resultant moment about point O is equal to the time derivative of angular momentum. This means that the change in angular momentum of a
particle is caused by a moment impulse!
This moment impulse is known as the angular impulse. The net torque equals the rate of change of the angular momentum! The net torque acting over a time interval is the angular impulse.
Vibrations: periodic motion of a body displaced from equilibrium
UNDAMPED FREE VIBRATION: When the motion is maintained by gravitational or elastic restoring forces (usually
spring force), vibration is called free vibration.
If there are no frictional or other similar forces that would resist the motion, the vibrations can continue indefinitely, the vibration is undamped
one degree of freedom (DOF): 1 DOF = one particle with mass m that can only move in 1 direction
In undamped free vibration, there is only one force acting on
the direction of movement: the restoring force F
−𝑘𝑥 = 𝑚a
If we move ma to the left and divide by mass, we get:
𝑎 + 𝑘/m = 0
Acceleration is the 2nd time derivative of displacement, so:
x'' + 𝑘/m = 0
This is usually written in form
This equation is a homogenous differential equation (DE) of 2nd order with constant coefficients. Later in math courses you will (hopefully) learn why, but at this point, let’s
just believe that the solution for it is
The solution can also be written in the form:
Here p is called circular frequency or natural frequency (often denoted by 𝜔 instead of p, but you’ll see a bit later why I favor p).
C is the amplitude (= maximum displacement) of our particle
The time period required to complete one cycle is:
We can also calculate the frequency of the vibration (number of cycles per time period):
VISCOUS DAMPED FREE VIBRATION: If there exists a viscous damping force in the system, the vibration becomes damped. our vibration will no longer continue forever; its amplitude will
decrease gradually towards zero
The damping force is proportional to the velocity of our particle, so our Newton’s 2 nd law –equation takes the form:
−𝑘𝑥 − 𝑐𝑣 = 𝑚a
Here c is the coefficient of viscous damping (unit kg/s).
If we move everything to the left again and remember that velocity is the 1st time derivative of displacement, we get the differential equation.
Again, we end up with a homogenous differential equation of 2nd order with constant coefficients. Now when there exists a damping term, the behavior of the system can be divided to three different cases – depending on the coefficient of damping
The critical coefficient of damping cc can be calculated from the equation:
When c is greater that cc, the system is called overdamped.
When c is smaller that cc, the system is called underdamped
When c is equal to cc, the system is called critically damped
For an overdamped system, the solution is:
For a critically damped system, the values for 𝜆1 and 𝜆2 are equal and therefore the solution simplifies to form:
Generally speaking, these 2 cases don’t actually represent vibration; in these cases the damping is so “strong” that the displacement x will reduce to zero without any oscillation.
Critically damped system will do this faster than an overdamped system
For an underdamped system, both 𝜆1 and 𝜆2 will be complex (the part under the square root will be negative)
Now the general solution takes the form:
In here, constant pd is called the damped natural frequency of the system:
The ratio c/cc is often referred to as damping factor. The time period of damped vibration can be written as:
UNDAMPED FORCED VIBRATION: When the motion is caused by an external periodic or intermittent force, the vibration is said to be forced vibration. External force has an angular frequency (now, this is denoted by 𝜔!)
The equation of motion will now be:
𝐹0 sin 𝜔𝑡 − 𝑘𝑥 = 𝑚a
Rearranging gives us:
𝑚𝑥''+ 𝑘𝑥 = 𝐹0 sin 𝜔t
This is a nonhomogenous differential equation of 2nd order with constant coefficients. The solution is a linear combination of these two:
+
If we take a closer look at the formula of particular solution:
we notice that if values of angular frequency ω and natural frequency p are really close to each other, the denominator will become very small, and therefore the amplitude will become very large
This is called resonance, and it happens when 𝜔 ≈ p. Must be avoided “at all costs”, because this leads to uncontrollable situations!
VISCOUS DAMPED FORCED VIBRATION: the forced
vibration case with damping
Now our differential equation takes the form:
Again, the solution for this differential equation is a linear combination of homogenous solution and particular solution:
The particular solution takes the form:
Amplitude C’ and phase difference 𝜙′ are given by following equations:
The free vibration part xc will gradually fade away due to damping, so ultimately the particular solution is the part that will remain (and hence, most interesting). This xp is therefore called “steady-state solution”
EQUILIBRIUM OF A RIGID BODY
RIGID BODY OR PARTICLE: Particles are small and only have mass; their size is considered negligible, so
they have no “dimensions”. Rigid body, on the other hand, has dimensions or in other words Acting points of forces become significant (not always on mass center)!
The equations for equilibrium are 2/3 the same as for particles:
FREE-BODY DIAGRAMS: In order to apply these equations, we must first find out the forces and moments
acting on our rigid body AND their locations
SUPPORT TYPES AND SUPPORT REACTIONS: Usually we know the external forces (and moments) acting on our rigid body. In most cases, equilibrium equations are used in order to find out the support reactions of our structure
Structures can be attached to their surroundings via several different support types (pin joints, fixed joints etc.) It is extremely important to know which support reactions will these support types exert on our structure
Different types of support reactions:
Usually our structures have as many unknown support reactions as we are able to solve using the equations of equilibrium (so, maximum of 3 in 2D cases). In reality our structures may have more supports and hence more support
reactions
These cases are referred as statically indeterminate.
If our structure has n support reactions and we have m equations, the difference (n – m) is referred to as degree of redundancy
Generally speaking, higher degree of redundancy gives our structure more reliability. On the other hand, calculation of support reactions becomes harder
This can be used to analyze knots, pulleys, connecting rings, and any small detail where all forces acting on it can be described as
acting on its center
STRUCTURAL ANALYSIS OF TRUSSES
A truss = structure composed of slender members joined
together at their end points. In calculation, joints of the truss are considered as pin joints
If the truss is welded or the joints are strengthened by gusset
plates, this does not reflect reality 100%.
Still, if the truss members are long compared to their other dimensions (mathematically: 𝐿 ≫ (ℎ, 𝑤, 𝑑)), this simplification provides a viable ground for analysis
If all loadings are applied at the joints, the truss can be
analyzed as a simple truss. As a result of the pin joint simplification, the truss members will only transmit normal forces – i.e. forces in axial direction
These rods can either be in tension or compression
Simple trusses consist of triangles. Most practical uses feature planar trusses, which are joined together by purlins or
floor beams in order to form a 3D-structure
Elementary statics gives us tools for evaluating the following things from zeroredundancy trusses: Support reaction and member forces
Essential to find dimesion of trusses members.
Mechanical II
If the truss has a degree of redundancy greater than zero, we will need additional
information about sizes & possible materials of members
Mechanical II and III
Tools for deformation analysis of trusses will be learned in further courses
Mechanical II and III
METHOD OF JOINTS
The truss member forces can be solved using the method of joints. In this method, we analyze each joint of the truss separately as a particle that
must be in equilibrium, so we apply
If there are n joints in the truss, we can form 2n equations. Draw a FBD of the whole truss and solve the support reactions – or at least write the equations of equilibrium for it (the truss itself is a rigid body!)
Solve the unknown forces (if possible) and move on to the next joint. Always draw unknown forces as tension forces – i.e. away from the joint! When done this way, you can see directly from the results whether the member is in tension (positive result) or compression (negative result)
In the end, present your solutions as positive values with a T or C in brackets (Tension or Compression)! For example: FAB = 1.47 kN (C)
ZERO-FORCE MEMBERS
It is rather common that our truss has members which don’t carry primary load. Therefore, the force in these members is considered to be zero.
These members are added to the structure in order to: Increase stability, Increase redundancy, Provide support in case of secondary loading or if the loading changes.
Two general rules can be used in order to spot zero-force members:
If only two members form a truss joint and no external load or support reaction is applied to the joint, the members must be zero-force members
If three members form a truss joint in which 2 of the members are collinear, the 3rd member is a zero-force member, if no external load or support reactionis applied to the joint
THE METHOD OF SECTIONS
This is based on the principle that if the whole body is in equilibrium, also any part of it must be in equilibrium.
Hence, we can “cut” any part of the truss and inspect it as a rigid body that must be in equilibrium
INTERNAL FORCES
When we move on to investigating beams, limiting our analysis only to normal forces will not be sufficient, because they are capable of carrying also shear force and moment. In general case, these three are called internal forces
These internal forces are not constant but take different values along the length axis of the beam. In design, it is essential to be able to find out the maximum values of internal
forces as well as their distribution
Internal forces can be solved at any point of our structure using the method of sections – just like we did with the trusses
Now when we are dealing with shear and moment, it is important to remember the sign convention rule for these forces:
- Shear force: if we use the left part of the cut, positive V is down
- Moment: left part of the cut = counter-clockwise
- For moment, also the “smiley face rule” can be used
METHOD OF SECTIONS
Usually, we need to solve the support reactions first.
If the question is just “calculate the internal forces at point X”, we can perform the cut at this point, draw FBD and simply solve our three unknowns using the three equations of equilibrium
BUT, if the question is to solve internal forces in the whole structure, we need to perform a much longer analysis in order to get the shear and moment equations (so, V and M as functions of length coordinate) and diagrams
In the latter case, we need to specify a length coordinate x starting from one end
of the beam and perform cuts after every discontinuity in the beam
Discontinuities are caused by:
➢ Point loads and point moments
➢ Start/end points of distributed loading
➢ Supports
SHEAR AND MOMENT EQUATIONS AND DIAGRAMS
As a result of the previous process, we get the equations (or actually functions, but “equations” is a more commonly used term) for shear and moment
From these, we can recognize the most critical parts of our structure and
calculate the maximum (as well as minimum) values for V and M
eg:
Using a FBD of the whole beam, we can solve the support reactions:
𝐹𝑦 = 0 ⇒ 𝐴𝑦 + 𝐶𝑦 − 15 − 25 = 0
𝑀𝐴 = 0, −80 − 15 ∙ 5 − 25 ∙ 7.5 + 𝐶𝑦 ∙ 10 = 0
10𝐶𝑦 = 80 + 75 + 187.5 ⇒ 𝑪𝒚 = 𝟑𝟒. 𝟐𝟓 𝒌N
𝐴𝑦 + 34.25 − 15 − 25 = 0 ⇒ 𝑨𝒚 = 𝟓. 𝟕𝟓 𝒌N
Now, let’s cut the beam first between A and B. We set a length coordinate x1 starting from left end of the beam. By applying the equations of equilibrium, we get functions for V and M:
0 ≤ x1 ≤ 5 m
𝐹𝑦 = 0 ⇒ 5.75 − 𝑉 = 0 ⇒ 𝑽 = 𝟓. 𝟕𝟓 𝒌N
𝑀 = 0 ⇒ 𝑀 − 5.75𝑥1 − 80 = 0
𝑴 = 𝟓. 𝟕𝟓𝒙𝟏 + 𝟖𝟎 𝒌𝑵
Then, let’s section the beam between B and C. We set a new coordinate x2 starting from left end of the beam. In a similar fashion, we get:
5 m ≤ x2 ≤ 10 m
𝐹𝑦 = 0 ⇒ 5.75 − 15 − 5 𝑥2 − 5 − 𝑉 = 0 ⇒ 𝑽 = (𝟏𝟓. 𝟕𝟓 − 𝟓𝒙𝟐) 𝒌N
𝑀 = 0 ⇒ 𝑀 − 5.75𝑥2 + 15 ∙ (𝑥2 − 5) + 5(𝑥2 − 5)∙(x2-5)/2 - 80 = 0
⇒ 𝑴 = −𝟐. 𝟓𝒙𝟐^𝟐 + 𝟏𝟓. 𝟕𝟓𝒙𝟐 + 𝟗𝟐. 𝟓 𝒌𝑵
Now we can draw the shear and moment diagrams and also calculate maximum values by substitution:
RELATIONS BETWEEN DISTR. LOAD, SHEAR AND MOMENT
GRAPHICAL METHOD
Usually, we are not that interested about the equations but more interested in the general shape of the diagrams and maximum values. The before mentioned relations can be used in order to sketch the shear and moment diagrams in a much easier fashion. Maximum values can be calculated using integration
click to edit
v = dS/dt
a = dv/dt = d2S/dt^2
a ds = v dv
Kinematics of a Rigid body
Motion of rigid bodies can be divided (in 2D) to three cases:
- Translation (rectilinear or curvilinear)
- Rotation about a fixed axis
- General plane motion (GPM)
TRANSLATION: We may use either a fixed coordinate system or a translating coordinate system in our calculations. In vector form, the position of point B can be given as
The time derivatives of this give us velocity and acceleration. If the rigid body is only in translation, vector 𝒓ത𝑩/𝑨 is constant (both in magnitude as well as in direction). • Therefore, taking the derivatives will yield us
ROTAT: ION ABOUT A FIXED AXIS: When a body is rotating about a fixed axis, any point P of the body travels a circular path subject to the fixed axis.
If we draw a perpendicular line from the fixed axis to point P, we can specify the angular motion of this line:
- Angular position is defined by angle θ
- Angular displacement is defined by 𝑑θ
- Angular velocity is the time derivative of angle
- Angular acceleration is the 2nd time derivative of angle
The angular motion of the line can be linked to point P. This is done by multiplication with radius r
• If we eliminate dt from previous differential equations, we get a 3rd equation:
If the angular acceleration is constant (𝛼 = 𝛼𝑐) , integration of these 3 gives us
When converted to motion of point P, we must also remember that the total acceleration is the Pythagorean sum of tangential and normal components:
GENERAL PLANE MOTION (VELOCITY): General plane motion of a rigid body is a combination of translation and rotation. In order to assess these separately, we must establish a relative-motion analysis
If we think back to the position vector again. If we took the derivative of this, in pure translation was constant hence disappeared in derivation. Now when our body is also rotating, that means that the direction of 𝒓ത𝑩/𝑨 is
changing, so its derivative will not be zero
Therefore, now our velocity will take the form
V B/A is the relative circular velocity of B subject to A
So, how to find out the relative velocity 𝒗𝑩/𝑨?. Think that momentarily the body is “pinned” at A. Since the motion is circular, the magnitude can be calculated by
The direction of 𝒗𝑩/𝑨 is perpendicular to 𝒓𝑩/�
In planar case, this vector equation breaks down to x- and y-components. As a result, we get magnitudes for velocity components:
The minus sign is caused by sign rule counterclockwise). Remember: here, also r has direction!
GENERAL PLANE MOTION (ACCELERATION): Since in general plane motion there is both translation and rotation, we must remember this when calculating acceleration
In rotation, the total acceleration is caused by tangential and normal components. Therefore, in general plane motion the acceleration will be
Here (𝒂𝑩/𝑨) 𝒕 is the relative tangential acceleration. Direction is perpendicular to 𝒓𝑩/�
Likewise, (𝒂𝑩/𝑨) 𝒏 is the relative normal acceleration. Direction is from B to A.
This is also a vector equation, which
yields us two scalar equations (x and y)!
Force and acceleration
MOMENT OF INERTIA: Mass m is the measure for a body’s resistance to acceleration. Likewise, there exists a measure for body’s resistance to angular acceleration; this measure is called moment of inertia or second moment of mass and its symbol is I
Mathematically speaking, this moment is calculated as an integral
Considering that 𝑑𝑚 = 𝜌 𝑑𝑉, this can also be presented in form
It is common that density 𝜌 is constant, so it can be factored out of the integral
We could naturally calculate moments of inertia for different-shaped bodies using integration, but in order to simplify calculations, there are formulas for the most commonly used ones
KINETIC EQUATIONS IN TRANSLATION: The body has acceleration in x- and y-directions, but no rotation. Therefore the equations of motion become. Note: aG
refers to the acceleration of mass center
KINETIC EQUATIONS IN TRANSLATION: No rotation, but the acceleration consists of tangential and normal components. Therefore the equations of motion become
Note: Here 𝜌 is radius of curvature, not density!
KINETIC EQUATIONS IN ROTATION ABOT A FIXED AXIS: If the rotation is accelerating or decelerating at angular acceleration α, the moment sum equals to this times moment of inertia. Hence, the equations of motion become
: Here 𝑟𝐺 refers to distance from O to G
KINETIC EQUATIONS IN GENERAL PLANE MOTION: In case of general plane motion, our equations will be
If the body is in GPM, it means that it rotates and translates at the same time. Hence, we often have to deal with unknown friction and support forces. Because we only have 3 equations with usually 3 unknown accelerations, this means that we need additional equations in order to solve the problem
In this case following tricks are used:
➢ Choose your coordinate system (x,y) in such a way that the other coordinate axis is in the direction of movement and the other is perpendicular (this will make the perpendicular acceleration component zero)
➢ If the movement is rolling with no slipping, use fourth equation that relates straight-line acceleration and angular acceleration:
➢ If the movement has slipping, use fourth equation that relates friction and support forces:
We can solve this kind of problem even if we don’t know whether the body will slip or not. First assume there is no slipping and use 𝑎𝐺 = 𝛼𝑟 as the 4th equation in order to
solve the group of equations
After this, perform a check:
If 𝐹𝜇 < 𝜇𝑘𝑁, the assumption stands and the solution is valid
If 𝐹𝜇 > 𝜇𝑘𝑁, the assumption was wrong (slipping occurs); change the 4th equation and rework the solution
click to edit
APPLYING THE ROTATIONAL EQUATION ABOUT POINT P: It is not compulsory to write the rotational equation subject to mass center G; it is also possible to write this equation about any point P in our body
This is done by fixing the origin of our coordinate system to P and using a kinetic diagram and calculating the “kinetic moments” of acceleration components
In rotation about fixed axis, the normal component of acceleration is directed towards the point of rotation O; so, if we take this moment sum subject to point O, the equation takes the form
Applying the rotational equation about some other point than G is not necessarily needed, but in some cases, it leads in much easier calculations
This is due to some forces getting cancelled out from the moment equation
Most usual case: support reactions of rotation point are unknown → if we apply the rotational equation subject to rotation point, we don’t need to solve the support reactions
Work and Energy
KINETIC ENERGY OF A RIGID BODY::
In pure translation (rectilinear or curvilinear), kinetic energy can be calculated the same way as for particles.
Subscript “G” is added to velocity to accentuate that this is the velocity of the mass center
If our rigid body is in pure rotation around its mass center, the velocity of its mass center is zero. However, the body still has kinetic energy due to rotation. The kinetic energy of this rotation can be calculated analogously
If our rigid body is in general plane motion, then its kinetic energy consists of translational and rotational energy. Due to the nature of energy, we can just calculate the total kinetic energy as the
sum of these two:
Note: if the rigid body is in pure rotation about a fixed axis passing through point O, then the kinetic energy can be expressed as
Because calculation of IO can be troublesome, in general I recommend using the GPM form
Exception #1: if we know the radius of gyration 𝑘𝑂, because this makes the
calculation of IO very simple
Exception #2: Bar attached at its end point O (because also this makes the
calculation of IO very simple)
WORK OF FORCE AND COUPLE MOMENT:
Work done by a force can be calculated as an integral of F over s
Work done by a couple moment can be calculated as an integral of M over theta
If positive work is done to the system, the total energy of the system increases!
Work of a force is positive if F and displacement s are in the same direction
➢ Work of a couple moment is positive if M and angle of rotation 𝜃 are in the same direction
PRINCIPLE OF WORK AND ENERGY: Same principle of work and energy that we already used for particles is also valid
for rigid bodies. We must only remember that for a rigid body the kinetic energy is a sum of translational and rotational energy
Good method: draw the rigid body in the beginning and in the end in order to see
the changes throughout the movement process
CONSERVATION OF ENERGY: Revision: if the work done by the force is independent of the path, the force issaid to be conservative.
If only conservative forces do work on our body, we can write the previous equation in a simpler form using kinetic (T) and potential (V) energies:
3D Force Vectors
EQUILIBRIUM OF A PARTICLE IN 3D
It is common that our structures are three-dimensional. If we can work out 2D cases, then adding one more dimension is actually not that difficult; procedure of particle equilibrium stays pretty much the same
➢ Draw FBD of the particle
➢ Mark all the forces acting on the particle
➢ Divide forces to components (now there are three: x, y and z)
➢ Set sum of components equal to zero in all directions
We end up with three equations, where we can solve three unknowns. We need at least 2 angles in order to specify the direction
Usually we give three:
➢ α = angle between force vector and x-axis
➢ β = angle between force vector and y-axis
➢ γ = angle between force vector and z-axis
Calculation of these angles can be troublesome in scalar form. In three-dimensional cases, it is often wisest to make the switch from scalar equations to vectors. Denote each force by a vector, whose i-, j- and k-components represent x-, y- and
z-components of the force (correspondingly)
If we know the direction vector v of the force, we can construct the force vector easily by first calculating the unit vector uv and then write
The Bolded letter are vector and the unbolded ones are scalar values
The sum of all forces is then set to equal a zero vector. This vector equation can then be transformed to 3 scalar equations. The division of forces to components is quicker to do via the unit vector than via trigonometry using different angles
If coordinates are A = 𝑥1, 𝑦1, 𝑧1 and B = 𝑥2, 𝑦2, 𝑧2 , the direction vector is
:
When we multiply the unit vector by the magnitude of force, we get our hands directly on the components:
If we want to know the direction angles, we can calculate them right from the
components using inverse cosine
Moment of a force in 3D
Earlier during Engineering Mechanics I we specified moment of a force about a point O as “force times shortest (perpendicular) distance from point to force vector” – or simply as force times moment arm
This was completely okay in 2D applications, but now when we move on to 3- dimensional structures, a vector definition makes a lot more sense. Moment of a force F about point O is the cross product of vectors r and F
• The magnitude of the cross product equals the area of the parallelogram, which has original vectors as sides.
The direction of the cross product we can get using the right hand rule
• A much quicker way is to can calculate the cross product using matrix notation as a determinant
Vector r is a position vector that goes from point O to any point which lies on the line of action of F
F is the forece of the vector
Notice: the i-component of the moment is not affected by Fx , the j-component of the moment is not affected by Fy etc. This is natural, because they produce no moment in the said direction; they are parallel.
Direction angles of the moment vector can be calculated using inverse cosines
IMPORTANT! The cross product is NOT
commutative!. 𝒂 × 𝒃 ≠ 𝒃 × a
On the other hand, it is anti-commutative
𝒂 × 𝒃 = −𝒃 × a
EG: Here don't forget how to get the componenet of the force.
RESULTANT MOMENT OF A SYSTEM OF FORCES: If a body is acted upon by several forces, we can calculate the moments for each force separately. The resultant moment is then the sum of all these moments
If two or more forces act on the same point, we can apply the distributive law of cross product in order to ease the calculation So, instead of calculating two cross products, we can first sum the forces and then calculate moment for the resulting force vector:
𝑴𝑶 = 𝒓 × 𝑭𝟏 + 𝒓 × 𝑭𝟐 = 𝒓 × 𝑭𝟏 + 𝑭𝟐 = 𝒓 × F
To find the angles of the moment vector, you have to first convert the Moment vector into a unit vector by dividing it by its magnitude and then using the cosin rule.
MOMENT OF A COUPLE USING CROSS PRODUCT: Moment of a couple can be calculated using the cross product as a resultant of
𝑴𝑶 = 𝒓𝑨 × (−𝑭) + 𝒓𝑩 × F
Anyhow, if we think that the minus sign belongs to rA
, this can be written as
𝑴𝑶 = (𝒓𝑩 − 𝒓𝑨) × 𝑭
Using the triangle rule of vector addition, we get
𝑴𝑶 = 𝒓 × F
Example:
Replace the two cples acting on the pipe coloumn by a resultatn couple moment.
MOMENT OF A FORCE ABOUT COORDINATE AXIS: If we calculate the moment of a force about a specific point, we will get a moment vector whose direction is perpendicular to the plane containing the force and the moment arm. This moment is not always the one that we want to use in calculations; we may
have to break it into parts
Different components of the moment have different effects on our structure:
➢ x-component of the moment is trying to turn our structure around x-axis
➢ y-component of the moment is trying to turn our structure around y-axis
➢ z-component of the moment is trying to turn our structure around z-axis
The structure usually has different “ability” to fight against the moment
depending on its direction
If we calculate the moment subject to point O, we will get a moment vector which acts on direction Ob; the magnitude of this moment has no clear physical meaning
The moment about y-axis does; it tries to turn the pipe around y-axis and hence creates torsion.
The moment about x-axis tries to turn the pipe around x-axis and hence creates bending
So, Mx is -8 Nm and My is 6 Nm
Hence, the moment that tries to twist the pipe is 6 Nm and the moment that tries to bend the pipe is 8 Nm
MOMENT OF A FORCE ABOUT AN ARBITRARY AXIS: What if the component that we are interested in is not in the direction of a coordinate axis – i.e. what if the pipe in the previous example would stick out of point O to some arbitrary direction, say 𝒂 = 𝑎𝑥𝒊 + 𝑎𝑦𝒋 + 𝑎𝑧𝒌 (instead of 𝒋)? In this case we need to calculate the projection of the moment vector 𝑴𝑶 on to
vector 𝒂 – let’s call it 𝑴a or actually we just need the magnitude Ma
Luckily, this can be calculated easily using the dot product:
Here 𝒖𝒂 is the unit vector in the direction of 𝒂
TRIPLE SCALAR PRODUCT: Using rules of mathematics, we can simplify our calculation by merging the dot
product into our cross product that we need in calculating the 𝑴o
Hence, we only need to calculate one 3x3-determinant and we get the magnitude of our moment subject to axis in the direction of 𝒂
The sign is important: if 𝑀𝑎 is positive, then the moment vector 𝑴𝒂 is in the same direction as 𝒂.
If 𝑀𝑎 is negative, it is in the opposing direction
If we want to know the vector form, we can get it by multiplying the unit vector
𝒖𝒂 by Ma
EQUILIBRIUM OF A RIGID BODY IN 3D
In the 3D-world our number of suppoert reactions increase by the following ones:
Support force in the z direction
support momment Mx
Suppoert momnet My
This means that we have a maximum of six support reactions or in other words six unknows. Normally not all six are present, In order to help in deciding whiich support reactions some osupport will develop think about the suport in the following way
If the support reaction stop translation in some direction, then it willl develop a support forece in said direction
If support restricts rotation subject to some axes then it will develop a support moment about said axis
Note: Many time we acn ingnor the couple moments of single supports, if there are other support s which are porperly aligned
Single connections theoretically are able to produce support couple moment, they are not very "good at it"
Moment carrying capacity is limited
Stress coaused by the mommnet often harms operations and shortens lifespan (especially bearings)
This is why single connection are often "shielded" from producin these moments by use of multiple properly alinged supports.
Multiple supports togethere develope adequate support reactions in order to maintain equilibrium
Journal bearing A, B and C would theoretically produce couple moments. Since they are used together properly aligned, together they will need to develop only support forces in two directions each
The red spots is where we would add journal bearings if we wanted to "kill" those moments
In 3D application, the formula for equilibirum is the same as for 2D
Cartesian vector from:
Scalar form:
CONSTRAINTS FOR A RIGID BODY: Sometimes there may occur problems, where the structure has more supports
than necessary in order to maintain equilibrium. These “additional” supports are called redundant supports and the number of them is often referred to as degree of redundancy
This kind of structures are called statically indeterminate. Very common in real-life structures. Additional supports provide additional safety margin (if designed right)
Calculation is problematic: this means that we have more unknowns than equations!
Solving this kind of problems is possible, but it requires to inlvolve the physical porperties of the body.
So, the solution also depends on material parameters and “inner”
measurements of the rigid body
Example: degree of redundancy = 2
It is also possible that we have as many unknowns as we have equations, but we can’t solve the support reactions. This is the case when our body is improperly constrained
Improper constraints occur when all support reactions in one direction intersect a common axis.
In calculation this shows as a case where we can write a moment equation where none of the unknown forces are present, hence the body will not be in equilibrium
This is a common error f.ex. in finite element analysis, where stresses and deformations are calculated numerically using a computer program: if your supports aka. constraints for the body are not sufficient, the calculation will fail and produce an error message
Improperly constrained case: rotation is not prevented about x-axis even though the amount of support reactions looks good
SOLUTION SUMMARY
click to edit
1) Draw a FBD of the structure (supports to support reactions!)
2) Find out the forces in vector form
3) Write the moment equilibrium subject to some point in the structure (tip:
choose a point that has as many unknown support forces as possible)
4) Find out suitable moment arm vectors
5) Calculate cross products and substitute results to the moment equilibrium
6) Convert the moment equilibrium to 3 scalar equations
7) Solve, if possible
8) Write the force equilibrium and convert to 3 scalar equations
9) Solve the remaining unknowns
Space Truss
A three dimensional truss structure is referred as a space truss. To calculate requires external loading applied at the joing. Joints are ball and socket type and weight of memebers is small and can be neglected.
Just in the 2D cases, these simplification don't align 100 % wiht reality, but they provide a good enough initial version for design
We can use the same two methods that we used in 2D. Method of joint and Method of sections.
If all of the truss must be determined for the method of joints. Start from a joint where there is one force you know and only three members. Sometime you might need to find the support forces before this. Then work your way up to the next joint that has 1-3 unknowns.
If you onlt need some of the foreces, it's often possible to ge tinto these quicker by using the method of sections.
Don't forget about symmetry and how they affect the direction of the vectors as well as how the change the final equation of equilibrium
Use this as a road map for the method of joints:
1) Chose a joint(preferably with max 3 unknown member forces, if possible) and draw FBD of the join(assume unknown forcess in tension)
2) Find out direction of vectors for the members
3) Write all forces in vector form
4) Sum of all force acting on the joint must be zero, so write the vector equation for force equilibirum
5) Group your vector components and convert the vector equation to three scalar equations
6) Solve the scalar equations and mark the solved forces somewhere
7) Go back to 1) and repeat until all forces have been solved
Another of doing this using momentum instead of equilibirum:
1) Find a suitable section that goes through 6 unknown member forces at maximum
2) Draw a FBD of one half of the section (you may choose which one)
3) Find out direction vectors for the unknown members & write forces in vector form
4) Sum of moments subject to any point must be zero, so choose a suitable point for moment equilibrium and write it in vector form. A joint that has 3 unknonw forces would be best
5) Find out required moment arm vectors and calculate cross products
6) Solve the moment equilibrium equation (breakdown to 3 scalar eqs etc.)
7) Sum of all forces acting on the section must be zero, so write the force equilibrium
8) Solve the remaining unknown forces (again breakdown to 3 scalar eqs)
Example:
An experienced truss designer knows that when the truss on the right is loaded as shown, the biggest (absolute) member force is developed in some of the diagonal members. Solve the diagonal member forces using method of sections.
INTERNAL LOADINGA IN 3D
Internal loadings: The member forces in trusses( SPACE TrUSS) are internal nomal forces. Trusses whihc consist of pinned joints have no other internal forces than the normal one, but in general case, there are more.
In 2D internal loadings are normal force, shear force and moment. All of these internal loadings can be calculated at the point of interest by the following proccedure:
1) Cute the structure at this point and draw a FBD
2) In FBD, mark all the external forces and at the cut point the internal loadings N, V, M
3) The "Cutaway" must be in equilibrium, so formulate equations of equilibrium
4) Solve unknown internal loadings
The only thing you have to take into account is the third dimesion, which means tha there can be up to 3 internal forces and 3 internal moments
Fx is the internal normal for in this example
Fx = N
Fy and Fz both act as shearing forces, so the total internal shear force is
My and Mz try to bend the structure, so together they produce the internal bending moment
Mx tries to twist the structure, so it is the internal torsional moment
Calculating the internal loadings are improtant as it also helps with calculating stresses in our structure
Also these formular are derived for the case where x-axis is the longitudinal axis of our structure. If the longitudinal axis is y or z, change the subscripts
Tips for calculation. Make the cure and draw FBD. Note: if our cut part infludes support reactions, it is good idea to solve support reactions (even for the whole structure if needed) first
Beding is rotation around the lateral axis
Torsion is rotation around longitudinal axis
At the cut, mark internal loading components by a force vector and a moment vector
Use equation of equilibrium inorder to solve for all six unknowns(start from vector equations)
Calculate the internal loadings N, V, M and T according to provide formulae(remember to check the longitudinal axis)
And check your units
Characteristics of Dry Friction
Friction = force of resistance acting on a body which prevents or retards slipping. Always acts tangent to the surface at points of contact, in the opposing direction of the existing or possible motion
There are two types of friction: fluid friction and dry friction
➢ In fluid friction, contacting surfaces are separated by a film of fluid (which can be gas or liquid)
➢ In dry friction, contacting surfaces are in direct contact
In this course, we will only analyze problems involving dry friction
CHARACTERISTICS OF DRY FRICTION: Dry friction is caused by irregularities on the surfaces of our contacting bodies.
The actual distribution of normal forces is as shown in picture,
because the resultant forces in these contacting points get smaller as we proceed
Therefore, when we replace the distributed normal force by its resultant force N, it will not act in the same point as the weight W, but at distance x towards the direction of motion
Same applies to resultant friction force F too, but this doesn’t make a differencebecause F acts along the plane
If the mass center is at height h measured from the ground, we can see why the block has two possible ways of motion when pulling force P is increased:
➢ If h is small or the surfaces are slippery, the frictional force resultant Fs may not be great enough to balance P (i.e. 𝑃 > 𝐹𝑠 ), so the block will slip
➢ If h is large or the surfaces are very rough, the moment caused by the frictional force (or initially by force P, since P causes Fs ) becomes so large that the resultant normal force can not create enough countermoment to balance the situation (i.e. 𝐹𝑠ℎ > 𝑁𝑥), so the block will tip over
If the block is slip-constrained and on the verge of sliding, the total resultant Rs (created by support resultant N and friction resultant Fs ) is in angle 𝜑s
This angle is called the angle of static friction, and it is defined by equation
The friction force is always no greater than it is needed in order to prevent motion, so if we plot friction force F as a function of external force P, we get the following graph
The maximum value Fs is called limiting static-frictional force, and it is defined via coefficient of static friction 𝜇s
When the block starts moving, the friction force will drop to a constant valueFk which is defined via coefficient of kinetic friction 𝜇k
There are three possible types of cases when calculating problems involving dry friction:
Equilibrium problems: total number of unknowns is equal to the number of available equilibrium equations
Friction forces can be determined via equilibrium equations
However, we must check that they satisfy the inequality 𝐹𝑠 ≤ 𝜇𝑠𝑁 (if they don’t, the solution is not valid and motion will occur)
Impending motion at all points: total number of unknowns is equal to the sum of available equilibrium equations and available frictional equations
We can solve friction forces as well as other unknowns
Solution will be “automatically” valid
Impending motion at some points: total number of unknowns is less than the sum of available equilibrium equations and available frictional equations (or conditional equations for tipping)
Several possibilities of motion; calculate which case will occur first
Wedge: A wedge is a simple “machine” which is used to transform applied forces into much larger forces. Give small displacements or adjustments to heavy loads
A wedge is subjected to five forces:
➢ Support forces from up and down (N2 and N1)
➢ Frictional forces along top plane and bottom plane (F2 and F1)
➢ Appplied Force P
If friction forces hold the block in place after insertion in such a way that an applied force is not needed (i.e. P = 0), the wedge is referred to as self-locking
FRICTIONAL FORCES ON SCCREWS: Screws are usually used as fasteners, but they can also be used to transmit power. Let’s investigate a square-threaded screw, which is the most common choice for power transmission (jacks, lathes, vises etc.)
If we cut one thread off, we get
If we are trying to lift our weight W upward, the moment needed for this is
click to edit
If surface of the screw is very slippery in such a way that 𝜃 > 𝜑, the weight W istrying to move down by itself. The moment needed to hold the weight in place is
If the surface of the screw is rough in such a way that 𝜃 < 𝜑, the weight stays in place. The moment needed to lower the weight down is
In this case, the scre is said to be self-locking
FRICTIONAL FORCES ON FLAT BELTS: Power transmission is often done by belts. Many attachments are also done by ropes relying on friction (ropes can be
approximated as flat belts)
Direction of impending motion is important, since the force on this side is
greater!
If we cut off a differential element of the belt and draw a FBD, we get:
Because 𝑑𝜃 is infinitesimal, we can approx. sin 𝑑𝜃/2 = 𝑑𝜃/2 and cos 𝑑𝜃/2 = 1. Term dT * 𝑑𝜃/2 is very small (product of two infinitesimal values), so it is = 0
If we substitute dN into the first equation, we get
Integration from T1 to T2 gives
Solving with respect to T2 gives
Left one is the “driving” force, right one is the smaller force (decreased by the friction)!
VITUAL WORK
Let’s recall the formulas for work done by a force and work done by a moment:
If we now introduce a virtual differenetially small displacement, the resulting work will also be virtual and differentially small:
The reason why this is called virtual is because the displacement doesn't actually exist.
Sign rule:
➢ If F (or M) and displacement are in the same direction, work is positive
➢ If F (or M) and displacement are in opposing directions, work is negative
PRINCIPLE OF VISTUAL WORK:
If a body is in equilibrium, the sum of the virtual work done by all forces and moments acting on the body is zero.
Idea: even if the displacement would happen, the forces and moments would maintain the body in equilibrium.
Because the displacement is small, it doesn’t change the forces/moments from the original case.
Example:A beam loaded at center by a point force. If a small virtual displacement y is applied at A, the only forces that do work are P and Ay (point B remains in place, so Bx and By do no work )
Total work done must be zero, so
Because d(theta) is not 0, it must be that
This example shows that the method works, this equation is still possible to get witht eh sum of y forces eqaul to zero. If we use the principle of virtual work to a particle or single rigid body, then all displacements will be linearly dependent on the virtual displacement.
This means that virtual desplacement gets factored out of the equation without leaving any traces behind.
The power of this principle lies in applying it to a sysytem of connected bodies. Equilibrium equations can't be used for whole mechanisms.
So analyzing equilibrium of a mechanism requires desmemebering its parts and writing equilibirum equations separatlely for each part.
It allows us to prevent this breakdown and porvides us a quicker way to solve the problem
There are two different way to apply this principle of virtual work
Virtual work for rigid bodies (one or several)
and virtual work for deformable bodies (one or severeal)
The second one will be studie in Mechanics III
VIRTUAL WORK FOR CONNECTED BODIES:
Note that her, we will stick with analysis of mechanism with only one degree of freedom. In this case all displacement can be related to just one displacement.
This means that all formula can be writen as a function of one coordinate s. This one coordinate is quite commonly angle (theta)
Steps to solve an equilibirum problem of a mechanism using virtual work:
- Draw FBD of the whole mechanism
- Sketch the defelected position of the system when it undergoes positive vitrual desplacement
- Express the coordinates of moving points measured from a fixed point asfucntions
- relate these coordinates to a single coordinate so that you can express all coordinates as a fucntion of on coordinate
- Differentiate thse coordinates in order to get the virtual displacements
- write the virtual work eaution
- Simplify and factor out the common displacement
- Solve the unknow force/moment/postion
example:
The spring is unstretched when 𝜃 = 30°. Determine the required pulling force P that is needed to maintain equilibrium in the scissors linkage when angle 𝜃 = 60°. Neglect the mass of the links.
Spirng force is a function of theta:
Coordinates of B and D as function of theta:
Differentiation of these Gives the virtual displacements:
Total virtual work is then:
Because total virtual work must be zero:
Simplify and factor out common displacement:
subtitute theta = 60 and solve for P:
POTENTIAL ENERGY CRITERION: Remember that hwen our systems had only conservatives foreces, then its total energy remained constant, so we could write
Now when we're working with systems that start from rest and will be in equilibrium in the end, their kinetic energies T1 and T2 will be zero. The amoutn of potential enrgy is alawys relative to the chosen datum, therefore we can write the total potential energy of the system as function V(s).
Eg:
If a differential small work is done to the system, it results in differentially small change to s and therfor ethe work is
The total differential work is
If we apply a small virtual displacement, it will result in virtual change in V:
If in this equation 3 we subtiture 2 and thne 1 we get
On the other hand, principle of virtual work requires that 𝛿𝑈 = 0, so
(𝛿𝑠 ≠ 0)
The result means that a frictionless system of rigid bodies is in equilibrium, when the first derivative of its potential function is zero.
So the steps are the following:
- Sketch the system in its initial state (and if it help, also in it finals state)
- Establish a datum through a fixed point
- Express gravitational potential energies as a function of distance from datum
- Express potential energies of springs as function of thier stretch
- Link distances to streches orvice versa in such a way that all potential energies are expressed using only one coordinate
- Formulate the total potential energy function
- Take the first derivative of V and set it equal to zero
- Solve in order to get the coordinate value when system is in equilibrium
NOTE: The equation V' = 0 might be hard to solve analytically. - Numerical solution is then preferred
Eg: The unifoorm lin shonw has a mass of 10 kg. If the spring is unstretched when 𝜃 = 0°, determine the angle 𝜃 when
the system is in equilibrium.
If we set the datum at center of B, the potential
energy function of the system will be:
Express s and y as functions of one coordinate (𝜃):
Substitute to V in order to express it as a function of 𝜃:
Differentiate 𝑉(𝜃):
Set the first derivative of V equal to zero and solve for 𝜃:
Solve for 𝜃:
Substitute the values:
STABILITY OF EQUILIBIRUM: The potential energy function derived in POTENTIAL ENERGY CRITERION section, can also be used to investigate the stability of the equilibrium.
Equilibrium has three types of stability:
- Neutral (if displaced, will move but remain in equilibrium)
- Unstable (if displaced, will no longer be in equilibrium)
- Stable equilibrium (if displaced, will return to equilibirum)
We can find out the equilibirum type simply by taking the 2nd derivative of our potential energy function (V''(s)). Then subtituting the value that results in equilirbium to it and checking the sign.
Remember if the value is positive, its concave up
and if the value is negative its concave donw.
Eg: If the spring AD has a stifness of 18kN/m and it is unstreached when when 𝜃 = 60°, determine the angle 𝜃𝑒𝑞 when the system is in equilibrium. The load has a mass of 1500 kg and its mass center is in G. Investigate also the stability of the equilibrium position. Ignore the mass of links AB and CD.
Datum at level AB, Potential energy function:
Express s and y as functions of one coordinate (theta)
Substitute to V in order to express it as a function of theta (note: H is a constant parameter, not a variable)
Differentiate V(theta)
For stability analysis, calculate the 2nd derivative:
click to edit