پرش به محتویات

زیرروال‌های مربوط به تغییرشکل الاستیک (uelastic.f90)

زیرروال‌هایی برای محاسبه ماتریس سختی الاستیک و به‌روزرسانی تنش در مسائل الاستیک و هایپرالاستیک ارائه می‌شوند. هنگام استفاده از رابطه ساختاری الاستیک یا هایپرالاستیک تعریف‌شده توسط کاربر، ابتدا باید !ELASTIC, TYPE=USER یا !HYPERELASTIC, TYPE=USER را در فایل ورودی تنظیم کرده و ثابت‌های ماده لازم را وارد کنید؛ سپس لازم است زیرروال‌های uElasticMatrix و uElasticUpdate را ایجاد کنید.

(1) زیرروال محاسبه ماتریس سختی الاستیک

subroutine uElasticMatrix( matl, strain, D )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: D(6,6)
  • matl: آرایه‌ای برای ذخیره ثابت‌های ماده (حداکثر 100)
  • strain: کرنش گرین–لاگرانژ
  • D: ماتریس الاستیک

(2) زیرروال محاسبه تنش

subroutine uElasticUpdate ( matl, strain, stress )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: stress(6)
  • matl: آرایه‌ای برای ذخیره ثابت‌های ماده (حداکثر 100)
  • strain: کرنش گرین–لاگرانژ
  • stress: تنش