For the only rigid body in question, choose point \(A \) as the pole, therefore the 3 Lagrangian parameters (\(n = 3 \)) that can be used for the kinematic description of the system are:
Before proceeding with the solution of the system, it is possible to classify the system by calculating the rank of the kinematic matrix using the following MATLAB® instructions. Listing4.4.3.
Since the \(\textrm{min} (m, n) = 3 == \text{rank} \ mat {A} \) and \(n == m \) then the system is kinematically determined. Finally, the solution of the linear system provides the values assumed by the Lagrangian parameters \(u_{1A} \text{,}\)\(u_{2A} \) and \(\varphi_{A} \) which, in this case, are identically null.
Therefore on the system of equations to be analyzed, the only effect is constituted by the modification of the vector of the known terms and there is no effect on the classification of the system which always remains kinematically determined (the matrix \(A \) has not changed). Therefore the system of equations becomes
The solution of the linear system provides the new values assumed by the Lagrangian parameters \(u_{1A} \text{,}\)\(u_{2A} \) and \(\varphi_{A} \text{.}\)
The proposed scheme does not present any particular novelty with respect to the previous scheme and the following MATLAB® instructions allow to obtain the desired solution.
$ % generic 2D displacement for a rigid body
rigidDispl = ...
@(u0, phi0, X0, X)...
[u0(1)-phi0*(X(2)-X0(2));...
u0(2)+phi0*(X(1)-X0(1))];
% geometric description of the beam
syms l;
A = [0; 0];
B = [l; 0];
C = [l; l/2];
% displacement description by using the point A as pole
POLE = A;
syms phiA;
phi0 = phiA;
u0 = sym('uA', [2 1]);
uA = rigidDispl(u0, phi0, POLE, A);
uB = rigidDispl(u0, phi0, POLE, B);
uC = rigidDispl(u0, phi0, POLE, C);
% constraint equations
eqns = [
uA(2) == 0,
uB(2) == 0,
uC(1) == 0
];
% kinematic matrix and vector of the assigned displacements
[A,d] = equationsToMatrix(eqns, [uA(1), uA(2), phiA]);
% degrees of freedom, n
% number of constrains, m
[m,n] = size(A);
% rank of A
r = rank(A);
% if the system is kinematically determined the solution is calculated
if and(r == min(m,n), m == n)
x = linsolve(A,d);
end
The rigid bodies under examination are two: for the body \(AB \) the point \(A \) is chosen as the pole and for the second body the point \(B'\) is chosen. Therefore the 6 Lagrangian parameters (\(n = 6 \)) adopted for the kinematic description of the system are:
Since \(\textrm{min}(m,n) = 6 == \textrm{rango}\mat{A}\) and \(n == m\) then the system is kinematically determined. The solution of the linear system provides the values assumed by the Lagrangian parameters \(u_{1A}\text{,}\)\(u_{2A}\text{,}\)\(\varphi_{A}\text{,}\)\(u_{1B'}\text{,}\)\(u_{2B'}\text{,}\)\(\varphi_{B'}\text{.}\)
Also in this case the solution of a kinematically determined system with respect to a vector of the known terms which is identically null still provides the trivial solution.
The graphical representation of the displacement field is obtained by evaluating the displacement of the end points of the straight parts of the beams, obtaining what follows.
$ % geometric description
L = 100;
A = [0; 0];
B = [2*L; 0];
C = [2*L; -L];
D = [4*L; -L];
beam1 = [A B];
beam2 = [B C D];
% function for the evaluation of the new position of a point
rigidT = @(u0, X0, phi0, X) ...
[X(1)+u0(1)-phi0*(X(2)-X0(2)); ...
X(2)+u0(2)+phi0*(X(1)-X0(1))];
% new configuration of the first beam (pole in A)
u0 = [L/4; 0];
phi0 = 1/4;
X0 = A;
TA = rigidT(u0, X0, phi0, A);
TB = rigidT(u0, X0, phi0, B);
beam1T = [TA TB];
% new configuration of the first beam (pole in D)
u0 = [0; 0];
phi0 = -1/4;
X0 = D;
TC = rigidT(u0, X0, phi0, C);
TD = rigidT(u0, X0, phi0, D);
beam2T = [TB TC TD];
% drawing
clf
x = beam1(1,:);
y = beam1(2,:);
line(x,y,'LineWidth',2,'Color','black')
x = beam2(1,:);
y = beam2(2,:);
line(x,y,'LineWidth',2,'Color','black')
x = beam1T(1,:);
y = beam1T(2,:);
line(x,y,'LineWidth',2,'Color','red')
x = beam2T(1,:);
y = beam2T(2,:);
line(x,y,'LineWidth',2,'Color','red')
xlim([0 4*L])
ylim([-L 0.5*L])
pbaspect([2.667 1 1])
For the two rigid bodies under examination, points B and D are chosen as poles, thus the kinematics is described through the following 6 Lagrangian parameters (\(n=6\)):
Always using the plane motion model represented by the equation (4.1.5) the constraint conditions can be expressed with respect to the chosen Lagrangian coordinates obtaining the following result,
In the present case \(\textrm{rango}\mat{A} == \textrm{min}(m,n) = 5\) and \(m < n\) therefore the system is a mechanism and characterized by infinite solutions. To describe all possible solutions, choose the rotation \(\varphi_{B}\) as parameter to be eliminated from the group of unknowns. In this way the system can be rewritten as follows
