Arruda-Boyce (/MAT/LAW92)

LAW92は、超弾性の挙動のモデル化に使用できるArruda-Boyce材料モデルを記述します。このモデルは、対角線方向に沿って8個のチェーンを擁する立方体の代表的なボリューム要素の材料の統計的な機構に基づいています。

チェーンの分子は、主ストレッチ空間内で立方体の対角線に沿って平均的に配置されています。

材料パラメータ

ひずみエネルギー密度関数は:(1)
W = μ i = 1 5 c i ( λ m ) 2 i 2 ( I ¯ 1 i 3 i ) W ( I ¯ 1 ) + 1 D ( J 2 1 2 + ln ( J ) ) U ( J ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGxbGaey ypa0ZaaGbaaeaacqaH8oqBdaaeWbqaamaalaaabaGaam4yamaaBaaa leaacaWGPbaabeaaaOqaaiaacIcacqaH7oaBdaWgaaWcbaGaamyBaa qabaGccaGGPaWaaWbaaSqabeaacaaIYaGaamyAaiabgkHiTiaaikda aaaaaaqaaiaadMgacqGH9aqpcaaIXaaabaGaaGynaaqdcqGHris5aO WaaeWaaeaaceWGjbGbaebadaqhaaWcbaGaaGymaaqaaiaadMgaaaGc cqGHsislcaaIZaWaaWbaaSqabeaacaWGPbaaaaGccaGLOaGaayzkaa aaleaacaWGxbWaaeWaaeaaceWGjbGbaebadaWgaaadbaGaaGymaaqa baaaliaawIcacaGLPaaaaOGaayjo+dGaey4kaSYaaGbaaeaadaWcaa qaaiaaigdaaeaacaWGebaaamaabmaabaWaaSaaaeaacaWGkbWaaWba aSqabeaacaaIYaaaaOGaeyOeI0IaaGymaaqaaiaaikdaaaGaey4kaS IaciiBaiaac6gaciGGOaGaamOsaiaacMcaaiaawIcacaGLPaaaaSqa aiaadwfadaqadaqaaiaadQeaaiaawIcacaGLPaaaaOGaayjo+daaaa@692A@

材料定数 c i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGJbWaaS baaSqaaiaadMgaaeqaaaaa@3861@ は:

c 1 = 1 2 , c 2 = 1 20 , c 3 = 11 1050 , c 4 = 19 7000 , c 5 = 519 673750 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yamaaBa aaleaacaaIXaaabeaakiabg2da9maalaaabaGaaGymaaqaaiaaikda aaGaaiilaiaaysW7caWGJbWaaSbaaSqaaiaaikdaaeqaaOGaeyypa0 ZaaSaaaeaacaaIXaaabaGaaGOmaiaaicdaaaGaaiilaiaaysW7caWG JbWaaSbaaSqaaiaaiodaaeqaaOGaeyypa0ZaaSaaaeaacaaIXaGaaG ymaaqaaiaaigdacaaIWaGaaGynaiaaicdaaaGaaiilaiaaysW7caWG JbWaaSbaaSqaaiaaisdaaeqaaOGaeyypa0ZaaSaaaeaacaaIXaGaaG yoaaqaaiaaiEdacaaIWaGaaGimaiaaicdaaaGaaGjbVlaacYcacaaM e8Uaam4yamaaBaaaleaacaaI1aaabeaakiabg2da9maalaaabaGaaG ynaiaaigdacaaI5aaabaGaaGOnaiaaiEdacaaIZaGaaG4naiaaiwda caaIWaaaaaaa@625B@
I ¯ 1 = λ ¯ 1 2 + λ ¯ 2 2 + λ ¯ 3 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaceWGjbGbae badaWgaaWcbaGaaGymaaqabaGccqGH9aqpcuaH7oaBgaqeamaaDaaa leaacaaIXaaabaGaaGOmaaaakiabgUcaRiqbeU7aSzaaraWaa0baaS qaaiaaikdaaeaacaaIYaaaaOGaey4kaSIafq4UdWMbaebadaqhaaWc baGaaG4maaqaaiaaikdaaaaaaa@4567@
1つ目のひずみ不変量
λ i
i t h MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyAamaaBa aaleaacaWG0bGaamiAaaqabaaaaa@38F6@ 主工学ストレッチ
LAW92を伴う材料は、次の2つの方法で定義することができます:
  • パラメータ入力

    せん断弾性係数、体積弾性係数、およびひずみストレッチ( μ , D , λ m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH8oqBca GGSaGaamiraiaacYcacqaH7oaBdaWgaaWcbaGaamyBaaqabaaaaa@3D10@

    明らかな物理的意味を有する上記の3つのパラメータのみ、材料に定義される必要があります。

    μ はゼロひずみにおけるせん断弾性係数。(2)
    D = 2 K MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGebGaey ypa0ZaaSaaaeaacaaIYaaabaGaam4saaaaaaa@39CA@
    ここで、
    K MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4saaaa@36C6@
    ゼロひずみにおける体積係数
    λ m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBda WgaaWcbaGaamyBaaqabaaaaa@3931@
    ストレッチの限界を定義
    ロッキングストレッチとも呼ばれ、引張りにおける硬化フェイズ(引張りにおけるロッキングひずみ)の開始を指定。デフォルト = 7.0


    図 1. ロッキングストレッチ λ m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBda WgaaWcbaGaamyBaaqabaaaaa@3931@
    パラメータ入力では、ポアソン比が次のように計算されます:(3)
    ν = 3 K 2 μ 6 K + 2 μ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH9oGBcq GH9aqpdaWcaaqaaiaaiodacaWGlbGaeyOeI0IaaGOmaiabeY7aTbqa aiaaiAdacaWGlbGaey4kaSIaaGOmaiabeY7aTbaaaaa@42FD@
  • 関数入力の使用時、ポアソン比 ν MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH9oGBaa a@3817@ Itypeを定義する必要があります。Itypeは、入力としてどのタイプの工学的応力ひずみ試験データを使用するかを定義します。


    • 図 2. Itype = 1: 単軸データ試験


    • 図 3. Itype = 2: 等2軸データ試験


    • 図 4. Itype = 3: 平面データ試験

ポアソン比と材料の非圧縮性

関数入力が定義されると、パラメータ μ , D , λ m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH8oqBca GGSaGaamiraiaacYcacqaH7oaBdaWgaaWcbaGaamyBaaqabaaaaa@3D10@ は無視され、Radiossは入力関数をフィッティングすることで材料定数を計算します。RadiossによるArruda-Boyceパラメータのフィッティングには、非線形最小2乗アルゴリズムが使用されます。カーブのフィッティングは、ポアソンの値が0.5に近い、すなわち材料が非圧縮性であるという仮定のもとで行われます。その他の超弾性材料モデルと同様、ポアソン比の値が0.5に近いと、体積弾性係数は大きくなり、時間ステップは小さくなります。非圧縮性とリーズナブルな時間ステップの良好なバランスのためには、ポアソン比の値は0.495が推奨されます。

材料フィッティングの情報は、Starter出力ファイル(*0000.out)で見つかります。


図 5. LAW92関数の例
フィッティングのエラーおよびフィッティングされた材料パラメータは、Starter出力ファイルに出力されます。


図 6.

粘性(速度)効果

粘性効果を含めるには、LAW92と共に/VISC/PRONYを使用する必要があります。

参考文献

1 Arruda, E. M. and Boyce, M. C., 1993, “A three-dimensional model for the large stretch behavior of rubber elastic materials”, J. Mech.Phys.Solids, 41(2), pp. 389–412