9 svar
455 visningar
econo behöver inte mer hjälp
econo 64 – Fd. Medlem
Postad: 23 sep 2020 16:38 Redigerad: 23 sep 2020 16:43

Vinkelhastighet uttryckt som vinkel theta

Hej! Jag har lite problem med att lösa denna typ av uppgift med hjälp av vektor produkt.

(Har löst uppgiften med hjälp av sinussatsen och deriveringsoperatorn ddt\frac{d}{dt}.)

 

Jag har försökt följande:

vB=ω0×OB=00-ω0×-rcosθrsinθ0=ω0rsinθω0rcosθ0\vec{v}_B=\vec{\omega}_0\times\overrightarrow{OB}=\begin{bmatrix}0\0\-\omega_0\end{bmatrix}\times\begin{bmatrix}-r\cos\theta\r\sin\theta\0\end{bmatrix}=\begin{bmatrix}\omega_0r\sin\theta\ \omega_0r\cos\theta\0\end{bmatrix}

 

Vektorn OA=-rcosθ-l2-r2sin2θ00\overrightarrow{OA}=\begin{bmatrix}-r\cos\theta-\sqrt{l^2-r^2\sin^2\theta}\0\0\end{bmatrix}

och därmed är vinkelhastigheten ωAB\boldsymbol{\omega}_{AB}

 

ωAB=(OB-OA)×vB\boldsymbol{\omega}_{AB}=(\overrightarrow{OB}-\overrightarrow{OA})\times\boldsymbol{v}_B

 

Matar in allt i MATLAB men får då inte rätt svar jfrt. med metoden av sinussatsen.

clear all; close all;
syms r theta l w0
W0=[0;0;-w0]
OB=[-r*cosd(theta);r*sind(theta);0]
vB=cross(W0,OB)
OA=[-r*cosd(theta)-(l^2-r^2*(sind(theta))^2)^(1/2);0;0]
AB=OB-OA
wAB=cross(AB,vB)

% Testar parametrar
r=0.5; l=1; w0=5;
fplot(subs(wAB(3)),[0,90])

 

Det rätta svaret i MATLAB:

clear all
syms theta
r = 0.5; l=1; w0=5 % Testar parametrar
omega_AB = r*w0*cosd(theta)/(l*sqrt(1-(r/l*sind(theta))^2))
fplot(omega_AB,[0,120])

econo 64 – Fd. Medlem
Postad: 23 sep 2020 17:29 Redigerad: 23 sep 2020 18:27

Jag har löst det! 

 

Jag har testat att ta bort x komposanten i vB\vec{v}_{B} så att vB=0ω0rcosθ0\vec{v}_{B}=\begin{bmatrix}0\\ \omega_0r\cos\theta\\0\end{bmatrix} där tanken är att det bara är y-komposanten som ger upphov till en vinkelhastighet vid γ\gamma. Jag har nu jämfört graferna och får följande:

Det skiljer sig lite från graferna, där den blåa grafen ska vara den rätta så att säga. Vet ej vad detta kan bero på!

JohanF 5660 – Moderator
Postad: 24 sep 2020 08:17

Jag tror inte du kan ta bort x-komposanten i vB, för den kommer att ge upphov till vinkelhastighet på ωAB också. (Enda tillfället när den inte gör det är när AB är parallell med x-axeln, dvs då täta är 0).

Jag tror att din ekvation för uträkning av ωAB ifrån den spatiellla hastigheten vBkan vara fel. Jag kan inte säga vad som är fel eftersom jag är för rostig på räkneregler för kryssprodukter. Men en enhetsanalys visar att om du multiplicerar en längd m med en hastighet ms, så får du inte en vinkelhastighet s-1.

econo 64 – Fd. Medlem
Postad: 24 sep 2020 10:01 Redigerad: 24 sep 2020 10:02
JohanF skrev:

Jag tror inte du kan ta bort x-komposanten i vB, för den kommer att ge upphov till vinkelhastighet på ωAB också. (Enda tillfället när den inte gör det är när AB är parallell med x-axeln, dvs då täta är 0).

Jag tror att din ekvation för uträkning av ωAB ifrån den spatiellla hastigheten vBkan vara fel. Jag kan inte säga vad som är fel eftersom jag är för rostig på räkneregler för kryssprodukter. Men en enhetsanalys visar att om du multiplicerar en längd m med en hastighet ms, så får du inte en vinkelhastighet s-1.

Vinkelhastigheten, ωAB\omega_{AB} är som störst då θ=0\theta=0, vi har en konstant vinkelhastighet ω0\omega_{0} vid B. Vinkelhastigheten ωAB\omega_{AB} är noll då θ=90°\theta=90^{\circ} och ωAB\omega_{AB} byter tecken då θ>90°\theta>90^{\circ}. Vinkelhastigheten ωAB\omega_{AB} är snarare noll då riktningsvektorn för vB\boldsymbol{v}_B är parallel med AB\overrightarrow{AB}

