本稿は ティモシェンコ梁でも モールの定理が使えることを紹介します。
本稿では 弾性曲線式を以下の解釈で扱っています。
今回は 1端ピン他端固定 で ティモシェンコ梁を再考します。
(8)~(11)式を 数式処理ソフト maxima で 解きます。
単純ばりで計算して、固定端モーメント分の曲げモーメントの荷重を補正して求めます。
つぎのプログラムを 400-7-2.bat で保存して実行すると
: 1端ピン他端固定のモールの定理(再考)
@echo off
path C:\maxima-5.47.0\bin;%path%
ruby -x %~f0
pause
goto:eof
#!ruby
f = open "| maxima --very-quiet", "w"
f.print <<~Maxima
C: w * L^2 / (8*(3*__g__+1))$
m : w * x * (L - x) / 2$
Q: C1 - integrate(m, x)$
M: C2 + integrate(Q, x)$
[Q, M]: ev([Q, M], solve([ev(M, x = 0) = 0, ev(M, x = L) = 0], [C1, C2]))$
Q: Q - C * (L^2 - 3 * x^2) / (6 * L)$
M: M - C * x * (L - x) * (L + x) / (6 * L)$
M: M + (m - ev(m, x=0)) * κ * EI / GA$
print("T =", string(factor(ratsubst(__g__, κ * EI / GA / L^2, Q))))$
print("V =", string(factor(ratsubst(__g__, κ * EI / GA / L^2, M))))$
quit()$
Maxima
f.close
__END__
計算結果は
となります。__g__ は λ です。固定端モーメント C はせん断変形の影響を考慮したものを使う必要があります。
曲げモーメントを荷重として 解きます。
たわみが 数式処理ソフト maxima の結果と微妙に違っています。中央部は一致するのですが 原因がわかりません。maxima の結果は 弾性曲線式を解いた結果と一致するので 正と考えています。
1端ピン他端固定でも ティモシェンコ梁でモールの定理が利用できました。
固定端モーメント C は せん断変形を考慮したものを使う必要があることに、前回は気づいていませんでした。
ティモシェンコ梁でも モールの定理は使えます。支持条件は片持ちや単純支持にかぎらず、両端固定でも1端ピン他端固定でも問題ありませでした。
あとがき の訂正と補足としてください。
ほっとしました。ひとには気を付けてと言いながら、自分が穴にはまりました。それでも、よかったよかった。大阪わんぱく博覧会。木を植えたらよかったのに。