s8zderim3.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    s8zderim3   ../engine/source/elements/solid/solide8z/s8zderim3.F
!||--- called by ------------------------------------------------------
!||    s8zforc3    ../engine/source/elements/solid/solide8z/s8zforc3.f
!||====================================================================
      SUBROUTINE s8zderim3(
     1   OFFG,    OFF,     VOL,     NGL,
     2   KSI,     ETA,     ZETA,    WI,
     3   PXC1,    PXC2,    PXC3,    PXC4,
     4   PYC1,    PYC2,    PYC3,    PYC4,
     5   PZC1,    PZC2,    PZC3,    PZC4,
     6   PX1H1,   PX1H2,   PX1H3,   PX1H4,
     7   PX2H1,   PX2H2,   PX2H3,   PX2H4,
     8   PX3H1,   PX3H2,   PX3H3,   PX3H4,
     9   PX4H1,   PX4H2,   PX4H3,   PX4H4,
     A   PX1,     PX2,     PX3,     PX4,
     B   PX5,     PX6,     PX7,     PX8,
     C   PY1,     PY2,     PY3,     PY4,
     D   PY5,     PY6,     PY7,     PY8,
     E   PZ1,     PZ2,     PZ3,     PZ4,
     F   PZ5,     PZ6,     PZ7,     PZ8,
     G   PXY1,    PXY2,    PXY3,    PXY4,
     H   PXY5,    PXY6,    PXY7,    PXY8,
     I   PYX1,    PYX2,    PYX3,    PYX4,
     J   PYX5,    PYX6,    PYX7,    PYX8,
     K   PXZ1,    PXZ2,    PXZ3,    PXZ4,
     L   PXZ5,    PXZ6,    PXZ7,    PXZ8,
     M   PZX1,    PZX2,    PZX3,    PZX4,
     N   PZX5,    PZX6,    PZX7,    PZX8,
     O   PYZ1,    PYZ2,    PYZ3,    PYZ4,
     P   PYZ5,    PYZ6,    PYZ7,    PYZ8,
     Q   PZY1,    PZY2,    PZY3,    PZY4,
     R   PZY5,    PZY6,    PZY7,    PZY8,
     S   BXY1,    BXY2,    BXY3,    BXY4,
     T   BXY5,    BXY6,    BXY7,    BXY8,
     U   BYX1,    BYX2,    BYX3,    BYX4,
     V   BYX5,    BYX6,    BYX7,    BYX8,
     W   BXZ1,    BXZ2,    BXZ3,    BXZ4,
     X   BXZ5,    BXZ6,    BXZ7,    BXZ8,
     Y   BZX1,    BZX2,    BZX3,    BZX4,
     Z   BZX5,    BZX6,    BZX7,    BZX8,
     1   BYZ1,    BYZ2,    BYZ3,    BYZ4,
     2   BYZ5,    BYZ6,    BYZ7,    BYZ8,
     3   BZY1,    BZY2,    BZY3,    BZY4,
     4   BZY5,    BZY6,    BZY7,    BZY8,
     5   CJ1,     CJ2,     CJ3,     CJ4,
     6   CJ5,     CJ6,     CJ7,     CJ8,
     7   CJ9,     HX,      HY,      HZ,
     8   JAC1,    JAC2,    JAC3,    JAC4,
     9   JAC5,    JAC6,    JACI1,   JACI2,
     A   JACI3,   JACI4,   JACI5,   JACI6,
     B   JACI7,   JACI8,   JACI9,   SMAX,
     C   DELTAX,  NU,      ICP,     IDEG,
     D   VOLDP,   NEL,     MTN,     ISMSTR,
     E   JHBE)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr17_c.inc"
#include      "impl1_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: MTN
      INTEGER, INTENT(IN) :: ISMSTR
      INTEGER, INTENT(IN) :: JHBE
      INTEGER ICP,IDEG(*)
C     REAL
      my_real
     .   OFF(*),OFFG(*),VOL(*),KSI,ETA,ZETA,WI,
     .   PXC1(*), PXC2(*), PXC3(*), PXC4(*),  
     .   PYC1(*), PYC2(*), PYC3(*), PYC4(*),  
     .   PZC1(*), PZC2(*), PZC3(*), PZC4(*),  
     .   PX1H1(*), PX1H2(*), PX1H3(*),PX1H4(*),  
     .   PX2H1(*), PX2H2(*), PX2H3(*),PX2H4(*),  
     .   PX3H1(*), PX3H2(*), PX3H3(*),PX3H4(*),  
     .   PX4H1(*), PX4H2(*), PX4H3(*),PX4H4(*),  
     .   PX1(*), PX2(*), PX3(*), PX4(*),
     .   PX5(*), PX6(*), PX7(*), PX8(*),
     .   PY1(*), PY2(*), PY3(*), PY4(*),
     .   PY5(*), PY6(*), PY7(*), PY8(*),
     .   PZ1(*), PZ2(*), PZ3(*), PZ4(*), 
     .   PZ5(*), PZ6(*), PZ7(*), PZ8(*),
     .   PXY1(*),PXY2(*),PXY3(*),PXY4(*),
     .   PXY5(*),PXY6(*),PXY7(*),PXY8(*),
     .   PYX1(*),PYX2(*),PYX3(*),PYX4(*),
     .   PYX5(*),PYX6(*),PYX7(*),PYX8(*),
     .   PXZ1(*),PXZ2(*),PXZ3(*),PXZ4(*),
     .   PXZ5(*),PXZ6(*),PXZ7(*),PXZ8(*),
     .   PZX1(*),PZX2(*),PZX3(*),PZX4(*),
     .   PZX5(*),PZX6(*),PZX7(*),PZX8(*),
     .   PYZ1(*),PYZ2(*),PYZ3(*),PYZ4(*),
     .   PYZ5(*),PYZ6(*),PYZ7(*),PYZ8(*),
     .   PZY1(*),PZY2(*),PZY3(*),PZY4(*),
     .   PZY5(*),PZY6(*),PZY7(*),PZY8(*),
     .   BXY1(*),BXY2(*),BXY3(*),BXY4(*),
     .   BXY5(*),BXY6(*),BXY7(*),BXY8(*),
     .   BYX1(*),BYX2(*),BYX3(*),BYX4(*),
     .   BYX5(*),BYX6(*),BYX7(*),BYX8(*),
     .   BXZ1(*),BXZ2(*),BXZ3(*),BXZ4(*),
     .   BXZ5(*),BXZ6(*),BXZ7(*),BXZ8(*),
     .   BZX1(*),BZX2(*),BZX3(*),BZX4(*),
     .   BZX5(*),BZX6(*),BZX7(*),BZX8(*),
     .   BYZ1(*),BYZ2(*),BYZ3(*),BYZ4(*),
     .   BYZ5(*),BYZ6(*),BYZ7(*),BYZ8(*),
     .   BZY1(*),BZY2(*),BZY3(*),BZY4(*),
     .   BZY5(*),BZY6(*),BZY7(*),BZY8(*),
     .   CJ1(*),CJ2(*),CJ3(*),
     .   CJ4(*),CJ5(*),CJ6(*),
     .   CJ7(*),CJ8(*),CJ9(*),
     .   HX(MVSIZ,4),  HY(MVSIZ,4), HZ(MVSIZ,4),   
     .   jac1(*),jac2(*),jac3(*),
     .   jac4(*),jac5(*),jac6(*),
     .   jaci1(mvsiz), jaci2(mvsiz), jaci3(mvsiz),
     .   jaci4(mvsiz), jaci5(mvsiz), jaci6(mvsiz),
     .   jaci7(mvsiz), jaci8(mvsiz), jaci9(mvsiz),
     .   smax(*),deltax(*),nu(*)
      DOUBLE PRECISION 
     .   VOLDP(MVSIZ)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NGL(*), I, J ,ICOR,ep
C     REAL
C                                                                     12
      my_real
     .   NU1(MVSIZ),  
     .   D1X(MVSIZ) , D2X(MVSIZ) , D3X(MVSIZ) , D4X(MVSIZ) ,
     .   D1Y(MVSIZ) , D2Y(MVSIZ) , D3Y(MVSIZ) , D4Y(MVSIZ) ,
     .   D1Z(MVSIZ) , D2Z(MVSIZ) , D3Z(MVSIZ) , D4Z(MVSIZ) ,
     .   XG1(MVSIZ), XG2(MVSIZ), XG3(MVSIZ), XG4(MVSIZ),  
     .   YG1(MVSIZ), YG2(MVSIZ), YG3(MVSIZ), YG4(MVSIZ),   
     .   ZG1(MVSIZ), ZG2(MVSIZ), ZG3(MVSIZ), ZG4(MVSIZ),   
     .   F1,F2,F3,XS,XAS,YS,YAS,ZS,ZAS,CS,CAS ,BXHI,BYHI,BZHI
C-----------------------------------------------
C------repeating but no memory
      DO I=1,nel
       jac1(i)=cj1(i)+hx(i,3)*eta+(hx(i,2)+hx(i,4)*eta)*zeta
       jac2(i)=cj2(i)+hy(i,3)*eta+(hy(i,2)+hy(i,4)*eta)*zeta
       jac3(i)=cj3(i)+hz(i,3)*eta+(hz(i,2)+hz(i,4)*eta)*zeta
C
       jac4(i)=cj4(i)+hx(i,1)*zeta+(hx(i,3)+hx(i,4)*zeta)*ksi
       jac5(i)=cj5(i)+hy(i,1)*zeta+(hy(i,3)+hy(i,4)*zeta)*ksi
       jac6(i)=cj6(i)+hz(i,1)*zeta+(hz(i,3)+hz(i,4)*zeta)*ksi
      ENDDO
C------ICP= 0 : Bindeman&Belytschko assumed strain
C------ICP= 1 : fast assumed volumic strain for the laws(elastoplastic)P=f(rho,rho0)
C------ICP= 11 : EAD nominal assumed strain for the incmpressible laws(hyperelastic)P/=f(rho,rho0)
C
      DO i=1,nel
       nu1(i) = nu(i)/(one - nu(i))
       off(i) = offg(i)
      ENDDO
       f1 = eta*zeta
       f2 = ksi*zeta
       f3 = ksi*eta
      DO i=1,nel
       d1x(i)=jaci3(i)*eta +jaci2(i)*zeta
       d2x(i)=jaci1(i)*zeta +jaci3(i)*ksi
       d3x(i)=jaci2(i)*ksi +jaci1(i)*eta
       d4x(i)=jaci1(i)*f1+jaci2(i)*f2+jaci3(i)*f3
C
       d1y(i)=jaci6(i)*eta +jaci5(i)*zeta
       d2y(i)=jaci4(i)*zeta +jaci6(i)*ksi
       d3y(i)=jaci5(i)*ksi +jaci4(i)*eta
       d4y(i)=jaci4(i)*f1+jaci5(i)*f2+jaci6(i)*f3
C
       d1z(i)=jaci9(i)*eta +jaci8(i)*zeta
       d2z(i)=jaci7(i)*zeta +jaci9(i)*ksi
       d3z(i)=jaci8(i)*ksi +jaci7(i)*eta
       d4z(i)=jaci7(i)*f1+jaci8(i)*f2+jaci9(i)*f3
      ENDDO
C-------Node 1,7-----
      DO i=1,nel
       xg1(i) = px1h1(i)*d1x(i)
       xg2(i) = px1h2(i)*d2x(i)
       xg3(i) = px1h3(i)*d3x(i)
       xg4(i) = px1h4(i)*d4x(i)
       yg1(i) = px1h1(i)*d1y(i)
       yg2(i) = px1h2(i)*d2y(i)
       yg3(i) = px1h3(i)*d3y(i)
       yg4(i) = px1h4(i)*d4y(i)
       zg1(i) = px1h1(i)*d1z(i)
       zg2(i) = px1h2(i)*d2z(i)
       zg3(i) = px1h3(i)*d3z(i)
       zg4(i) = px1h4(i)*d4z(i)
      ENDDO
C
      DO i=1,nel
       xs =d1x(i)+d2x(i)+d3x(i)
       xas =d4x(i)+xg1(i)+xg2(i)+xg3(i)+xg4(i)-pxc1(i)
       px1(i)=xs-xas
       px7(i)=xs+xas
       ys =d1y(i)+d2y(i)+d3y(i)
       yas =d4y(i)+yg1(i)+yg2(i)+yg3(i)+yg4(i)-pyc1(i)
       py1(i)=ys-yas
       py7(i)=ys+yas
       zs =d1z(i)+d2z(i)+d3z(i)
       zas =d4z(i)+zg1(i)+zg2(i)+zg3(i)+zg4(i)-pzc1(i)
       pz1(i)=zs-zas
       pz7(i)=zs+zas
C
       ys =d1y(i)+d2y(i)
       yas =yg1(i)+yg2(i)-pyc1(i)
       pxy1(i)=ys-yas
       pxy7(i)=ys+yas
       xs =d1x(i)+d2x(i)
       xas =xg1(i)+xg2(i)-pxc1(i)
       pyx1(i)=xs-xas
       pyx7(i)=xs+xas
       zs =d1z(i)+d3z(i)
       zas =zg1(i)+zg3(i)-pzc1(i)
       pxz1(i)=zs-zas
       pxz7(i)=zs+zas
       xs =d1x(i)+d3x(i)
       xas =xg1(i)+xg3(i)-pxc1(i)
       pzx1(i)=xs-xas
       pzx7(i)=xs+xas
       zs =d2z(i)+d3z(i)
       zas =zg2(i)+zg3(i)-pzc1(i)
       pyz1(i)=zs-zas
       pyz7(i)=zs+zas
       ys =d2y(i)+d3y(i)
       yas =yg2(i)+yg3(i)-pyc1(i)
       pzy1(i)=ys-yas
       pzy7(i)=ys+yas
      ENDDO
C  
      IF (icp /= 1 .AND. icp /= 11) THEN
       DO i=1,nel
       cas =-nu(i)*(d4x(i)+xg1(i)+xg4(i))
       xas =-nu1(i)*xg3(i)+cas
       bxy1(i)=-xas
       bxy7(i)=xas
       xas =-nu1(i)*xg2(i)+cas
       bxz1(i)=-xas
       bxz7(i)=xas
       cas =-nu(i)*(d4y(i)+yg2(i)+yg4(i))
       yas =-nu1(i)*yg3(i)+cas
       byx1(i)=-yas
       byx7(i)=yas
       yas =-nu1(i)*yg1(i)+cas
       byz1(i)=-yas
       byz7(i)=yas
       cas =-nu(i)*(d4z(i)+zg3(i)+zg4(i))
       zas =-nu1(i)*zg2(i)+cas
       bzx1(i)=-zas
       bzx7(i)=zas
       zas =-nu1(i)*zg1(i)+cas
       bzy1(i)=-zas
       bzy7(i)=zas
       ENDDO
      ENDIF
C-------Node 2,8-----
      DO i=1,nel
       xg1(i) = px2h1(i)*d1x(i)
       xg2(i) = px2h2(i)*d2x(i)
       xg3(i) = px2h3(i)*d3x(i)
       xg4(i) = px2h4(i)*d4x(i)
       yg1(i) = px2h1(i)*d1y(i)
       yg2(i) = px2h2(i)*d2y(i)
       yg3(i) = px2h3(i)*d3y(i)
       yg4(i) = px2h4(i)*d4y(i)
       zg1(i) = px2h1(i)*d1z(i)
       zg2(i) = px2h2(i)*d2z(i)
       zg3(i) = px2h3(i)*d3z(i)
       zg4(i) = px2h4(i)*d4z(i)
      ENDDO
C
      DO i=1,nel
       xs =d1x(i)-d2x(i)-d3x(i)
       xas =-d4x(i)+xg1(i)+xg2(i)+xg3(i)+xg4(i)-pxc2(i)
       px2(i)=xs-xas
       px8(i)=xs+xas
       ys =d1y(i)-d2y(i)-d3y(i)
       yas =-d4y(i)+yg1(i)+yg2(i)+yg3(i)+yg4(i)-pyc2(i)
       py2(i)=ys-yas
       py8(i)=ys+yas
       zs =d1z(i)-d2z(i)-d3z(i)
       zas =-d4z(i)+zg1(i)+zg2(i)+zg3(i)+zg4(i)-pzc2(i)
       pz2(i)=zs-zas
       pz8(i)=zs+zas
C
       ys =d1y(i)-d2y(i)
       yas =yg1(i)+yg2(i)-pyc2(i)
       pxy2(i)=ys-yas
       pxy8(i)=ys+yas
       xs =d1x(i)-d2x(i)
       xas =xg1(i)+xg2(i)-pxc2(i)
       pyx2(i)=xs-xas
       pyx8(i)=xs+xas
       zs =d1z(i)-d3z(i)
       zas =zg1(i)+zg3(i)-pzc2(i)
       pxz2(i)=zs-zas
       pxz8(i)=zs+zas
       xs =d1x(i)-d3x(i)
       xas =xg1(i)+xg3(i)-pxc2(i)
       pzx2(i)=xs-xas
       pzx8(i)=xs+xas
       zs =-d2z(i)-d3z(i)
       zas =zg2(i)+zg3(i)-pzc2(i)
       pyz2(i)=zs-zas
       pyz8(i)=zs+zas
       ys =-d2y(i)-d3y(i)
       yas =yg2(i)+yg3(i)-pyc2(i)
       pzy2(i)=ys-yas
       pzy8(i)=ys+yas
      ENDDO
C
      IF (icp /= 1 .AND. icp /= 11) THEN
      DO i=1,nel
       cas =-nu(i)*(-d4x(i)+xg1(i)+xg4(i))
       xas =-nu1(i)*xg3(i)+cas
       bxy2(i)=-xas
       bxy8(i)=xas
       xas =-nu1(i)*xg2(i)+cas
       bxz2(i)=-xas
       bxz8(i)=xas
       cas =-nu(i)*(-d4y(i)+yg2(i)+yg4(i))
       yas =-nu1(i)*yg3(i)+cas
       byx2(i)=-yas
       byx8(i)=yas
       yas =-nu1(i)*yg1(i)+cas
       byz2(i)=-yas
       byz8(i)=yas
       cas =-nu(i)*(-d4z(i)+zg3(i)+zg4(i))
       zas =-nu1(i)*zg2(i)+cas
       bzx2(i)=-zas
       bzx8(i)=zas
       zas =-nu1(i)*zg1(i)+cas
       bzy2(i)=-zas
       bzy8(i)=zas
      ENDDO
      ENDIF
C
C-------Node 3,5-----
      DO i=1,nel
       xg1(i) = px3h1(i)*d1x(i)
       xg2(i) = px3h2(i)*d2x(i)
       xg3(i) = px3h3(i)*d3x(i)
       xg4(i) = px3h4(i)*d4x(i)
       yg1(i) = px3h1(i)*d1y(i)
       yg2(i) = px3h2(i)*d2y(i)
       yg3(i) = px3h3(i)*d3y(i)
       yg4(i) = px3h4(i)*d4y(i)
       zg1(i) = px3h1(i)*d1z(i)
       zg2(i) = px3h2(i)*d2z(i)
       zg3(i) = px3h3(i)*d3z(i)
       zg4(i) = px3h4(i)*d4z(i)
      ENDDO
