SolidMechanicsPDEComponent

SolidMechanicsPDEComponent[vars,pars]

変数 vars,パラメータ pars の固体力学の偏微分方程式の項を与える.

詳細

  • SolidMechanicsPDEComponentは,固体力学解析のための偏微分方程式を返す.
  • SolidMechanicsPDEComponentは,偏微分方程式の部分として使われる微分演算子の和を返す.
  • SolidMechanicsPDEComponentは,加えられた荷重と拘束を受ける物体の結果の変位をモデル化する.
  • SolidMechanicsPDEComponentは,定常,時間依存,パラメトリック,周波数応答,固有モードの各分析のための偏微分方程式の成分を作成する.
  • SolidMechanicsPDEComponentは,従属変数としての変位 (単位:[TemplateBox[{InterpretationBox[, 1], "m", meters, "Meters"}, QuantityTF]]),独立変数 (単位:[TemplateBox[{InterpretationBox[, 1], "m", meters, "Meters"}, QuantityTF]]),時間変数 (単位:秒 [TemplateBox[{InterpretationBox[, 1], "s", seconds, "Seconds"}, QuantityTF]]),角周波数 (単位:ラジアン毎秒)で固体力学現象をモデル化する.
  • SolidMechanicsPDEComponentは2Dおよび3Dで偏微分方程式成分を作成する.
  • 定常変数 varsvars={{u[x1,,xn],v[x1,,xn],},{x1,,xn}}である.
  • 時間依存変数または固有モード変数の varsvars={{u[t,x1,,xn],v[t,x1,,xn],},t,{x1,,xn}}である.
  • 周波数依存変数 varsvars={{u[x1,,xn],v[x1,,xn],},ω,{x1,,xn}}である.
  • 異なる解析タイプのための方程式SolidMechanicsPDEComponentvars の形式によって生成される.
  • 質量密度 [TemplateBox[{InterpretationBox[, 1], {"kg", , "/", , {"m", ^, 3}}, kilograms per meter cubed, {{(, "Kilograms", )}, /, {(, {"Meters", ^, 3}, )}}}, QuantityTF]],応力テンソル [TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]],歪みテンソル ,変位ベクトル [TemplateBox[{InterpretationBox[, 1], "m", meters, "Meters"}, QuantityTF]],物体荷重ベクトル [TemplateBox[{InterpretationBox[, 1], {"N", , "/", , {"m", ^, 3}}, newtons per meter cubed, {{(, "Newtons", )}, /, {(, {"Meters", ^, 3}, )}}}, QuantityTF]]または [TemplateBox[{InterpretationBox[, 1], {"m", , "/", , {"s", ^, 2}}, meters per second squared, {{(, "Meters", )}, /, {(, {"Seconds", ^, 2}, )}}}, QuantityTF]]の固体力学の偏微分方程式SolidMechanicsPDEComponentの定常平衡方程式は以下に基づく.
  • 固体力学モデルSolidMechanicsPDEComponentの時間依存平衡方程式は以下に基づく.
  • 固有値が の固体力学モデルSolidMechanicsPDEComponentの固有周波数方程式は以下に基づく.
  • 角周波数 の固体力学モデルSolidMechanicsPDEComponentの周波数応答方程式は以下に基づく.
  • 固体力学モデル項の単位は力密度(単位:)である.
  • 次のパラメータ pars を与えることができる.
  • パラメータデフォルトシンボル
    "AnalysisType"Automaticなし
    "BodyLoad"0
  • ,体積力密度(単位:[TemplateBox[{InterpretationBox[, 1], {"N", , "/", , {"m", ^, 3}}, newtons per meter cubed, {{(, "Newtons", )}, /, {(, {"Meters", ^, 3}, )}}}, QuantityTF]]
  • "BodyLoadValue"0
  • ,体加速度(単位:[TemplateBox[{InterpretationBox[, 1], {"m", , "/", , {"s", ^, 2}}, meters per second squared, {{(, "Meters", )}, /, {(, {"Seconds", ^, 2}, )}}}, QuantityTF]]
  • "MassDensity"-
  • ,密度(単位:[TemplateBox[{InterpretationBox[, 1], {"kg", , "/", , {"m", ^, 3}}, kilograms per meter cubed, {{(, "Kilograms", )}, /, {(, {"Meters", ^, 3}, )}}}, QuantityTF]]
  • "Material"-なし
    "MaterialSymmetry""Isotropic"なし
    "SolidMechanicsMaterialModel""LinearElastic"なし
    "SolidMechanicsModelForm""Solid"なし
  • "AnalysisType"Automaticなら,生成されるモデルは vars の形式に依存する.
  • 固有周波数解析については,"AnalysisType""Eigenmode"に設定し,時間依存変数 tvars を使う必要がある.
  • "Material"が指定されているなら,材料定数は材料データから抽出される.それ以外の場合は関連する材料パラメータを指定しなければならない.
  • 線形弾性等方性材料については,任意の2つの係数が使用できる.
  • パラメータ名
    "BulkModulus"
    "LameParameter"
    "PoissonRatio"
    "PWaveModulus"
    "ShearModulus"
    "YoungModulus"
  • パラメータ"SolidMechanicsModelForm"は,"Solid""PlaneStress"あるいは"PlaneStrain"でよい.
  • "PlaneStress"モデルと"PlaneStress"モデルについては,"Thickness"パラメータを定義しなければならない.
  • デフォルトの材料モデルは線形弾性等方性材料モデルである.
  • 線形弾性で小さな変形があるモデルについて,次の材料対称性が使用可能である.
  • 材料対称性名
    "Isotropic"
    "Orthotropic"
    "TransverselyIsotropic"
    "Anisotropic"
  • 次の等方性可縮超弾性の大きく変形する材料モデルが利用できる.
  • 材料モデル名
    "StVenantKirchhoff"
    "NeoHookean"
  • 次のほぼ等方性非可縮超弾性の大きく変形する材料モデルが利用できる.
  • 材料モデル名
    "ArrudaBoyce"
    "Gent"
    "MooneyRivlin"
    "NeoHookean"
    "Yeoh"
  • 運動方程式は,微量な小さい変形歪み測定値を使用する.
  • 非線形材料法則の場合は,運動方程式は変形勾配 に基づいてグリーン・ラグランジュ(GreenLagrange)歪みの測定値を使う.ここで, は恒等行列である.
  • 弾性行列 ,初期応力 ,初期歪み ,熱歪み の線形弾性材料モデルの構成方程式は以下で与えられる.
  • 垂直歪み成分 とせん断歪み成分 は次の次数のフォークト(Voigt)記法を使う.
  • 垂直応力成分 とせん断応力成分 は次の次数のフォークト(Voigt)記法を使う.
  • "PlaneStrain"モデルについて 方向の歪みが0であると仮定する.
  • "PlaneStress"モデルについて 方向の応力が0であると仮定する.
  • "RegionSymmetry""Axisymmetric"は切頂円筒座標系を使い, 方向に の変位があってであると仮定する.である点に注意のこと.
  • ヤング率 ,ポアソン比 の線形弾性等方性材料についての構成方程式は以下で与えられる.
  • 熱歪み は,熱膨張係数 (単位:[1TemplateBox[{InterpretationBox[, 1], {"/", , "K"}, reciprocal kelvins, {1, /, {(, "Kelvins", )}}}, QuantityTF]]),熱歪み温度(単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]]),熱歪み参照温度 (単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]])で加えることができる.
  • "LinearElastic" "Isotropic"材料モデルには,次のサブパラメータが使用できる.
  • パラメータデフォルトシンボル
    "InitialStrain"0,初期歪み
    "InitialStress"0,初期応力(単位:[TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]]
    "PoissonRatio"Automatic,ポアソン比
    "ThermalExpansion"0,熱膨張係数(単位:[1TemplateBox[{InterpretationBox[, 1], {"/", , "K"}, reciprocal kelvins, {1, /, {(, "Kelvins", )}}}, QuantityTF]]
    "ThermalStrainTemperature"0,温度(単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]]
    "ThermalStrainReferenceTemperature"0,温度(単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]]
    "YoungModulus"Automatic,ヤング率(単位:[TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]]
  • コンプライアンス行列 を含む線形弾性直交異方性材料モデルの構成式は以下で与えられる.
  • 弾性行列はコンプライアンス行列 の逆行列である.
  • せん断弾性率が の線形弾性直交異方性コンプライアンス行列 は以下で与えられる.
  • 線形弾性直交異方性材料モデルについては,熱弾性係数 は方向に依存する.
  • "LinearElastic" "Orthotropic"材料モデルには次のサブパラメータを使うことができる.
  • パラメータデフォルトシンボル
    "InitialStrain"0,初期歪み
    "InitialStress"0,初期応力(単位:[TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]])
    "PoissonRatio"-, , , , , ポアソン比
    "ShearModulus"-, , 剛性率(単位:[TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]])
    "ThermalExpansion"0, , 熱膨張係数(単位:[TemplateBox[{{InterpretationBox[, 1], 1}, {"/", , "K"}, reciprocal kelvins, {1, /, {(, "Kelvins", )}}}, QuantityTF]])
    "ThermalStrainTemperature"0,温度(単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]]
    "ThermalStrainReferenceTemperature"0,温度(単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]]
    "YoungModulus"-, , ,ヤング率(単位:[TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]])
  • ポアソン(Poisson)比,剛性率,ヤング(Young)率は指標付きの形式変数として指定される.
  • 異方性材料モデルについては,完全弾性行列 を指定しなければならない.
  • 代りにコンプライアンス行列 を指定してもよい.
  • 線形弾性異方性材料モデルについては,熱膨張 の係数は方向に依存する.
  • "LinearElastic" "Anisotropic"材料モードについては次のサブパラメータを使うことができる.
  • パラメータデフォルトシンボル
    "ComplianceMatrix"-,コンプライアンス行列
    "ElasticityMatrix"-,弾性行列
    "InitialStrain"0,初期歪み
    "InitialStress"0,初期応力(単位:[TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]]
    "ThermalExpansion"0熱膨張係数 , , , , ,(単位: [TemplateBox[{{InterpretationBox[, 1], 1}, {"/", , "K"}, reciprocal kelvins, {1, /, {(, "Kelvins", )}}}, QuantityTF]]
    "ThermalStrainTemperature"0,温度(単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]]
    "ThermalStrainReferenceTemperature"0,温度(単位:[TemplateBox[{InterpretationBox[, 1], "K", kelvins, "Kelvins"}, QuantityTF]]
  • 圧縮性およびほぼ非圧縮性の超弾性材料モデルが使用可能である.使用可能な場合は,ほぼ非圧縮性のモデルがデフォルトである.
  • "ArrudaBoyce"材料モデルはエネルギー密度関数 に基づいている.ただし, はネットワーク中のポリマー鎖の数, は単一の鎖内のセグメントの数,はボルツマン定数, は絶対温度, は鎖のストレッチ,L はランジュバン(Langevin)関数である.
  • "Gent"材料モデルはエネルギー密度関数 に基づいている.ただし,は最初の歪み不変量, はせん断弾性率,は極限値である.
  • "NeoHookean"材料モデルはエネルギー密度関数 に基づいている.ただし, は第1ラメ(Lamé)定数, は第2定数, は右コーシー(Cauchy)・グリーン(Green)テンソル, は変形勾配である.
  • 圧縮可能な"StVenantKirchhoff"材料モデルはエネルギー密度関数 に基づいている.ただし, は第1ラメ定数, は第2定数, はグリーン・ラグランジュ歪みである.
  • "MooneyRivlin"材料モデルはエネルギー密度関数 に基づいている.ただし,は材料係数,は第1および第2の等重歪み不変量である.
  • "Yeoh"材料モデルはエネルギー密度関数 に基づいている.ただし,はモデル定数で は最初の不変量である.
  • 可能な場合はモデルの"Compressibility""NearlyIncompressibile"または"Compressibile"として指定できる.
  • ほぼ圧縮不可能性は,歪みエネルギー密度関数 に加えられた静水圧 として実装されている. は材料の体積弾性率である.
  • すべての超弾性モデルで,平面歪みと平面応力のモデル形式が使用できる.
  • 超弾性材料モデルはどれも標準の強化材料モデルを利用して繊維強化材料などの横方向等方性材料がモデル化できる.
  • SolidMechanicsPDEComponent"SIBase"単位を使う.幾何学配置は偏微分方程式と同じ単位でなければならない.
  • SolidMechanicsPDEComponent,keypi,pivi,]として連想 pars で指定されているパラメータ に依存するなら,パラメータ で置き換えられる.

