Sources Esope | Cast3M

loadsf
C LOADSF    SOURCE    AM        11/05/26    21:15:58     6982
      subroutine LOADSF (da, force, q, yield)
C     Routine permettant le calcul de la valeur de la surface de charge
 
        IMPLICIT INTEGER(I-N)
        IMPLICIT REAL*8(A-H,O-Z)
        DIMENSION da(20), Xk9(4), q(5), force(5)
*
        a=da(8)
        b=da(9)
        c=da(10)
        d=da(11)
        e=da(12)
        f=da(13)
 
        gamma=1.
        if (force(1).le.0) then
        Xnum1=0
        Xnum2=0
        else
        Xnum1=((force(1))**c)*((gamma-force(1))**d)
        Xnum2=((force(1))**(c-e))*((gamma-force(1))**(d-f))
        endif
        Xk1 = force(2)/a
        Xk2 = force(3)*Xnum2/b
        Xk3 = force(4)/a
        Xk4 = force(5)*Xnum2/b
        Xk5 = q(1)*Xnum1
        Xk6 = q(2)*Xnum1
        Xk7 = q(3)*Xnum1
        Xk8 = q(4)*Xnum1
 
        Xk9(1)=Xk1-Xk5
        Xk9(2)=Xk2-Xk6
        Xk9(3)=Xk3-Xk7
        Xk9(4)=Xk4-Xk8
 
        call SCALPR(4,Xk9,Xk9,Xk10)
        yield=Xk10-q(5)**2*Xnum1**2
 
        return
        end