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