Υπορουτίνες για ελαστική παραμόρφωση (uelastic.f90)¶

Παρέχονται υπορουτίνες για τον υπολογισμό του ελαστικού πίνακα δυσκαμψίας και την ενημέρωση της τάσης σε ελαστικά και υπερελαστικά προβλήματα. Για να χρησιμοποιήσετε καταστατική σχέση ελαστικότητας ή υπερελαστικότητας που ορίζεται από τον χρήστη, ορίστε πρώτα !ELASTIC, TYPE=USER ή !HYPERELASTIC, TYPE=USER στο αρχείο εισόδου και εισαγάγετε τις απαιτούμενες σταθερές υλικού, και κατόπιν υλοποιήστε τις υπορουτίνες uElasticMatrix και uElasticUpdate.

(1) Υπορουτίνα για τον υπολογισμό του ελαστικού πίνακα δυσκαμψίας¶

subroutine uElasticMatrix( matl, strain, D )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: D(6,6)

• matl: Πίνακας που αποθηκεύει τις σταθερές υλικού (έως 100)
• strain: Παραμόρφωση Green-Lagrange
• D: Ελαστικός πίνακας

(2) Υπορουτίνα για τον υπολογισμό της τάσης¶

subroutine uElasticUpdate ( matl, strain, stress )
    REAL(KIND=kreal), INTENT(IN) :: matl(:)
    REAL(KIND=kreal), INTENT(IN) :: strain(6)
    REAL(KIND=kreal), INTENT(OUT) :: stress(6)

• matl: Πίνακας που αποθηκεύει τις σταθερές υλικού (έως 100)
• strain: Παραμόρφωση Green-Lagrange
• stress: Τάση