Télécharger loadsf.eso

Retour à la liste

Numérotation des lignes :

  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.  
  38. call SCALPR(4,Xk9,Xk9,Xk10)
  39. yield=Xk10-q(5)**2*Xnum1**2
  40.  
  41. return
  42. end
  43.  
  44.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales