Télécharger la plaquette Cast3M


fr en

Télécharger eigglo.eso

Retour à la liste

Numérotation des lignes : 

eigglo

  1. C EIGGLO    SOURCE    CB215821  17/07/21    21:15:08     9513           
  2. c---------------------------------------------------------------------
  3. c
  4.       SUBROUTINE EIGGLO (NDIME, NSTR1, NTENS,
  5.      .                   betam, prosi, sgten, soxyz)
  6. c
  7. c=====================================================================
  8. c                                                                    =
  9. c   This subroutine performs rotation of tension stress matrix,      =
  10. c   from eigen to global system (optimized).                         =
  11. c                                                                    =
  12. c    Input variable:                                                 =
  13. c                                                                    =
  14. c       ntens                  no. of tension componentes            =
  15. c       betam                  rotation matrix                       =
  16. c       sgten (i)  (i=1,3)     principal tension stress tensor       =
  17. c                                                                    =
  18. c    Output variable:                                                =
  19. c                                                                    =
  20. c       sgten                  tension stress tensor contribution    =
  21. c                                                                    =
  22. c   Author: Rui Faria                           Version: 96.09.18    =
  23. c                                               --------             =
  24. c=====================================================================
  25.       IMPLICIT INTEGER(I-N)
  26.       integer    ndime, nstr1, ntens
  27.       real*8     sgten (nstr1),
  28.      .           betam (ndime,ndime), prosi (ndime,ndime),
  29.      .           soxyz (ndime,ndime)
  30. c
  31.       integer    istrs, jstrs
  32.       real*8     provi
  33. c
  34.       if (ntens .gt. 0)  then
  35.           do 10 istrs = 1, ndime
  36.              do 20 jstrs = 1, ndime
  37.                 prosi (istrs,jstrs) = betam (jstrs,istrs) *
  38.      .                                sgten (jstrs)
  39.    20        continue
  40.    10     continue
  41. c
  42.           call proma1 (soxyz, prosi, betam, ndime, ndime, ndime)
  43. c
  44.           if (nstr1 .eq. 6)  then
  45.               sgten (1) = soxyz (1,1)
  46.               sgten (2) = soxyz (2,2)
  47.               sgten (3) = soxyz (3,3)
  48.               sgten (4) = soxyz (1,2)
  49.               sgten (5) = soxyz (1,3)
  50.               sgten (6) = soxyz (2,3)
  51.           else
  52.               sgten (1) = soxyz (1,1)
  53.               sgten (2) = soxyz (2,2)
  54.               provi     = sgten (3)
  55.               sgten (3) = provi
  56.               sgten (4) = soxyz (1,2)
  57.           end if
  58.       else
  59.           call zeror  (sgten, nstr1)
  60.       end if
  61. c
  62.       return
  63.       end
  64.  
  65.  
  66.  
  67.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales