0
0

Last updated at Posted at 2024-01-27

・演習問題1

# 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
``````

(テンプレート)

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

#### 参考文献

＞等分布荷重を受ける単純支持はり
＞JSME p85

＞単純支持ばり 7)
＞構造力学公式集 p136

0
0
0

Register as a new user and use Qiita more conveniently

1. You get articles that match your needs
2. You can efficiently read back useful information
3. You can use dark theme
What you can do with signing up
0
0