Search Site | Contact Details | FAQ

ADAS Subroutine eiqip

       SUBROUTINE EIQIP(EI,EIJ,EM,Z,PHI,SC,WI,WJ,R,EIQ,FLAG)            
C
       IMPLICIT REAL*8 (A-H,O-Z)                                      
C
C-----------------------------------------------------------------------
C
C  ****************** FORTRAN77 SUBROUTINE: EIQIP **********************
C
C PURPOSE UNKNOWN
C
C NOTES: THIS ROUTINE IS NOT YET PROPERLY ANNOTATED
C
C UNIX-IDL PORT:
C
C VERSION: 1.1                          DATE: 16-1-96
C MODIFIED: TIM HAMMOND (TESSELLA SUPPORT SERVICES PLC)
C               - FIRST VERSION
C
C-----------------------------------------------------------------------
C
       FLAG=1.0                                                       
       EF=EI-EIJ                                                     
       T=EF/EI                                                      
       C=1.0                                                       
       IF(T-0.2)1,4,4                                             
    1  IF(Z-0.01)2,3,3                                           
    2  C=DSQRT(4.0*EF/EIJ)                                      
       GO TO 21                                                
    3  IF(EM-1.5)21,21,20                                     
   20     C=DABS(EIJ)/(EI*DSQRT(EF)+EF*DSQRT(EI))            
       C=Z*DSQRT(EM)*(1.1056/DSQRT(EIJ)-C)                  
       C=DEXP(3.142*C)                                     
   21     EI=1.25*EIJ                                     
       EF=0.25*EIJ                                       
       GO TO 8                                          
    4  IF(T-5.0)8,8,5                                  
    5  IF(Z-0.01)6,6,7                                
    6  C=DSQRT(-4.0*EI/EIJ)                          
       GO TO 23                                     
    7  IF(EM-1.5)23,23,22                          
   22     C=DABS(EIJ)/(EI*DSQRT(EF)+EF*DSQRT(EI)) 
       C=Z*DSQRT(EM)*(1.1056/DSQRT(-EIJ)-C)        
       C=DEXP(3.142*C)                              
   23     EI=-0.25*EIJ                               
       EF=-1.25*EIJ                                   
    8  TI=DSQRT(EI)                                    
       TF=DSQRT(EF)                                     
       TIF=DABS(EIJ)/(TI+TF)                             
       XI=Z*TIF/(TI*TF)                                   
       D=TIF*R                                             
       E=TI*TF                                              
       IF(EM-1.5)25,24,24                                    
   24  EM2=DSQRT(EM)                                          
          TIF=EM2*TIF                                          
          XI=-EM2*XI                                            
          D=EM2*D                                                
   25  CONTINUE                                                   
       T=Z+E*R                                                     
       T1=XIP(XI,D)/(T*T)                                           
       T2=4.0*PHI*E*SC*EM                                            
       IF(WI-WJ)10,10,9                                               
    9  T2=T2*WI/WJ                                                     
   10  P=T2*T1                                                          
       EIQW=8.0*C*PHI*YIP(XI,D)*EM                  
       IF(P-0.5)15,15,11                             
   11  R1=R                                           
   12  A=R1                                            
       FLAG=0.0                                         
       VA=0.5-P                                          
       R1=R1+R1                                           
       T=Z+E*R1                                            
       D=TIF*R1                                             
       P=T2*XIP(XI,D)/(T*T)                                  
       IF(P-0.5)13,14,12                                      
   13  B=R1                                                    
       VB=0.5-P                                                 
       CALL ZERO1(A,B,VA,VB,0.01D0,R1,XI,Z,E,TIF,T2)             
   14  T=Z+E*R1                                                   
       D=TIF*R1                                                    
       T1=XIP(XI,D)/(T*T)                                           
   15  EIQ=8.0*C*PHI*(YIP(XI,D)+0.5*(T*T-Z*Z)*T1)*EM                 
       EIQ=DMIN1(EIQ,EIQW)                                            
       RETURN                                                          
       END                                                               
      REAL*8              EI,          EIJ,         EIQ,         EM
      REAL*8              FLAG,        PHI,         R,           SC
      REAL*8              WI,          WJ,          Z
© Copyright 1995-2018 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk