Télécharger la plaquette Cast3M


fr en

Télécharger partit.eso

Retour à la liste

Numérotation des lignes : 

partit

  1. C PARTIT    SOURCE    CB215821  17/07/21    21:15:22     9513           
  2. c---------------------------------------------------------------------
  3. c
  4.       SUBROUTINE PARTIT (NDIME, NSTR1, NTENS, NCOMP,
  5.      .                   betam, sgeff, sgten, sgcom)
  6. c
  7. c=====================================================================
  8. c                                                                    =
  9. c   This subroutine determines principal stresses and directions and =
  10. c   performs split of stress tensor into tensile and compressive     =
  11. c   contributions (eigen system).                                    =
  12. c                  ------------                                      =
  13. c   Variables:                                                       =
  14. c                                                                    =
  15. c       ntens                  no. of tension componentes            =
  16. c       ncomp                   "   " compre.      "                 =
  17. c                                                                    =
  18. c       betam                  rotation matrix                       =
  19. c       sgeff                  effect. stress tensor to be splitted  =
  20. c       sgten                  tension stress tensor contribution    =
  21. c       sgcom                  compre.    "      "        "          =
  22. c                                                                    =
  23. c       stpri                  eigen stresses                        =
  24. c                                                                    =
  25. c   Author: Rui Faria                           Version: 93.08.25    =
  26. c                                               --------             =
  27. c=====================================================================
  28.       IMPLICIT INTEGER(I-N)
  29.       integer   ndime, nstr1, ntens, ncomp
  30.       real*8    sgeff (nstr1), sgten (nstr1), sgcom (nstr1),
  31.      .          betam (ndime,ndime)
  32. c
  33.       integer   istrs
  34.       real*8    eigen, stpri (3)
  35. c=====----------------------------------------------------------------
  36. c  determination of principal stresses and directions                |
  37. c=====----------------------------------------------------------------
  38. C+CCd
  39.       if (nstr1 .eq. 6)  then
  40.           call pstres (6, sgeff, stpri)
  41.           call pdirec (3, nstr1, sgeff, stpri, betam)
  42.       else
  43.           call pdirec (2, nstr1, sgeff, stpri, betam)
  44.       end if
  45. C CCd
  46. c=====----------------------------------------------------------------
  47. c  partition of stress tensor                                        |
  48. c=====----------------------------------------------------------------
  49.       ntens = 0
  50.       ncomp = 0
  51. c
  52.       do 10 istrs = 1, 3
  53.          eigen = stpri (istrs)
  54.          if (eigen .ge. 0.0 d0)  then
  55.              ntens = ntens + 1
  56.              sgten (istrs) = eigen
  57.              sgcom (istrs) = 0.0 d0
  58.          else
  59.              ncomp = ncomp + 1
  60.              sgten (istrs) = 0.0 d0
  61.              sgcom (istrs) = eigen
  62.          end if
  63.    10 continue
  64. c
  65.       return
  66.       end
  67.  
  68.  
  69.  
  70.  

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