Det är väl frågan om man kan få ut vinkelhastigheten mha ω=r×v\omega=\boldsymbol{r}\times\boldsymbol{v}.

Jroth 1191 – Fd. Medlem
Postad: 24 sep 2020 10:02 Redigerad: 24 sep 2020 10:23

Vi kan uttrycka hastigheten för punkten B på två olika sätt

vB=vA+ωAB×rBA\mathbf{v}_B=\mathbf{v}_A+\mathbf{\omega}_{AB}\times \mathbf{r}_{BA}

vB=ω0×rB0\mathbf{v}_B=\mathbf{\omega}_{0}\times \mathbf{r}_{B0}

Där ω0=-ω0z^\mathbf{\omega}_0=-\omega_0\hat{z}, ωAB=ωABz^\mathbf{\omega}_{AB}=\omega_{AB}\hat{z}, rB0=(-rcos(θ),rsin(θ),0)\mathbf{r}_{B0}=(-r\cos(\theta), r\sin(\theta),0) och rBA=(lcos(φ),lsin(φ),0)\mathbf{r}_{BA}=(l\cos(\varphi), l\sin(\varphi),0).

Eftersom hastigheten för A i y-led är 0 blir ekvationen vB=vB\mathbf{v}_B=\mathbf{v}_B särskilt enkel i y-led

rω0cos(θ)=lωABcos(φ)r\omega_0\cos(\theta)=l\omega_{AB}\cos(\varphi)

ωAB=rω0cos(θ)lcos(φ)\omega_{AB}=\frac{r\omega_0\cos(\theta)}{l\cos(\varphi)}

Ett samband mellan vinkeln φ\varphi och θ\theta ges av sinussatsen

rsin(θ)=lsin(φ)r\sin(\theta)=l\sin(\varphi)

cos(φ)=1-r2sin2(θ)l2\cos(\varphi)=\sqrt{1-\frac{r^2\sin^2(\theta)}{l^2}}

Alltså blir svaret på uppiften

ωAB=rω0cos(θ)l1-r2sin2(θ)l2\omega_{AB}=\frac{r\omega_0\cos(\theta)}{l\sqrt{1-\frac{r^2\sin^2(\theta)}{l^2}}}

econo 64 – Fd. Medlem
Postad: 24 sep 2020 10:30 Redigerad: 24 sep 2020 10:35

Riktigt snyggt!

 

Jag antar då att ett uttryck för vinkelaccelerationen är om man deriverar ekvationen map tiden

rω0cosθ=lωABcosφr\omega_0\cos\theta=l\omega_{AB}\cos\varphi?

 

Eller skulle det vara bökigare att kryssa igen enligt

a=ω×(ω×r)+α×r\boldsymbol{a}=\omega\times(\omega\times\boldsymbol{r})+\alpha\times\boldsymbol{r}?

JohanF 5660 – Moderator
Postad: 24 sep 2020 11:09 Redigerad: 24 sep 2020 11:11
econo skrev:
JohanF skrev:

Jag tror inte du kan ta bort x-komposanten i vB, för den kommer att ge upphov till vinkelhastighet på ωAB också. (Enda tillfället när den inte gör det är när AB är parallell med x-axeln, dvs då täta är 0).

Jag tror att din ekvation för uträkning av ωAB ifrån den spatiellla hastigheten vBkan vara fel. Jag kan inte säga vad som är fel eftersom jag är för rostig på räkneregler för kryssprodukter. Men en enhetsanalys visar att om du multiplicerar en längd m med en hastighet ms, så får du inte en vinkelhastighet s-1.

Vinkelhastigheten, ωAB\omega_{AB} är som störst då θ=0\theta=0, vi har en konstant vinkelhastighet ω0\omega_{0} vid B. Vinkelhastigheten ωAB\omega_{AB} är noll då θ=90°\theta=90^{\circ} och ωAB\omega_{AB} byter tecken då θ>90°\theta>90^{\circ}. Vinkelhastigheten ωAB\omega_{AB} är snarare noll då riktningsvektorn för vB\boldsymbol{v}_B är parallel med AB\overrightarrow{AB}

Det är väl frågan om man kan få ut vinkelhastigheten mha ω=r×v\omega=\boldsymbol{r}\times\boldsymbol{v}.

Ja. Det är det som är frågan. Jag tror inte du kan beräkna vinkelhastigheten sådär, eftersom enhetsanalysen visar att det blir fel.

SaintVenant 3956
Postad: 25 sep 2020 00:28
econo skrev:

Riktigt snyggt!

 

Jag antar då att ett uttryck för vinkelaccelerationen är om man deriverar ekvationen map tiden

rω0cosθ=lωABcosφr\omega_0\cos\theta=l\omega_{AB}\cos\varphi?

