Télécharger la plaquette Cast3M


fr en

Télécharger b2coq2.eso

Retour à la liste

Numérotation des lignes : 

b2coq2

  1. C B2COQ2    SOURCE    CHAT      05/01/12    21:31:09     5004
  2.       SUBROUTINE B2COQ2(XSTRS,NSTRS,xddl,nbddl,xel,nbnn,t,p,igau)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8 (A-H,O-Z)
  5.  
  6.       parameter(xzero=0.d0,demi=0.5d0,un=1.d0,deux=2.d0,trois=3.d0,
  7.      +          quatre=4.d0,six=6.d0)
  8.  
  9.       dimension xstrs(*),xddl(*),xel(3,*),t(*),p(*)
  10.       dimension ddlloc(6),z(3,3)
  11.  
  12.  
  13. -INC PPARAM
  14. -INC CCOPTIO
  15.  
  16.       D=SQRT((XEL(1,2)-XEL(1,1))**2+(XEL(2,2)-XEL(2,1))**2)
  17.       SP=(XEL(1,2)-XEL(1,1))/D
  18.       CP=(XEL(2,2)-XEL(2,1))/D
  19.  
  20.       z(1,1)=sp
  21.       z(1,2)=cp
  22.       z(1,3)=xzero
  23.       z(2,1)=-cp
  24.       z(2,2)=sp
  25.       z(2,3)=xzero
  26.       z(3,1)=xzero
  27.       z(3,2)=xzero
  28.       z(3,3)=un
  29.  
  30.       do 100 i=1,3
  31.          ddlloc(i  )=xzero
  32.          ddlloc(i+3)=xzero
  33.          do 101 j=1,3
  34.             ddlloc(i  )=ddlloc(i  )+z(i,j)*xddl(j  )
  35.             ddlloc(i+3)=ddlloc(i+3)+z(i,j)*xddl(j+3)
  36.  101     continue
  37.  100  continue
  38.  
  39.       s=t(igau)*demi+demi
  40.  
  41.       dudx = xstrs(1)
  42.       terme1=six*s*(s-un)
  43.       dwdx = terme1*(ddlloc(2)-ddlloc(5))/D
  44.       dwdx = dwdx + ddlloc(3)*(trois*s*s-quatre*s+un)
  45.      +            + ddlloc(6)*(trois*s*s-deux*s     )
  46.  
  47. c ... on modifie ep_ss dans tous les cas :
  48. c         PLAN GENE       PLAN CONT      PLAN DEFO
  49.       if(ifour.eq.-3.or.ifour.eq.-2.or.ifour.eq.-1.or.
  50. c            AXIS                FOURIER N° 0
  51.      &   ifour.eq. 0.or.(ifour.eq.1.and.nifour.eq.0)) then
  52.          xstrs(1)=xstrs(1)+demi*(dudx*dudx+dwdx*dwdx)
  53.       endif
  54.  
  55. c ... on modifie ep_theta_theta ou ep_zz en modes :
  56. c          AXIS        PLAN GENE            FOURIER N° 0
  57.       if(ifour.eq.0.or.ifour.eq.-3.or.(ifour.eq.1.and.nifour.eq.0)) then
  58.          xstrs(2) = xstrs(2) + demi*xstrs(2)*xstrs(2)
  59.       endif
  60.  
  61. c ... on ne modifie pas de cisaillement ep_s_theta en mode Fourier 0,
  62. c     car rien ne depend de theta, or ce terme contient les derivées par
  63. c     rapport à theta ...
  64.  
  65.       RETURN
  66.       END
  67.  
  68.  
  69.  
  70.  

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