九州大学 工学府 機械系専攻 2015年8月実施 機械力学 問3
Author
祭音Myyura (refined with GPT 6 Astra)
Description
図のように、両端が固定された、密度 ρ \rho ρ 、縦弾性係数 E E E 、長さ l l l の弾性丸棒がある。弾性丸棒の長手方向に左端を原点とする座標 x x x をとり、縦振動を考える。断面積 A ( x ) A(x) A ( x ) は x x x の関数であり、座標 x x x 、時刻 t t t における弾性丸棒の長手方向変位を u ( x , t ) u(x,t) u ( x , t ) とする。以下の問いに答えよ。ただし、c = E / ρ c=\sqrt{E/\rho} c = E / ρ を使用してよい。(50点)
(1) 座標 x x x における弾性丸棒の縦振動の運動方程式を導出せよ。
(2) A ( x ) = A 0 e − 2 γ x A(x)=A_0e^{-2\gamma x} A ( x ) = A 0 e − 2 γ x として (1) の運動方程式を書き直せ。ただし、A 0 , γ A_0,\gamma A 0 , γ は定数とする。
(3) u ( x , t ) = U ( x ) sin ω t u(x,t)=U(x)\sin\omega t u ( x , t ) = U ( x ) sin ω t と仮定したとき、U ( x ) U(x) U ( x ) が満たすべき微分方程式を求めよ。また、境界条件を満たし得る U ( x ) U(x) U ( x ) の一般解を求めよ。
(4) r r r 次の固有角振動数 ω r \omega_r ω r を求めよ。
(5) 断面積が一様な場合と比較して、固有角振動数の大小関係を求めよ。
题目描述
如图,有一密度为 ρ \rho ρ 、纵向弹性模量为 E E E 、长度为 l l l 、两端固定的弹性圆杆。以左端为原点,沿杆轴向取坐标 x x x ,考虑纵振动。截面积 A ( x ) A(x) A ( x ) 随 x x x 变化,坐标 x x x 、时刻 t t t 处的轴向位移为 u ( x , t ) u(x,t) u ( x , t ) 。回答下列问题,可使用 c = E / ρ c=\sqrt{E/\rho} c = E / ρ 。(50分)
(1) 推导弹性圆杆在坐标 x x x 处纵振动的运动方程。
(2) 令 A ( x ) = A 0 e − 2 γ x A(x)=A_0e^{-2\gamma x} A ( x ) = A 0 e − 2 γ x ,改写 (1) 的运动方程,其中 A 0 , γ A_0,\gamma A 0 , γ 为常数。
(3) 假设 u ( x , t ) = U ( x ) sin ω t u(x,t)=U(x)\sin\omega t u ( x , t ) = U ( x ) sin ω t ,求 U ( x ) U(x) U ( x ) 应满足的微分方程,并求能够满足边界条件的 U ( x ) U(x) U ( x ) 通解。
(4) 求第 r r r 阶固有角频率 ω r \omega_r ω r 。
(5) 与截面积均匀的情况比较,判断固有角频率的大小关系。
Kai
(1)
軸力は N = E A ( x ) u x N=EA(x)u_x N = E A ( x ) u x 。微小要素の運動量収支から
ρ A ( x ) u t t = ∂ ∂ x ( E A ( x ) u x ) . \boxed{\rho A(x)u_{tt}=\frac{\partial}{\partial x}\left(EA(x)u_x\right).} ρ A ( x ) u tt = ∂ x ∂ ( E A ( x ) u x ) .
(2)
A ′ / A = − 2 γ A'/A=-2\gamma A ′ / A = − 2 γ なので
u t t = c 2 ( u x x − 2 γ u x ) . \boxed{u_{tt}=c^2(u_{xx}-2\gamma u_x).} u tt = c 2 ( u xx − 2 γ u x ) .
(3)
U ′ ′ − 2 γ U ′ + ω 2 c 2 U = 0. U''-2\gamma U'+\frac{\omega^2}{c^2}U=0. U ′′ − 2 γ U ′ + c 2 ω 2 U = 0.
U = e γ x Y U=e^{\gamma x}Y U = e γ x Y とおくと Y ′ ′ + ( ω 2 / c 2 − γ 2 ) Y = 0 Y''+(\omega^2/c^2-\gamma^2)Y=0 Y ′′ + ( ω 2 / c 2 − γ 2 ) Y = 0 。両端で Y = 0 Y=0 Y = 0 の非自明解があるためには ω 2 / c 2 > γ 2 \omega^2/c^2>\gamma^2 ω 2 / c 2 > γ 2 が必要である。このとき
U ( x ) = e γ x ( C cos β x + D sin β x ) , β = ω 2 / c 2 − γ 2 . \boxed{\begin{gathered}
U(x)=e^{\gamma x}(C\cos\beta x+D\sin\beta x),\\
\beta=\sqrt{\omega^2/c^2-\gamma^2}.
\end{gathered}} U ( x ) = e γ x ( C cos β x + D sin β x ) , β = ω 2 / c 2 − γ 2 .
ω 2 / c 2 = γ 2 \omega^2/c^2=\gamma^2 ω 2 / c 2 = γ 2 の場合は Y = C + D x Y=C+Dx Y = C + D x 、小さい場合は双曲線関数となり、どちらも両端条件を課すとゼロ解しか残らない。
(4)
U ( 0 ) = 0 U(0)=0 U ( 0 ) = 0 より C = 0 C=0 C = 0 、U ( l ) = 0 U(l)=0 U ( l ) = 0 より β l = r π \beta l=r\pi βl = r π (r = 1 , 2 , … r=1,2,\ldots r = 1 , 2 , … )。したがって
ω r = c γ 2 + ( r π l ) 2 , U r ( x ) ∝ e γ x sin r π x l . \boxed{\begin{gathered}
\omega_r=c\sqrt{\gamma^2+\left(\frac{r\pi}{l}\right)^2},\\
U_r(x)\propto e^{\gamma x}\sin\frac{r\pi x}{l}.
\end{gathered}} ω r = c γ 2 + ( l r π ) 2 , U r ( x ) ∝ e γ x sin l r π x .
(5)
一様断面では ω r ( 0 ) = c r π / l \omega_r^{(0)}=cr\pi/l ω r ( 0 ) = cr π / l 。よって
ω r > ω r ( 0 ) ( γ ≠ 0 ) , \boxed{\omega_r>\omega_r^{(0)}\quad(\gamma\ne0),} ω r > ω r ( 0 ) ( γ = 0 ) ,
γ = 0 \gamma=0 γ = 0 では等しい。r = 0 r=0 r = 0 は非自明な固有モードを与えない。