例題

すべて開くすべて閉じる

  (3)

固体力学の偏微分方程式モデルを定義する:

固体力学の偏微分方程式モデルを記号的に定義する:

固体力学の時間依存偏微分方程式モデルを記号的に定義する:

スコープ  (19)

基本的な例題  (3)

固体力学の偏微分方程式モデルをアクティベートする:

定常固体力学モデルを特定の材料で定義する:

材料の値を指定してモデルを定義する:

定常解析  (2)

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

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

たわみを可視化する:

熱膨張を考慮する固体力学の偏微分方程式を記号的に定義する:

定常平面応力解析  (3)

左端が固定されて右端が強制的に動かされる長方形鋼板の変位を計算する.領域,変数,パラメータを設定する:

方程式を解く:

変位を可視化する:

左端が固定され右端に圧力が加えられた長方形鋼板の変位を計算する.領域,変数,パラメータを設定する:

方程式を解く:

変位を可視化する:

底辺が固定され他の3辺に圧力が加えられた長方形鋼板の変位を計算する.領域,変数,パラメータを設定する:

方程式を解く:

変位を可視化する:

歪みを計算する:

方向の垂直歪みがほぼ0であることを確かめる:

方向の垂直歪みを可視化する:

方向の垂直歪みがほぼ0であることを確かめる:

