⑧3次元弾性体の応力とひずみとは？「youtube 【初めて学ぶ人のための材料力学の教室】演習問題1/2」をsympyでやってみたい。

・演習問題1

youtubeです。

途中4:38から (0:00〜10:57)

＜再生リスト
https://www.youtube.com/@zairiki/playlists

sympyで

・ver0.1

# ver0.1
from sympy import *
from decimal import Decimal, ROUND_HALF_UP, ROUND_HALF_EVEN
var('E,ν' ,real=True)
var('σx,σy,σz'    ,real=True)
var('εx,εy,εz,εv',real=True)
var('τxy,τyz,τyz' ,real=True)
var('α ,ΔT'        ,real=True)
rep={E    :200,ν :0.3,        \
     σx  :1  ,σy:1  ,σz :0, \
     #τxy:0  ,τyz:1,τzx:0,  \
     τxy:0  ,τyz:1  ,        \
     α   :0  ,ΔT:0           \
     }
εx=    σx/E-ν*σy/E-ν*σz/E+α*ΔT  ;print("# εx=",εx)
εz=-ν*σx/E-ν*σy/E+   σz/E+α*ΔT  ;print("# εz=",εz)
# print("# εx =",εx.subs(rep)*10**3)     
# print("# εz =",εz.subs(rep)*10**3)     
# print("# εx =",round(εx.subs(rep)*10**3,0))     
# print("# εz =",round(εz.subs(rep)*10**3,0))     
# 
print("# εx =",Decimal(str(εx.subs(rep)*10**3)).quantize(Decimal('0.1'), ROUND_HALF_UP))
print("# εz =",Decimal(str(εz.subs(rep)*10**3)).quantize(Decimal('0.1'), ROUND_HALF_UP))
# εx= ΔT*α - ν*σy/E - ν*σz/E + σx/E
# εz= ΔT*α - ν*σx/E - ν*σy/E + σz/E
# εx = 3.5
# εz = -3.0

いつもの? sympyの実行環境と 参考のおすすめです。

(テンプレート)

いつもと違うおすすめです。

小数を任意の桁数で四捨五入　＜　

参考文献

＞等分布荷重を受ける単純支持はり
＞JSME p85
　
＞単純支持ばり 7)
＞構造力学公式集 p136