Då gör du det du redan gjort på ett mer omständligt sätt. 

Eller skulle det vara bökigare att kryssa igen enligt

a=ω×(ω×r)+α×r\boldsymbol{a}=\omega\times(\omega\times\boldsymbol{r})+\alpha\times\boldsymbol{r}?

Du kan ställa upp accelerationen enligt tillvägagångssättet de går igenom i avsnitt 5/6. Det finns flera sample problem där som visar hur man gör.

SaintVenant 3956
Postad: 25 sep 2020 01:36 Redigerad: 25 sep 2020 01:41

En exempellösning:

aB=aO+αBO×rOB+ωBO×ωBO×rOB

Här har vi att aO=αBO=0 eftersom punkten O är fix och vinkelhastigheten ωBO\boldsymbol{\omega}_{BO} är konstant. Vi får:

aB=ω0k×ω0k×-rcosθi+rsinθj=rω02cosθi-sinθj

Vi ställer nu upp ett uttryck för aA\boldsymbol{a}_{A}:

aA=aB+αAB×rBA+ωAB×ωAB×rBA

Vi får:

 aAi=rω02cosθi-sinθj+αABk×-lcosφi-lsinφj         +ωABk×ωABk×-lcosφi-lsinφj       =rω02cosθi-sinθj+lαABsinφi-cosφj+lωAB2cosφi+sinφj

Detta ger om du jämför i\boldsymbol{i}- och j\boldsymbol{j}-termer:

aA=rω02cosθ+lαABsinφ+lωAB2cosφ0=-rω02sinθ-lαABcosφ+lωAB2sinφ

Den andra ekvationen ger oss:

αAB=lωAB2sinφ-rω02sinθlcosφ\alpha_{AB}=\dfrac{l\omega_{AB}^{2}\sin\varphi-r\omega_{0}^{2}\sin\theta}{l\cos\varphi}

Med det som Jroth tagit fram kan du ta fram ett uttryck som enbart beror på ursprungsvariablerna θ\theta, ω0\omega_{0}, rr och ll.

econo 64 – Fd. Medlem
Postad: 11 okt 2020 18:24 Redigerad: 11 okt 2020 18:32

Om man inte ska inblanda sinussats eller dylikt kan uppgiften lösas med vektorer enligt:

 

OA=-rcosθ-l2-r2sinθ00,OB=-rcosθrsinθ0\vec{OA}=\begin{bmatrix}-r\cos\theta-\sqrt{l^2-r^2\sin\theta}\\ 0\\ 0\end{bmatrix}, \overrightarrow{OB}=\begin{bmatrix}-r\cos\theta\\ r\sin\theta\\ 0\end{bmatrix}.

BA=OA-OB\overrightarrow{BA}=\overrightarrow{OA}-\overrightarrow{OB}

Vi har att

vx00\begin{bmatrix}v_x\\ 0 \\0 \end{bmatrix}=00ω0×-rcosθrsinθ0+00-ωAB×-l2-r2sin2θ-rsinθ0\begin{bmatrix}0\\ 0\\ \omega_0\end{bmatrix}\times\begin{bmatrix}-r\cos\theta\\ r\sin\theta\\ 0\end{bmatrix}+\begin{bmatrix}0\\ 0\\ -\omega_{AB}\end{bmatrix}\times\begin{bmatrix}-\sqrt{l^2-r^2\sin^2\theta}\\ -r\sin\theta\\ 0\end{bmatrix}

 

Detta ger


vAx00=-rω0sinθ-rω0sinθ0+-rωABsinθωABl2-r2sin20.\begin{bmatrix}v_{Ax}\\ 0\\ 0\end{bmatrix} = \begin{bmatrix}-r\omega_0\sin\theta\\ -r\omega_0\sin\theta\\ 0\end{bmatrix}+\begin{bmatrix}-r\omega_{AB}\sin\theta\\ \omega_{AB}\sqrt{l^2-r^2\sin^2}\\ 0\end{bmatrix}.


Lösning av systemet ger vxv_x resp. ωAB\omega_{AB} uttryckt i θ\theta så att

 

ωAB=ω0rcosθl2-r2sin2θ\omega_{AB}=\frac{\omega_0r\cos\theta}{\sqrt{l^2-r^2\sin^2\theta}}

 

Uttrycket för αAB\alpha_{AB} fås mha av aB=aO+αBO×rOB+ωBO×ωBO×rOB\boldsymbol{a}_{B}=\boldsymbol{a}_{O}+\boldsymbol{\alpha}_{\mathrm{BO}} \times \mathbf{r}_{\mathrm{OB}}+\boldsymbol{\omega}_{\mathrm{BO}} \times\left(\boldsymbol{\omega}_{\mathrm{BO}} \times \mathbf{r}_{\mathrm{OB}}\right)

 

resten är historia...

Svara
Close