方向の垂直歪みを可視化する:

せん断歪みを求める:

- 方向の垂直歪みを可視化する:

定常平面歪み解析  (1)

2D平面歪みの偏微分方程式モデルを定義する:

時間依存解析  (3)

特定の材料のための時間依存固体力学モデルを定義する:

梁にかかる時間依存の力のシミュレーションを行う.領域,変数,材料を設定する:

偏微分方程式モデルを設定する:

梁の右端にかかる時間依存の力を作成する:

梁の左端を固定する:

初期条件と初期速度条件を0に設定する:

時間依存偏微分方程式を解き,進捗状況をモニターする:

参照点における 変位を経時的に可視化する:

レイリー(Rayleigh)減衰モデルを考慮して,スプーンの断面にかかる時間依存の力のシミュレーションを行う.領域,変数,材料を設定する:

偏微分方程式を設定し,過程をモニターしながらこれを解く:

クエリ点における時間の経過に伴う の変位を可視化する:

変位を可視化する:

固有モード解析  (2)

固有モードの固体力学モデルを特定の材料で定義する:

鉄のブラケットの固有モードを計算する.領域,変数,材料を設定する:

固有値と固有モードを計算する:

7番目から10番目までの固有モードを可視化する:

超弾性材料モデル  (2)

左側が固定され,右側に圧力が加えられたゴム片の変数とモデルパラメータを設定する.Neo-Hookean超弾性材料モデルが使用される:

