Télécharger pmpb.procedur

Retour à la liste

Numérotation des lignes :

  1. * PMPB PROCEDUR PASCAL 20/11/13 21:15:13 10778
  2. DEBPROC PMPB TSIG*chpoint LIG*maillage grapl*logique
  3. nnc*entier nne*entier ;
  4. *************************************************************************
  5. * Analyse des contraintes sur le segment d appui LIG ( SEG2 )suivant
  6. * specification des codes
  7. * ENTREES : TSIG chpoint des contraintes projetees sur LIG
  8. * (TSIG = PROI LIG SIG )
  9. * LIG segment d'appui
  10. * GRAPL Logique si VRAI --> traces des decompositions
  11. * SORTIES :
  12. * TABLES contenant pour chacune des composantes des contraintes
  13. * caractéristiques sur les extrémites de segments d'appui
  14. *
  15. ************************************************************************
  16. P1 = lig point init ;
  17. P2 = lig point fina ;
  18. lig1 = chan 'POI1' lig ;
  19. Listcomp = extr Tsig comp ;
  20. nn = lig nbno ;
  21. idi = vale DIME ;
  22. si (ega idi 2);
  23. x y = coor lig ;
  24. z = x * 0. ;
  25. sinon ;
  26. x y z = coor lig ;
  27. z = z - ( coor 3 P1) ;
  28. finsi ;
  29. x = x - ( coor 1 P1);
  30. y = y - ( coor 2 P1);
  31. xabs = ((x * x) + ( y * y) + ( z * z)) ** .5 ;
  32. ep = maxi xabs ;
  33. SAUT LIGNE ;
  34. MESS ' ETAT ' NNE ' COUPE' NNC ': CONTRAINTES LINEARISEES ' ;
  35. SAUT LIGNE ;
  36. nbcomp = dime Listcomp ;
  37. compm = Table ;
  38. compb = Table ;
  39. compmb = Table ;
  40. MESS ' ----------------------------------------------'
  41. '----------' ;
  42. MESS ' ! MEMBRANE ! MEM+FLEX ! MEM+FLEX ! TOTALES '
  43. ' ! TOTALES !' ;
  44. MESS ' ! ! INT. ! EXT. ! INT. '
  45. ' ! EXT. !' ;
  46. TIRET=
  47. '-----------------------------------------------------------------';
  48. MESS TIRET;
  49. *
  50. REPETER BOU1 NBCOMP ;
  51. CCC = extr Listcomp &bou1 ;
  52. ev1 = evol roug chpo TSIG CCC lig ;
  53. rabs = extr ev1 absc ;
  54. * ep = maxi rabs ;
  55. *
  56. *------------------------ contrainte moyenne membrane
  57. *
  58. pm =((INTG ev1 ) ) / ep ;
  59. compm.&bou1 = manu chpo lig 1 SCAL pm ;
  60.  
  61. *
  62. *----------------------- calcul du moment de flexion
  63. *
  64. xl = rabs - ( prog nn * (ep / 2 )) ;
  65. z = manu chpo lig1 1 SCAL xl ;
  66. sigzz= ((exco TSIG CCC) * z ) nomc SCAL ;
  67. evm = evol bleu chpo sigzz SCAL lig ;
  68. mm = (INTG evm ) ;
  69. *----------------------------- contrainte flexion
  70. pb = mm * 6. / ( ep * ep) ;
  71. a =( pb * 2. ) / ep ;
  72. compb.&bou1 = ((xabs * a) - pb ) ;
  73. *
  74. compmb.&bou1 = compm.&bou1 + compb.&bou1 ;
  75. MEMBR = extr compm.&bou1 SCAL P1 ;
  76. MEFLRI =(extr compmb.&bou1 SCAL P1 ) ;
  77. MEFLRE =(extr compmb.&bou1 SCAL P2 ) ;
  78. FLEXI = (extr compb.&bou1 SCAL P1 ) ;
  79. FLEXE = (extr compb.&bou1 SCAL P2 ) ;
  80. TOTIR = EXTR TSIG CCC P1 ;
  81. TOTER = EXTR TSIG CCC P2 ;
  82.  
  83.  
  84. si GRAPL ;
  85. *
  86. 'TITRE' 'Seg_app' nnc 'Etat' nne CCC ' Membrane ' membr
  87. 'Flexion max ' (MEFLRI - membr) ;
  88. evpm = evol vert chpo compm.&bou1 SCAL lig ;
  89. evpmb = evol turq chpo compmb.&bou1 SCAL lig ;
  90. evsignl = ev1 - evpmb ;
  91. evt = ev1 et evpm et evpmb et evsignl ;
  92. *
  93. TUB = TABLE ;TIB = TABLE ;
  94. TUB.1 = 'MARQ ETOI'; TIB.1 = 'Totale ' ;
  95. TUB.2 = 'MARQ CROI'; TIB.2 = 'Pm ' ;
  96. TUB.3 = 'MARQ TRIA'; TIB.3 = 'Pm+Pb ' ;
  97. TUB.4 = 'MARQ CARR'; TIB.4 = 'NonLin ' ;
  98. TUB.'TITRE' = TIB ;
  99. *
  100. 'DESS' evt lege tub xbor 0. ep AXES ;
  101. finsi ;
  102.  
  103. AA = CHAINE FORMAT '(F9.2)' ' ! ' CCC ' ! 'MEMBR ' !' MEFLRI
  104. ' !' MEFLRE ' !' TOTIR ' !' TOTER ' !';
  105. MESS AA;
  106. MESS TIRET;
  107. SI ( EGA &BOU1 1) ;
  108. P_MEMB = PROG MEMBR ;
  109. P_FLEXI = PROG FLEXI ;
  110. P_MPFI = PROG MEFLRI ;
  111. P_MPFE = PROG MEFLRE;
  112. P_FLEXE = PROG FLEXE ;
  113. P_TOTI = PROG TOTIR ;
  114. P_TOTE = PROG TOTER ;
  115. SINON ;
  116. P_MEMB = P_MEMB et (PROG MEMBR ) ;
  117. P_FLEXI = P_FLEXI et (PROG FLEXI ) ;
  118. P_MPFI = P_MPFI et (PROG MEFLRI ) ;
  119. P_MPFE = P_MPFE et (PROG MEFLRE ) ;
  120. P_FLEXE = P_FLEXE et (PROG FLEXE ) ;
  121. P_TOTI = P_TOTI et (PROG TOTIR ) ;
  122. P_TOTE = P_TOTE et (PROG TOTER ) ;
  123. FINSI ;
  124.  
  125. FIN BOU1 ;
  126. *
  127. TABV = TABLE ;
  128. TABV . 1 = P_FLEXI ;
  129. TABV . 2 = P_MPFI ;
  130. TABV . 3 = P_TOTI ;
  131. *
  132. TABV . 4 = P_FLEXE ;
  133. TABV . 5 = P_MPFE ;
  134. TABV . 6 = P_TOTE ;
  135. *
  136. * Contraintes équivalentes - Critère de Von Mises en peau
  137. * -----------------------------------------------
  138. * - Critère de Tresca en peau
  139. * --------------------------------------------
  140. *
  141. tabana = table ;
  142. tabana.1 = P_MEMB ;
  143. tabana.2 = P_MPFI;
  144. tabana.3 = P_MPFE ;
  145. tabana.4 = P_TOTI ;
  146. tabana.5 = P_TOTE ;
  147. REPETER BTB 5 ;
  148. prcour = tabana.&btb ;
  149. S1 = EXTR prcour 1 ;S2 = EXTR prcour 2 ;S3 = EXTR prcour 3 ;
  150. S4 = EXTR prcour 4 ;S5= 0. ;S6 = 0. ;
  151. SI ( EGA NBCOMP 6 ) ;
  152. S5 = EXTR prcour 5 ;S6 = EXTR prcour 6 ;
  153. FINSI ;
  154. *
  155. TRMM VMMM = @RCCMTRV S1 S2 S3 S4 S5 S6 ;
  156. si ( ega &btb 1 ) ;
  157. bbb = prog TRMM ;
  158. VVV = prog VMMM ;
  159. sinon ;
  160. bbb = bbb et (prog TRMM ) ;
  161. VVV = VVV et(prog VMMM ) ;
  162. finsi ;
  163. 'FIN' BTB ;
  164. TRESMI = 'EXTR' BBB 1 ; MEMISI = 'EXTR' VVV 1 ;
  165. TREMFI = 'EXTR' BBB 2 ; BEMISI = 'EXTR' VVV 2 ;
  166. TREMFE = 'EXTR' BBB 3 ; BEMISE = 'EXTR' VVV 3 ;
  167. TTREI = 'EXTR' BBB 4 ; VMV1 = 'EXTR' VVV 4 ;
  168. TTREE = 'EXTR' BBB 5 ; VMV2 = 'EXTR' VVV 5 ;
  169. MESS TIRET;
  170. AA = CHAINE FORMAT '(F9.2)' ' ! ' TRESCA ' ! 'TRESMI ' !' TREMFI
  171. ' !' TREMFE ' !' TTREI ' !' TTREE ' !';
  172. MESS AA;
  173. MESS TIRET;
  174. AA = CHAINE FORMAT '(F9.2)' ' ! ' VMISES ' ! 'MEMISI ' !' BEMISI
  175. ' !'BEMISE ' !' VMV1 ' !' VMV2 ' !';
  176. MESS AA;
  177. MESS TIRET;
  178. 'FINPROC' TABV;
  179.  

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