સ્થિતિસ્થાપક-પ્લાસ્ટિક વિરૂપતા સંબંધિત સબરૂટિન (uyield.f90)¶
સ્થિતિસ્થાપક-પ્લાસ્ટિક કઠોરતા મેટ્રિક્સ અને તાણનું રિટર્ન મૅપિંગ ગણવા માટે સબરૂટિન આપવામાં આવ્યા છે. વપરાશકર્તા-વ્યાખ્યાયિત યીલ્ડ વિધેય વાપરતી વખતે પહેલાં ઇનપુટ ફાઇલમાં !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) તાણનું રિટર્ન મૅપિંગ ગણતું સબરૂટિન¶
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: સ્થિતિસ્થાપક વિરૂપતા ધારણા પરથી મળેલું ટ્રાયલ તાણ, દ્વિતીય Piola-Kirchhoff તાણplstrain: વર્તમાન સબસ્ટેપની શરૂઆતનું પ્લાસ્ટિક વિકૃતિistat: પૂર્ણાંક સ્થિતિ ચલfstat: વાસ્તવિક સ્થિતિ ચલોનું એરેtemp: તાપમાનhdflag: બધા ઘટકો (0), ફક્ત વિચલનાત્મક ઘટકો (1), અથવા ફક્ત ઘનફળીય ઘટકો (2) ગણવા