小さい字は補足説明や余談なので、読み飛ばしてもいいです。
第5回で、 M(z)=−EIv″(z)という式を導いた。 静定梁で梁の曲げモーメントがわかるなら、 この式を使えば、 それを2回積分してたわみが求まるので、これはなかなか便利な式だ。 では、不静定梁で、力のつりあいから曲げモーメントが求まらない場合は どうしたらいいだろう。 今回は、もう少し強力?な式を導いてみる。
任意の境界条件で任意の外力を受けてつりあっている梁の微小部分を切り取って、
この微小部分の力のつりあいを考えてみる。
図では単純梁にしているが、別に不静定梁でも構わない。
外力は、鉛直方向には、q(z)で表される分布外力が作用し、
軸方向には、p(z)で表されるz方向分布外力が作用しているものとする。
微小部分の長さはdzとし、
この微小範囲では鉛直方向の分布荷重外力q(z)や
軸方向の分布荷重外力p(z)は等分布していると見なせることにしておく。
切り取られたことによってできた両端の切断面には、
この切り取られた部分だけでつりあいが成り立つような
内力としての断面力が
図心(z軸)に作用していることにする。
この微小部分は外力q(z)とp(z)を受けているので、
左端の断面力と右端の断面力は大きさが同じ(で向きが逆)にはならない。
左端での軸力をN(z), せん断力をS(z), 曲げモーメントをM(z)とすると、
そこからdzだけずれた右端でのこれらの変化量は、
「瞬間の変化率」
dN(z)dz,dS(z)dz,dM(z)dz
にdzをかけて、
それぞれ、
dN(z)dzdz,dS(z)dzdz,dM(z)dzdz
と表せるから、右端の軸力、せん断力、曲げモーメントは、それぞれ、
N(z)+dN(z)dzdz,S(z)+dS(z)dzdz,M(z)+dM(z)dzdz
と表せる。
(参考:
(応力のつりあい)
これらを用いてこの微小部分の力のつりあいを考えると
∑↓=−S(z)+q(z)dz+(S(z)+dS(z)dzdz)=0
∑→=−N(z)+p(z)dz+(N(z)+dN(z)dzdz)=0
\sum_{A}\circlearrowleft=-M(z)-q(z)dz\frac{dz}{2}-(S(z)+\frac{dS(z)}{dz}dz)dz
+(M(z)+\frac{dM(z)}{dz}dz)
=0
となる。
分布荷重q(z)は長方形分布とみなしているから、
荷重の合計はq(z)dzを長方形の重心に作用させると、
点Aからのモーメントの腕の長さは\frac{dz}{2}。
3つめのモーメントのつりあいの式に関して、
q(z)dz\frac{dz}{2}や(\frac{dS(z)}{dz}dz)dz
はdのつく微小量の2次項なので微小と見なして無視して整理する。
もう少し厳密に言うと、2次項は、
分子のdの数ひく分母のdの数が2つ以上の項ということ。
すると、
-S(z)dz+\frac{dM(z)}{dz}dz=0
となるから、両辺をdzで割って整理すると、以下の関係が得られる。
S(z)=
\frac{dM(z)}{dz}
(重要)
これは、構造力学Iで
モーメント図を描くときにさんざん使ってきた S(z)=M'(z)の式である。
つまり、梁モデルのせん断力は曲げモーメントの微分で与えられるということになる。
また、
y方向のつりあい式を同様に整理すると:
\frac{dS(z)}{dz}+q(z)=0
となるが、これに上の
S(z)= \frac{dM(z)}{dz}
を代入すると、曲げに関するつりあい式は、以下のようになる。
\frac{d^{2}M(z)}{dz^{2}}+q(z)=0
となる。微分を'で表すなら、以下のように書いてもよい。
M''(z)+q(z)=0
また、z方向のつりあい式を整理すると、
軸力のつりあい式は以下のようになる。
\frac{dN(z)}{dz}+p(z)=0
微分を'で書けば、以下のように書いてもよい。
N'(z)+p(z)=0
梁モデルの仮定では、
初等梁にはせん断変形がないことになっていた。
だから、初等梁でゼロでないひずみは、
軸方向の直ひずみ\varepsilon_{zz}だけで、
ゼロでない応力は
軸方向の直応力\sigma_{zz}だけであった。
だから、この軸方向の直応力\sigma_{zz}を面積積分して合計することによって
軸力Nを定義したり、
この軸方向の直応力\sigma_{zz}に「腕の長さ」をかけながら
面積積分して合計する
ことによって
曲げモーメントM
を定義することはできたが、
断面に作用するせん断応力\sigma_{yz}はない(無視できる)ことになっていたので、
せん断力Sを面積積分で定義することはできなかった。
それなのに梁の微小部分の分布外力と断面力の力のつりあいを考えると、
(断面に作用するせん断応力\sigma_{yz}は無視されているにもかかわらず)
せん断力が、曲げモーメントの微分として
S(z)=\frac{dM(z)}{dz}
と与えられてしまう。これは
梁モデルの仮定
(でやや強引な仮定をしたことによる)
初等梁理論のちょっと変な点である。
そもそも、断面のせん断応力\sigma_{yz}が無視できるからといって、
断面のせん断力Sもないことにしてしまうと、
力のつりあいは曲げモーメントが一定の場合しかなりたたなくなってしまう。
実はベルヌーイ・オイラーの仮定が
(それほど細長くない梁にも一般的に)成り立つ変形は、
梁が円弧のように変形する場合で、これは曲げモーメントが一定
の状態で「純曲げ」と呼ばれる状態である。
しかし、初等梁理論では純曲げ状態でなくても
ベルヌーイ・オイラーの仮定が
成り立つことにしてしまったために、
「断面のせん断応力\sigma_{yz}は無視できても、
断面力としてのせん断力Sはあることにしないとモーメントのつりあいが
満たされない」
という矛盾点を抱える理論なので、その点に注意しておく必要がある。
例えば、図のように微小部分の左端と右端とで曲げモーメントが異なる場合、
A点まわりのモーメントのつりあいは、\frac{dM}{dz}dz-Sdz=0となり、
つまり
S(z)=\frac{dM(z)}{dz}
となる訳だが、S=0にしてしまうと、
モーメントの変化率\frac{dM(z)}{dz}も0になってしまうので、
純曲げしか表せなくなってしまう。
つまり、梁のたわみは円弧でしか表せなくなる。
いくら細長くてせん断変形が無視できる梁だからといって、
断面力は純曲げ状態にしかならないというモデルでは、
精度が悪すぎて実用的ではないだろう。
そこで、(せん断応力\sigma_{yz}は無視できたとしても)
モーメントのつりあいを満たすために必要なせん断力Sを
モーメントの変化率として
S(z)=\frac{dM(z)}{dz}
と定義してやることで、細長い梁に対する実用的な
予測が行えるモデルとなっているのである。
さて、初等梁の式を導く際には、せん断応力\sigma_{yz}を無視しているが、
実際には、細長い梁でもせん断応力はそれなりにあるのである。
初等梁の
直応力\sigma_{zz}は、三角形分布しているが、
この直応力とつりあうようなせん断応力を求めると、
せん断応力\sigma_{yz}は、図のような放物線分布をしていて、
y=0の位置で最大となる。
ちなみに、長方形断面の場合、断面の最大せん断応力は、
\sigma_{yz}^{max}=\frac{3S}{2A}となる。
こうしたせん断応力に関する詳しい解説をもっと読みたい人は、
このテキストのネタ本でもある
岩熊哲夫・小山茂『鬆徒労苦衷有迷禍荷苦痛--
計算機による構造解析の基礎としての構造力学を独習する』
を参照してほしい。
さて、上で導いた
軸力に関するつりあい式
N'(z)+p(z)=0
と
曲げに関するつりあい式
M''(z)+q(z)=0
に、
軸力-図心変位関係
N(z)=EAw'(z)
と
曲げモーメント-たわみ関係
M(z)=
-EIv''(z)
をそれぞれ代入すると、以下のような
のような図心変位v(たわみ)とw(軸方向変位)の微分方程式が導かれる。
EAw''(z)+p(z)=0
-EIv''''(z)+q(z)=0
(最重要)
これらを梁の支配(微分)方程式とか、
特にたわみに関する式をたわみの微分方程式とか言う。
上の式でzの関数になっているのは、
変位(wやv)の微分と分布外力(pやq)だけだから、
外力(pやq)の式がわかっていれば、
モーメントの式M(z)がわからなくても、
これらの式を積分して、
静定梁のたわみをM=-EIv''を積分して解いた
のと同じ要領で境界条件と連続条件
(と後述するが、集中荷重の場合は載荷点の微小部分の外力と断面力のつりあい条件)を
使って積分定数を決定すれば、
不静定梁でもたわみや軸方向変位を求めることができる。
これについては、次回やるので楽しみにしておいてほしい。
梁の微分方程式を積分して、 たわみや断面力などを求める際、 たわみvの微分(v', v'', v'''など)がたわみ角や断面力と どのように関係付けられているかということを知っていないと、 境界条件や連続条件などを適切に利用することができない。 そこで、以下にたわみの微分とたわみ角や断面力との関係式をまとめておく。
図の不静定梁のたわみが
v(z)=\frac{q}{48EI}(2z^{4}-5\ell z^{3}+3\ell^{2}z^{2})
で与えられるとき、
せん断力S(z), 曲げモーメントM(z),
反力V_\text{A}, M_\text{A}, V_\text{B}を求めよ。
答えはここ。
図の不静定梁のたわみが
v(z)=\frac{q}{24EI}(z^{4}-2\ell z^{3}+\ell^{2}z^{2})
で与えられるとき、
せん断力S(z), 曲げモーメントM(z)を求めよ。
答えはここ。
梁の鉛直方向のつりあい式が
S'(z)+q(z)=0
で与えられるとき、せん断力と曲げモーメントの関係
S(z)=M'(z)と
曲げモーメントとたわみの関係M(z)=-EIv''(z)を用いて、
たわみの微分方程式を導け。
答えはここ。
図の不静定梁のたわみが
v(z)=\frac{q}{48EI}(2z^{4}-3\ell z^{3}+\ell^{3}z)
で与えられるとき、
せん断力S(z), 曲げモーメントM(z),
反力V_\text{A}, V_\text{B}, M_\text{B}を求めよ。
答えはここ。
求め方は、例題1の求め方と同様。
2022年度小テスト:
小テスト221214
2021年度小テスト:
問
(問題文の\frac{(3)}{(4)}, \frac{(5)}{(6)}とすべきところが、
すべて\frac{(1)}{(2)}になってました。
すいません。適宜、答案を見て判断します)
解答
2020年度小テスト:
問1-4,
解答
メモ: