C Rutina 092
C---------------------------------------------------------DER_SOL_FUND1---
C     
C     SUBRUTINA DER_SOL_FUND1
C                                     
C     PROBLEMA DE ELASTICIDAD AXISIMETRICA ISOTROPA
C     CALCULO DE LA SOLUCION FUNDAMENTAL.
C
      SUBROUTINE DER_SOL_FUND1(X,Y,PMAT,Aab,Bab,Dab,nr,nz,COD)
C
	USE MSIMSLMD
C
      IMPLICIT NONE
C
	INCLUDE '_param.for'
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C     VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C X() : Coordenadas del punto de colocacion
      REAL*8 X(2)
C Y() : Coordenadas del punto de integracion
      REAL*8 Y(2)
C PMAT() : Propiedades del material
      REAL*8 PMAT(2)
C Aab,Bab,Dab : Nucleos en el calculo de tensiones en el dominio
C		   Dab(1,1)=Drrr, Bab(1,1)=Brrr, Aab(1)=Arrr 
C		   Dab(2,1)=Drzr, Bab(2,1)=Brzr, Aab(2)=Arzr 
C		   Dab(3,1)=Dzrr, Bab(3,1)=Bzrr, Aab(3)=Azrr 
C		   Dab(4,1)=Dzzr, Bab(4,1)=Bzzr, Aab(4)=Azzr 
C		   Dab(5,1)=Door, Bab(5,1)=Boor, Aab(5)=Aoor 
C		   Dab(1,2)=Drrz, Bab(1,1)=Brrz
C		   Dab(2,2)=Drzz, Bab(2,1)=Brzz
C		   Dab(3,2)=Dzrz, Bab(3,1)=Bzrz
C		   Dab(4,2)=Dzzz, Bab(4,1)=Bzzz
C		   Dab(5,2)=Dooz, Bab(5,1)=Booz
      REAL*8 Aab(5),Bab(5,2),Dab(5,2)
C Nr,Nz : Componentes del vector normal al elemento
      REAL*8 Nr,Nz                                     
C COD: 0 Regular
C      1 Singularidad en el origen
C      2 Singularidad en el extremo
      INTEGER COD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C     VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C Yr,Yz,XR,XZ : Coordenadas
      REAL*8 Yr,Yz,XR,XZ
C a,b,b2,c,c2,d,s2,k : Variables auxiliares
      REAL*8 a,b,b2,c,c2,d,s,s2,k,m2,c4,c6,Yr2,Yr4,XR2,XR4
C G,nu : Propiedades del material
      REAL*8 G,nu
C FCA,FCB,FCD : Factor comun que afecta a todos los terminos
      REAL*8 FCA,FCB,FCD
C Kk,Ek : Integral eliptica completa de primer y segundo orden
      REAL*8 Kk,Ek
C Ck,Ce,Ckr,Ckz,Cer,Cez : Coeficientes que afectan a Kk y Ek
      REAL*8 Ck,Ce,Ckr,Ckz,Cer,Cez
C PI : Numero Pi
      REAL*8 PI
C Tr,Tz : Componentes del vector tangente al elemento
      REAL*8 Tr,Tz
C I,J: Contadores
	INTEGER I,J
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C     VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C                                    
      IF(COD.NE.0)THEN
	  WRITE(F_W,09510) 1
        WRITE (*,09510) 1 
	  STOP
      END IF
C
      DO i= 1 , 5
        Aab(i)=0.d0
        DO j=1 , 2
          Bab(i,j)=0.d0
	    Dab(i,j)=0.d0
        END DO
      END DO      
C
      Yr=Y(1)
      Yz=Y(2)
      XR=X(1)
      XZ=X(2)
	Tr=-Nz
	Tz=Nr
      G=PMAT(1)/2.d0/(1.d0+PMAT(2))
      nu=PMAT(2)
      PI=DACOS(-1.d0)
C
      b=Yr-XR
      b2=b**2.d0
      c=Yz-XZ
      c2=c**2.d0
      s2=b**2.d0+c**2.d0
      s=DSQRT(s2)
      a=4.d0*XR*Yr+s2
      d=2.d0*XR*Yr+s2
      k=2.d0*DSQRT(XR*Yr)/DSQRT(a)
	m2=s2/a
C
C     INTEGRALES ELIPTICAS
C
      IF(XR.NE.0) THEN
xxx	  Kk=DELK(k*k)
	  Kk=DELRF(0.d0,m2,1.d0)
xxx	  Ek=DELE(k*k)
	  Ek=DELRF(0.d0,m2,1.d0)-k*k/3*DELRD(0.d0,m2,1.d0)
      ELSE
        Kk=PI/2.d0
        Ek=PI/2.d0
      END IF
C
C     Factor comun
C
      IF(XR.NE.0) THEN
        FCA=G/2.d0/PI/(1.d0-nu)/XR/DSQRT(a)
	  FCB=G/2.d0/PI/(1.d0-nu)/DSQRT(a)
	  FCD=1.d0/4.d0/PI/(1.d0-nu)/DSQRT(a)
 	ELSE
	  FCA=G/(4.d0*(1.d0-nu)*s*s2)
	  FCB=G/(4.d0*(1.d0-nu)*s*s2*s2)
	  FCD=1.d0/(8.d0*(1.d0-nu)*s*s2*s2)
      END IF
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C																		C
C				  NUCLEOS   Aabi                                        C
C																		C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C     Arrr
C
      IF(XR.NE.0) THEN
        Ckr=c*(XR/yr-2.d0*(3.d0-2.d0*nu)*(1.d0+s2/2.d0/XR/Yr))
        Ckz=2.d0*(2.d0*Yr+s2/XR)
        Cer=(5.d0-4.d0*nu)*c+2.d0*(3.d0-2.d0*nu)*c*(1.d0+s2/2.d0/XR/Yr)+
	@	  b*c/Yr+2.d0*XR*b*c/s2-4.d0*nu*XR*Yr*c/s2
        Cez=-(XR+7.d0*Yr+2.d0*s2/XR+b)
C
        Aab(1)=FCA*((Ckr*Tr+Ckz*Tz)*Kk + (Cer*Tr+Cez*Tz)*Ek)
C
	ELSE
C
        Aab(1)=FCA*(Yr*(-c*(1.d0+2.d0*nu)*Tr+2.d0*Yr*Tz))
C
      END IF
C
C     Arzr
C
      IF(XR.NE.0) THEN
        Ckr=(1.d0-2.d0*nu)*Yr+(XR*XR-c2)/Yr-
	@      2.d0*XR*(3.d0-2.d0*nu)*(1.d0+s2/2.d0/XR/Yr)
        Ckz=2.d0*nu*c
        Cer=2.d0*Yr*nu*(1.d0+2.d0*XR*b/s2)+
	@      2.d0*XR*(5.d0-4.d0*nu+c2/s2-b2/s2)+(3.d0-2.d0*nu)*s2/Yr +
	@      (c2-b2)/Yr
        Cez=-2.d0*nu*(2.d0*XR*Yr*c/s2+c)
C      
        Aab(2)=FCA/2.d0*((Ckr*Tr+Ckz*Tz)*Kk + (Cer*Tr+Cez*Tz)*Ek)
C
      END IF
C
C     Azrr
C
	Aab(3)=Aab(2)
C
C     Azzr 
C
      IF(XR.NE.0) THEN
        Ckr=-c/Yr
        Ckz=2.d0*nu
        Cer=-c/Yr*(2.d0*Yr*b/s2-1.d0)
        Cez=2.d0*nu*(2.d0*Yr*b/s2-1.d0)
C
	  Aab(4)=FCA*XR*((Ckr*Tr+Ckz*Tz)*Kk + (Cer*Tr+Cez*Tz)*Ek)
C
      ELSE
C
        Aab(4)=FCA*( 2.d0*Yr*(-c*Tr + 2.d0*nu*Yr*Tz) )
C
      END IF
C
C     Aoor 
C
      IF(XR.NE.0) THEN
