Télécharger la plaquette Cast3M


fr en

Télécharger majw3d.eso

Retour à la liste

Numérotation des lignes : 

majw3d

  1. C MAJW3D    SOURCE    FD218221  24/02/07    21:15:19     11834          
  2.       subroutine majw3d(epspt60,epspt6,wplt06,wplt6,wpltx06,wpltx6,
  3.      # wpl3,vwpl33,vwpl33t,wplx3,vwplx33,vwplx33t,XE3D,NBNMAX3D,
  4.      # NBNB3D,IDIMB3D,Method_N,LcH,Method_H,err1)
  5.  
  6. c     mise a jour du tenseur des ouvertures: Sellier mai 2015
  7.  
  8. c     tables de dimension fixe pour resolution des sytemes lineaires 
  9.       implicit real*8 (a-h,o-z)
  10.       implicit integer (i-n)
  11.  
  12. c     variables externes
  13.       integer NB_HELM
  14.       real*8 LCH
  15.       logical Method_H,Method_N
  16.       real*8 epspt60(6),epspt6(6)
  17.       real*8 wplt06(6),wplt6(6),wpltx06(6),wpltx6(6)
  18.       real*8 wpl3(3),vwpl33(3,3),vwpl33t(3,3)
  19.       real*8 wplx3(3),vwplx33(3,3),vwplx33t(3,3)
  20.       integer NBNMAX3D,NBNB3D,idimb3d,err1
  21.       real*8 xe3d(3,NBNMAX3D)
  22.  
  23. c     variables locales
  24.       real*8 depspt33(3,3),depspt3(3)
  25.       real*8 vdepspt33(3,3),vdepspt33t(3,3)
  26.       real*8 long3(3),dwp6(6),dw6(6)
  27.       real*8 wplt33(3,3)
  28.       real*8 wpltx061(6),wpltx61(6)
  29.       real*8 wpltx33(3,3),wplt61(6),aux6(6),aux3(3)       
  30.  
  31. c     *** actualisation ouvertures de fissures actuelles  **************
  32.  
  33. c     base principale de l increment d ouverture
  34.       call bsfs3d(epspt6,epspt60,depspt3,
  35.      # vdepspt33,vdepspt33t)
  36. c     calcul des tailles dans les directions principales des increments
  37. c     choix de la methode pour le calcul de la taille de la zone de traction
  38. c     pour une raison de stabilité dans la methode  de calcul de la taille
  39. c     on verifie que la methode de calcul est la meme dans les trois directions
  40.       err1=0
  41. c      print*,'dans majw3d av tail3d',Method_N,Method_H,LcH
  42.       call tail3d(long3,vdepspt33,Method_N,Method_H,
  43.      # XE3D,NBNMAX3D,NBNB3D,IDIMB3D,LcH,err1)
  44. c      print*,'dans majw3d ap tail3d',Method_N,Method_H,LcH
  45. C     do i=1,3
  46. C          write(*,'(a9,i2,1x,a2,3(1x,e10.3),a4,e10.3)') 
  47. C      #  'Direction', i ,'(',(vdepspt33(j,i),j=1,3),'),L=',long3(i)
  48. C     end do
  49.       if(err1.eq.1) then
  50.           print*,'Pb de calcul de taille dans majw3d'
  51.           return
  52.       end if
  53.  
  54. c     mise a jour des ouvertures actuelles et maximales
  55.       call wfis3d(depspt3,long3,vdepspt33t,wplt06,wplt6,
  56.      # vwpl33,vwpl33t,wpl3,wpltx06,wpltx6,vwplx33,vwplx33t,wplx3)
  57.  
  58. c     ******************************************************************  
  59.  
  60.       return
  61.       end      
  62.  
  63.  
  64.  
  65.  

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