Esope source files

   1 : C LOADSF    SOURCE    AM        11/05/26    21:15:58     6982           
   2 :       subroutine LOADSF (da, force, q, yield)
   3 : C     Routine permettant le calcul de la valeur de la surface de charge
   4 : 
   5 :         IMPLICIT INTEGER(I-N)
   6 :         IMPLICIT REAL*8(A-H,O-Z)
   7 :         DIMENSION da(20), Xk9(4), q(5), force(5)
   8 : *
   9 :         a=da(8)
  10 :         b=da(9)
  11 :         c=da(10)
  12 :         d=da(11)
  13 :         e=da(12)
  14 :         f=da(13)
  15 : 
  16 :         gamma=1.
  17 :         if (force(1).le.0) then
  18 :         Xnum1=0
  19 :         Xnum2=0
  20 :         else
  21 :         Xnum1=((force(1))**c)*((gamma-force(1))**d)
  22 :         Xnum2=((force(1))**(c-e))*((gamma-force(1))**(d-f))
  23 :         endif
  24 :         Xk1 = force(2)/a
  25 :         Xk2 = force(3)*Xnum2/b
  26 :         Xk3 = force(4)/a
  27 :         Xk4 = force(5)*Xnum2/b
  28 :         Xk5 = q(1)*Xnum1
  29 :         Xk6 = q(2)*Xnum1
  30 :         Xk7 = q(3)*Xnum1
  31 :         Xk8 = q(4)*Xnum1
  32 : 
  33 :         Xk9(1)=Xk1-Xk5
  34 :         Xk9(2)=Xk2-Xk6
  35 :         Xk9(3)=Xk3-Xk7
  36 :         Xk9(4)=Xk4-Xk8
  37 : 
  37 :         call SCALPR(4,Xk9,Xk9,Xk10)
  38 :         yield=Xk10-q(5)**2*Xnum1**2
  39 : 
  40 :         return
  41 :         end
  42 :