प्रत्यास्थ-प्लास्टिक विरूपणसँग सम्बन्धित सबरुटिन (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: प्रत्यास्थ विरूपण मानेर प्राप्त trial stress, अर्थात् 2nd Piola-Kirchhoff प्रतिबल
• plstrain: हालको सबस्टेप सुरु हुने समयको प्लास्टिक विकृति
• istat: पूर्णाङ्क अवस्था चल
• fstat: वास्तविक-संख्या अवस्था चलहरूको एरे
• temp: तापक्रम
• hdflag: सबै अवयव (0), केवल विचलन अवयव (1), वा केवल आयतनी अवयव (2) गणना गर्छ