SolidMechanicsStress

✖
`SolidMechanicsStress`

[実験的]

✖

SolidMechanicsStress[vars,pars,strain]

固体力学内部応力を，変数 vars，パラメータ pars，全歪み strain で与える．

✖

SolidMechanicsStress[vars,pars,strain,displacement]

非線形材料規則に対する固体力学応力を与える．

詳細

SolidMechanicsStressは，与えられた全歪みに対する力学応力を，変位の従属変数，，（単位：），独立変数（単位：），時間変数（単位：）で返す．

応力は断面積あたりの抵抗力である．
応力の単位はである．
SolidMechanicsStressはSolidMechanicsPDEComponentと同じ変数 vars 指定を使う．
SolidMechanicsStressはSolidMechanicsPDEComponentと同じパラメータpars s指定を使う．
SymmetrizedArrayは strain を指定する．
通常，strain はSolidMechanicsStrainの結果である．
弾性材料については，SolidMechanicsStressは弾性歪みを使って応力を計算する．

SolidMechanicsStressは熱歪みまたは初期歪みのような非弾性歪みを与えられた全歪みから引く．
SolidMechanicsStressは応力のSymmetrizedArrayを以下の形式で返す．

は垂直応力を，はせん断応力を表す．
SolidMechanicsStressは，内部応力を計算する．

初期応力は，SolidMechanicsPDEComponentについて"InitialStress"パラメータで指定できる．
SolidMechanicsStressは，非線形の材料については，全歪みと変位を使って応力を計算する．

例題

すべて開くすべて閉じる

例 (1)基本的な使用例

歪みによる応力を計算する：

In[1]:=1

✖

https://wolfram.com/xid/0tp02rczjhk79i-x185m3

Out[1]=1

応力を可視化する：

In[2]:=2

✖

https://wolfram.com/xid/0tp02rczjhk79i-4jx412

スコープ (4)標準的な使用例のスコープの概要

定常解析 (1)

端が固定され先端に力が加わったスプーンの変位を計算する．変数とパラメータを設定する：

In[1]:=1

✖

https://wolfram.com/xid/0tp02rczjhk79i-me42fu

偏微分方程式と幾何学配置を設定する：

In[2]:=2

✖

https://wolfram.com/xid/0tp02rczjhk79i-52jbt6

変位の応力を計算する：

In[3]:=3

✖

https://wolfram.com/xid/0tp02rczjhk79i-074nfv

Out[3]=3

歪みの応力を計算する：

In[4]:=4

✖

https://wolfram.com/xid/0tp02rczjhk79i-0qzv5l

Out[4]=4

応力を可視化する：

In[5]:=5

✖

https://wolfram.com/xid/0tp02rczjhk79i-5pkawf

Out[5]=5

定常平面応力解析 (1)

平面応力のケースを拡張モデルとして計算する．こうすると，面外歪みが計算でき，面外応力が0であることが検証できる．左端が肯定され右側が強制的に変位された矩形の鋼板がある．領域，変数，パラメータを設定する．変数には3方向すべてが含まれている：

In[1]:=1

✖

https://wolfram.com/xid/0tp02rczjhk79i-n3hln

3つの準従属変数で方程式を解く：

In[4]:=4

✖

https://wolfram.com/xid/0tp02rczjhk79i-9scgt4

Out[4]=4

変位のリストには3つの出力変数が含まれるようになった．主変数についての変位を可視化する：

In[5]:=5

✖

https://wolfram.com/xid/0tp02rczjhk79i-jbpbkr

Out[5]=5

歪みを計算する：

In[6]:=6

✖

https://wolfram.com/xid/0tp02rczjhk79i-72dans

Out[6]=6

歪みは3×3配列である．面外歪みを可視化する：

In[7]:=7

✖

https://wolfram.com/xid/0tp02rczjhk79i-cuaa5q

Out[7]=7

歪みからの応力を計算する：

In[8]:=8

✖

https://wolfram.com/xid/0tp02rczjhk79i-uxc2y4

Out[8]=8

応力は3×3配列である．平面応力条件を検証する：

In[9]:=9

✖

https://wolfram.com/xid/0tp02rczjhk79i-45ogn

Out[9]=9

定常平面歪み解析 (1)

平面歪みのケースを拡張モデルとして計算する．こうすることで，面外応力が計算でき，面外歪みが0であることが検証できる．領域，変数，パラメータを設定する．変数には3方向すべてが含まれる：

In[1]:=1

✖

https://wolfram.com/xid/0tp02rczjhk79i-odgagf

固体力学PDE成分を設定する：

In[4]:=4

✖

https://wolfram.com/xid/0tp02rczjhk79i-lmscmx

Out[4]=4

PDEを設定する：

In[5]:=5

✖

https://wolfram.com/xid/0tp02rczjhk79i-rz669l

PDEを解く：

In[7]:=7

✖

https://wolfram.com/xid/0tp02rczjhk79i-molvgw

Out[7]=7

変位による歪みを計算する：

In[8]:=8

✖

https://wolfram.com/xid/0tp02rczjhk79i-8pbvhd

Out[8]=8

歪みは3×3配列である．平面歪み条件を検証する：

In[9]:=9

✖

https://wolfram.com/xid/0tp02rczjhk79i-027b51

Out[9]=9

歪みからの応力を計算する：

In[10]:=10

✖

https://wolfram.com/xid/0tp02rczjhk79i-il12ny

Out[10]=10

応力は3×3配列である．面外応力を可視化する：

In[11]:=11

✖

https://wolfram.com/xid/0tp02rczjhk79i-181jjq

Out[11]=11

定常超弾性平面応力分析 (1)

左側が固定され右側に力が加えられた矩形のゴム板の変位を計算する．領域，変数，パラメータを設定する：

In[1]:=1

✖

https://wolfram.com/xid/0tp02rczjhk79i-dug6s7

方程式を解く：

In[2]:=2

✖

https://wolfram.com/xid/0tp02rczjhk79i-1sk6i6

変位を可視化する：

In[3]:=3

✖

https://wolfram.com/xid/0tp02rczjhk79i-mx96p7

Out[3]=3

変位による歪みを計算する：

In[4]:=4

✖

https://wolfram.com/xid/0tp02rczjhk79i-fn57y7

Out[4]=4

歪みと変位による応力を計算する：

In[5]:=5

✖

https://wolfram.com/xid/0tp02rczjhk79i-uuyxu1

Out[5]=5

歪みを可視化する：

In[6]:=6

✖

https://wolfram.com/xid/0tp02rczjhk79i-pj66wz

Out[6]=6

考えられる問題 (1)よく起る問題と予期しない動作

軸対称の場合は，転移が空間座標の関数であることが重要である．その理由を説明するために，変数，パラメータ，空間の関数ではない転移関数を作成する：

In[1]:=1

✖

https://wolfram.com/xid/0tp02rczjhk79i-i804qn

Out[3]=3

生成されたInterpolatingFunctionは空間座標およびの関数ではない点に注意のこと．この状態で歪みと応力を計算すると，応力テンソルにやはり空間座標を持たない何らかの補間関数が現れる：

In[4]:=4

✖

https://wolfram.com/xid/0tp02rczjhk79i-ecka0q

Out[5]=5

これは，軸対称の場合の式では，歪みの成分がとして計算されるためである．ここで，は最初の従属変数，は円筒座標系の半径方向である．導出は軸対称モデルのセクションで説明されている．歪みから応力が計算されるとき，の値は最初に指定された変位関数 (この場合は usol) から取得される．その関数が独立変数とに依存しない場合，それらは出力に表示されない．