跳到主要内容

名古屋大学 工学研究科 機械航空系 2017年8月実施 専門部門 問題3

Author​

祭音Myyura (co-authored with GPT 6 Astra)

Description​

問題3.

(1) ばねとダンパが接続された質量 mm の台車を図1に示す。ばねとダンパの質量や,台車と地面との間の摩擦は,無視できる。k1,k2k_1,k_2 はばね定数,c1,c2c_1,c_2 は粘性減衰係数である。u1,y1,y2u_1,y_1,y_2 は変位であり,台車が静止しているとき変位は全て0とする。

図1 ばねとダンパに接続された台車

1) u1u_1 と y2y_2 を0に固定したとき,台車の運動方程式を示せ。

2) 問1)で求めた運動方程式の特性方程式から特性根を示せ。また,減衰比は1未満と仮定し,特性根から自由振動の減衰固有角振動数と減衰比を求めよ。

3) u1u_1 と y2y_2 を0に固定しないとき,台車の運動方程式を示せ。

4) y2y_2 を0に固定しないとき,u1=Asin⁡(ωt)u_1=A\sin(\omega t) を加えた。ここで,AA と ω\omega は正の定数,tt は時刻である。y1y_1 の強制振動解の振幅を示せ。

(2) 伸縮とねじりの振動が同時に生じるばね振り子を図2に示す。釣り合いの位置からみて,伸縮の変位を yy で表し,ねじり角を θ\theta で表す。このとき運動方程式は

{my¨+kyy+kdθ=0,Iθ¨+kdy+kθθ=0\begin{cases} m\ddot y+k_yy+k_d\theta=0,\\ I\ddot\theta+k_dy+k_\theta\theta=0 \end{cases}

で表されるとする。ここで,mm はおもりの質量,II はねじり振動軸まわりのおもりの慣性モーメント,kyk_y は伸縮のばね定数,kθk_\theta はねじりのばね定数を表す。ばねは伸縮するとねじれるため,kdk_d は伸縮とねじりとの連成を表す定数である。なお kd2<kykθk_d^2<k_yk_\theta であり,また kyk_y と kθk_\theta は

kym=kθI=ωn(ωnは定数)\sqrt{\frac{k_y}{m}}=\sqrt{\frac{k_\theta}{I}}=\omega_n \qquad(\omega_n\text{は定数})

を満たすとする。

図2 伸縮とねじりが連成するばね振り子

1) 与えられた運動方程式をもとに特性方程式を示せ。

2) ばね振り子の固有角振動数を全て求めよ。解は m,I,kd,ωnm,I,k_d,\omega_n を用いて表せ。

3) 伸縮の自由振動を観察したところ,振幅が徐々に小さくなっていき,振幅がほぼ0になると続いて徐々に大きくなり,ある程度大きくなると再び徐々に小さくなるという現象が見られた。うなりと考えられるこの現象の周期 TT を求めよ。

4) kd2≪kykθk_d^2\ll k_yk_\theta の場合,T≫2πωnT\gg\dfrac{2\pi}{\omega_n} となることを示せ。

题目描述​

问题3.

(1) 图1为连接了弹簧和阻尼器、质量为 mm 的小车。忽略弹簧和阻尼器的质量及小车与地面间的摩擦。k1,k2k_1,k_2 为弹簧常数,c1,c2c_1,c_2 为黏性阻尼系数。u1,y1,y2u_1,y_1,y_2 为位移,令小车静止时所有位移均为0。

1) 当 u1,y2u_1,y_2 固定为0时,写出小车的运动方程。

2) 由第1)问运动方程的特征方程求特征根。假定阻尼比小于1,由特征根求自由振动的有阻尼固有角频率和阻尼比。

3) 当 u1,y2u_1,y_2 不固定为0时,写出小车的运动方程。

4) 当 y2y_2 不固定为0时,施加 u1=Asin⁡(ωt)u_1=A\sin(\omega t)。其中 A,ωA,\omega 为正常数,tt 为时间。给出 y1y_1 的受迫振动解的振幅。

(2) 图2为同时发生伸缩与扭转振动的弹簧摆。以平衡位置为基准,用 yy 表示伸缩位移,θ\theta 表示扭转角。设运动方程为

{my¨+kyy+kdθ=0,Iθ¨+kdy+kθθ=0.\begin{cases} m\ddot y+k_yy+k_d\theta=0,\\ I\ddot\theta+k_dy+k_\theta\theta=0. \end{cases}

其中 mm 为重物质量,II 为重物绕扭振轴的转动惯量,kyk_y 为伸缩弹簧常数,kθk_\theta 为扭转弹簧常数。弹簧伸缩时也会扭转,kdk_d 为表示伸缩与扭转耦合的常数。设 kd2<kykθk_d^2<k_yk_\theta,且 ky,kθk_y,k_\theta 满足

kym=kθI=ωn,ωn为常数.\sqrt{\frac{k_y}{m}}=\sqrt{\frac{k_\theta}{I}}=\omega_n, \qquad\omega_n\text{为常数}.

1) 根据给定运动方程写出特征方程。

