応力

提供: testwiki
ナビゲーションに移動 検索に移動

テンプレート:出典の明記 テンプレート:連続体力学 テンプレート:物理量 応力(おうりょく、ストレステンプレート:Lang-en-short)とは、物体[注 1]の内部に生じる力の大きさや作用方向を表現するために用いられる物理量である。物体の変形破壊などに対する負担の大きさを検討するのに用いられる。

この物理量には応力ベクトル (テンプレート:En) と応力テンソル (テンプレート:En) の2つがあり、単に「応力」といえば応力テンソルのことを指すことが多い。応力テンソルは座標系などを特別に断らない限り、主に2階の混合テンソルおよび混合ベクトルとして扱われる(混合テンソルについてはテンソル積#テンソル空間を参照)。応力ベクトルと応力テンソルは、ともに連続体内部に定義した微小面積に作用する単位面積あたりの力として定義される。そのため、それらの単位は、SIではPa (N/m2)、重力単位系ではkgf/mm2で、圧力と同じである。

異なる定義

応力という物理量は、分野によって全く異なる使われ方がなされている。即ち、土木・建築分野においては連続体内部の面にかかる力(単位:ニュートン(N))のことを応力と呼び、その単位断面積当たりの力を「応力度(stress intensity)(単位:N/m2 = Pa)」と呼んでいる[1][2][3]

応力の定義の違い
物理量 計量法、物理学、材料工学、機械工学など 土木・建築分野
力(単位:N ) 応力
単位断面積当たりの力(単位:N/m2 = Pa) 応力 応力度

以下では、計量法体系の定義[4]にあるとおり、応力を「単位断面積当たりの力」の意味で用いる。

応力ベクトル

応力ベクトルとは、物体表面あるいは物体内に仮想的な微小面を考えたとき、その微小面に作用する単位面積あたりの力であり、ベクトル(1階のテンソル)で表される。後述する応力テンソルの説明にあるように、応力テンソルσの各成分の第1の下添字は「応力成分を考えている微小面の法線の向き」を、第2の下添字は「考えている微小面に作用する力の向き」をそれぞれ表している。このことから明らかなように、微小面の単位法線ベクトルを n とすると、その微小面での応力ベクトル t は次のように与えられる。

𝒕=σT𝒏

この式はコーシーの式テンプレート:Sfnと呼ばれる。例えば、3次元デカルト座標系 (x , y , z ) において、単位法線ベクトルを 𝒏=(nx,ny,nz)=(cosα,cosβ,cosγ) と表す[注 2]と、応力ベクトルの成分 tx,ty,tz は次のようになる。

(txtytz)=(σxxnx+σyxny+σzxnzσxynx+σyyny+σzynzσxznx+σyzny+σzznz)

応力テンソル

応力テンソルは、応力ベクトルの定め方の違いから、真応力テンソル・コーシー応力テンソル、公称応力テンソル・第1パイオラ・キルヒホッフ応力テンソル、第2パイオラ・キルヒホッフ応力テンソルの3種類が定義されておりいずれも(行列の形式で記述できる)2階のテンソルとなる。ただし、これらの応力テンソルに違いが生じるのは有限変形理論に基づいて物体の運動を記述した場合であり、材料力学応用力学で多用されている微小変位・微小変形の仮定の下では、これらの応力テンソルはすべて真応力テンソルに一致する。

真応力テンソル(微小変形理論における応力テンソル)を σ で表すものとすると、その成分は座標軸を テンプレート:Math と定めた3次元デカルト座標の下では、

σ=(σxxσxyσxzσyxσyyσyzσzxσzyσzz), or, σ=(σ11σ12σ13σ21σ22σ23σ31σ32σ33) or, σ=σij(𝒆i𝒆j)

のように表される。テンプレート:Math 等は座標軸 テンプレート:Math 方向の基底ベクトルである。このとき、各成分の第1の下添字は「応力成分を考えている微小面の法線の向き」を、第2の下添字は「考えている微小面に作用する力の向き」をそれぞれ表している。例えば、σxy とは、法線の方向がx 軸の向きに一致する微小面において考えている、y 軸方向の力の成分を意味する。そのため、応力テンソルの成分には、微小面の法線と力の作用方向が一致する垂直応力 (テンプレート:En) 成分と、一致しない(異なっている)せん断応力 (テンプレート:En) 成分の2種類に分類することができる。

垂直応力とせん断応力

上に示した3次元デカルト座標系における応力テンソルの成分について考えた場合、垂直応力σxx,σyy,σzz の3成分となる。垂直応力は、力の作用面と力の作用方向とが直交し、作用面を引っ張る方向に作用した場合には引張応力 (テンプレート:En)、作用面を押し込む方向に作用した場合には圧縮応力 (テンプレート:En) と呼ばれる。材料力学応用力学構造力学などにおいては、引張応力が正の垂直応力となるように応力テンソルを定義するのが一般的であるが、地盤工学土質力学)においては圧縮応力が正の垂直応力となるように力の正の向きを定義することもある。