C
        Ckr=(c**6*(3 - 2*Nu) + c**4*((9 - 6*Nu)*Yr**2 + (9 -5*Nu)*XR**2)
	@	 + c**2*((9 - 6*Nu)*Yr**4 - (-6 + Nu)*Yr**2*XR**2
	@	 + (9 - 4*Nu)*XR**4) - (Yr**2 - XR**2)*((-3 + 2*Nu)*Yr**4
	@	 - 2*Nu*Yr**2*XR**2 - (-3 + Nu)*XR**4))*c/Yr/XR/s2/a
	  Ckz=(-2*c**6 - 2*c**4*(3*Yr**2 + 2*XR**2)
	@	 - (Yr**2 - XR**2)**2*(2*Yr**2 + Nu*XR**2)
	@	 - c**2*(6*Yr**4 + Nu*Yr**2*XR**2 + (2 + Nu)*XR**4))/XR/s2/a
	  Cer=(c**8*(-3 + 2*Nu) + c**6*(4*(-3 + 2*Nu)*Yr**2
	@	  + (-12 + 7*Nu)*XR**2) + c**4*(6*(-3 + 2*Nu)*Yr**4
	@	  + 4*(-4 + Nu)*Yr**2*XR**2 + 9*(-2 + Nu)*XR**4)
	@	  + (Yr**2 - XR**2)*((-3 + 2*Nu)*Yr**6
	@	  + (5 - 8*Nu)*Yr**4*XR**2 - (5 + Nu)*Yr**2*XR**4
	@	  - (-3 + Nu)*XR**6) + c**2*(4*(-3 + 2*Nu)*Yr**6
	@	  + (4 - 13*Nu)*Yr**4*XR**2 - 4*(-1 + Nu)*Yr**2*XR**4
	@	  + (-12 + 5*Nu)*XR**6))*c/Yr/XR/s2/s2/a
	  Cez=(2*c**8 + 2*c**6*(4*Yr**2 + 3*XR**2)
	@	  + (Yr**2 - XR**2)**2*(2*Yr**4 + (-2 + Nu)*Yr**2*XR**2
	@	  + Nu*XR**4) + c**4*(12*Yr**4 + (6 + Nu)*Yr**2*XR**2
	@	  + (6 + Nu)*XR**4) + 2*c**2*(4*Yr**6 + (-3 + Nu)*Yr**4*XR**2
	@	  - 2*(1 + 3*Nu)*Yr**2*XR**4 + (1 + Nu)*XR**6))/XR/s2/s2/a
C
	  Aab(5)=FCA*((Ckr*Tr+Ckz*Tz)*Kk + (Cer*Tr+Cez*Tz)*Ek)
C
	ELSE
C
	  Aab(5)=Aab(1)
C
      END IF
C
C     Aorr,Aror,Aabz = 0
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C																		C
C				  NUCLEOS   Babi                                        C
C																		C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C     Brrr
C
      IF(XR.NE.0) THEN
        Ck=c/2.d0/XR*(1.d0/XR/Yr*(s2*(13.d0/2.d0-4.d0*nu)+
     @     d/2.d0/a*(b2-c2)-3.d0*b2)+7.d0-8.d0*nu+2.d0*d*b2/a/s2)
        Ce=c*s2/2.d0/XR/XR/Yr*(4.d0*nu-3.d0-3.d0*c2/s2)- 
	@     c*b2/2.d0/XR/XR/Yr+ 4.d0*Yr*c/a*(1.d0-2.d0*b2/s2)- 
     @     c/XR*(b2/s2+3.d0*c2/s2+7.d0-8.d0*nu)- 
     @     2.d0*Yr*c/s2*(1.d0+8.d0*Yr*XR*b2/a/s2)
C
        Bab(1,1)=FCB*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Bab(1,1)=FCB*( -c*Yr*(c**2.d0*(1.d0+2.d0*nu)
	@		  +2.d0*(2.d0+nu)*Yr**2.d0) )
C
      END IF
C
C     Brzr
C
      IF(XR.NE.0) THEN
        Ck=(13.d0/2.d0-4.d0*nu)*(1.d0+s2/2.d0/XR/Yr)- 3.d0*XR/2.d0/Yr- 
	@     3.d0*Yr/XR-
	@     (2.d0*XR*b+s2)/4.d0/XR/Yr/a*(2.d0*Yr*b/s-s)**2.d0+
	@     c2/2.d0/XR/Yr*(3.d0-d/a)+ b*d/XR/a*c2/s2
        Ce=s2/2.d0/XR/Yr*(4.d0*nu-3.d0-4.d0*XR*XR/a)- c2/XR/Yr+ 
	@     3.d0*(b2-c2)/2.d0/XR/Yr+ b*(12.d0*XR/a-1.d0/XR)+ 
	@     b2/s2*(Yr/XR-24.d0*XR*Yr/a+6.d0)- 
	@     c2/s2*(2.d0+Yr/XR-8.d0*XR*Yr/a)- 7.d0+	8.d0*nu +
	@     16.d0*XR*Yr*Yr*b/a/s2*(b2/s2-c2/s2)
