九州大学 工学府 機械系専攻 2015年度 機械力学 問2
Author
祭音Myyura (refined with GPT 6 Astra)
Description
長さ l l l 、質量を無視できる剛体棒の左端から a a a (0 ≤ a ≤ l 0\le a\le l 0 ≤ a ≤ l )の位置に質点 m m m を取り付ける。左端をばね定数 3 k 3k 3 k 、右端を k k k のばねで支持する。平衡時は棒が水平であり、平衡からの微小振動では重力の影響を考えなくてよい。ニュートンの第二法則によって 運動方程式を求め、以下に答えよ。(50 点)
質点の鉛直上向き変位を y y y 、棒の反時計回り角変位を θ \theta θ とする。a = l / 2 a=l/2 a = l /2 のとき、固有角振動数と自由振動時の y , θ y,\theta y , θ の関係を求めよ。
固有角振動数を最大にする a a a とその最大値、およびそのときの y , θ y,\theta y , θ の関係を求めよ。
質量 m m m を m / 2 m/2 m /2 ずつ棒の両端に分けて取り付けた図2の系について、固有角振動数と固有モードを求めよ。この場合は棒の中心 G の上向き変位を y y y 、反時計回り角変位を θ \theta θ とする。
题目描述
无质量刚杆长 l l l ,距左端 a a a 处装质量 m m m ,左、右两端分别用刚度 3 k , k 3k,k 3 k , k 的弹簧支撑,静平衡时杆水平。位移均相对静平衡,采用牛顿第二定律。以质点向上位移 y y y 和杆逆时针转角 θ \theta θ 为坐标,求 a = l / 2 a=l/2 a = l /2 时的固有频率及位移关系,再求使频率最大的 a a a 和相应关系。最后把质量等分为两端各 m / 2 m/2 m /2 ,改用杆中心的 y y y 与转角 θ \theta θ ,求频率及模态。
Kai
(1)
両端の上向き変位は y L = y − a θ y_L=y-a\theta y L = y − a θ 、y R = y + ( l − a ) θ y_R=y+(l-a)\theta y R = y + ( l − a ) θ 。質点の並進運動と質点まわりのモーメントの式は
m y ¨ = − 3 k y L − k y R = − 4 k y + k ( 4 a − l ) θ , m\ddot y=-3ky_L-ky_R=-4ky+k(4a-l)\theta, m y ¨ = − 3 k y L − k y R = − 4 k y + k ( 4 a − l ) θ ,
0 = 3 k a y L − k ( l − a ) y R = k ( 4 a − l ) y − k { 3 a 2 + ( l − a ) 2 } θ . 0=3ka\,y_L-k(l-a)y_R
=k(4a-l)y-k\{3a^2+(l-a)^2\}\theta. 0 = 3 ka y L − k ( l − a ) y R = k ( 4 a − l ) y − k { 3 a 2 + ( l − a ) 2 } θ .
第2式の左辺がゼロなのは、棒に質量がなく、質点まわりの慣性モーメントもゼロだからである。a = l / 2 a=l/2 a = l /2 では第2式から θ = y / l \theta=y/l θ = y / l 、第1式から m y ¨ + 3 k y = 0 m\ddot y+3ky=0 m y ¨ + 3 k y = 0 。したがって
ω = 3 k m , y = l θ . \boxed{\omega=\sqrt{\frac{3k}{m}},\qquad y=l\theta.} ω = m 3 k , y = lθ .
(2)
D ( a ) = 3 a 2 + ( l − a ) 2 D(a)=3a^2+(l-a)^2 D ( a ) = 3 a 2 + ( l − a ) 2 とおくと
θ = 4 a − l D ( a ) y , m y ¨ + 3 k l 2 D ( a ) y = 0 , \theta=\frac{4a-l}{D(a)}y,\qquad
m\ddot y+\frac{3kl^2}{D(a)}y=0, θ = D ( a ) 4 a − l y , m y ¨ + D ( a ) 3 k l 2 y = 0 ,
D ( a ) = 4 ( a − l 4 ) 2 + 3 l 2 4 . D(a)=4\left(a-\frac l4\right)^2+\frac{3l^2}{4}. D ( a ) = 4 ( a − 4 l ) 2 + 4 3 l 2 .
よって
a = l 4 , ω max = 2 k m , θ = 0. \boxed{a=\frac l4,\qquad\omega_{\max}=2\sqrt{\frac{k}{m}},\qquad\theta=0.} a = 4 l , ω m a x = 2 m k , θ = 0.
このとき棒は水平のまま並進する。
(3)
重心まわりの慣性モーメントは I G = m l 2 / 4 I_G=ml^2/4 I G = m l 2 /4 。ニュートンの並進・回転の式をまとめると
( m 0 0 m l 2 / 4 ) ( y ¨ θ ¨ ) + ( 4 k − k l − k l k l 2 ) ( y θ ) = 0 . \begin{pmatrix}m&0\\0&ml^2/4\end{pmatrix}
\begin{pmatrix}\ddot y\\\ddot\theta\end{pmatrix}
+\begin{pmatrix}4k&-kl\\-kl&kl^2\end{pmatrix}
\begin{pmatrix}y\\\theta\end{pmatrix}=\boldsymbol0. ( m 0 0 m l 2 /4 ) ( y ¨ θ ¨ ) + ( 4 k − k l − k l k l 2 ) ( y θ ) = 0 .
特性方程式は
( 4 k − m ω 2 ) 2 − 4 k 2 = 0. (4k-m\omega^2)^2-4k^2=0. ( 4 k − m ω 2 ) 2 − 4 k 2 = 0.
したがって、任意の非零定数倍を許して
ω 1 = 2 k m , ϕ 1 = ( 1 2 / l ) , ω 2 = 6 k m , ϕ 2 = ( 1 − 2 / l ) . \boxed{\omega_1=\sqrt{\frac{2k}{m}},\quad
\boldsymbol\phi_1=\begin{pmatrix}1\\2/l\end{pmatrix},\qquad
\omega_2=\sqrt{\frac{6k}{m}},\quad
\boldsymbol\phi_2=\begin{pmatrix}1\\-2/l\end{pmatrix}.} ω 1 = m 2 k , ϕ 1 = ( 1 2/ l ) , ω 2 = m 6 k , ϕ 2 = ( 1 − 2/ l ) .
第1モードでは左端が、第2モードでは右端が静止する。