【構造力学】たわみ曲線の微分方程式、厳密解調べてみた！

今回は私が高専で学習している構造力学の微分方程式について紹介します．

証明は参考文献にございます．
多くの書籍では，供試体の傾き$\frac{\mathrm{d}y}{\mathrm{d}x}$を$0$とみなして簡単な形で微分方程式を解いていますが，これを厳密に解きます．

ちなみに鋼などの弾性係数が大きいものだと近似を使ってもとても良い精度のたわみ量が算出できます．

曲げモーメントが$x$の関数のとき

$\frac{\mathrm{d}y}{\mathrm{d}x}=p$とします．
$$\begin{array}{rcl}\frac{\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{x}^2}}{(1+(\frac{\mathrm{d}y}{\mathrm{d}x}{)}^2{)}^{3/2}}& =& -\frac{M(x)}{EI}\\ \frac{\frac{\mathrm{d}p}{\mathrm{d}x}}{(1+p^2{)}^{3/2}}& =& -\frac{M(x)}{EI}\\ \int \frac{1}{(1+p^2{)}^{3/2}}\frac{\mathrm{d}p}{\mathrm{d}x}\mathrm{d}x& =& -\int \frac{M(x)}{EI}\mathrm{d}x\\ \frac{p}{\sqrt{1+p^2}}& =& -\int \frac{M(x)}{EI}\mathrm{d}x\\ \text{(1)}& p=\frac{\mathrm{d}y}{\mathrm{d}x}& =& -\frac{\int \frac{M(x)}{EI}\mathrm{d}x}{\sqrt{1-{(\int \frac{M(x)}{EI}\mathrm{d}x)}^2}}y& =& -\int \frac{\int \frac{M(x)}{EI}\mathrm{d}x}{\sqrt{1-{(\int \frac{M(x)}{EI}\mathrm{d}x)}^2}}\mathrm{d}x\end{array}$$
境界条件などを考え積分定数を決定して解ける．

具体例を考えて解いてみます．

曲げモーメントが$y$の関数のとき

先ほどと同様に，$\frac{\mathrm{d}y}{\mathrm{d}x}=p$とします．
$$\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{x}^2}=\frac{\mathrm{d}p}{\mathrm{d}x}=\frac{\mathrm{d}p}{\mathrm{d}y}\frac{\mathrm{d}y}{\mathrm{d}x}=p\frac{\mathrm{d}p}{\mathrm{d}y}$$を用いて，
$$\begin{array}{rcl}\frac{p\frac{\mathrm{d}p}{\mathrm{d}y}}{(1+p^2{)}^{3/2}}& =& -\frac{M(y)}{EI}\\ \int \frac{p}{(1+p^2{)}^{3/2}}\frac{\mathrm{d}p}{\mathrm{d}y}\mathrm{d}y& =& -\int \frac{M(y)}{EI}\mathrm{d}y\\ \frac{1}{\sqrt{1+p^2}}& =& \int \frac{M(y)}{EI}\mathrm{d}y\\ \text{(2)}& p=\frac{\mathrm{d}y}{\mathrm{d}x}& =& \sqrt{\frac{1}{{(\int \frac{M(y)}{EI}\mathrm{d}y)}^2}-1}\int \frac{\int \frac{M(y)}{EI}\mathrm{d}y}{\sqrt{1-{(\int \frac{M(y)}{EI}\mathrm{d}y)}^2}}\frac{\mathrm{d}y}{\mathrm{d}x}\mathrm{d}x& =& \int \mathrm{d}x\\ x& =& \int \frac{\int \frac{M(y)}{EI}\mathrm{d}y}{\sqrt{1-{(\int \frac{M(y)}{EI}\mathrm{d}y)}^2}}\mathrm{d}y\end{array}$$

座屈$$M(y)=Py$$とかで使えます．