අන්තර්ගතය වෙත යන්න

ප්‍රත්‍යාස්ථ-ප්ලාස්ටික් විරූපණයට අදාළ උපරූටීන (uyield.f90)

ප්‍රත්‍යාස්ථ-ප්ලාස්ටික් දෘඪතා න්‍යාසය සහ ආතතියේ ආපසු-නිරූපණය (return mapping) ක්‍රමය ගණනය කිරීමට උපරූටීනය සපයයි. පරිශීලක-නිර්වචිත අවනති ශ්‍රිතය භාවිත කරන විට, පළමුව ආදාන ගොනුවේ !PLASTIC, YIELD=USER සකසා අවශ්‍ය ද්‍රව්‍ය නියතයන් ඇතුළත් කළ යුතු අතර, ඊළඟට uElastoPlasticNumStatus ශ්‍රිතය සහ uElastoPlasticMatrix, uBackwardEuler උපරූටීනය නිර්මාණය කළ යුතුය.

(1) සැබෑ සංඛ්‍යා තත්ත්ව විචල්‍ය ගණන ආපසු ලබාදෙන ශ්‍රිතය

integer function uElastoPlasticNumStatus( matl )
    real(kind=kreal),   intent(in)  :: matl(:)
  • matl: ද්‍රව්‍ය නියත සුරකින අරාවක් (1-100: පද්ධතියෙන් අර්ථ දක්වන ද්‍රව්‍ය නියත, 101-200: පරිශීලක-නිර්වචිත ද්‍රව්‍ය නියත)

(2) ප්‍රත්‍යාස්ථ-ප්ලාස්ටික් දෘඪතා න්‍යාසය ගණනය කරන උපරූටීනය

  subroutine uElastoPlasticMatrix( matl, stress, istat, fstat, plstrain, D, temp, hdflag )
    real(kind=kreal),   intent(in)  :: matl(:)
    real(kind=kreal),   intent(in)  :: stress(6)
    integer(kind=kint), intent(in)  :: istat
    real(kind=kreal),   intent(in)  :: fstat(:)
    real(kind=kreal),   intent(in)  :: plstrain
    real(kind=kreal),   intent(out) :: D(:,:)
    real(kind=kreal),   intent(in)  :: temp
    integer(kind=kint), intent(in)  :: hdflag
  • matl: ද්‍රව්‍ය නියත සුරකින අරාවක් (1-100: පද්ධතියෙන් අර්ථ දක්වන ද්‍රව්‍ය නියත, 101-200: පරිශීලක-නිර්වචිත ද්‍රව්‍ය නියත)
  • stress: 2nd Piola-Kirchhoff ආතතිය
  • istat: පූර්ණ සංඛ්‍යා තත්ත්ව විචල්‍යය
  • fstat: සැබෑ සංඛ්‍යා තත්ත්ව විචල්‍ය අරාවක්
  • plstrain: වත්මන් උපපියවර ආරම්භයේ ප්ලාස්ටික් වික්‍රියාව
  • D: ප්‍රත්‍යාස්ථ-ප්ලාස්ටික් න්‍යාසය
  • temp: උෂ්ණත්වය
  • hdflag: සියලු සංරචක (0), අපගමන සංරචක පමණක් (1), පරිමාමිතික සංරචක පමණක් (2) ගණනය කරයි

(3) ආතතියේ ආපසු-නිරූපණය (return mapping) ගණනය කරන උපරූටීනය

  subroutine uBackwardEuler( matl, stress, plstrain, istat, fstat, temp, hdflag )
    real(kind=kreal),   intent(in)    :: matl(:)
    real(kind=kreal),   intent(inout) :: stress(6)
    real(kind=kreal),   intent(in)    :: plstrain
    integer(kind=kint), intent(inout) :: istat
    real(kind=kreal),   intent(inout) :: fstat(:)
    real(kind=kreal),   intent(in)    :: temp
    integer(kind=kint), intent(in)    :: hdflag
  • matl: ද්‍රව්‍ය නියත සුරකින අරාවක් (1-100: පද්ධතියෙන් අර්ථ දක්වන ද්‍රව්‍ය නියත, 101-200: පරිශීලක-නිර්වචිත ද්‍රව්‍ය නියත)
  • stress: පරීක්ෂණ ආතතිය, ප්‍රත්‍යාස්ථ විරූපණය උපකල්පනය කර ලබාගත් 2nd Piola-Kirchhoff ආතතිය
  • plstrain: වත්මන් උපපියවර ආරම්භයේ ප්ලාස්ටික් වික්‍රියාව
  • istat: පූර්ණ සංඛ්‍යා තත්ත්ව විචල්‍යය
  • fstat: සැබෑ සංඛ්‍යා තත්ත්ව විචල්‍ය අරාවක්
  • temp: උෂ්ණත්වය
  • hdflag: සියලු සංරචක (0), අපගමන සංරචක පමණක් (1), පරිමාමිතික සංරචක පමණක් (2) ගණනය කරයි