ಸ್ಥಿತಿಸ್ಥಾಪಕ-ಪ್ಲಾಸ್ಟಿಕ್ ವಿರೂಪಣೆಗೆ ಸಂಬಂಧಿಸಿದ ಉಪಕ್ರಮಗಳು (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: 2ನೇ 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: ಸ್ಥಿತಿಸ್ಥಾಪಕ ವಿರೂಪಣವನ್ನು ಊಹಿಸಿ ಪಡೆದ ಪರೀಕ್ಷಾ ಒತ್ತಡ, ಅಂದರೆ 2ನೇ Piola-Kirchhoff ಒತ್ತಡplstrain: ಪ್ರಸ್ತುತ ಉಪಹಂತದ ಆರಂಭದಲ್ಲಿನ ಪ್ಲಾಸ್ಟಿಕ್ ವಿಕೃತಿistat: ಪೂರ್ಣಾಂಕ ಸ್ಥಿತಿ ಚರfstat: ನಿಜ ಸಂಖ್ಯೆಯ ಸ್ಥಿತಿ ಚರಗಳ ಅರೆtemp: ತಾಪಮಾನhdflag: ಎಲ್ಲಾ ಘಟಕಗಳು (0), ವಿಚಲನ ಘಟಕಗಳು ಮಾತ್ರ (1), ಅಥವಾ ಆಯತನ ಘಟಕಗಳು ಮಾತ್ರ (2) ಲೆಕ್ಕಿಸಲಾಗುತ್ತದೆ