לדלג לתוכן

תת־שגרות לעיוות אלסטי (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: מאמץ