Skip to content

इलॅस्टोप्लास्टिक विकृतीशी संबंधित सबरूटीन (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: द्वितीय 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: लवचिक विकृती गृहित धरून मिळालेला trial प्रतिबल, म्हणजे द्वितीय Piola-Kirchhoff प्रतिबल
  • plstrain: वर्तमान सबस्टेपच्या सुरुवातीची प्लास्टिक विकृती
  • istat: पूर्णांक स्थिती चल
  • fstat: वास्तविक-संख्या स्थिती चलांचा अॅरे
  • temp: तापमान
  • hdflag: सर्व घटक (0), फक्त विचलनात्मक घटक (1), किंवा फक्त आयतनात्मक घटक (2) यांची गणना करते