形状を設定する:

方程式を設定して変位について解く:

変形された物体を物体のもとの形状上に可視化する:

Neo-Hookean超弾性材料モデルの変数とモデルパラメータを設定する:

境界条件がある偏微分方程式モデルを作成する:

実際の変位を表示する:

歪みを計算する:

応力を計算する:

繊維強化材料モデル  (1)

左側が固定され,右側に圧力が加えられた繊維強化ゴム片の変数とモデルパラメータを設定する.Neo-Hookean超弾性材料モデルが使用される:

形状を設定する:S

形状における繊維強化を可視化する:

方程式を設定して変位について解く:

変形された物体を物体のもとの形状の上に可視化する:

歪みを計算する:

応力を計算する:

vonMises応力を計算する:

vonMises応力を可視化する:

マルチマテリアルモデル  (1)

コンプライアンス行列は行列として指定する必要がある.これはマルチマテリアルモデルにも当てはまる.この例ではその方法を示す.2つのコンプライアンス行列を指定する:

マルチマテリアルコンプライアンス行列を作成するマテリアル1は の幾何形状の値に使用され,マテリアル2はその他すべての場合に使用される:

コンプライアンス行列が行列である点に注意のこと:

特定の成分を見る:

パラメータを設定してマルチマテリアルコンプライアンス行列を利用する:

