1. Synteza strukturalna i geometryczna mechanizmu
1. 1. Budowa łańcucha kinematycznego
1. 2. Ruchliwość i klasa mechanizmu
Ruchliwość mechanizmu:
Ilość członów: 3
Liczba par kinematycznych klasy 4: 0
Liczba par kinematycznych klasy 5: 4 [(0,1),(1,2),(2,3),(3,0)]
w=3*n-p4-2*p5
w=3*3-0-2*4
w=1
Klasa mechanizmu:
Po odłączeniu członu napędzającego 1 pozostałe człony 2 i 3 tworzą grupę strukturalną.
Ruchliwość grupy strukturalnej po połączeniu jej członów ruchomych
z podstawą n=2, p5=3 (0,2),(2,3),(3,0)
wgr=3*n-2*p5=0
Grupa strukturalna jest grupą klasy 2 postaci 1.
Analizowany mechanizm składa się z członu napędzającego 1 i grupy strukturalnej klasy 2, jest zatem mechanizmem klasy 2.
Nazwa strukturalna mechanizmu: mechanizm suwakowo – korbowy
Nazwa funkcjonalna mechanizmu: mechanizm chwytaka
Ograniczenia geometryczne:
|AB|<|BD|
Smax=$\sqrt{{|AE|}^{2} + {|CE|}^{2}}$
Smin=|AB|+|BD|
Model mechanizmu w programie SAM:
Przyjęte dane do obliczeń kinematycznych:
s0=330[mm]=0,033[m]
Δs=400[mm]=0,04[m]
V1=200[mm/s]
sA=s0+V1*t
Czas ruchu członu: t=0,5[s]
Położenie członu napędzającego: sA=430[mm]=0,043[m]
lAE= 220[mm]=0,022[m]
lEB= 200[mm]=0,02[m]
lEC= 300[mm]=0,03[m]
lBD= 620[mm]=0,062[m]
2. Analiza kinematyczna mechanizmu
2.1. Analiza kinematyczna mechanizmu metodą grafoanalityczną
Obliczenie prędkości punktu B:
VB=VA+VBA
VB | DB
VA ||AD
VA =V1
VBA | AB
Obliczenie prędkości punktu C:
VC=VA+VCA
VA=V1
VA ||AD
VCA=ω1*|CA|=170,05[mm/s] $\omega_{1} = \frac{V_{\text{BA}}}{|AB|}$=0,4571[rad/s]
VCA | AC
Obliczenie przyspieszenia punktu B:
aB=aA+aBA
aA=0
aBA=anBA+atBA
aB= anBA+atBA
aB= anB+atB
anB+atB = anBA+atBA
anBA || AB
atBA | AB
anB || BD
atB | BD
$$\omega_{2} = \frac{VB}{|BD|} = 0,1210\lbrack rad/s\rbrack$$
anBA = ω12|AB| = 59, 7550[mm/s2]
anB = ω22|BD| = 9, 07742[mm/s2]
Obliczenie przyspieszenia punktu C:
aC=aA+aCA
aA=0
aCA=anCA+atCA
aC= aB+anCB+atCB
aB+anCB+atCB = anCA+atCA
anCB || CB
atCB | CB
anCA || AC
atCA | AC
anCA = ω12|AC| = 61, 308[mm/s2]
anCB = ω12|BC| = 59, 1855[mm/s2]
2.2. Analiza kinematyczna mechanizmu metoda analityczną
Dane: l1, l2, l3(t), φ3 = π
Szukane: VB, aB, φ1
l1+l2+ l3=0
$$\left\{ \begin{matrix}
\ l_{1}\text{cos\ }\varphi_{1} + l_{2}\text{cos\ }\varphi_{2} - \ l_{3} = 0 \\
l_{1}\text{sin\ }\varphi_{1} + l_{2}\text{sin\ }\varphi_{2} = 0 \\
\end{matrix} \right.\ $$
$$\text{sin\ }\varphi_{2} = - \frac{l_{1}}{l_{2}}\sin\varphi_{1}$$
$$\text{cos\ }\varphi_{2} = \sqrt{1 - \sin^{2}\varphi_{2}} = \sqrt{1 - \frac{l_{1}^{2}}{l_{2}^{2}}\sin^{2}\varphi_{1}}$$
$$l_{1}\text{cos\ }\varphi_{1} - l_{2}\sqrt{1 - \frac{l_{1}^{2}}{l_{2}^{2}}\sin^{2}\varphi_{1}}{- \ l}_{3} = 0$$
Po podniesienu do kwadratu:
$$l_{2}^{2}\left( 1 - \frac{l_{1}^{2}}{l_{2}^{2}}\sin^{2}\varphi_{1} \right) = \ l_{3}^{2} - \ 2l_{1}l_{2}\cos\varphi_{1} + l_{1}^{2}\cos^{2}\varphi_{1}$$
Po zastosowaniu wzoru (jedynka trygonometryczna):
l22 − l12 + l12cos2φ1 = l32 − 2l1l2cosφ1 + l12cos2φ1
i ostatecznie:
$$\cos\varphi_{1} = \frac{l_{2}^{2} - l_{1}^{2} - l_{3}^{2}}{- 2l_{1}l_{3}}$$
po podstawieniu: t=0,5[s], l1 = 420[mm], l2=620[mm], l3 = 610[mm]
cosφ1=0,3203 φ1= 71,319[⁰]
Obliczanie φ2
$$\cos\varphi_{2} = \ \frac{l_{3} - l_{1}\cos\varphi_{1}}{l_{2}}$$
cosφ2=0,7669 φ2=309,924[⁰]
Wyznaczenie prędkości kątowych:
−ω1l1sinφ1−ω2l2sinφ2 − VA = 0
Obliczenie ω1(obracamy układ współrzędnych o kąt φ2)
−ω1l1sin(φ1− φ2) − ω2l2sin(φ2−φ2)−VAcosφ2 = 0
$$\omega_{1} = \frac{- V_{A}\cos\varphi_{2}}{l_{1}{sin(\varphi}_{1}{- \ \varphi}_{2})}$$
ω1 = 0, 4277[rad/s]
Obliczenie ω2(obracamy układ współrzędnych o kąt φ1)
−ω1l1sin(φ1− φ1) − ω2l2sin(φ2−φ1)−VAcosφ1 = 0
$$\omega_{2} = \frac{- V_{A}\cos\varphi_{1}}{l_{2}{sin(\varphi}_{2}{- \ \varphi}_{1})}$$
ω20, 4277[rad/s]
Obliczenie przyśpieszeń kątowych:
−ε1l1sinφ1−ω12l1cosφ1 − ε2l2sinφ2 − ω22l2cosφ2 = 0
Obliczenie ε1(obracamy układ współrzędnych o kąt φ2)
−ε1l1sin(φ1 − φ2)−ω12l1cos(φ1− φ2)−ε2l2sin(φ2 − φ2)−ω22l2cos(φ2 − φ2)=0
$$\varepsilon_{1} = \frac{\omega_{2}^{2}l_{2} - \omega_{1}^{2}l_{1}\cos{(\varphi}_{1}{- \ \varphi}_{2})}{l_{1}\sin{(\varphi}_{1}{- \ \varphi}_{2})}$$
ε1 = 0, 02936[rad/s2]
Obliczenie ε2(obracamy układ współrzędnych o kąt φ1)
−ε1l1sin(φ1 − φ1)−ω12l1cos(φ1− φ1)−ε2l2sin(φ2 − φ1)−ω22l2cos(φ2 − φ1)=0
$$\varepsilon_{2} = \frac{\omega_{1}^{2}l_{1} - \omega_{2}^{2}l_{2}\cos{(\varphi}_{2}{- \ \varphi}_{1})}{l_{2}\sin{(\varphi}_{2}{- \ \varphi}_{1})}$$
ε2 = 0, 10893[rad/s2]
Obliczenie VB i aB
$${\dot{\varphi}}_{1} = \omega_{1}\ ,\ {\ddot{\varphi}}_{1} = \varepsilon_{1}$$
rB=SA+l1
rBX=SA+l1 *cosφ1
rBY=l1 *sinφ1
$$V_{\text{BX}} = \ \frac{Dr_{\text{BX}}}{\text{Dt}} = V_{A} - \sin\varphi_{1}$$
$$V_{\text{BY}} = \ \frac{Dr_{\text{BY}}}{\text{Dt}} = l_{1}{\dot{\varphi}}_{1}\cos\varphi_{1}$$
$$V_{B} = \ \sqrt{V_{\text{BX}}^{2} + V_{\text{BY}}^{2}\ } = \sqrt{l_{1}^{2}{\dot{\varphi}}_{1}^{2}\cos^{2}\varphi_{1} + \left( V_{A} - \sin\varphi_{1} \right)^{2}}$$
VB = 76, 0127[mm/s]
$$a_{\text{BX}} = \ \frac{{D^{2}r}_{\text{BX}}}{Dt^{2}} = 0 - l_{1}{\ddot{\varphi}}_{1}\sin\varphi_{1} - \ l_{1}{\dot{\varphi}}_{1}^{2}\cos\varphi_{1} = {- l_{1}{\ddot{\varphi}}_{1}\sin\varphi_{1} - l}_{1}{\dot{\varphi}}_{1}^{2}\cos\varphi_{1}$$
$$a_{\text{BY}} = \ \frac{D^{2}r_{\text{BY}}}{Dt^{2}} = \ l_{1}{\ddot{\varphi}}_{1}\cos\varphi_{1} - \ l_{1}{\dot{\varphi}}_{1}^{2}\sin\varphi_{1}$$
$$a_{B} = \ \sqrt{a_{\text{BX}}^{2} + a_{\text{BY}}^{2}} = \sqrt{\left( - l_{1}{\ddot{\varphi}}_{1}\sin\varphi_{1} - \ l_{1}{\dot{\varphi}}_{1}^{2}\sin\varphi_{1} \right)^{2} + \left( \ l_{1}{\ddot{\varphi}}_{1}\cos\varphi_{1} - \ l_{1}{\dot{\varphi}}_{1}^{2}\sin\varphi_{1} \right)^{2}\text{\ \ }}$$
aB = 65, 9931[mm/s2]
$$l_{4} = \sqrt{\frac{l_{1}^{2}}{4} + l_{\text{EC}}^{2}}$$
φ4 = φ1 + α
rC=SA+l1
rCX=SA+l4 *cosφ4
rCY=l4 *sinφ4
$${\dot{\text{\ φ}}}_{4} = \omega_{1}{,\ \ \ddot{\varphi}}_{4} = \varepsilon_{1}\ $$
$$V_{\text{CX}} = \ \frac{Dr_{\text{CX}}}{\text{Dt}} = V_{A} - \sin\varphi_{4}$$
$$V_{\text{CY}} = \ \frac{Dr_{\text{CY}}}{\text{Dt}} = l_{4}{\dot{\varphi}}_{4}\cos\varphi_{4}$$
$$V_{C} = \ \sqrt{V_{\text{CX}}^{2} + V_{\text{CY}}^{2}\ } = \sqrt{l_{4}^{2}{\dot{\varphi}}_{4}^{2}\cos^{2}\varphi_{4} + \left( V_{A} - \sin\varphi_{4} \right)^{2}}$$
VC = 114, 0691[mm/s]
$$a_{\text{CX}} = \ \frac{{D^{2}r}_{\text{CX}}}{Dt^{2}} = 0 - l_{4}{\ddot{\varphi}}_{4}\sin\varphi_{4} - \ l_{4}{\dot{\varphi}}_{4}^{2}\cos\varphi_{4}$$
$$a_{\text{CY}} = \ \frac{D^{2}r_{\text{CY}}}{Dt^{2}} = \ l_{4}{\ddot{\varphi}}_{4}\cos\varphi_{4} - \ l_{4}{\dot{\varphi}}_{4}^{2}\sin\varphi_{4}$$
$$a_{C} = \ \sqrt{a_{\text{CX}}^{2} + a_{\text{CY}}^{2}} = \sqrt{\left( l_{4}{\ddot{\varphi}}_{4}\sin\varphi_{4} - \ l_{4}{\dot{\varphi}}_{4}^{2}\sin\varphi_{4} \right)^{2} + \left( \ l_{4}{\ddot{\varphi}}_{4}\cos\varphi_{4} - l_{4}{\dot{\varphi}}_{4}^{2}\cos\varphi_{4} \right)^{2}\text{\ \ }}$$
aC = 58, 2587[mm/s2]
2.3. Analiza kinematyczna mechanizmu w programie SAM
aC = 57, 57886[mm/s2]
VC = 113, 72608[mm/s]
aB = 65, 44174[mm/s2]
VB = 75, 44050[mm/s]
ε1 = 0, 02462[rad/s2]
ε2 = 0, 10473[rad/s2]
ω1 = 0, 39168[rad/s]
ω2 = −0, 12202[rad/s]
2.4 Porównanie wyników
Lp | Parametr | SAM | Metoda grafoanalityczna | Metoda analityczna |
---|---|---|---|---|
1 | ω1 |
0, 39168 |
0,4571 | 0, 4277 |
2 | ω2 |
−0, 12202 |
0, 1210 |
−0, 1210 |
3 | ε1 |
0, 02462 |
0, 03134 |
0, 02936 |
4 | ε2 |
0, 10473 |
0, 11207 |
0, 10893 |
5 | VB |
75, 44050 |
76,5 | 76, 0127 |
6 | aB |
65, 44174 |
66,5 | 65, 9931 |
7 | VC |
113, 72608 |
115,1 | 114, 0691 |
8 | aC |
57, 57886 |
59,2 | 58, 2587 |
3. Analiza kinetostatyczna mechanizmu
Obliczenie sił ciężkości, sił bezwładności i momentów od sił bezwładności oraz przyjęcie zewnętrznych sił i momentów oporu
Moment bezwładności:
m2=2[kg]
JS2=$m_{2}*\frac{l_{2}^{2}}{12}$
JS2=0.0294[kg*m2]
P2=2[N]
M3=1[Nm]
B2=-m2aS2
B2= m2aS2=2*0,03272=0,065[N]
MB2=- JS2ε1
MB2= JS2ε1=0,0294*0,1047=0,0031[Nm]
G2=m2*g=19,62[N]
$\sum_{}^{}{P_{i} = 0\ ,\ \ \ R_{12}^{n}} + R_{12}^{t} + B_{2} + G_{2} + P_{2} + R_{30}^{n} + R_{30}^{t} = 0$
$$\sum_{}^{}{M_{iB(2)} = 0,\ \ {- R}_{12}^{t}\left| \text{AB} \right| - B_{2}h_{2} + G_{2}h_{3} + P_{2}h_{4} + M_{B2} = 0\ \ }$$
$$R_{12}^{t} = \frac{{- B}_{2}h_{2} + G_{2}h_{3} + P_{2}h_{4} + M_{B2}}{|AB|}$$
R12t = 6, 5061[N]
$$\sum_{}^{}{M_{iB(3)} = 0,\ \ R_{30}^{t}\left| \text{BD} \right| + M_{3} = 0\ \ }$$
$$R_{03}^{t} = \frac{- M_{3}}{|BD|}$$
R03t = 1, 6129[N]
R12n + R12t + B2 + G2 + P2 + R30n + R30t = 0
Równanie równowagi sił zewnętrznych i reakcji działających na człon napędzający
PR1 + R01 + R12 = 0
Po odczytaniu:
PR1 = 2[N]
R12 = 17, 5[N]
Analiza kinetostatyczna metodą mocy chwilowych.
PR1VA + B2VS2 + G2VS2 + MB2ω1 + P2VC + M3ω2 = 0
PR1VA + B2VS2cos(70, 50)+G2VS2cos(115, 40)+MB2ω1 + P2VCcos(33, 50)+M3ω2 = 0
$$P_{R1} = \frac{{- \ B}_{2}V_{S2}\cos\left( {70,5}^{0} \right) - G_{2}V_{S2}\cos\left( {115,4}^{0} \right) - M_{B2}\omega_{1} - P_{2}V_{C}\cos\left( {33,5}^{0} \right) - M_{3}\omega_{2}}{V_{A}}$$
PR1 = 2, 2367[N]
Analiza kinetostatyczna w programie SAM
Porównanie Wyników:
Metoda grafoanalityczna | Metoda mocy chwilowych | Obliczenie w programie SAM |
---|---|---|
2[N] | 2,2367[N] | 2,158511[N] |