-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathMacchiaMarco_Assignment_1_solution.m
105 lines (72 loc) · 1.94 KB
/
MacchiaMarco_Assignment_1_solution.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
clear
clc
addpath("matlab_functions\");
ex1_tau = zeros(2,6);
ex2_tau_and_force = zeros(2,5);
%% Ex 1.1
angle1 = 90;
angle2 = -90;
F_p1 = [0 0 0];
F_p2 = [0 0 0];
F_p3 = [0 0 0];
M_p1 = [0 0 0];
M_p2 = [0 0 0];
M_p3 = [0 0 0];
g = [0 -9.81 0];
ex1_tau(:,1) = compute_tau1(angle1,angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 1.2
angle1 = 0;
angle2 = 90;
ex1_tau(:,2) = compute_tau1(angle1,angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 1.3.1
angle1 = 30;
angle2 = 60;
F_p1 = [-0.7 -0.5 0];
g = [0 0 0];
ex1_tau(:,3) = compute_tau1(angle1,angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 1.3.2
F_p1 = [0 0 0];
F_p2 = [-0.7 -0.5 0];
ex1_tau(:,4) = compute_tau1(angle1,angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 1.4
F_p2 = [0 0 0];
F_p3 = [1.5 -0.3 0];
M_p1 = [0 0 1.2];
ex1_tau(:,5) = compute_tau1(angle1,angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 1.5
angle2 = -60;
F_p2 = [1.2 -0.2 0];
F_p3 = [-0.4 1.2 0];
M_p1 = [0 0 0];
g = [0 -9.81 0];
ex1_tau(:,6) = compute_tau1(angle1,angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 2.1
angle2 = 45;
F_p1 = [0 0 0];
F_p2 = [0 0 0];
F_p3 = [0 0 0];
M_p1 = [0 0 0];
M_p2 = [0 0 0];
M_p3 = [0 0 0];
g = [0 -9.81 0];
ex2_tau_and_force(:,1) = compute_tau2(-90+angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 2.2
angle2 = 90;
ex2_tau_and_force(:,2) = compute_tau2(-90+angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 2.3
angle2 = 45;
F_p1 = [-0.8 -0.8 0];
g = [0 0 0];
ex2_tau_and_force(:,3) = compute_tau2(-90+angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 2.4
F_p1 = [0 0 0];
F_p2 = [-0.8 -0.2 0];
M_p1 = [0 0 0.5 ];
ex2_tau_and_force(:,4) = compute_tau2(-90+angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
%% Ex 2.5
F_p2 = [1 -0.4 0];
F_p1 = [0.5 -0.6 0];
M_p1 = [0 0 0];
g = [0 -9.81 0];
ex2_tau_and_force(:,5) = compute_tau2(-90+angle2,g,F_p1,F_p2,F_p3,M_p1,M_p2,M_p3);
clearvars -except ex1_tau ex2_tau_and_force