C
        Bab(2,1)=FCB/2.d0*( Ck*Kk + Ce*Ek )
C
      END IF
C
C     Bzrr
C
	Bab(3,1)=Bab(2,1)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CC    Bzzr
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
      IF(XR.NE.0) THEN
        Ck=c/2.d0/Yr*(3.d0-1.d0/a*(2.d0*Yr*b/s-s)**2.d0)
        Ce=c/Yr*(2.d0*d/a-3.d0)+ 2.d0*b*c/s2*(3.d0-4.d0*d/a)+ 
	@     2.d0*Yr*c/s2*(-1.d0+4.d0*d/a*b2/s2)
C
        Bab(4,1)=FCB*( Ck*Kk + Ce*Ek )
C
      ELSE
C
        Bab(4,1)=FCB*( -2.d0*c**3.d0*Yr+4.d0*c*Yr**3.d0 )
C
      END IF
C
C     Boor
C
      IF(XR.NE.0) THEN
        Ck=c/XR*(s2/XR/Yr*(2.d0*nu-3.d0)+Yr/XR+2.d0*(2.d0*nu-3.d0))
        Ce=c*s2/XR/XR/Yr*(3.d0-2.d0*nu)- b*c/XR/XR+
	@  c/XR*(11.d0-8.d0*nu)- 2.d0*Yr*b*c/XR/s2- 4.d0*nu*Yr*c/s2
C
C	Los nucleos de Sladek estan mal y hace falta esta correccion
C
	  Ck=Ck-Nu*c/Yr
	  Ce=Ce+Nu*c*d/Yr/s2
C
        Bab(5,1)=FCB*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Bab(5,1)=Bab(1,1)
C
      END IF
C
C     Borr,Bror = 0
C
C     Brrz
C
      IF(XR.NE.0) THEN
        Ck=s2/XR/XR*(2.d0*nu-1.d0+d*c2/2.d0/a/s2)- 3.d0/2.d0*c2/XR/XR+ 
	@     d*b*c2/a/XR/s2+ 2.d0*(2.d0*nu-1.d0)*Yr/XR+ 1.d0
        Ce=s2/XR/XR*(1.d0-2.d0*nu+c2/s2)+ b/XR+ 
	@     Yr/XR*(3.d0-8.d0*nu+2.d0*c2/s2)+ c2/s2*(1.d0-8.d0*Yr*Yr/a)+ 
     @     2*Yr*b/s2*(1.d0-8.d0*XR*Yr*c2/a/s2)
C
        Bab(1,2)=FCB*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Bab(1,2)=FCB*( Yr**2.d0*(2.d0*c**2.d0*(-1.d0+nu)
	@		  +(1.d0+2.d0*nu)*Yr**2.d0) )
C
      END IF
C
C     Brzz
C
      IF(XR.NE.0) THEN
        Ck=c*( (c2-b2)/2.d0/XR/a- 2*Yr*b2/a/s2+ d*c2/a/XR/s2- 
	@     3.d0/2.d0/XR )
        Ce=c*b2/XR/s2*(2.d0*d/a-3.d0)+ 
	@     c/XR*(3.d0-(2.d0*d/a+1.d0)*c2/s2)+ 
     @     8.d0*Yr*c/s2*(d/a*b2/s2-2.d0*XR*Yr/a*c2/s2)
C
        Bab(2,2)=FCB/2.d0*( Ck*Kk + Ce*Ek )
C
      END IF
C
C     Bzrz
C
	Bab(3,2)=Bab(2,2)
C
C     Bzzz
C
      IF(XR.NE.0) THEN
        Ck=1.d0-c2/s2*(2.d0*Yr*b-s2)/a
        Ce=c2/s2*(3.d0-4.d0*d/a)- 1.d0+ 
	@     2.d0*Yr*b/s2*(1.d0+4.d0*d/a*c2/s2)
C
        Bab(4,2)=FCB*( Ck*Kk + Ce*Ek )
C
      ELSE
C
        Bab(4,2)=FCB*( 2.d0*Yr*(4.d0*c**2.d0*Yr+Yr**3.d0) )
C
      END IF
