|
1
|
---------------------------------------------------------------------
|
|
2
|
>>>> algorithme (1): dynamique_relaxation_dynam >>>>
|
|
3
|
---------------------------------------------------------------------
|
|
4
|
modification de lambda 2 --> 4 (delta= 2)
|
|
5
|
recalcul de la pseudo-masse (type: 0)
|
|
6
|
modification de lambda 4 --> 2 (delta= -2)
|
|
7
|
recalcul de la pseudo-masse (type: 0)
|
|
8
|
======================================================================
|
|
9
|
INCREMENT DE CHARGE : 3 intensite 1 t= 3 dt= 1
|
|
10
|
======================================================================
|
|
11
|
modification de lambda 2 --> 4 (delta= 2)
|
|
12
|
*** pb sur le calcul de la contrainte: le resultat est soit infini soit un nan mail: 1 ele: 161 npti: 1
|
|
13
|
-- deformation: ---
|
|
14
|
nb noeud concerne = 2 numInteg= 1
|
|
15
|
-- information concernant les noeuds:
|
|
16
|
|
|
17
|
--noeud --
|
|
18
|
Numero : 3, du maillage: 1
|
|
19
|
Coordonnee(s) initiale(s) du noeud :
|
|
20
|
{ 0.05 0.1 0 }
|
|
21
|
Coordonnee(s) a l'instant t :
|
|
22
|
{ 0.05 0.1 0 }
|
|
23
|
Coordonnee(s) a l'instant t+dt :
|
|
24
|
{ 0.05 0.016246 0 }
|
|
25
|
Degre(s) de liberte lie(s) au noeud :
|
|
26
|
Taille du tableau : 12 , a 0, a t ( et t+dt ?)
|
|
27
|
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
|
|
28
|
ddl_etendu: pression_ext 6.43154e-06 Sigma_principaleI 2.19636e+06 Sigma_principaleII 371003
|
|
29
|
ddl_quelconque: NN_SURF_t0 Coordonnee dim= 3 1 0 3.2051034673266428e-09
|
|
30
|
NN_SURF_t Coordonnee dim= 3 0 0 0
|
|
31
|
XI_ITER_0 Coordonnee dim= 3 0 0 0
|
|
32
|
MASSE_RELAX_DYN Coordonnee dim= 3 43709.702790363539 43709.702790363539 43709.702790363539
|
|
33
|
FORCE_GENE_INT_t Coordonnee dim= 3 2.8561967312964032 -0.30755123671937262 2.0825030946375294
|
|
34
|
FORCE_GENE_EXT_t Coordonnee dim= 3 3.2426955974737797e-05 1.6640258543983685e-08 -0.0012009033822452167
|
|
35
|
CONTRAINTE_COURANTE Grandeur_TenseurHH= Tenseur3HH 1404198.0309479809 416455.35061452439 746709.5868425098 196701.67378671496 143459.65782369769 1023976.6492183005
|
|
36
|
DIRECTION_PLI Tableau_de_Coordonnee_nb= 2 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0
|
|
37
|
-- fin noeud --
|
|
38
|
|
|
39
|
--noeud --
|
|
40
|
Numero : 6, du maillage: 1
|
|
41
|
Coordonnee(s) initiale(s) du noeud :
|
|
42
|
{ 0.04999999948449272 0.1312770006688897 0.16083951359128573 }
|
|
43
|
Coordonnee(s) a l'instant t :
|
|
44
|
{ 0.18029820338444699 0.1312770006688897 0.095018213959489908 }
|
|
45
|
Coordonnee(s) a l'instant t+dt :
|
|
46
|
{ 0.18029820338444699 0.1312770006688897 0.095018213959489908 }
|
|
47
|
Degre(s) de liberte lie(s) au noeud :
|
|
48
|
Taille du tableau : 12 , a 0, a t ( et t+dt ?)
|
|
49
|
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
|
|
50
|
ddl_etendu: reaction_normale 0.21051508199761262 reaction_tangentielle 0 pression_ext 0.070957322566171002 Sigma_principaleI 2002143.7504815042 Sigma_principaleII 312639.31521304499
|
|
51
|
ddl_quelconque: NN_SURF_t0 Coordonnee dim= 3 1 -4.1308696573974632e-18 3.20510346106244e-09
|
|
52
|
NN_SURF_t Coordonnee dim= 3 0 0 0
|
|
53
|
XI_ITER_0 Coordonnee dim= 3 0 0 0
|
|
54
|
MASSE_RELAX_DYN Coordonnee dim= 3 74168.972671200245 74168.972671200245 74168.972671200245
|
|
55
|
FORCE_GENE_INT_t Coordonnee dim= 3 0.00037300470074796976 -1.5553315842905366 0.0033047357077766737
|
|
56
|
FORCE_GENE_EXT_t Coordonnee dim= 3 0.00023794466577214474 2.7225891886925192e-07 -0.0027984221919892431
|
|
57
|
CONTRAINTE_COURANTE Grandeur_TenseurHH= Tenseur3HH 1272759.2485056848 363827.50003149424 678196.31715737004 231298.48361126054 169066.09899349156 929074.44778083672
|
|
58
|
DIRECTION_PLI Tableau_de_Coordonnee_nb= 2 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0 Grandeur_coordonnee= Coordonnee dim= 3 0 0 0
|
|
59
|
-- fin noeud --
|
|
60
|
|
|
61
|
-- information concernant la metrique actuelle au pti --
|
|
62
|
-- info elements de base de la metrique --
|
|
63
|
giB_0:
|
|
64
|
BaseB dimension: 3 nombre de vecteurs : 1
|
|
65
|
les vecteurs de la base: Tableau :taille= 1 CoordonneeB dim= 3 -2.5775364134728029e-10 0.015638500334444846 0.080419756795642866
|
|
66
|
giB_t:
|
|
67
|
BaseB dimension: 3 nombre de vecteurs : 1
|
|
68
|
les vecteurs de la base: Tableau :taille= 1 CoordonneeB dim= 3 0.065149101692223488 0.015638500334444846 0.047509106979744954
|
|
69
|
giB_tdt:
|
|
70
|
BaseB dimension: 3 nombre de vecteurs : 1
|
|
71
|
les vecteurs de la base: Tableau :taille= 1 CoordonneeB dim= 3 0.065149101692223488 0.057515507934444848 0.047509106979744954
|
|
72
|
giH_0:
|
|
73
|
BaseH dimension: 3 nombre de vecteurs : 1
|
|
74
|
les vecteurs de la base: Tableau :taille= 1 CoordonneeH dim= 3 -3.8402485477638015e-08 2.329966237708371 11.981667945861759
|
|
75
|
giH_t:
|
|
76
|
BaseH dimension: 3 nombre de vecteurs : 1
|
|
77
|
les vecteurs de la base: Tableau :taille= 1 CoordonneeH dim= 3 9.6573223195752895 2.3181599500480456 7.0424725330226154
|
|
78
|
giH_tdt:
|
|
79
|
BaseH dimension: 3 nombre de vecteurs : 1
|
|
80
|
les vecteurs de la base: Tableau :taille= 1 CoordonneeH dim= 3 6.6413926021425871 5.86321313698007 4.8431463125966889
|
|
81
|
gijBB_0: Tenseur1BB 0.0067118999757807786
|
|
82
|
gijBB_t: Tenseur1BB 0.0067460833900269603
|
|
83
|
gijBB_tdt: Tenseur1BB 0.0098095543502737161
|
|
84
|
gijHH_0: Tenseur1HH 148.98910943375202
|
|
85
|
gijHH_t: Tenseur1HH 148.23415931655177
|
|
86
|
gijHH_tdt: Tenseur1HH 101.94143019066885
|
|
87
|
jacobien_0= 0.081926186142043611 jacobien_t= 0.082134544437934026 jacobien_tdt= 0.099043194366264839
|
|
88
|
-- fin info elements de metrique --
|
|
89
|
recalcul de la pseudo-masse (type: 0)
|
|
90
|
contact: reaction ==> F_N max 0.12812998616710702 F_T max 0, maxi gap_N : -0.079645315450655854 gap_T : 3.4099255362731258
|
|
91
|
contact: reaction ==> F_N max 0.12812998616710702 F_T max 0, maxi gap_N : -0.079645315450655854 gap_T : 3.4099255362731258
|
|
92
|
IT:1
|
|
93
|
(assemblage 3) reaction maxi : 103.00775058637555, du ddl GAMMA1, du noeud/maillage 3/1
|
|
94
|
|
|
95
|
contact: reaction ==> F_N max 0.12810277588007843 F_T max 0, maxi gap_N : -0.079628401596540013 gap_T : 3.4099255503801582
|
|
96
|
contact: reaction ==> F_N max 0.12810277588007843 F_T max 0, maxi gap_N : -0.079628401596540013 gap_T : 3.4099255503801582
|
|
97
|
IT:2
|
|
98
|
(assemblage 3) reaction maxi : 94.463911917036128, du ddl GAMMA1, du noeud/maillage 3/1
|
|
99
|
|
|
100
|
max puissance exterieure = 0.44926640958824149
|
|
101
|
max puissance interieurs = 94.463892229083754
|
|
102
|
max des reactions = 94.463911917036128
|
|
103
|
max du residu total = 80.976360004044096
|
|
104
|
E_cinetique = 0.33084625453978095
|
|
105
|
E_ext = -37.743241925932296
|
|
106
|
E_int = 3.7737815484582344
|
|
107
|
**** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.85722005748779906
|
|
108
|
|
|
109
|
contact: reaction ==> F_N max 0.12806937538712265 F_T max 0, maxi gap_N : -0.079607639924138038 gap_T : 3.4099256421823014
|
|
110
|
contact: reaction ==> F_N max 0.12806937538712265 F_T max 0, maxi gap_N : -0.079607639924138038 gap_T : 3.4099256421823014
|
|
111
|
IT:3
|
|
112
|
(assemblage 3) reaction maxi : 82.065065001501893, du ddl GAMMA1, du noeud/maillage 3/1
|
|
113
|
|
|
114
|
max puissance exterieure = 0.44925855199480347
|
|
115
|
max puissance interieurs = 82.065046185933895
|
|
116
|
max des reactions = 82.065065001501893
|
|
117
|
max du residu total = 63.668750150473286
|
|
118
|
E_cinetique = 0.79295682981549376
|
|
119
|
E_ext = -37.743540921107723
|
|
120
|
E_int = 3.3008198626390186
|
|
121
|
**** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.77583256833231129
|
|
122
|
|
|
123
|
contact: reaction ==> F_N max 0.12803595446781763 F_T max 0, maxi gap_N : -0.079586865554762645 gap_T : 3.4099258703657584
|
|
124
|
contact: reaction ==> F_N max 0.12803595446781763 F_T max 0, maxi gap_N : -0.079586865554762645 gap_T : 3.4099258703657584
|
|
125
|
IT:4
|
|
126
|
(assemblage 3) reaction maxi : 66.287471288261287, du ddl GAMMA1, du noeud/maillage 3/1
|
|
127
|
|
|
128
|
max puissance exterieure = 0.44924874593371783
|
|
129
|
max puissance interieurs = 66.287453522887347
|
|
130
|
max des reactions = 66.287471288261287
|
|
131
|
max du residu total = 48.838310268157336
|
|
132
|
E_cinetique = 1.2636553829606565
|
|
133
|
E_ext = -37.743964858221503
|
|
134
|
E_int = 2.8214074127945432
|
|
135
|
**** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.73676532411044038
|
|
136
|
|
|
137
|
contact: reaction ==> F_N max 0.12800614180684994 F_T max 0, maxi gap_N : -0.079568334070773433 gap_T : 3.4099262695403176
|
|
138
|
contact: reaction ==> F_N max 0.12800614180684994 F_T max 0, maxi gap_N : -0.079568334070773433 gap_T : 3.4099262695403176
|
|
139
|
IT:5
|
|
140
|
(assemblage 3) reaction maxi : 47.927912057458485, du ddl GAMMA1, du noeud/maillage 3/1
|
|
141
|
|
|
142
|
max puissance exterieure = 0.44923748923692552
|
|
143
|
max puissance interieurs = 55.805264712597918
|
|
144
|
max des reactions = 47.927912057458485
|
|
145
|
max du residu total = 55.752481746300795
|
|
146
|
E_cinetique = 1.603201422287055
|
|
147
|
E_ext = -37.744483922117247
|
|
148
|
E_int = 2.4780236307317005
|
|
149
|
**** E_cinetique/E_statique_ET_Res/Reac_et_Fext --> 0.99905415794425567
|