一方、せん断応力は、力の作用面の法線の向きと力の作用方向とが一致しない応力成分であり、σxy,σyx,σyz,σzy,σzx,σxzの6つが該当する。なお、微小変形の力学においては、せん断応力を記号τで表すことがある。この場合の応力テンソルの表記は以下のようになる。

σ=(σxτxyτxzτyxσyτyzτzxτzyσz)

応力テンソルの対称性

応力を定義している物体内でモーメントのつりあい条件(角運動量保存則)を満たすものと仮定する[注 3]と、応力テンソル(真応力テンソル)は対称テンソルとなる[注 4]。すなわち、

σ=σT

が成り立つ。例えば、上に示した3次元デカルト座標系での成分については、

σxy=σyx,σyz=σzy,σzx=σxz

が成り立ち、応力テンソルσの独立な成分は6成分となることがわかる。

この性質のため、固体物性CAEなどの分野では、独立な6成分を並べてベクトルとする表記がしばしば用いられる。これをフォークト表記 (Voigt notation)というテンプレート:Sfn

σ=(σxxσyyσzzσxyσyzσzx)(σxσyσzτxyτyzτzx)

任意座標系への応力の変換

応力テンソルの座標変換

真応力はテンソル量であり、座標系によってその成分は変化することとなる。以下のように座標系を変換する。x=Ax応力テンソルの座標系変換式は以下で表される。

σ=AσAT

ここで、 σは変換前の座標系における応力テンソル、σ' は変換後の座標系における応力テンソル、A回転行列ATA転置行列である。各成分で表すと以下の通りである。

