इलॅस्टोप्लास्टिक विकृतीशी संबंधित सबरूटीन (uyield.f90)¶
इलॅस्टोप्लास्टिक कडकपणा मॅट्रिक्स आणि प्रतिबलाचे return mapping गणना करण्यासाठी सबरूटीन उपलब्ध करून देते. वापरकर्ता-परिभाषित यील्ड फलन वापरताना प्रथम इनपुट फाइलमध्ये !PLASTIC, YIELD=USER सेट करून आवश्यक सामग्री स्थिरांक इनपुट करा, आणि त्यानंतर uElastoPlasticNumStatus फंक्शन व uElastoPlasticMatrix आणि uBackwardEuler ही सबरूटीन तयार करा.
(1) वास्तविक-संख्या स्थिती चलांची संख्या परत करणारे फंक्शन¶
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) यांची गणना करते