C
C     Booz
C
      IF(XR.NE.0) THEN
        Ck=s2/XR/XR*(1.d0-2.d0*nu) + c2/XR/XR+ 2.d0*nu+ 
	@     2.d0*Yr/XR*(1.d0-2.d0*nu)
        Ce=2.d0*nu/XR*b+ s2/XR/XR*(2.d0*nu-1.d0-c2/s2)+ 
	@     2.d0*Yr/XR*(3.d0*nu-2.d0-c2/s2)+ 4.d0*nu*Yr*b/s2
C
        Bab(5,2)=FCB*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Bab(5,2)=Bab(1,2)
C
      END IF
C
C     Borz,Broz = 0
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C																		C
C				  NUCLEOS   Dabi                                        C
C																		C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C     Drrr
C
      IF(XR.NE.0) THEN
        Ck=3.d0*b/2.d0/XR*(b/XR-1.d0)+ s2/XR/XR*(4.d0*nu-17.d0/4.d0)-
	@      1.d0/a*(b/s+s/2.d0/XR)**2.d0*(2.d0*Yr*b-s2) + 
	@     (1.d0-2.d0*nu)*(1.d0-4.d0*Yr/XR)
        Ce=2.d0*Yr*b/s2*(1.d0-2.d0*nu+8.d0*XR*Yr/a*b2/s2)+ 
	@     b2/s2*(1.d0-8.d0*XR*Yr/a)+ 2.d0*Yr/XR*(5.d0-7.d0*nu)+ 
	@     2.d0*b/XR*(2.d0-nu-2.d0*XR*Yr/a)- 
	@     2.d0*b*b2/XR/a-b2/XR/XR*(1.d0+4.d0*XR*XR/a)+ 
	@     2.d0*s2/XR/XR*(2.d0-2.d0*nu+XR*XR/a)
C
        Dab(1,1)=FCD*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Dab(1,1)=FCD*( 3.d0*Yr**4.d0 )
C
      END IF
C
C     Drzr
C
      IF(XR.NE.0) THEN
        Ck=c*( 1.d0/a*((c2-b2)/2.d0/XR-2.d0*Yr*b2/s2)+ 
	@     1.d0/XR*(2.d0*nu-5.d0/2.d0) )
        Ce=c/XR*(1.d0-2.d0*nu+3.d0*c2/s2+4.d0*d/a*b2/s2-2.d0*d/a)+ 
	@     2.d0*Yr*c/s2*(1.d0-2.d0*nu+4.d0*d/a*b2/s2)
C
        Dab(2,1)=FCD*( Ck*Kk + Ce*Ek )
C
      END IF
C
C     Dzrr
C
	Dab(3,1)=Dab(2,1)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CC    Dzzr
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
      IF(XR.NE.0) THEN
        Ck=c2/s2*(s2-2.d0*Yr*b)/a+ 2.d0*nu- 1.d0
        Ce=1.d0- 2.d0*nu+ (3.d0-4.d0*d/a)*c2/s2+ 
	@     2.d0*Yr*b/s2*(2.d0*nu-1.d0+4.d0*d/a*c2/s2)
C
        Dab(4,1)=FCD*( Ck*Kk + Ce*Ek )
C
      ELSE
C
        Dab(4,1)=FCD*( 2.d0*Yr**2.d0*(2.d0*c**2.d0*(1.d0+nu)
	@		  +(-1.d0+2.d0*nu)*Yr**2.d0) )
C
      END IF
C
C     Door
C
      IF(XR.NE.0) THEN
        Ck=4.d0*(1.d0-nu)*s2/XR/XR+ 8.d0*(1.d0-nu)*Yr/XR- Yr*Yr/XR/XR+ 
	@     2.d0*nu-1.d0
        Ce=4.d0*s2/XR/XR*(nu-1.d0)+ b/XR*(Yr/XR+2.d0*nu-1.d0)+ 
	@     14.d0*(nu-1.d0)*Yr/XR+ 2.d0*Yr*b/s2*(Yr/XR+2.d0*nu-1.d0)
C
        Dab(5,1)=FCD*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Dab(5,1)=Dab(1,1)
C
      END IF