C
      DO i=1,nel
       xs =-d1x(i)-d2x(i)+d3x(i)
       xas =d4x(i)+xg1(i)+xg2(i)+xg3(i)+xg4(i)-pxc3(i)
       px3(i)=xs-xas
       px5(i)=xs+xas
       ys =-d1y(i)-d2y(i)+d3y(i)
       yas =d4y(i)+yg1(i)+yg2(i)+yg3(i)+yg4(i)-pyc3(i)
       py3(i)=ys-yas
       py5(i)=ys+yas
       zs =-d1z(i)-d2z(i)+d3z(i)
       zas =d4z(i)+zg1(i)+zg2(i)+zg3(i)+zg4(i)-pzc3(i)
       pz3(i)=zs-zas
       pz5(i)=zs+zas
C
       ys =-d1y(i)-d2y(i)
       yas =yg1(i)+yg2(i)-pyc3(i)
       pxy3(i)=ys-yas
       pxy5(i)=ys+yas
       xs =-d1x(i)-d2x(i)
       xas =xg1(i)+xg2(i)-pxc3(i)
       pyx3(i)=xs-xas
       pyx5(i)=xs+xas
       zs =-d1z(i)+d3z(i)
       zas =zg1(i)+zg3(i)-pzc3(i)
       pxz3(i)=zs-zas
       pxz5(i)=zs+zas
       xs =-d1x(i)+d3x(i)
       xas =xg1(i)+xg3(i)-pxc3(i)
       pzx3(i)=xs-xas
       pzx5(i)=xs+xas
       zs =-d2z(i)+d3z(i)
       zas =zg2(i)+zg3(i)-pzc3(i)
       pyz3(i)=zs-zas
       pyz5(i)=zs+zas
       ys =-d2y(i)+d3y(i)
       yas =yg2(i)+yg3(i)-pyc3(i)
       pzy3(i)=ys-yas
       pzy5(i)=ys+yas
      ENDDO
C  
      IF (icp /= 1 .AND. icp /= 11) THEN
      DO i=1,nel
       cas =-nu(i)*(d4x(i)+xg1(i)+xg4(i))
       xas =-nu1(i)*xg3(i)+cas
       bxy3(i)=-xas
       bxy5(i)=xas
       xas =-nu1(i)*xg2(i)+cas
       bxz3(i)=-xas
       bxz5(i)=xas
       cas =-nu(i)*(d4y(i)+yg2(i)+yg4(i))
       yas =-nu1(i)*yg3(i)+cas
       byx3(i)=-yas
       byx5(i)=yas
       yas =-nu1(i)*yg1(i)+cas
       byz3(i)=-yas
       byz5(i)=yas
       cas =-nu(i)*(d4z(i)+zg3(i)+zg4(i))
       zas =-nu1(i)*zg2(i)+cas
       bzx3(i)=-zas
       bzx5(i)=zas
       zas =-nu1(i)*zg1(i)+cas
       bzy3(i)=-zas
       bzy5(i)=zas
      ENDDO
      ENDIF
C-------Node 4,6-----
      DO i=1,nel
       xg1(i) = px4h1(i)*d1x(i)
       xg2(i) = px4h2(i)*d2x(i)
       xg3(i) = px4h3(i)*d3x(i)
       xg4(i) = px4h4(i)*d4x(i)
       yg1(i) = px4h1(i)*d1y(i)
       yg2(i) = px4h2(i)*d2y(i)
       yg3(i) = px4h3(i)*d3y(i)
       yg4(i) = px4h4(i)*d4y(i)
       zg1(i) = px4h1(i)*d1z(i)
       zg2(i) = px4h2(i)*d2z(i)
       zg3(i) = px4h3(i)*d3z(i)
       zg4(i) = px4h4(i)*d4z(i)
      ENDDO
C
      DO i=1,nel
       xs =-d1x(i)+d2x(i)-d3x(i)
       xas =-d4x(i)+xg1(i)+xg2(i)+xg3(i)+xg4(i)-pxc4(i)
       px4(i)=xs-xas
       px6(i)=xs+xas
       ys =-d1y(i)+d2y(i)-d3y(i)
       yas =-d4y(i)+yg1(i)+yg2(i)+yg3(i)+yg4(i)-pyc4(i)
       py4(i)=ys-yas
       py6(i)=ys+yas
       zs =-d1z(i)+d2z(i)-d3z(i)
       zas =-d4z(i)+zg1(i)+zg2(i)+zg3(i)+zg4(i)-pzc4(i)
       pz4(i)=zs-zas
       pz6(i)=zs+zas
C
       ys =-d1y(i)+d2y(i)
       yas =yg1(i)+yg2(i)-pyc4(i)
       pxy4(i)=ys-yas
       pxy6(i)=ys+yas
       xs =-d1x(i)+d2x(i)
       xas =xg1(i)+xg2(i)-pxc4(i)
       pyx4(i)=xs-xas
       pyx6(i)=xs+xas
       zs =-d1z(i)-d3z(i)
       zas =zg1(i)+zg3(i)-pzc4(i)
       pxz4(i)=zs-zas
       pxz6(i)=zs+zas
       xs =-d1x(i)-d3x(i)
       xas =xg1(i)+xg3(i)-pxc4(i)
       pzx4(i)=xs-xas
       pzx6(i)=xs+xas
       zs =d2z(i)-d3z(i)
       zas =zg2(i)+zg3(i)-pzc4(i)
       pyz4(i)=zs-zas
       pyz6(i)=zs+zas
       ys =d2y(i)-d3y(i)
       yas =yg2(i)+yg3(i)-pyc4(i)
       pzy4(i)=ys-yas
       pzy6(i)=ys+yas
      ENDDO
C
      IF (icp /= 1 .AND. icp /= 11) THEN
       DO i=1,nel
        cas =-nu(i)*(-d4x(i)+xg1(i)+xg4(i))
        xas =-nu1(i)*xg3(i)+cas
        bxy4(i)=-xas
        bxy6(i)=xas
        xas =-nu1(i)*xg2(i)+cas
        bxz4(i)=-xas
        bxz6(i)=xas
        cas =-nu(i)*(-d4y(i)+yg2(i)+yg4(i))
        yas =-nu1(i)*yg3(i)+cas
        byx4(i)=-yas
        byx6(i)=yas
        yas =-nu1(i)*yg1(i)+cas
        byz4(i)=-yas
        byz6(i)=yas
        cas =-nu(i)*(-d4z(i)+zg3(i)+zg4(i))
        zas =-nu1(i)*zg2(i)+cas
        bzx4(i)=-zas
        bzx6(i)=zas
        zas =-nu1(i)*zg1(i)+cas
        bzy4(i)=-zas
        bzy6(i)=zas
       ENDDO
C
       DO i=1,nel
        cs =-nu(i)*d1x(i)
        xs =-nu1(i)*d3x(i)
        bxy1(i)=bxy1(i)+cs+xs
        bxy2(i)=bxy2(i)+cs-xs
        bxy3(i)=bxy3(i)-cs+xs
        bxy4(i)=bxy4(i)-cs-xs
        bxy5(i)=bxy5(i)-cs+xs
        bxy6(i)=bxy6(i)-cs-xs
        bxy7(i)=bxy7(i)+cs+xs
        bxy8(i)=bxy8(i)+cs-xs
        xs =-nu1(i)*d2x(i)
        bxz1(i)=bxz1(i)+cs+xs
        bxz2(i)=bxz2(i)+cs-xs
        bxz3(i)=bxz3(i)-cs-xs
        bxz4(i)=bxz4(i)-cs+xs
        bxz5(i)=bxz5(i)-cs-xs
        bxz6(i)=bxz6(i)-cs+xs
        bxz7(i)=bxz7(i)+cs+xs
        bxz8(i)=bxz8(i)+cs-xs
       ENDDO
       DO i=1,nel
        cs =-nu(i)*d2y(i)
        ys =-nu1(i)*d3y(i)
        byx1(i)=byx1(i)+cs+ys
        byx2(i)=byx2(i)-cs-ys
        byx3(i)=byx3(i)-cs+ys
        byx4(i)=byx4(i)+cs-ys
        byx5(i)=byx5(i)-cs+ys
        byx6(i)=byx6(i)+cs-ys
        byx7(i)=byx7(i)+cs+ys
        byx8(i)=byx8(i)-cs-ys
        ys =-nu1(i)*d1y(i)
        byz1(i)=byz1(i)+cs+ys
        byz2(i)=byz2(i)-cs+ys
        byz3(i)=byz3(i)-cs-ys
        byz4(i)=byz4(i)+cs-ys
        byz5(i)=byz5(i)-cs-ys
        byz6(i)=byz6(i)+cs-ys
        byz7(i)=byz7(i)+cs+ys
        byz8(i)=byz8(i)-cs+ys
       ENDDO
       DO i=1,nel
        cs =-nu(i)*d3z(i)
        zs =-nu1(i)*d2z(i)
        bzx1(i)=bzx1(i)+cs+zs
        bzx2(i)=bzx2(i)-cs-zs
        bzx3(i)=bzx3(i)+cs-zs
        bzx4(i)=bzx4(i)-cs+zs
        bzx5(i)=bzx5(i)+cs-zs
        bzx6(i)=bzx6(i)-cs+zs
        bzx7(i)=bzx7(i)+cs+zs
        bzx8(i)=bzx8(i)-cs-zs
        zs =-nu1(i)*d1z(i)
        bzy1(i)=bzy1(i)+cs+zs
        bzy2(i)=bzy2(i)-cs+zs
        bzy3(i)=bzy3(i)+cs-zs
        bzy4(i)=bzy4(i)-cs-zs
        bzy5(i)=bzy5(i)+cs-zs
        bzy6(i)=bzy6(i)-cs-zs
        bzy7(i)=bzy7(i)+cs+zs
        bzy8(i)=bzy8(i)-cs+zs
       ENDDO
      ENDIF
      DO i=1,nel
        off(i) = offg(i)
        vol(i) = voldp(i)
      ENDDO
      IF (idts6==0) THEN
       DO i=1,nel
        deltax(i)=min(deltax(i),hundred28*vol(i)*smax(i)/wi)
       ENDDO
      ELSE
      DO i=1,nel
       IF (ideg(i)==0) deltax(i)=min(deltax(i),hundred28*vol(i)*smax(i)/wi)
      ENDDO
      END IF
