Error calculos cinematicos
Publicado por iker (1 intervención) el 25/05/2017 09:48:14
Buenas, el programa mostrado a continuación muestra una serie de errores si alguno es capaz de resolverlas sería de gran utilidad. Gracias de antemano
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
clear all
close all
%aldagaiak
syms x1 x2 x3 x4 x5 x6 %x1=F x2=Na x3=Nb x4=Nc x5=Ne x6=Nd
%parametroak
syms nuk m1 g b h L m2 tita VerlA A w
syms sinBeta cosBeta Ro atanB
sinBeta = (A*cos(tita)*w)/(sqrt(1+A*A*(cos(tita))*(cos(tita))*w*w));
cosBeta = 1/(sqrt(1+A*A*(cos(tita))*(cos(tita))*w*w));
Ro = (1+A*A*A*(cos(tita))*(cos(tita))*(cos(tita))*w*w*w)/(A*sin(tita)*w*w);
atanB = -((VerlA/cosBeta)*(VerlA/cosBeta)*sinBeta)/(Ro*cosBeta);
eqn1 = x1-nuk*x2-nuk*x3-sinBeta*x4-nuk*x4*cosBeta;
eqn2 = x3-x2+sinBeta*nuk*x4+cosBeta*x4-m1*g;
eqn3 = x3*b-nuk*x3*h-m1*g*(b/2);
eqn4 = x5-x6-cosBeta*nuk*x4+sinBeta*x4;
eqn5 = nuk*x5+nuk*x6-cosBeta*x4-sinBeta*nuk*x4-m2*g-m2*atanB;
eqn6 = -L*sinBeta*x4+L*cosBeta*nuk*x4-x5*(L-(L*(h/2)*cosBeta));
[x1, x2, x3, x4, x5, x6] = solve(eqn1,eqn2,eqn3,eqn4,eqn5,eqn6);
soluzio=subs([x1 x2 x3 x4 x5 x6],[nuk,m1,g,b,h,L,m2,A,w],[0.57,0.25,9.8,2*pi,0.5,0.2,0.5,20,742.6]);
figure
ezplot(soluzio(1),[0,2*pi])
grid on
title('V_b');
xlabel('tita (rad)');
ylabel('V_b (m/s)');
figure
ezplot(soluzio(2),[0,2*pi])
grid on
title('w_2');
xlabel('tita (rad)');
ylabel('w_2 (rad/s)');
figure
ezplot(soluzio(3),[0,2*pi])
grid on
title('a_b');
xlabel('tita (rad)');
ylabel('a_b (m/s^2)');
figure
ezplot(soluzio(4),[0,2*pi])
grid on
title('alfa_2');
xlabel('tita (rad)');
ylabel('x1 (rad/s^2)');
figure
ezplot(soluzio(5),[0,2*pi])
grid on
title('V_b');
xlabel('tita (rad)');
ylabel('V_b (m/s)');
figure
ezplot(soluzio(6),[0,2*pi])
grid on
title('V_b');
xlabel('tita (rad)');
ylabel('V_b (m/s)');
Valora esta pregunta
0