C
C     Dorr,Dror = 0
C
C     Drrz
C
      IF(XR.NE.0) THEN
        Ck=c/2.d0*Yr/XR/XR*(3.d0-1.d0/a*(2.d0*XR*b/s+s)**2.d0)
        Ce=4.d0*Yr*c/a*(2.d0*b2/s2-1.d0)- Yr*c/XR/XR- 2.d0*Yr*b*c/XR/s2+ 
	@     2.d0*Yr*c/s2*(2.d0*nu-1.d0+8.d0*XR*Yr/a*b2/s2)
C
        Dab(1,2)=FCD*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Dab(1,2)=FCD*( c*Yr*(c**2.d0*(-2.d0+4.d0*nu)
	@		  +(1.d0+4.d0*nu)*Yr**2.d0) )
C
      END IF
C
C     Drzz
C
      IF(XR.NE.0) THEN
        Ck=Yr/XR*(2.d0*nu- 1.d0 -c2/s2*(2.d0*XR*b+s2)/a)
        Ce=Yr/XR*(1.d0- 2.d0*nu+ c2/s2*(4.d0*d/a-3.d0)+ 
	@     2.d0*XR*b/s2*(1.d0-2.d0*nu+4.d0*d/a*c2/s2))
C
        Dab(2,2)=FCD*( Ck*Kk + Ce*Ek )
C
      END IF
C
C     Dzrz
C
	Dab(3,2)=Dab(2,2)
C
C     Dzzz
C
      IF(XR.NE.0) THEN
        Ck=-2.d0*Yr*c/a*c2/s2
        Ce=2.d0*Yr*c/s2*(1.d0-2.d0*nu+4.d0*d/a*c2/s2)
C
        Dab(4,2)=FCD*( Ck*Kk + Ce*Ek )
C
      ELSE
C
        Dab(4,2)=FCD*( -2.d0*c*Yr*(2.d0*c**2.d0*(-2.d0+nu)
	@		  +(-1.d0+2.d0*nu)*Yr**2.d0) )
C
      END IF
C
C     Dooz
C
      IF(XR.NE.0) THEN
        Ck=-Yr*c/XR/XR
        Ce=Yr*c/XR/XR+ 2.d0*Yr*c/s2*(2.d0*nu-1.d0+Yr/XR)
C
        Dab(5,2)=FCD*( Ck*Kk + Ce*Ek )
C
	ELSE
C
        Dab(5,2)=Dab(1,2)
C
      END IF
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
CC    Dorz,Droz = 0
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C     FORMATOS DE ESCRITURA. 
C                  
09510 FORMAT(/,2x,'ERROR 09510 - DER_SOL_FUND1 - Error Programa (#',
     @       i1,'#)') 
C
C     FIN DE LA RUTINA 
C        
      RETURN
C
      END
C<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="es">
	<head>
            <title>e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.</title>

                <meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>
<meta name="description" content="Buscador de proyectos ETSI" />
<meta name="author" content="Daniel Callejo"/>
<meta name="author" content="Federico L&oacute;pez"/>
<meta name="author" content="Juli&aacute;n P&eacute;rez"/>
<meta name="author" content="Jos&eacute; Mar&iacute;a Vidal Vidal"/>
<base href="/bibing/proyectos/" /><link rel="shortcut icon" href="./favicon.ico" />
<link rel="stylesheet" href="./publico/css/estructura.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/bootstrap.css" type="text/css" media="screen"/>
<link rel="stylesheet" href="./publico/css/enlaces-capas.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/pestanas.css" type="text/css"  /> 
<link rel="stylesheet" href="./publico/css/formularios.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/lightbox.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/autocompleter.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/resultados.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/milista.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/explorador.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/ultimos.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/contactar.css" type="text/css"  />
<!-- <link rel="stylesheet" href="./css/print.css" type="text/css" media="print" /> -->
<script src="./publico/javascript/jquery-1.9.1.min.js" type="text/javascript"></script>
<script src="./publico/javascript/bootstrap.js" type="text/javascript"></script>

<script src="./publico/javascript/mootools-1.2.4-core.js" type="text/javascript"></script>
<script src="./publico/javascript/mootools-1-1.2.4.2-more.js" type="text/javascript"></script>
<script src="./publico/javascript/targetbl.js" type="text/javascript"></script>
<script src="./publico/javascript/mensajes.js" type="text/javascript"></script>
<script src="./publico/javascript/form.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Observer.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.Request.js" type="text/javascript"></script>
<script src="./publico/javascript/domready.js" type="text/javascript"></script>