C
      IF ((impl_s >0.AND.jhbe /= 14.AND.jhbe < 20).OR.
     .    ((ismstr==10.OR.ismstr==12).AND.
     .    (icp==1.OR.icp==11.OR.mtn==99.OR.mtn==1).AND.jhbe==14)) THEN
       DO i=1,nel
        pxy1(i)=py1(i)
        pxy2(i)=py2(i)
        pxy3(i)=py3(i)
        pxy4(i)=py4(i)
        pxy5(i)=py5(i)
        pxy6(i)=py6(i)
        pxy7(i)=py7(i)
        pxy8(i)=py8(i)
C
        pyx1(i)=px1(i)
        pyx2(i)=px2(i)
        pyx3(i)=px3(i)
        pyx4(i)=px4(i)
        pyx5(i)=px5(i)
        pyx6(i)=px6(i)
        pyx7(i)=px7(i)
        pyx8(i)=px8(i)
       ENDDO
       DO i=1,nel
        pyz1(i)=pz1(i)
        pyz2(i)=pz2(i)
        pyz3(i)=pz3(i)
        pyz4(i)=pz4(i)
        pyz5(i)=pz5(i)
        pyz6(i)=pz6(i)
        pyz7(i)=pz7(i)
        pyz8(i)=pz8(i)
C
        pzy1(i)=py1(i)
        pzy2(i)=py2(i)
        pzy3(i)=py3(i)
        pzy4(i)=py4(i)
        pzy5(i)=py5(i)
        pzy6(i)=py6(i)
        pzy7(i)=py7(i)
        pzy8(i)=py8(i)
       ENDDO
       DO i=1,nel
        pxz1(i)=pz1(i)
        pxz2(i)=pz2(i)
        pxz3(i)=pz3(i)
        pxz4(i)=pz4(i)
        pxz5(i)=pz5(i)
        pxz6(i)=pz6(i)
        pxz7(i)=pz7(i)
        pxz8(i)=pz8(i)
C
        pzx1(i)=px1(i)
        pzx2(i)=px2(i)
        pzx3(i)=px3(i)
        pzx4(i)=px4(i)
        pzx5(i)=px5(i)
        pzx6(i)=px6(i)
        pzx7(i)=px7(i)
        pzx8(i)=px8(i)
       ENDDO
      ENDIF 
      IF ((ismstr==10.OR.ismstr==12).AND.mtn==1.AND.jhbe==14.AND.
     .    (icp/=1.AND.icp/=11)) THEN
       DO i=1,nel
        bxy1(i)=zero
        bxz1(i)=zero
        byx1(i)=zero
        byz1(i)=zero
        bzx1(i)=zero
        bzy1(i)=zero
        bxy2(i)=zero
        bxz2(i)=zero
        byx2(i)=zero
        byz2(i)=zero
        bzx2(i)=zero
        bzy2(i)=zero
        bxy3(i)=zero
        bxz3(i)=zero
        byx3(i)=zero
        byz3(i)=zero
        bzx3(i)=zero
        bzy3(i)=zero
        bxy4(i)=zero
        bxz4(i)=zero
        byx4(i)=zero
        byz4(i)=zero
        bzx4(i)=zero
        bzy4(i)=zero
        bxy5(i)=zero
        bxz5(i)=zero
        byx5(i)=zero
        byz5(i)=zero
        bzx5(i)=zero
        bzy5(i)=zero
        bxy6(i)=zero
        bxz6(i)=zero
        byx6(i)=zero
        byz6(i)=zero
        bzx6(i)=zero
        bzy6(i)=zero
        bxy7(i)=zero
        bxz7(i)=zero
        byx7(i)=zero
        byz7(i)=zero
        bzx7(i)=zero
        bzy7(i)=zero
        bxy8(i)=zero
        bxz8(i)=zero
        byx8(i)=zero
        byz8(i)=zero
        bzx8(i)=zero
        bzy8(i)=zero
       ENDDO
      ENDIF 
