--------------------------------------------------------------------- >>>> algorithme (1): dynamique_relaxation_dynam >>>> --------------------------------------------------------------------- modification de lambda 2 --> 4 (delta= 2) recalcul de la pseudo-masse (type: 0) modification de lambda 4 --> 2 (delta= -2) recalcul de la pseudo-masse (type: 0) ====================================================================== INCREMENT DE CHARGE : 3 intensite 1 t= 3 dt= 1 ====================================================================== modification de lambda 2 --> 4 (delta= 2) *** pb sur le calcul de la contrainte: le resultat est soit infini soit un nan mail: 1 ele: 161 npti: 1 -- deformation: --- nb noeud concerne = 2 numInteg= 1 -- information concernant les noeuds: --noeud -- Numero : 3, du maillage: 1 Coordonnee(s) initiale(s) du noeud : { 0.05 0.1 0 } Coordonnee(s) a l'instant t : { 0.05 0.1 0 } Coordonnee(s) a l'instant t+dt : { 0.05 0.016246 0 } Degre(s) de liberte lie(s) au noeud : Taille du tableau : 12 , a 0, a t ( et t+dt ?) 0.05 X1 = 0.05 val_fixe, 0.05 0.1 X2 = 0.1 val_fixe, 0.016246 0 X3 = 0 val_fixe, 0 0 V1 = 0 val_fixe, 0 0 V2 = 0 val_fixe, -0.083754 0 V3 = 0 val_fixe, 0 0 GAMMA1 = 0 val_fixe, 0 0 GAMMA2 = 0 val_fixe, 0 0 GAMMA3 = 0 val_fixe, 0 0 R_X1 = 2.85623 lisible_sur_fixe, 0 0 R_X2 = -0.307551 lisible_sur_fixe, 0 0 R_X3 = 2.0813 lisible_sur_fixe, 0 ddl_etendu: pression_ext 6.43154e-06 Sigma_principaleI 2.19636e+06 Sigma_principaleII 371003 ddl_quelconque: NN_SURF_t0 Coordonnee dim= 3 1 0 3.2051034673266428e-09 NN_SURF_t Coordonnee dim= 3 0 0 0 XI_ITER_0 Coordonnee dim= 3 0 0 0 MASSE_RELAX_DYN Coordonnee dim= 3 43709.702790363539 43709.702790363539 43709.702790363539 FORCE_GENE_INT_t Coordonnee dim= 3 2.8561967312964032 -0.30755123671937262 2.0825030946375294 FORCE_GENE_EXT_t Coordonnee dim= 3 3.2426955974737797e-05 1.6640258543983685e-08 -0.0012009033822452167 CONTRAINTE_COURANTE Grandeur_TenseurHH= Tenseur3HH 1404198.0309479809 416455.35061452439 746709.5868425098 196701.67378671496 143459.65782369769 1023976.6492183005 DIRECTION_PLI Tableau_de_Coordonnee_nb= 2 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0 -- fin noeud -- --noeud -- Numero : 6, du maillage: 1 Coordonnee(s) initiale(s) du noeud : { 0.04999999948449272 0.1312770006688897 0.16083951359128573 } Coordonnee(s) a l'instant t : { 0.18029820338444699 0.1312770006688897 0.095018213959489908 } Coordonnee(s) a l'instant t+dt : { 0.18029820338444699 0.1312770006688897 0.095018213959489908 } Degre(s) de liberte lie(s) au noeud : Taille du tableau : 12 , a 0, a t ( et t+dt ?) 0.04999999948449272 X1 = 0.18029820338444699 val_libre, 0.18029820338444699 0.1312770006688897 X2 = 0.1312770006688897 val_libre, 0.1312770006688897 0.16083951359128573 X3 = 0.095018213959489908 val_libre, 0.095018213959489908 0 V1 = 0 val_libre, 0 0 V2 = 0 val_libre, 0 0 V3 = 0 val_libre, 0 0 GAMMA1 = 0 val_libre, 0 0 GAMMA2 = 0 val_libre, 0 0 GAMMA3 = 0 val_libre, 0 0 R_X1 = 0.076923207059187509 lisible_fixe, 0.065052684781512451 0 R_X2 = -1.5553313120316177 lisible_sur_fixe, -0.20021175777455918 0 R_X3 = 0 lisible_fixe, 0 ddl_etendu: reaction_normale 0.21051508199761262 reaction_tangentielle 0 pression_ext 0.070957322566171002 Sigma_principaleI 2002143.7504815042 Sigma_principaleII 312639.31521304499 ddl_quelconque: NN_SURF_t0 Coordonnee dim= 3 1 -4.1308696573974632e-18 3.20510346106244e-09 NN_SURF_t Coordonnee dim= 3 0 0 0 XI_ITER_0 Coordonnee dim= 3 0 0 0 MASSE_RELAX_DYN Coordonnee dim= 3 74168.972671200245 74168.972671200245 74168.972671200245 FORCE_GENE_INT_t Coordonnee dim= 3 0.00037300470074796976 -1.5553315842905366 0.0033047357077766737 FORCE_GENE_EXT_t Coordonnee dim= 3 0.00023794466577214474 2.7225891886925192e-07 -0.0027984221919892431 CONTRAINTE_COURANTE Grandeur_TenseurHH= Tenseur3HH 1272759.2485056848 363827.50003149424 678196.31715737004 231298.48361126054 169066.09899349156 929074.44778083672 DIRECTION_PLI Tableau_de_Coordonnee_nb= 2 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0 -- fin noeud -- -- information concernant la metrique actuelle au pti -- -- info elements de base de la metrique -- giB_0: BaseB dimension: 3 nombre de vecteurs : 1 les vecteurs de la base: Tableau :taille= 1 CoordonneeB dim= 3 -2.5775364134728029e-10 0.015638500334444846 0.080419756795642866 giB_t: BaseB dimension: 3 nombre de vecteurs : 1 les vecteurs de la base: Tableau :taille= 1 CoordonneeB dim= 3 0.065149101692223488 0.015638500334444846 0.047509106979744954 giB_tdt: BaseB dimension: 3 nombre de vecteurs : 1 les vecteurs de la base: Tableau :taille= 1 CoordonneeB dim= 3 0.065149101692223488 0.057515507934444848 0.047509106979744954 giH_0: BaseH dimension: 3 nombre de vecteurs : 1 les vecteurs de la base: Tableau :taille= 1 CoordonneeH dim= 3 -3.8402485477638015e-08 2.329966237708371 11.981667945861759 giH_t: BaseH dimension: 3 nombre de vecteurs : 1 les vecteurs de la base: Tableau :taille= 1 CoordonneeH dim= 3 9.6573223195752895 2.3181599500480456 7.0424725330226154 giH_tdt: BaseH dimension: 3 nombre de vecteurs : 1 les vecteurs de la base: Tableau :taille= 1 CoordonneeH dim= 3 6.6413926021425871 5.86321313698007 4.8431463125966889 gijBB_0: Tenseur1BB 0.0067118999757807786 gijBB_t: Tenseur1BB 0.0067460833900269603 gijBB_tdt: Tenseur1BB 0.0098095543502737161 gijHH_0: Tenseur1HH 148.98910943375202 gijHH_t: Tenseur1HH 148.23415931655177 gijHH_tdt: Tenseur1HH 101.94143019066885 jacobien_0= 0.081926186142043611 jacobien_t= 0.082134544437934026 jacobien_tdt= 0.099043194366264839 -- fin info elements de metrique -- recalcul de la pseudo-masse (type: 0) contact: reaction ==> F_N max 0.12812998616710702 F_T max 0, maxi gap_N : -0.079645315450655854 gap_T : 3.4099255362731258 contact: reaction ==> F_N max 0.12812998616710702 F_T max 0, maxi gap_N : -0.079645315450655854 gap_T : 3.4099255362731258 IT:1 (assemblage 3) reaction maxi : 103.00775058637555, du ddl GAMMA1, du noeud/maillage 3/1 contact: reaction ==> F_N max 0.12810277588007843 F_T max 0, maxi gap_N : -0.079628401596540013 gap_T : 3.4099255503801582 contact: reaction ==> F_N max 0.12810277588007843 F_T max 0, maxi gap_N : -0.079628401596540013 gap_T : 3.4099255503801582 IT:2 (assemblage 3) reaction maxi : 94.463911917036128, du ddl GAMMA1, du noeud/maillage 3/1 max puissance exterieure = 0.44926640958824149 max puissance interieurs = 94.463892229083754 max des reactions = 94.463911917036128 max du residu total = 80.976360004044096 E_cinetique = 0.33084625453978095 E_ext = -37.743241925932296 E_int = 3.7737815484582344 **** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.85722005748779906 contact: reaction ==> F_N max 0.12806937538712265 F_T max 0, maxi gap_N : -0.079607639924138038 gap_T : 3.4099256421823014 contact: reaction ==> F_N max 0.12806937538712265 F_T max 0, maxi gap_N : -0.079607639924138038 gap_T : 3.4099256421823014 IT:3 (assemblage 3) reaction maxi : 82.065065001501893, du ddl GAMMA1, du noeud/maillage 3/1 max puissance exterieure = 0.44925855199480347 max puissance interieurs = 82.065046185933895 max des reactions = 82.065065001501893 max du residu total = 63.668750150473286 E_cinetique = 0.79295682981549376 E_ext = -37.743540921107723 E_int = 3.3008198626390186 **** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.77583256833231129 contact: reaction ==> F_N max 0.12803595446781763 F_T max 0, maxi gap_N : -0.079586865554762645 gap_T : 3.4099258703657584 contact: reaction ==> F_N max 0.12803595446781763 F_T max 0, maxi gap_N : -0.079586865554762645 gap_T : 3.4099258703657584 IT:4 (assemblage 3) reaction maxi : 66.287471288261287, du ddl GAMMA1, du noeud/maillage 3/1 max puissance exterieure = 0.44924874593371783 max puissance interieurs = 66.287453522887347 max des reactions = 66.287471288261287 max du residu total = 48.838310268157336 E_cinetique = 1.2636553829606565 E_ext = -37.743964858221503 E_int = 2.8214074127945432 **** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.73676532411044038 contact: reaction ==> F_N max 0.12800614180684994 F_T max 0, maxi gap_N : -0.079568334070773433 gap_T : 3.4099262695403176 contact: reaction ==> F_N max 0.12800614180684994 F_T max 0, maxi gap_N : -0.079568334070773433 gap_T : 3.4099262695403176 IT:5 (assemblage 3) reaction maxi : 47.927912057458485, du ddl GAMMA1, du noeud/maillage 3/1 max puissance exterieure = 0.44923748923692552 max puissance interieurs = 55.805264712597918 max des reactions = 47.927912057458485 max du residu total = 55.752481746300795 E_cinetique = 1.603201422287055 E_ext = -37.744483922117247 E_int = 2.4780236307317005 **** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.99905415794425567