<!-- Google tag (gtag.js) -->
<script async src="https://www.googletagmanager.com/gtag/js?id=G-8T1TSNC1QY"></script>
<script>
  window.dataLayer = window.dataLayer || [];
  function gtag(){dataLayer.push(arguments);}
  gtag('js', new Date());

  gtag('config', 'G-8T1TSNC1QY');
</script>



<script type="text/javascript">
function seleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=1
}
function deseleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=0
} 
</script>
	</head>
    <body style="min-width: 960px; background-color: #e2e2e2; height: 100%">	
		<a name="arriba"></a>
<div>
<div id="header">
    <table id="Tabla_01" height="105" border="0" cellpadding="0" cellspacing="0">
    <tr>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_01.jpg" width="3" height="105" alt="">
        </td>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_02.jpg" width="100" height="105" alt="">
        </td>
        <td width="100%" align="right"  style="min-width: 613px; background-repeat: no-repeat; background-position: top right; font-size: medium; font-weight: bold; color: white;">
            <table width="100%" height="105" border="0">
                <tr>
                    <td style=" vertical-align: top; width: 547px">
                        <!-- <a style="text-decoration: none" href="./"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a> -->
                        <a style="text-decoration: none" href="/bibing/proyectos/"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a>
                     </td>
                     <td  background="./publico/imagenes/cabecera/banner_new2_04.jpg" style="background-repeat: no-repeat; background-color: #c6d4e8;"><br><br><br>
                         &nbsp;
                        </td>
                    </tr>
             </table>
        </td>
    </tr>
    <tr style="border-top: 1px solid #ffffff;">
        <td colspan="3" id="menu_sup" style="font-size: small; text-align: right; background-color: #9c1f2f; vertical-align: bottom">
             
					<a href="/bibing/proyectos/" class="enlace " style=" color: #ffffff">Inicio&nbsp;&nbsp;|</a>
					<a href="/bibing/proyectos/acerca/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Acerca de&nbsp;&nbsp;|</a>
					
					<a href="/bibing/proyectos/contactar/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Contactar&nbsp;&nbsp;</a>
			 <!--       <a href="/bibing/proyectos/ayuda/" style=" color: #ffffff" class="enlace  ">&nbsp;&nbsp;Ayuda&nbsp;&nbsp;|</a> -->
					<!-- <form id="idioma" class=""  method="post" action="">
					<div class="inline" style="color: white">
					<input type="submit" class="tipo_letra " style=" color: #ffffff"name="idioma" value="Castellanossss" />&nbsp;&nbsp;
					</div>	
					</form> -->
			        </td>

    </tr>
    </table>
</div>
<br><br>
<div class="colmask rightmenu">
<div class="panel panel-default" style="margin:20px">
            <div class="panel-heading"><h4>ANÁLISIS DEL PROBLEMA TERMOELÁSTICO DE CONTACTO MEDIANTE EL MÉTODO DE LOS ELEMENTOS DE CONTORNO CON MALLAS CONFORMES Y NO CONFORMES</h4><h5 style="color:grey"><span class="glyphicon glyphicon-user""></span> : <b>MORALES SIERRA, DIÓGENES</b><br><span class="glyphicon glyphicon-book""></span>  :<b> Ingenier&iacute;a Industrial</b></h5></div><div class="panel-body"><div class="contenido_proyecto"><div class="titulo2_explorador">Contenido del proyecto:</div><div class="ruta"><a href="abreproy/4093">Directorio ra&iacute;z</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="abreproy/4093/direccion/RUTINAS%252F">RUTINAS</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="abreproy/4093/direccion/RUTINAS%252F17-2-2-dsf1.for%252F">17-2-2-dsf1.for</a></div><div class="listing"></div></div></div></div>
</div>
</div>

        
	</body>
 
 
<!-- Inicio Google Analytics -->
<!--
	<script type="text/javascript">
		var gaJsHost = (("https:" == document.location.protocol) ? "https://ssl." : "http://www.");
		document.write(unescape("%3Cscript src='" + gaJsHost + "google-analytics.com/ga.js' type='text/javascript'%3E%3C/script%3E"));
	</script>
	
	<script type="text/javascript">
		try {
		var pageTracker = _gat._getTracker("UA-7789620-1");
		pageTracker._trackPageview();
	} 
		catch(err) {}
	</script>
-->
<!-- Fin Google Analytics -->

</html>