2) 求弹簧摆的全部固有角频率,用 m,I,kd,ωnm,I,k_d,\omega_n 表示。

3) 观察伸缩自由振动时,发现振幅逐渐减小至接近0后又逐渐增大,增大到一定程度后又逐渐减小。将此现象视为拍振,求其周期 TT。

4) 当 kd2≪kykθk_d^2\ll k_yk_\theta 时,证明 T≫2πωnT\gg\dfrac{2\pi}{\omega_n}。

Kai​

(1)​

1), 2)​

my¨1+c1y˙1+(k1+k2)y1=0.\boxed{m\ddot y_1+c_1\dot y_1+(k_1+k_2)y_1=0}.

特性根と振動諸量は

s=−c12m±i4m(k1+k2)−c122m,\boxed{s=-\frac{c_1}{2m}\pm i\frac{\sqrt{4m(k_1+k_2)-c_1^2}}{2m}},
ωd=k1+k2m−c124m2,ζ=c12m(k1+k2).\boxed{\omega_d=\sqrt{\frac{k_1+k_2}{m}-\frac{c_1^2}{4m^2}}},\qquad \boxed{\zeta=\frac{c_1}{2\sqrt{m(k_1+k_2)}}}.

3)​

台車および質量のない接続点の力のつり合いから

my¨1+c1y˙1+(k1+k2)y1−k2y2=k1u1+c1u˙1,\boxed{m\ddot y_1+c_1\dot y_1+(k_1+k_2)y_1-k_2y_2=k_1u_1+c_1\dot u_1},
c2y˙2+k2y2=k2y1.\boxed{c_2\dot y_2+k_2y_2=k_2y_1}.

4)​

s=iωs=i\omega として後式から Y2=k2Y1/(k2+c2s)Y_2=k_2Y_1/(k_2+c_2s) を消去すると

Y1U1=(k1+c1s)(k2+c2s)(ms2+c1s+k1+k2)(k2+c2s)−k22.\frac{Y_1}{U_1}= \frac{(k_1+c_1s)(k_2+c_2s)}{(ms^2+c_1s+k_1+k_2)(k_2+c_2s)-k_2^2}.

従って定常振幅は

A(k12+ω2c12)(k22+ω2c22){k2(k1−mω2)−ω2c1c2}2+ω2{c2(k1+k2−mω2)+c1k2}2.\boxed{ A\frac{\sqrt{(k_1^2+\omega^2c_1^2)(k_2^2+\omega^2c_2^2)}}{ \sqrt{\begin{gathered}\{k_2(k_1-m\omega^2)-\omega^2c_1c_2\}^2\\ +\omega^2\{c_2(k_1+k_2-m\omega^2)+c_1k_2\}^2\end{gathered}}} }.

(2)​

1), 2)​

調和解を代入した係数行列の行列式をゼロにして

(ky−mω2)(kθ−Iω2)−kd2=0.\boxed{(k_y-m\omega^2)(k_\theta-I\omega^2)-k_d^2=0}.

δ=∣kd∣/mI\delta=|k_d|/\sqrt{mI} とおくと、大小順に

ω−=ωn2−δ,ω+=ωn2+δ.\boxed{\omega_-=\sqrt{\omega_n^2-\delta},\qquad \omega_+=\sqrt{\omega_n^2+\delta}}.

kd2<kykθk_d^2<k_yk_\theta によりいずれも実数である。kdk_d の符号はモードの位相関係を変えるが、固有振動数の集合は変えない。

3)​

両モードが励起される場合のうなり周期は

T=2πω+−ω−(kd≠0).\boxed{T=\frac{2\pi}{\omega_+-\omega_-}} \qquad(k_d\ne0).

例えば cos⁡ω−t+cos⁡ω+t=2cos⁡[(ω+−ω−)t/2]cos⁡[(ω++ω−)t/2]\cos\omega_-t+\cos\omega_+t=2\cos[(\omega_+-\omega_-)t/2]\cos[(\omega_++\omega_-)t/2] であり、振幅の絶対値が繰り返す周期が上式である。kd=0k_d=0 では両周波数が一致し、有限周期のうなりは生じない。

4)​

弱い連成では ω+−ω−≃δ/ωn\omega_+-\omega_-\simeq\delta/\omega_n なので

T≃2πωnmI∣kd∣,T2π/ωn≃kykθ∣kd∣≫1.\boxed{T\simeq\frac{2\pi\omega_n\sqrt{mI}}{|k_d|}},\qquad \boxed{\frac{T}{2\pi/\omega_n}\simeq\frac{\sqrt{k_yk_\theta}}{|k_d|}\gg1}.

Reference​