C-----------for the imcompressible laws
C-------|BXX BYX BZX |
C-------|BXY BYY BZY |  same than H8A
C-------|BXZ BYZ BZZ |
      IF (icp==11) THEN
       DO i=1,nel
        bxhi=third*(px1(i)-pxc1(i))
        byhi=third*(py1(i)-pyc1(i))
        bzhi=third*(pz1(i)-pzc1(i))
        px1(i)=pxc1(i)+two*bxhi
        py1(i)=pyc1(i)+two*byhi
        pz1(i)=pzc1(i)+two*bzhi
        bxy1(i)=-bxhi
        bxz1(i)=-bxhi
        byx1(i)=-byhi
        byz1(i)=-byhi
        bzx1(i)=-bzhi
        bzy1(i)=-bzhi
       ENDDO
       DO i=1,nel
        bxhi=third*(px2(i)-pxc2(i))
        byhi=third*(py2(i)-pyc2(i))
        bzhi=third*(pz2(i)-pzc2(i))
        px2(i)=pxc2(i)+two*bxhi
        py2(i)=pyc2(i)+two*byhi
        pz2(i)=pzc2(i)+two*bzhi
        bxy2(i)=-bxhi
        bxz2(i)=-bxhi
        byx2(i)=-byhi
        byz2(i)=-byhi
        bzx2(i)=-bzhi
        bzy2(i)=-bzhi
       ENDDO
       DO i=1,nel
        bxhi=third*(px3(i)-pxc3(i))
        byhi=third*(py3(i)-pyc3(i))
        bzhi=third*(pz3(i)-pzc3(i))
        px3(i)=pxc3(i)+two*bxhi
        py3(i)=pyc3(i)+two*byhi
        pz3(i)=pzc3(i)+two*bzhi
        bxy3(i)=-bxhi
        bxz3(i)=-bxhi
        byx3(i)=-byhi
        byz3(i)=-byhi
        bzx3(i)=-bzhi
        bzy3(i)=-bzhi
       ENDDO
       DO i=1,nel
        bxhi=third*(px4(i)-pxc4(i))
        byhi=third*(py4(i)-pyc4(i))
        bzhi=third*(pz4(i)-pzc4(i))
        px4(i)=pxc4(i)+two*bxhi
        py4(i)=pyc4(i)+two*byhi
        pz4(i)=pzc4(i)+two*bzhi
        bxy4(i)=-bxhi
        bxz4(i)=-bxhi
        byx4(i)=-byhi
        byz4(i)=-byhi
        bzx4(i)=-bzhi
        bzy4(i)=-bzhi
       ENDDO
       DO i=1,nel
        bxhi=third*(px5(i)+pxc3(i))
        byhi=third*(py5(i)+pyc3(i))
        bzhi=third*(pz5(i)+pzc3(i))
        px5(i)=-pxc3(i)+two*bxhi
        py5(i)=-pyc3(i)+two*byhi
        pz5(i)=-pzc3(i)+two*bzhi
        bxy5(i)=-bxhi
        bxz5(i)=-bxhi
        byx5(i)=-byhi
        byz5(i)=-byhi
        bzx5(i)=-bzhi
        bzy5(i)=-bzhi
       ENDDO
       DO i=1,nel
        bxhi=third*(px6(i)+pxc4(i))
        byhi=third*(py6(i)+pyc4(i))
        bzhi=third*(pz6(i)+pzc4(i))
        px6(i)=-pxc4(i)+two*bxhi
        py6(i)=-pyc4(i)+two*byhi
        pz6(i)=-pzc4(i)+two*bzhi
        bxy6(i)=-bxhi
        bxz6(i)=-bxhi
        byx6(i)=-byhi
        byz6(i)=-byhi
        bzx6(i)=-bzhi
        bzy6(i)=-bzhi
       ENDDO
       DO i=1,nel
        bxhi=third*(px7(i)+pxc1(i))
        byhi=third*(py7(i)+pyc1(i))
        bzhi=third*(pz7(i)+pzc1(i))
        px7(i)=-pxc1(i)+two*bxhi
        py7(i)=-pyc1(i)+two*byhi
        pz7(i)=-pzc1(i)+two*bzhi
        bxy7(i)=-bxhi
        bxz7(i)=-bxhi
        byx7(i)=-byhi
        byz7(i)=-byhi
        bzx7(i)=-bzhi
        bzy7(i)=-bzhi
       ENDDO
       DO i=1,nel
        bxhi=third*(px8(i)+pxc2(i))
        byhi=third*(py8(i)+pyc2(i))
        bzhi=third*(pz8(i)+pzc2(i))
        px8(i)=-pxc2(i)+two*bxhi
        py8(i)=-pyc2(i)+two*byhi
        pz8(i)=-pzc2(i)+two*bzhi
        bxy8(i)=-bxhi
        bxz8(i)=-bxhi
        byx8(i)=-byhi
        byz8(i)=-byhi
        bzx8(i)=-bzhi
        bzy8(i)=-bzhi
       ENDDO
      ENDIF
C      
      RETURN
      END