Skip to content

સ્થિતિસ્થાપક-પ્લાસ્ટિક વિરૂપતા સંબંધિત સબરૂટિન (uyield.f90)

સ્થિતિસ્થાપક-પ્લાસ્ટિક કઠોરતા મેટ્રિક્સ અને તાણનું રિટર્ન મૅપિંગ ગણવા માટે સબરૂટિન આપવામાં આવ્યા છે. વપરાશકર્તા-વ્યાખ્યાયિત યીલ્ડ વિધેય વાપરતી વખતે પહેલાં ઇનપુટ ફાઇલમાં !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) તાણનું રિટર્ન મૅપિંગ ગણતું સબરૂટિન

  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) ગણવા