(σ'11σ'12σ'13σ'21σ'22σ'23σ'31σ'32σ'33)=(a11a12a13a21a22a23a31a32a33)(σ11σ12σ13σ21σ22σ23σ31σ32σ33)(a11a21a31a12a22a32a13a23a33)

ここで、aij は2つの座標間の方向余弦で、各座標軸とは下記の表のような関係となる。

x1 x2 x3
x'1 a11 a12 a13
x'2 a21 a22 a23
x'3 a31 a32 a33

上式を展開すると、3次元応力状態での各応力の変換式は以下のようになる[5]

σ11=a112σ11+a122σ22+a132σ33+2a11a12σ12+2a11a13σ13+2a12a13σ23
σ22=a212σ11+a222σ22+a232σ33+2a21a22σ12+2a21a23σ13+2a22a23σ23
σ33=a312σ11+a322σ22+a332σ33+2a31a32σ12+2a31a33σ13+2a32a33σ23
σ12=a11a21σ11+a12a22σ22+a13a23σ33+(a11a22+a12a21)σ12+(a12a23+a13a22)σ23+(a11a23+a13a21)σ13
σ23=a21a31σ11+a22a32σ22+a23a33σ33+(a21a32+a22a31)σ12+(a22a33+a23a32)σ23+(a21a33+a23a31)σ13
σ13=a11a31σ11+a12a32σ22+a13a33σ33+(a11a32+a12a31)σ12+(a12a33+a13a32)σ23+(a11a33+a13a31)σ13

平面応力状態での応力変換式は以下の通りである[6]

σ11=a112σ11+a122σ22+2a11a12σ12
σ22=a212σ11+a222σ22+2a21a22σ12
σ12=a11a21σ11+a12a22σ22+(a11a22+a12a21)σ12

ここで座標軸間の角度θを用いて上式を書き直した場合は以下の通りである。

σ11=σ11cos2θ+σ22sin2θ+2σ12cosθsinθ
σ22=σ11sin2θ+σ22cos2θ2σ12cosθsinθ
σ12=(σ11σ22)cosθsinθ+σ12(cos2θsin2θ)
x1 x2
x'1 a11(=cosθ) a12(=sinθ)
x'2 a21(=sinθ) a22(=cosθ)

この変換を図示する方法として、モールの応力円が知られている。

主応力

せん断応力成分がゼロとなるように座標系を取ったときの垂直応力を主応力 (テンプレート:En) と呼ぶ。その座標系の基底ベクトルを応力テンソルの主軸あるいは主応力軸と呼ぶ。さらに主軸に垂直な面を主面あるいは主応力面と呼ぶ[7]。各点での主軸の方向(主方向)を連ねていくと、物体の中には互いに直交する曲線群を描くことができる。これを主応力線という[8]。なお、真応力テンソル(コーシー応力テンソル)は対称テンソルであるため、ある応力状態を表す真応力テンソルに対して、せん断応力が見掛け上現れず主応力のみが垂直応力として現れる主軸が必ず一組存在する。

せん断応力がゼロとなるときの垂直応力が主応力であるが、同時に主応力はあらゆる座標系の中で垂直応力が最大、最小となる値を示している[9]。3つの主応力をσ1 ≥ σ2 ≥ σ3 の関係となるようにとったとき、最大の主応力σ1最大主応力 、最小となる主応力σ3最小主応力 、これら2つに直交する主応力σ2中間主応力 と呼び、ある座標系での応力状態(σx,σy,σz,τxy,τyz,τzx)が与えられているとき、主応力は以下の関係から求められる[9]

det(σλI)=|(σxλ)τxyτzxτxy(σyλ)τyzτzxτyz(σzλ)|=0

上式を展開したλに関する3次方程式の根が主応力となる。実際に上式を展開すると、

λ3J1λ2+J2λJ3=0,J1=σx+σy+σz,J2=σxσy+σyσz+σzσxτxy2τyz2τzx2,J3=σxσyσz+2τxyτyzτzxσxτyz2σyτzx2σzτxy2

となる。一方、上式の根はσ1、σ2、σ3となるので、上式は以下のようも書き表せる。

0=(λσ1)(λσ2)(λσ3)=λ3(σ1+σ2+σ3)λ2+(σ1σ2+σ2σ3+σ3σ1)λ(σ1σ2σ3)

以上の2式を等値すれば、

J1=σ1+σ2+σ3=σx+σy+σz=tr(σ),J2=σ1σ2+σ2σ3+σ3σ1=σxσy+σyσz+σzσxτxy2τyz2τzx2,J3=σ1σ2σ3=σxσyσz+2τxyτyzτzxσxτyz2σyτzx2σzτxy2=det(σ)

を得る。J1J2J3は、ある応力状態において座標系に関わらず常に一定値となるので応力不変量(テンプレート:En)と総称される。それぞれ第一次応力不変量第二次応力不変量第三次応力不変量と呼ぶ[9]。第一次応力普遍量、第三次応力不変量は、それぞれ応力テンソルの行列式に等しい。応力不変量は以下のように表されることもある[10]

テンプレート:Math

平面応力状態における主応力

2次元における一般的な応力状態
2次元における主応力面

平面応力状態では テンプレート:Math が 0 なので、主応力は以下の関係から求められる[9]

|(σxλ)τxyτxy(σyλ)|=0

上式を展開するとλに関する2次方程式が得られ、これを解くと、平面応力状態での主応力 テンプレート:Math は次のようになる。

σ1,σ2=(σx+σy)±(σxσy)2+4τxy22

主軸の方向は次のようになる。

θ1=tanσ1σxτxy , θ2=tanσ2σxτxy

ここでθは、x 軸とσ1、σ2の主軸がなす角度である。

主せん断応力

あらゆる座標系の中で最大となるせん断応力を主せん断応力または最大せん断応力と呼ぶ。主せん断応力が働く面は、主軸に対して45°あるいは135°傾いた面となる。主せん断応力τ1、τ2、τ3 は、主応力σ1、σ2、σ3 より次式で求まる[9]

τ1=|σ2σ3|2,τ2=|σ3σ1|2,τ3=|σ1σ2|2

一般的に、主応力とは異なり、主せん断応力が働く面にはせん断応力だけでなく垂直応力も働く。

平衡方程式

外力F を受けて静的な釣り合い状態にある物体内部の任意の点では、その応力σは次の平衡方程式あるいはつりあい方程式を満たすテンプレート:Sfn

σxxx+τyxy+τzxz+Fx=0,τxyx+σyyy+τzyz+Fy=0,τxzx+τyzy+σzzz+Fz=0

あるいは次のような書き方もされる。

σij,j+Fi=0,(i=1,2,3),divσ+𝑭=0

応力場σが平衡方程式と、表面力規定境界∂Rt における境界条件(コーシーの式)

𝒕=σT𝒏,atRt

を満たすとき、その応力場σを静的に許容な場というテンプレート:Sfn

パイオラ・キルヒホッフ応力テンソル

真応力(コーシー応力)テンソルσと変形勾配テンソルF を用いて定義される次のテンソルをパイオラ・キルヒホッフ応力テンソル(Piola-Kirchhoff stress tensor)というテンプレート:Sfn

第1パイオラ・キルヒホッフ応力テンソル
Π=det(F)σ(F1)T
第2パイオラ・キルヒホッフ応力テンソル
K=F1Π=det(F)F1σ(F1)T

真応力に関するコーシーの式は上述のとおり現配置での応力ベクトルt と法線ベクトルn で表されるが、パイオラ・キルヒホッフ応力テンソルを用いても類似の関係式が成り立つ。

𝒕0=Π𝒏0,𝒕^=K𝒏0

ここで、

  • 𝒏0:基準配置の微小面の法線ベクトル
  • 𝒕0:現配置の微小面に作用している力を、基準配置の微小面の面積で割って定義される応力ベクトル
  • 𝒕^:現配置の微小面に作用している力を基準配置で求めなおし、それを基準配置の微小面の面積で割って定義される応力ベクトル

である。

仮想仕事の原理を適用する際には、これらの応力テンソルと共役な関係にあるひずみテンソルは以下のようになる。

  • コーシー応力 - アルマンシーひずみ
  • 第1パイオラ・キルヒホッフ応力 - 変形勾配
  • 第2パイオラ・キルヒホッフ応力 - グリーンひずみ

偏差応力

偏差応力(テンプレート:En) は、応力テンソルからその等方成分を差し引いたものとして定義される。物体に等方的な圧縮・引張り以外のせん断変形が生じた場合に、偏差応力が発生する。偏差応力 dev[σ] は次のように定義される。

dev[σ]:=σ+pI

ここでI は2階の単位テンソル、

p:=13tr[σ]=σkk3

は非決定(静水圧)応力であり、平均応力(3つの主応力の平均値)のマイナスに等しい。テンプレート:Math は平均応力テンソルと呼ばれる[11]

偏差応力の固有値 テンプレート:Math は、元の応力テンソルの固有値(主応力)と次の関係がある[12]

si=σi+p,i=1,2,3

偏差応力の主軸は元の応力テンソルの主軸と一致する。

材料の降伏と等価応力

上記にあるとおり、応力は3次元的なテンソルである。一般の応力について材料の特性値を調べるのは困難であるため、降伏に対して等価とみなせる1軸応力に対応するスカラー量である等価応力に換算すると便利である。等価応力は材料の降伏する条件に応じて以下のようなものがある。

最大主応力説

ある点で最大主応力 σ1 が材料の降伏を決定するというのが最大主応力説である。すなわち、

σ1σY

が降伏の条件である。ここで σY は材料の降伏応力である。最大主応力説はガラスなどの脆性材料で良く当てはまる[13]

せん断ひずみエネルギー説

単位体積あたりのせん断ひずみエネルギーが限界を越えると、材料が破壊されるという説である。ともいう。全ひずみエネルギーから静ひずみエネルギーを差し引いたせん断ひずみエネルギー U を評価基準とする。

U=1+ν6E((σ1σ2)2+(σ2σ3)2+(σ3σ1)2)

ここで、νはポアソン比Eヤング率である。

せん断ひずみエネルギーに比例する相当応力をMisesの相当応力σMisesとよび、主応力を用いて以下の式で表される。

σMises2=(σ1σ2)2+(σ2σ3)2+(σ3σ1)22

降伏条件は以下の通り。

σMisesσY

せん断ひずみエネルギー説は鋼材などの延性材料に比較的良く当てはまる[13]

最大せん断応力説

延性材料が降伏するときすべりが観察されることに着目し、最大せん断応力が降伏を決定するという説を最大せん断応力説、またはトレスカの応力説と呼ぶ。このときに用いられる相当応力をトレスカ応力とよび、最大せん断応力を記号 τmax 、トレスカ応力をσTrescaで表すと、主応力とは次式に示す関係がある。

τmax=12(max|σiσj|),σTresca=2τmax

降伏条件は以下の通り。

σTrescaσY

最大せん断応力説も延性材料に当てはまることが多い[13]。また、σTresca ≥ σ1, σTresca ≥ σMises であり、上記2説に対して安全側であることから評価基準として利用されることがある。

残留応力

残留応力とは、外力が作用していない物体の内部に生じている応力である。残留応力は機械的または熱的な原因で物体に不均一に弾塑性変形が生じることにより発生する[14]テンプレート:Main

脚注

テンプレート:脚注ヘルプ

注釈

テンプレート:Reflist

出典

テンプレート:Reflist

参考文献

テンプレート:Refbegin

テンプレート:Refend

関連項目

外部リンク

テンプレート:Commonscat テンプレート:Wiktionary

テンプレート:Normdaten


引用エラー: 「注」という名前のグループの <ref> タグがありますが、対応する <references group="注"/> タグが見つかりません

  1. 平凡社大百科事典、応力の項、「応力の大きさは単位面積に作用する内力のおおきさにより定義され、これを応力度あるいは応力強さともいうが、一般には応力度のことを単に応力と呼び・・・」第2巻、p.831、平凡社、1984年11月2日 初版
  2. テンプレート:Cite webウェイバックマシンより)
  3. テンプレート:Cite web
  4. 計量単位令 別表第一 項番23、応力、「一平方メートルにつき一ニュートンの応力」
  5. 引用エラー: 無効な <ref> タグです。「弾性力学_8-9」という名前の注釈に対するテキストが指定されていません
  6. 引用エラー: 無効な <ref> タグです。「弾性力学_5-6」という名前の注釈に対するテキストが指定されていません
  7. 引用エラー: 無効な <ref> タグです。「機械工学辞典_567-568」という名前の注釈に対するテキストが指定されていません
  8. テンプレート:Cite
  9. 9.0 9.1 9.2 9.3 9.4 引用エラー: 無効な <ref> タグです。「弾性力学_10-15」という名前の注釈に対するテキストが指定されていません
  10. テンプレート:Cite
  11. テンプレート:Cite
  12. テンプレート:Cite
  13. 13.0 13.1 13.2 引用エラー: 無効な <ref> タグです。「材料強度_9-12」という名前の注釈に対するテキストが指定されていません
  14. テンプレート:Cite