東京大学 工学系研究科 機械工学専攻 機械工学第2部 2025年8月実施 問題 1
Author
Andrea
Description
図 1-1 に示すような両端において単純支持された,高さ h h h の矩形横断面を有する長さ 2 L 2L 2 L のはり AB が,中央の点 D において長さ L L L ,断面積 S S S の棒 CD にピン結合されている構造を考える.また,棒 CD は点 C において水平な床面とピン結合されている.はり AB は水平方向,棒 CD は鉛直方向を向いている.はり AB は降伏応力 σ Y \sigma_Y σ Y の弾塑性体,棒 CD は線膨張係数 α \alpha α の線形弾性体である.ヤング率,断面二次モーメントは,はり AB と棒 CD で共通で,それぞれ E E E , I I I である.初期状態では,はり AB と棒 CD には応力は加わっていない.はり AB,棒 CD の自重は無視してよい.
いま,棒 CD のみ温度が Δ T ( > 0 ) \Delta T (>0) Δ T ( > 0 ) 上昇した.以下の設問に答えよ.
(1) 点 D において,はり AB と棒 CD 間が互いに及ぼしあっている力を R D R_D R D とする.このとき,はり AB の変形は,図 1-2 に示すようにはりの中央に集中荷重 R D R_D R D が加わる単純支持はりの変形として考えることができる.点 D におけるはり AB のたわみ δ A B \delta_{AB} δ A B を R D , L , E , I R_D, L, E, I R D , L , E , I を用いて表せ.
(2) 棒 CD の伸び δ C D \delta_{CD} δ C D を R D , L , E , S , α , Δ T R_D, L, E, S, \alpha, \Delta T R D , L , E , S , α , Δ T を用いて表せ.
(3) R D R_D R D を L , E , I , S , α , Δ T L, E, I, S, \alpha, \Delta T L , E , I , S , α , Δ T を用いて表せ.
(4) Δ T \Delta T Δ T が Δ T Y \Delta T_Y Δ T Y に達すると,はり AB は降伏を始める.Δ T Y \Delta T_Y Δ T Y を求めよ.ただし,本設問(4)において棒 CD は座屈しないものとする. 5. 棒 CD が座屈する最低の温度上昇量を Δ T c r \Delta T_{cr} Δ T cr とする.Δ T c r < Δ T Y \Delta T_{cr} < \Delta T_Y Δ T cr < Δ T Y となるための条件を,L , h , E , σ Y L, h, E, \sigma_Y L , h , E , σ Y を用いて表せ. ただし,両端が回転軸に支持されている長さ L L L ,ヤング率 E E E ,断面二次モーメント I I I の棒の臨界座屈荷重 P c r P_{cr} P cr は
P c r = π 2 L 2 E I P_{cr} = \frac{\pi^2}{L^2}EI P cr = L 2 π 2 E I
で与えられるものとする.
题目描述
长度 2 L 2L 2 L 、矩形截面高度 h h h 的简支梁 AB,在跨中 D 与竖直杆 CD 铰接;杆长 L L L 、截面积 S S S ,下端 C 与地面铰接。梁为屈服应力 σ Y \sigma_Y σ Y 的弹塑性体,杆为线膨胀系数 α \alpha α 的线弹性体;两者杨氏模量均为 E E E 、截面二次矩均为 I I I 。初始无应力,忽略自重,仅使杆 CD 升温 Δ T > 0 \Delta T>0 Δ T > 0 。
设 D 点梁、杆之间的相互作用力为 R D R_D R D 。把 AB 看成跨中受集中力的简支梁,用 R D , L , E , I R_D,L,E,I R D , L , E , I 表示 D 点挠度 δ A B \delta_{AB} δ A B 。
综合自由热膨胀与轴向弹性压缩,用 R D , L , E , S , α , Δ T R_D,L,E,S,\alpha,\Delta T R D , L , E , S , α , Δ T 表示杆的净伸长 δ C D \delta_{CD} δ C D 。
用 D 点位移协调条件求 R D R_D R D 关于结构参数和温升的表达式。
假定杆不座屈,求梁最外缘首次达到 σ Y \sigma_Y σ Y 时的温升 Δ T Y \Delta T_Y Δ T Y 。
杆两端均为铰支,其 Euler 临界力为 P c r = π 2 E I / L 2 P_{\rm cr}=\pi^2EI/L^2 P cr = π 2 E I / L 2 。求杆在梁屈服前先发生座屈,即 Δ T c r < Δ T Y \Delta T_{\rm cr}<\Delta T_Y Δ T cr < Δ T Y ,所需的 L , h , E , σ Y L,h,E,\sigma_Y L , h , E , σ Y 条件。
梁挠度:积分 Euler–Bernoulli 挠曲线方程,并用支座位移和跨中对称条件求集中荷载下的中央挠度。
热变形协调:把杆的热伸长与受压力产生的弹性缩短叠加,再令连接点的杆伸长等于梁挠度。
屈服:由跨中最大弯矩和 σ max = M c / I \sigma_{\max}=Mc/I σ m a x = M c / I 求梁开始屈服对应的相互作用力及温升。
Euler 屈曲:比较铰支杆临界压力与梁屈服载荷,得到哪种失效模式先发生的判据。
Kai
(1) 梁 AB 的挠度推导
设梁 AB 左端为原点 x = 0 x=0 x = 0 ,梁长 2 L 2L 2 L 。
由于结构对称,支座反力为 R A = R D / 2 R_A = R_D/2 R A = R D /2 (方向向下)。
任意位置 x x x (0 ≤ x ≤ L 0 \le x \le L 0 ≤ x ≤ L ) 的弯矩为 M ( x ) = − R D 2 x M(x) = - \frac{R_D}{2}x M ( x ) = − 2 R D x 。
根据挠曲线微分方程 E I y ′ ′ = − M ( x ) EI y'' = -M(x) E I y ′′ = − M ( x ) :
E I d 2 y d x 2 = R D 2 x EI \frac{d^2y}{dx^2} = \frac{R_D}{2}x E I d x 2 d 2 y = 2 R D x
积分两次并利用以下边界条件确定积分常数:
对称性 :在梁中点 (x = L x=L x = L ),转角为 0 ⇒ y ′ ( L ) = 0 \Rightarrow y'(L) = 0 ⇒ y ′ ( L ) = 0
支持条件 :在支座处 (x = 0 x=0 x = 0 ),位移为 0 ⇒ y ( 0 ) = 0 \Rightarrow y(0) = 0 ⇒ y ( 0 ) = 0
解得挠度方程:
E I y ( x ) = R D 12 x 3 − R D L 2 4 x EI y(x) = \frac{R_D}{12}x^3 - \frac{R_D L^2}{4}x E I y ( x ) = 12 R D x 3 − 4 R D L 2 x
当 x = L x=L x = L 时,即为 D 点挠度 δ A B \delta_{AB} δ A B (取绝对值):
δ A B = ∣ y ( L ) ∣ = ∣ R D L 3 12 E I − R D L 3 4 E I ∣ = R D L 3 6 E I \delta_{AB} = |y(L)| = \left| \frac{R_D L^3}{12EI} - \frac{R_D L^3}{4EI} \right| = \frac{R_D L^3}{6EI} δ A B = ∣ y ( L ) ∣ = 12 E I R D L 3 − 4 E I R D L 3 = 6 E I R D L 3
(2) 棒 CD 的伸长量 δ C D \delta_{CD} δ C D
应用叠加原理 ,棒 CD 的总变形由以下两部分组成:
热膨胀 (向上伸长):
由温度升高 Δ T \Delta T Δ T 引起:δ T = α L Δ T \delta_T = \alpha L \Delta T δ T = αL Δ T
弹性压缩 (向下收缩):
由梁的反力 R D R_D R D (压力)引起:δ F = R D L E S \delta_F = \frac{R_D L}{ES} δ F = ES R D L
因此,棒 CD 的总伸长量为:
δ C D = δ T − δ F = α L Δ T − R D L E S \delta_{CD} = \delta_T - \delta_F = \alpha L \Delta T - \frac{R_D L}{ES} δ C D = δ T − δ F = αL Δ T − ES R D L
(3) 求相互作用力 R D R_D R D
利用几何协调条件 :在连接点 D,梁的向上位移必须等于棒的向上伸长量。
δ A B = δ C D \delta_{AB} = \delta_{CD} δ A B = δ C D
代入 (1) 和 (2) 的结果:
R D L 3 6 E I = α L Δ T − R D L E S \frac{R_D L^3}{6EI} = \alpha L \Delta T - \frac{R_D L}{ES} 6 E I R D L 3 = αL Δ T − ES R D L
我们需要解出 R D R_D R D 。首先消去等式两边的 L L L ,并将含有 R D R_D R D 的项移到左边:
R D L 2 6 E I + R D E S = α Δ T \frac{R_D L^2}{6EI} + \frac{R_D}{ES} = \alpha \Delta T 6 E I R D L 2 + ES R D = α Δ T
提取公因式 R D R_D R D :
R D ( L 2 6 E I + 1 E S ) = α Δ T R_D \left( \frac{L^2}{6EI} + \frac{1}{ES} \right) = \alpha \Delta T R D ( 6 E I L 2 + ES 1 ) = α Δ T
通分括号内的项(公分母为 6 E I S 6EIS 6 E I S ):
R D ( S L 2 + 6 I 6 E I S ) = α Δ T R_D \left( \frac{S L^2 + 6I}{6EIS} \right) = \alpha \Delta T R D ( 6 E I S S L 2 + 6 I ) = α Δ T
解得最终结果:
R D = 6 E I S α Δ T S L 2 + 6 I R_D = \frac{6EIS \alpha \Delta T}{S L^2 + 6I} R D = S L 2 + 6 I 6 E I S α Δ T
(4) 求梁 AB 开始降伏时的温度 Δ T Y \Delta T_Y Δ T Y
梁 AB 发生降伏的条件是最大弯曲正应力达到降伏应力 σ Y \sigma_Y σ Y 。
对于跨度为 2 L 2L 2 L 、中点受力 R D R_D R D 的简支梁,最大弯矩 M m a x M_{max} M ma x 发生在跨中:
M m a x = R D ⋅ ( 2 L ) 4 = R D L 2 M_{max} = \frac{R_D \cdot (2L)}{4} = \frac{R_D L}{2} M ma x = 4 R D ⋅ ( 2 L ) = 2 R D L
最大弯曲应力 σ m a x \sigma_{max} σ ma x 发生在截面最外缘(距离中性轴 y = h / 2 y = h/2 y = h /2 处):
σ m a x = M m a x ⋅ h 2 I = R D L 2 ⋅ h 2 I = R D L h 4 I \sigma_{max} = \frac{M_{max} \cdot \frac{h}{2}}{I} = \frac{\frac{R_D L}{2} \cdot \frac{h}{2}}{I} = \frac{R_D L h}{4I} σ ma x = I M ma x ⋅ 2 h = I 2 R D L ⋅ 2 h = 4 I R D L h
令 σ m a x = σ Y \sigma_{max} = \sigma_Y σ ma x = σ Y ,解出对应的临界力 R D , Y R_{D,Y} R D , Y :
R D , Y L h 4 I = σ Y ⟹ R D , Y = 4 I σ Y L h \frac{R_{D,Y} L h}{4I} = \sigma_Y \implies R_{D,Y} = \frac{4I \sigma_Y}{L h} 4 I R D , Y L h = σ Y ⟹ R D , Y = L h 4 I σ Y
将 R D , Y R_{D,Y} R D , Y 代入第 (3) 问求得的 R D R_D R D 与 Δ T \Delta T Δ T 的关系式中:
6 E I S α Δ T Y S L 2 + 6 I = 4 I σ Y L h \frac{6EIS \alpha \Delta T_Y}{S L^2 + 6I} = \frac{4I \sigma_Y}{L h} S L 2 + 6 I 6 E I S α Δ T Y = L h 4 I σ Y
消去两边的 I I I ,解出 Δ T Y \Delta T_Y Δ T Y :
Δ T Y = 4 σ Y ( S L 2 + 6 I ) 6 E S α L h = 2 σ Y ( S L 2 + 6 I ) 3 E S α L h \Delta T_Y = \frac{4 \sigma_Y (S L^2 + 6I)}{6 E S \alpha L h} = \frac{2 \sigma_Y (S L^2 + 6I)}{3 E S \alpha L h} Δ T Y = 6 ES αL h 4 σ Y ( S L 2 + 6 I ) = 3 ES αL h 2 σ Y ( S L 2 + 6 I )
(5) 棒 CD 发生座屈的条件
题目给定棒 CD 的临界座屈载荷为:
P c r = π 2 E I L 2 P_{cr} = \frac{\pi^2 EI}{L^2} P cr = L 2 π 2 E I
我们要求的是“棒 CD 在梁 AB 降伏之前发生座屈”(即 Δ T c r < Δ T Y \Delta T_{cr} < \Delta T_Y Δ T cr < Δ T Y )。
由于系统是线性的,内力 R D R_D R D 与温升 Δ T \Delta T Δ T 成正比。因此,这一条件等价于:使棒座屈所需的力,必须小于使梁降伏所需的力 。
P c r < R D , Y P_{cr} < R_{D,Y} P cr < R D , Y
代入前面求出的表达式(注意 R D , Y = 4 I σ Y L h R_{D,Y} = \frac{4I \sigma_Y}{L h} R D , Y = L h 4 I σ Y ):
π 2 E I L 2 < 4 I σ Y L h \frac{\pi^2 EI}{L^2} < \frac{4I \sigma_Y}{L h} L 2 π 2 E I < L h 4 I σ Y
化简步骤 :
消去两边的 I I I (材料性质相同,截面二次矩相同)。
两边同乘 L L L (消去分母中的一个 L L L )。
π 2 E L < 4 σ Y h \frac{\pi^2 E}{L} < \frac{4 \sigma_Y}{h} L π 2 E < h 4 σ Y
整理得到最终条件:
σ Y > π 2 E h 4 L \sigma_Y > \frac{\pi^2 E h}{4 L} σ Y > 4 L π 2 E h