固体力学 PDE成分を作成する:

非軸整列材料モデルl  (1)

軸が揃っていない異方性弾性材料の場合,方向行列を指定して材料が調整できる.ここに示されているように,これは熱膨張にも適用される.変数とパラメータを次のように設定する:

偏微分方程式を解く:

変位を可視化する:

アプリケーション  (1)

地質工学  (1)

地質工学のアプリケーションで土をモデリングする際は,土の深さでヤング率を変えることができる.この例は,位置に依存するヤング率について調べる.縦横がそれぞれ100メートルの長方形の固体スラブを使う:

次に,変数とパラメータを設定する.この時点でヤング率の記号 を使う:

偏微分方程式を設定する:

上部に部分的に負の 方向に作用する力がある:

左側と右側には 方向の制約がある.土は 方向には自由に動くことができる.底では,土は 方向に動くことはできるが 方向には動けない.以下で,土が動かない「より固い」地面に立っている場合をモデル化する:

定ヤング率を設定する:

ヤング率 baseYoungModulusに設定し,細分化されたメッシュ上で偏微分方程式を解く:

移動を可視化する:

von Mises応力を計算するヘルパー関数を作成する:

定ヤング率に対するvon Mises応力を計算する:

von Mises応力の等高線プロットを作成する:

次に,深さ に依存するヤング率を作成する:

関数を可視化する:

を変数と深さのヤング率に置き換えて,細分化メッシュ上で偏微分方程式を解く:

移動を可視化する:

von Mises応力を計算する:

von Mises応力を可視化する:

2つのモデルのvon Mises応力の差をプロットする:

2つのモデルのvon Mises応力の差を,プロット範囲全体を使ってプロットする:

Wolfram Research (2021), SolidMechanicsPDEComponent, Wolfram言語関数, https://reference.wolfram.com/language/ref/SolidMechanicsPDEComponent.html (2024年に更新).

テキスト

Wolfram Research (2021), SolidMechanicsPDEComponent, Wolfram言語関数, https://reference.wolfram.com/language/ref/SolidMechanicsPDEComponent.html (2024年に更新).

CMS

Wolfram Language. 2021. "SolidMechanicsPDEComponent." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2024. https://reference.wolfram.com/language/ref/SolidMechanicsPDEComponent.html.

APA

Wolfram Language. (2021). SolidMechanicsPDEComponent. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/SolidMechanicsPDEComponent.html

BibTeX

@misc{reference.wolfram_2024_solidmechanicspdecomponent, author="Wolfram Research", title="{SolidMechanicsPDEComponent}", year="2024", howpublished="\url{https://reference.wolfram.com/language/ref/SolidMechanicsPDEComponent.html}", note=[Accessed: 24-November-2024 ]}

BibLaTeX

@online{reference.wolfram_2024_solidmechanicspdecomponent, organization={Wolfram Research}, title={SolidMechanicsPDEComponent}, year={2024}, url={https://reference.wolfram.com/language/ref/SolidMechanicsPDEComponent.html}, note=[Accessed: 24-November-2024 ]}