$$$$ MRCFRAM1 * MRCFRAM1 PROCEDUR BP208322 20/05/05 21:15:02 10601 DEBPROC MRCFRAM1 TYP_CAL2*'MOT ' TAB_IN2*'TABLE' TOL_2*'FLOTTANT'; * * -------------------------------------------------------- * Procedure FRCFRAM1 * Procedure appelee par MRCFRAME pour le calcul de la marge sismique * d'un element frame tyep POUTEAU COURT * * Developpers: * Alberto FRAU * DEN/DANS/DM2S/SEMT/EMSI * Nicolas ILE * DEN/DANS/DM2S/SEMT/EMSI * Giulia DI STEFANO * * * LAM1 = MRCFRAM1 MOT1 TAB1 TOL1; * * * INPUT: * * MOT1 [MOT] afin de specifier le type de calcul: * 'ELLIPSOIDE' calcul des marges avec l'enveloppe * ellipsiodale * 'RECTANGLE' calcul des marges avec l'enveloppe * prismatique * * TAB1 type TABLE pour definir les donnees suivantes: * TAB1.'BASE_Y1' dimension section direction Y en m [FLOTTANT] * TAB1.'BASE_Z1' dimension section direction Z en m [FLOTTANT] * TAB1.'LONG_Y1' longuer ly [FLOTTANT] * TAB1.'LONG_Z1' longuer lz [FLOTTANT] * TAB1.'ENROBA1' enrobage en m [FLOTTANT] * TAB1.'ACADRY1' aire de l'armature transversale * direction y en m2 [FLOTTANT] * TAB1.'ACADRZ1' aire de l'armature transversale * direction z en m2 [FLOTTANT] * TAB1.'S_CADR1' espacement des cadres en * m [FLOTTANT] * TAB1.'LI_PHI1' LISTREEL des diam barres d'acier en mm * TAB1.'PY_PHI1' LISTREEL de la position y en m de * chaque barre par rapport au * baricentre de la section * TAB1.'PZ_PHI1' LISTREEL de la position z en m de * chaque barre par rapport au * baricentre de la section * TAB1.'E_ACIE1' Module Young acier [FLOTTANT] * TAB1.'EP_BET1' deformation ultime beton [FLOTTANT] * TAB1.'EP_ACI1' deformation ultime acier [FLOTTANT] * TAB1.'OM_SY1' omega_s_y [FLOTTANT] * TAB1.'OM_SZ1' omega_s_z [FLOTTANT] * TAB1.'OM_WZ1' omega_w_y [FLOTTANT] * TAB1.'OM_WY1' omega_w_z [FLOTTANT] * TAB1.'LAMB_Y1' lambda_y [FLOTTANT] * TAB1.'LAMB_Z1' lambda_y [FLOTTANT] * TAB1.'CHI_Z1' chi_z [FLOTTANT] * TAB1.'CHI_Y1' chi_y [FLOTTANT] * TAB1.'EFFX0' EFFX statique [FLOTTANT] * TAB1.'EFFY0' EFFY statique [FLOTTANT] * TAB1.'EFFZ0' EFFZ statique [FLOTTANT] * TAB1.'MOMX0' MOMX statique [FLOTTANT] * TAB1.'MOMY0' MOMY statique [FLOTTANT] * TAB1.'MOMZ0' MOMZ statique [FLOTTANT] * TAB1.'NNNN1' (EFFX)2 - matrice X seisme [FLOTTANT] * TAB1.'NNVY1' Cov EFFX et EFFY -matrice X seisme [FLOTTANT] * TAB1.'NNVZ1' Cov EFFX et EFFZ -matrice X seisme [FLOTTANT] * TAB1.'NNMT1' Cov EFFX et MOMX -matrice X seisme [FLOTTANT] * TAB1.'NNMY1' Cov EFFX et MOMY -matrice X seisme [FLOTTANT] * TAB1.'NNMZ1' Cov EFFX et MOMZ -matrice X seisme [FLOTTANT] * TAB1.'VYVY1' (EFFY)2 - matrice X seisme [FLOTTANT] * TAB1.'VYVZ1' Cov EFFY et EFFZ -matrice X seisme [FLOTTANT] * TAB1.'VYMT1' Cov EFFY et MOMX -matrice X seisme [FLOTTANT] * TAB1.'VYMY1' Cov EFFY et MOMY -matrice X seisme [FLOTTANT] * TAB1.'VYMZ1' Cov EFFY et MOMZ -matrice X seisme [FLOTTANT] * TAB1.'VZVZ1' (EFFZ)2 - matrice X seisme [FLOTTANT] * TAB1.'VZMT1' Cov EFFZ et MOMX -matrice X seisme [FLOTTANT] * TAB1.'VZMY1' Cov EFFZ et MOMY -matrice X seisme [FLOTTANT] * TAB1.'VZMZ1' Cov EFFZ et MOMZ -matrice X seisme [FLOTTANT] * TAB1.'MTMT1' (MOMX)2 - matrice X seisme [FLOTTANT] * TAB1.'MTMY1' Cov MOMX et MOMY -matrice X seisme [FLOTTANT] * TAB1.'MTMZ1' Cov MOMX et MOMZ -matrice X seisme [FLOTTANT] * TAB1.'MYMY1' (MOMY)2 - matrice X seisme [FLOTTANT] * TAB1.'MYMZ1' Cov MOMY et MOMZ -matrice X seisme [FLOTTANT] * TAB1.'MZMZ1' (MOMZ)2 - matrice X seisme [FLOTTANT] * TAB1.'FCD1' fcd du beton [FLOTTANT] * TAB1.'FSD1' fsd du acier [FLOTTANT] * * TOL1 Tolerence du calcul (1.e-4) * * * OUTPUT: * * * LAM1 Marge sismique * * * * -------------------------------------------------------- * * * Input Phase BASE_Y1 = TAB_IN2.'BASE_Y1'; BASE_Z1 = TAB_IN2.'BASE_Z1'; LONG_Y1 = TAB_IN2.'LONG_Y1'; LONG_Z1 = TAB_IN2.'LONG_Z1'; ENROBA1 = TAB_IN2.'ENROBA1'; ACADRY1 = TAB_IN2.'ACADRY1'; ACADRZ1 = TAB_IN2.'ACADRZ1'; S_CADR1 = TAB_IN2.'S_CADR1'; LI_PHI1 = TAB_IN2.'LI_PHI1'; PY_PHI1 = TAB_IN2.'PY_PHI1'; PZ_PHI1 = TAB_IN2.'PZ_PHI1'; E_ACIE1 = TAB_IN2.'E_ACIE1'; EP_BET1 = TAB_IN2.'EP_BET1'; EP_ACI1 = TAB_IN2.'EP_ACI1'; OM_SY1 = TAB_IN2.'OM_SY1'; OM_SZ1 = TAB_IN2.'OM_SZ1'; OM_WZ1 = TAB_IN2.'OM_WZ1'; OM_WY1 = TAB_IN2.'OM_WY1'; LAMB_Y1 = TAB_IN2.'LAMB_Y1'; LAMB_Z1 = TAB_IN2.'LAMB_Z1'; CHI_Z1 = TAB_IN2.'CHI_Z1'; CHI_Y1 = TAB_IN2.'CHI_Y1'; EFFX0 = TAB_IN2.'EFFX0'; EFFY0 = TAB_IN2.'EFFY0'; EFFZ0 = TAB_IN2.'EFFZ0'; MOMX0 = TAB_IN2.'MOMX0'; MOMY0 = TAB_IN2.'MOMY0'; MOMZ0 = TAB_IN2.'MOMZ0'; NNNN1 = TAB_IN2.'NNNN1'; NNVY1 = TAB_IN2.'NNVY1'; NNVZ1 = TAB_IN2.'NNVZ1'; NNMT1 = TAB_IN2.'NNMT1'; NNMY1 = TAB_IN2.'NNMY1'; NNMZ1 = TAB_IN2.'NNMZ1'; VYVY1 = TAB_IN2.'VYVY1'; VYVZ1 = TAB_IN2.'VYVZ1'; VYMT1 = TAB_IN2.'VYMT1'; VYMY1 = TAB_IN2.'VYMY1'; VYMZ1 = TAB_IN2.'VYMZ1'; VZVZ1 = TAB_IN2.'VZVZ1'; VZMT1 = TAB_IN2.'VZMT1'; VZMY1 = TAB_IN2.'VZMY1'; VZMZ1 = TAB_IN2.'VZMZ1'; MTMT1 = TAB_IN2.'MTMT1'; MTMY1 = TAB_IN2.'MTMY1'; MTMZ1 = TAB_IN2.'MTMZ1'; MYMY1 = TAB_IN2.'MYMY1'; MYMZ1 = TAB_IN2.'MYMZ1'; MZMZ1 = TAB_IN2.'MZMZ1'; FCD1 = TAB_IN2.'FCD1'; FSD1 = TAB_IN2.'FSD1'; * ENR1 = ENROBA1; * * limites compression et traction NLIMC1 = (-1.0 - OM_SY1)*(((BASE_Y1)*(BASE_Z1))*(FCD1)); NLIMT1 = (OM_SY1)*(((BASE_Y1)*(BASE_Z1))*(FCD1)); * * Static Point TINI1 = TABLE; TINI1.'TYPE' = CHAINE 'P_COURT'; TINI1.'NN' = EFFX0; TINI1.'VY' = EFFY0; TINI1.'VZ' = EFFZ0; TINI1.'MT' = MOMX0; TINI1.'MY' = MOMY0; TINI1.'MZ' = MOMZ0; TINI1.'BY' = BASE_Y1; TINI1.'BZ' = BASE_Z1; TINI1.'FCD' = FCD1; TINI1.'FSD' = FSD1; TINI1.'WSY' = OM_SY1; TINI1.'WSZ' = OM_SZ1; TINI1.'WWY' = OM_WY1; TINI1.'WWZ' = OM_WZ1; TINI1.'LY' = LAMB_Y1; TINI1.'LZ' = LAMB_Z1; TINI1.'XIY' = CHI_Y1; TINI1.'XIZ' = CHI_Z1; TINI1.'PHI_LON' = LI_PHI1; TINI1.'PHI_Y' = PY_PHI1; TINI1.'PHI_Z' = PZ_PHI1; TINI1.'EPSI_BET' = EP_BET1; TINI1.'EPSI_ACI' = EP_ACI1; TINI1.'YOUNG_ACIER' = E_ACIE1; * * 'SI' ('EGA' TYP_CAL1 'ELLIPSOIDE'); * *------------------------------------------------------------- * Methode Ellipse * * Test about the static state stress T_ST1 = G_ULTIFR TINI1; * 'SI' (T_ST1 '' NLIMT1)); RBVAL2 = RBVAL1; RBVAL1 = ((NLIMT1) - (EXTR EF_ST0 'FZ' MAI1))/ (EXTR BJ 'FZ' MAI1); 'FINSI'; 'FINSI'; 'SI' ((((TINIC1.'NN') '>EG' NLIMC1)) 'ET' (((TINIC1.'NN') 'EG' 0.0)); RBVAL2 = RBVAL1; RBVAL1 = ((NLIMT1) - (EXTR EF_ST0 'FZ' MAI1) - ((NLIMT1)*(TOL_1)))/ (EXTR BJ 'FZ' MAI1); 'FINSI'; 'FINSI'; 'FINSI'; 'SINON'; * otherwise RBVAL2 = RBVAL1; RBVAL1 = RBVAL1 - ((F_VAL1)/(D_VAL1)); 'FINSI'; 'FINSI'; I_I1 = I_I1 + 1; 'FIN' IND1; * computation of lambda LAMBJ = RBJ/SBJ; TAB_ITE. I_M1 . 'LAMBDA' = LAMBJ; TAB_ITE. I_M1 . 'RB' = RBJ; XCOUB1 = (((TAB_ITE. I_M1 .'LAMBDA')*(TAB_ITE. I_M1 .'LAMBDA'))* (MATX1)); * amplified stress state using the lambda coef XALP_L = (XCOUB1*ALP_K)/(('XTMX' XCOUB1 ALP_K)**(0.5)); XALP_L = ('CHAN' ATTRIBUT XALP_L NATURE DISCRET) + EF_ST0; * Dh*ex ; Dh*ey ; Dh*ez DE_HX = 'MANU' 'CHPO' MAI1 3 'FX' H_VAL2 'FY' 0.0 'FZ' 0.0; DE_HY = 'MANU' 'CHPO' MAI1 3 'FX' 0.0 'FY' H_VAL2 'FZ' 0.0; DE_HZ = 'MANU' 'CHPO' MAI1 3 'FX' 0.0 'FY' 0.0 'FZ' H_VAL2; * x0*lam + Dh*ex ; x0*lam + Dh*ey ; x0*lam + Dh*ez VV_HXP = XALP_L + DE_HX; VV_HYP = XALP_L + DE_HY; VV_HZP = XALP_L + DE_HZ; * x0*lam + Dh*ex ; x0*lam + Dh*ey ; x0*lam + Dh*ez VV_HXN = XALP_L - DE_HX; VV_HYN = XALP_L - DE_HY; VV_HZN = XALP_L - DE_HZ; * x0*lam VV_0 = XALP_L; * VV_HXP1 = 'EXTR' VV_HXP 'FX' MAI1; VV_HXP2 = 'EXTR' VV_HXP 'FY' MAI1; VV_HXP3 = 'EXTR' VV_HXP 'FZ' MAI1; VV_HYP1 = 'EXTR' VV_HYP 'FX' MAI1; VV_HYP2 = 'EXTR' VV_HYP 'FY' MAI1; VV_HYP3 = 'EXTR' VV_HYP 'FZ' MAI1; VV_HZP1 = 'EXTR' VV_HZP 'FX' MAI1; VV_HZP2 = 'EXTR' VV_HZP 'FY' MAI1; VV_HZP3 = 'EXTR' VV_HZP 'FZ' MAI1; VV_HXN1 = 'EXTR' VV_HXN 'FX' MAI1; VV_HXN2 = 'EXTR' VV_HXN 'FY' MAI1; VV_HXN3 = 'EXTR' VV_HXN 'FZ' MAI1; VV_HYN1 = 'EXTR' VV_HYN 'FX' MAI1; VV_HYN2 = 'EXTR' VV_HYN 'FY' MAI1; VV_HYN3 = 'EXTR' VV_HYN 'FZ' MAI1; VV_HZN1 = 'EXTR' VV_HZN 'FX' MAI1; VV_HZN2 = 'EXTR' VV_HZN 'FY' MAI1; VV_HZN3 = 'EXTR' VV_HZN 'FZ' MAI1; * VV_01 = 'EXTR' VV_0 'FX' MAI1; VV_02 = 'EXTR' VV_0 'FY' MAI1; VV_03 = 'EXTR' VV_0 'FZ' MAI1; * 'SI' ((VV_HZP3) '>' (NLIMT1)); TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HXP3; TINIA1.'VY' = VV_HXP1; TINIA1.'VZ' = VV_HXP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HXN3; TINIB1.'VY' = VV_HXN1; TINIB1.'VZ' = VV_HXN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_X1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((2.0)*(H_VAL2)); * TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HYP3; TINIA1.'VY' = VV_HYP1; TINIA1.'VZ' = VV_HYP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HYN3; TINIB1.'VY' = VV_HYN1; TINIB1.'VZ' = VV_HYN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_Y1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((2.0)*(H_VAL2)); * TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_03; TINIA1.'VY' = VV_01; TINIA1.'VZ' = VV_02; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HZN3; TINIB1.'VY' = VV_HZN1; TINIB1.'VZ' = VV_HZN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_Z1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((1.0)*(H_VAL2)); 'SINON'; 'SI' ((VV_HZN3) '<' (NLIMC1)); TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HXP3; TINIA1.'VY' = VV_HXP1; TINIA1.'VZ' = VV_HXP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HXN3; TINIB1.'VY' = VV_HXN1; TINIB1.'VZ' = VV_HXN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_X1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((2.0)*(H_VAL2)); * TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HYP3; TINIA1.'VY' = VV_HYP1; TINIA1.'VZ' = VV_HYP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HYN3; TINIB1.'VY' = VV_HYN1; TINIB1.'VZ' = VV_HYN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_Y1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((2.0)*(H_VAL2)); * TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HZP3; TINIA1.'VY' = VV_HZP1; TINIA1.'VZ' = VV_HZP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_03; TINIB1.'VY' = VV_01; TINIB1.'VZ' = VV_02; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_Z1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((1.0)*(H_VAL2)); 'SINON'; TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HXP3; TINIA1.'VY' = VV_HXP1; TINIA1.'VZ' = VV_HXP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HXN3; TINIB1.'VY' = VV_HXN1; TINIB1.'VZ' = VV_HXN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_X1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((2.0)*(H_VAL2)); * TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HYP3; TINIA1.'VY' = VV_HYP1; TINIA1.'VZ' = VV_HYP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HYN3; TINIB1.'VY' = VV_HYN1; TINIB1.'VZ' = VV_HYN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_Y1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((2.0)*(H_VAL2)); * TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = VV_HZP3; TINIA1.'VY' = VV_HZP1; TINIA1.'VZ' = VV_HZP2; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = VV_HZN3; TINIB1.'VY' = VV_HZN1; TINIB1.'VZ' = VV_HZN2; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; VAL_Z1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((2.0)*(H_VAL2)); 'FINSI'; 'FINSI'; * * costruction n_limit DELTG1 = 'MANU' 'CHPO' MAI1 'UX' ((-1.0)*(VAL_X1)) 'UY' ((-1.0)*(VAL_Y1)) 'UZ' ((-1.0)*(VAL_Z1)); DELTG1 = DELTG1/((XTX DELTG1)**(0.5)); * * variable for the convergence test ID1_CON = 0; * check on the n_limit 'SI' ((('XTX' DELTG1)**(0.5)) < 1.E-7); * if |n_limit|=0 rotate the chosen direction by 10 degrees MATR1 = 'VIDE' 'RIGIDITE'/'RIGIDITE'; MATR1 = MATR1 'ET' ('MANU' 'RIGIDITE' MAI1 ('MOTS' 'UX' 'UY' 'UZ') 'QUEL' ('PROG' (COS(10.0)) ((-1.0)*(SIN(10.0))) 0.0 (SIN(10.0)) (COS(10.0)) 0.0 0.0 0.0 1.0)); AAA1 = MATR1*ALP_K; DELTG1 = MANU CHPO MAI1 'UX' (EXTR AAA1 MAI1 'FX') 'UY' (EXTR AAA1 MAI1 'FY') 'UZ' (EXTR AAA1 MAI1 'FZ'); 'FINSI'; * convergence test VET_TEST = (TAB_ITE. I_M1 . 'ALPA') - (TAB_ITE. (I_M1 - 1) . 'ALPA'); * 'SI' ((('XTX' VET_TEST)**(0.5)) '<' TOL_2); * deduction of the final value of Lambda 'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0); 'SI' ('NEG' ('TYPE' LAM_FIN1) 'FLOTTANT'); LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA'); 'SINON'; 'SI' (LAM_FIN1 '>' (TAB_ITE. I_M1 . 'LAMBDA')); LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA'); 'FINSI'; 'FINSI'; 'FINSI'; 'QUIT' IND2; 'SINON'; 'SI' (I_M1 'EGA' 200); ID1_CON = 1; 'SI' ((TAB_ITE. I_M1 . 'LAMBDA') '>' 0.0); 'SI' ('NEG' ('TYPE' LAM_FIN1) 'FLOTTANT'); LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA'); 'SINON'; 'SI' (LAM_FIN1 '>' (TAB_ITE. I_M1 . 'LAMBDA')); LAM_FIN1 = (TAB_ITE. I_M1 . 'LAMBDA'); 'FINSI'; 'FINSI'; 'FINSI'; 'SINON'; * otherwise updating the n_limit ALP_K = DELTG1; 'FINSI'; 'FINSI'; * Display * MESS '-----------------------------'; * LIST F_VAL1; * LIST I_I1; * LIST LAMBJ; * MESS '*'; * LIST NLIMT1; * LIST NLIMC1; * LIST TINIC1.'VY'; * LIST TINIC1.'VZ'; * LIST TINIC1.'NN'; * MESS '*'; * LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UX' MAI1); * LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UY' MAI1); * LIST ('EXTR' (TAB_ITE. I_M1 . 'ALPA') 'UZ' MAI1); * LIST ('EXTR' DELTG1 'UX' MAI1); * LIST ('EXTR' DELTG1 'UY' MAI1); * LIST ('EXTR' DELTG1 'UZ' MAI1); * MESS '-----------------------------'; I_M1 = I_M1 + 1; 'FIN' IND2; I_K1 = I_K1 + 1; 'FIN' IND3; * * screen message 'SI' (ID1_CON 'EGA' 0); TIT1 = 'CHAINE' ' Element ' I_ELE1 ' Lambda=' LAM_FIN1; 'MESS' TIT1; 'FINSI'; 'SI' (ID1_CON 'EGA' 1); TIT1 = 'CHAINE' ' Element ' I_ELE1 ' Lambda=' LAM_FIN1 '.....Attention non convergence'; 'FINSI'; 'FINSI'; 'SINON'; * *------------------------------------------------------------- * Methode Rectangle * * Test about the static state stress T_ST1 = G_ULTIFR TINI1; * ID1_CON = 0; 'SI' (T_ST1 'EG' NLIMT1); TINIA1 = 'COPI' TINI1; TINIB1 = 'COPI' TINI1; TINIA1.'NN' = Z1_0; TINIA1.'VY' = X1_0; TINIA1.'VZ' = Y1_0; TINIA1.'MT' = 0.0; TINIA1.'MY' = 0.0; TINIA1.'MZ' = 0.0; TINIB1.'NN' = Z1_DN; TINIB1.'VY' = X1_DN; TINIB1.'VZ' = Y1_DN; TINIB1.'MT' = 0.0; TINIB1.'MY' = 0.0; TINIB1.'MZ' = 0.0; D_VAL1 = ((G_ULTIFR TINIA1) - (G_ULTIFR TINIB1))/ ((1.0)*(H_VAL3)); 'SINON'; 'SI' (Z1_DN 'EG' (NLIMT1)); LAMB1 = LAMB0; LAMB0 = ((NLIMT1) - EFFX0 - ((NLIMT1)*(TOL_2)))/ (((NNNN1)**(0.5))*('EXTR' LL_LAMZ II1)); 'SINON'; 'SI' ((Z1_0) '