deplacement test debug vers la batterie : Test_R_dilatation_thermique_LOI_ADDITIVE_MELANGES; modif tests Test_R_dilatation_thermique_LOI_DES_MELANGES (la temperature passe desormais de 20°C a 5°C au lieu de 0°C pour avoir une contribution non nulle de la loi 0°C); modif test Test_R_traction_biaxiale_BLOCAGE_RELATIF_1 : modif d un mot-cle dans le README (BLOCAGE_RELATIF_ au lieu de BLOQUAGE_RELATIF_)
This commit is contained in:
troufflard
parent
6f955a2324
commit
6e0d2c7ec3
29 changed files with 16 additions and 1565 deletions
|
|
@ -8,14 +8,14 @@ Julien Troufflard (troufflard@univ-ubs.fr)
|
|||
Mots-cles
|
||||
------------------------------------------------------
|
||||
Traction biaxiale
|
||||
BLOQUAGE_RELATIF_
|
||||
BLOCAGE_RELATIF_
|
||||
TEMPS_MINI
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------
|
||||
Vérification du fonctionnement de la combinaison BLOQUAGE_RELATIF_ + TEMPS_MINI dans une condition de déplacement imposé
|
||||
Vérification du fonctionnement de la combinaison BLOCAGE_RELATIF_ + TEMPS_MINI dans une condition de déplacement imposé
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
|
|
@ -26,7 +26,7 @@ Calcul de traction biaxiale sur un cube unitaire.
|
|||
La gestion du chargement se fait via des courbes de charge :
|
||||
1) Entre le temps 0 et 1, c est de la traction uniaxiale classique selon X car la condition sur UY de la face Y=1 est inactive jusqu au temps 1 inclus (TEMPS_MINI= 1).
|
||||
|
||||
2) Entre le temps 1 et 2, la courbe dpl_UX permet de figer la face X=1 dans sa position actuelle (X=1.02). La condition sur la face Y=1 devient active, c'est-à-dire un déplacement UY=0.02 (piloté par la courbe dpl_UY_relatif) depuis la position atteinte au temps 1 (BLOQUAGE_RELATIF_).
|
||||
2) Entre le temps 1 et 2, la courbe dpl_UX permet de figer la face X=1 dans sa position actuelle (X=1.02). La condition sur la face Y=1 devient active, c'est-à-dire un déplacement UY=0.02 (piloté par la courbe dpl_UY_relatif) depuis la position atteinte au temps 1 (BLOCAGE_RELATIF_).
|
||||
|
||||
|
||||
-------------------------------------------------------------
|
||||
|
|
|
|||
|
|
@ -23,24 +23,11 @@ But du test
|
|||
------------------------------------------------------
|
||||
test de la combinaison "dilatation_thermique + LOI_ADDITIVE_EN_SIGMA + LOI_DES_MELANGES_EN_SIGMA" avec une loi des mélanges gouvernée par la température
|
||||
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
**
|
||||
** ----- remarques a supprimer une fois le bug resolu -----
|
||||
**
|
||||
** ce test renvoie "Segmentation fault".
|
||||
** une fois resolu, il faudra creer un .maple.ref. Mais avant ça, il ne faut pas oublier
|
||||
** le probleme sur la sortie contrainte_individuelle_a_chaque_loi_a_t
|
||||
** vu dans le test en attente debug : Tests_en_attente_debug/Test_R_sortie_contrainte_individuelle_a_chaque_loi_a_t_1
|
||||
**
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------
|
||||
traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE). La température initialement à 20°C évolue linéairement en fonction du temps pour atteindre 0°C à la fin du calcul
|
||||
traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE). La température initialement à 20°C évolue linéairement en fonction du temps pour atteindre 5°C à la fin du calcul.
|
||||
|
||||
Loi de comportement :
|
||||
- loi additive en contrainte consituée d'une loi ISOLEAS + un mélange de 2 lois élastiques ISOELAS
|
||||
|
|
@ -73,8 +73,8 @@ N_arriere UY
|
|||
N_bas UZ
|
||||
N_gauche 'UX=0.1'
|
||||
|
||||
#temperature initialement a 20 passe a 0 (delta = -20)
|
||||
N_tout 'TEMP= -20'
|
||||
#temperature initialement a 20 passe a 5 (delta = -15)
|
||||
N_tout 'TEMP= -15'
|
||||
|
||||
|
||||
initialisation
|
||||
|
|
@ -1,6 +1,6 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.768 #
|
||||
# Visualisation elements finis : Herezh++ V6.774 #
|
||||
# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
|
@ -72,11 +72,12 @@
|
|||
# noeud_8 [2]X [3]Y [4]Z [5]X1 [6]X2 [7]X3 [8]TEMP ;
|
||||
#
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 4 20 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
#1 4 30 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_1 pt_integ_1: [10]X [11]Y [12]Z [13] EPS11 [14] EPS22 [15] EPS33 [16] SIG11
|
||||
# [17]X [18]Y [19]Z
|
||||
# VOLUME_PTI TYPE_SIMPLE SCALAIRE 1 [20...20]
|
||||
# contrainte_individuelle_a_chaque_loi_a_t TABLEAU_T TENSEURHH 12 [21...32]
|
||||
# contrainte_individuelle_a_chaque_loi_a_t TABLEAU_T TENSEURHH 24 [21...44]
|
||||
# contrainte_individuelle_a_t_sans_proportion TABLEAU_T TENSEURHH 12 [45...56]
|
||||
#
|
||||
#
|
||||
#
|
||||
|
|
@ -101,4 +102,4 @@
|
|||
# puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs
|
||||
# etc. pour tous les points d'integration - element
|
||||
|
||||
1.000000000000e+00 1.100000000000e+00 9.426378600697e-01 9.426378600697e-01 1.100000000000e+00 9.426378600697e-01 9.426378600697e-01 0.000000000000e+00 1.000000000000e+00 2.324573519457e-01 7.434350411647e-01 7.434350411647e-01 1.067768595041e-01 -4.270431673240e-02 -4.270431673240e-02 3.203856676774e+02 2.324573519457e-01 7.434350411647e-01 7.434350411647e-01 1.221778435951e-01 1.067952225591e+02 1.217001211349e-13 2.295381876671e-02 -5.817482815057e-16 3.040041898616e-29 2.295381876663e-02 2.135904451183e+02 2.434002422697e-13 4.590763753342e-02 -1.163496563011e-15 6.080083797233e-29 4.590763753326e-02
|
||||
1.000000000000e+00 1.100000000000e+00 9.616930670808e-01 9.616930670808e-01 1.100000000000e+00 9.616930670808e-01 9.616930670808e-01 5.000000000000e+00 1.000000000000e+00 2.324573519457e-01 7.584634091188e-01 7.584634091188e-01 1.017768595041e-01 -2.562613172673e-02 -2.562613172673e-02 3.894272479161e+02 2.324573519457e-01 7.584634091188e-01 7.584634091188e-01 1.271673638498e-01 9.826584910042e+01 9.920918663753e-14 -1.755505201854e+01 -1.287619828688e-14 -1.134770635631e-13 -1.755505201854e+01 2.911613988157e+02 2.343499056149e-13 1.757164200580e+01 -3.041589146612e-14 -2.680531918140e-13 1.757164200580e+01 1.235608245011e+01 1.136771930222e-14 -9.150624697747e-01 -1.475397720371e-15 -1.300258019993e-14 -9.150624697747e-01 1.813028017614e+02 1.416592713045e-13 1.592391891368e+01 -1.838572543847e-14 -1.620321532607e-13 1.592391891368e+01 1.976973192017e+02 1.818835088355e-13 -1.464099951640e+01 -2.360636352594e-14 -2.080412831989e-13 -1.464099951640e+01 3.223160920203e+02 2.518387045414e-13 2.830918917987e+01 -3.268573411284e-14 -2.880571613524e-13 2.830918917987e+01
|
||||
|
|
@ -26,7 +26,7 @@ test de la combinaison "dilatation_thermique + LOI_DES_MELANGES_EN_SIGMA" avec u
|
|||
------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------
|
||||
traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE). La température initialement à 20°C évolue linéairement en fonction du temps pour atteindre 0°C à la fin du calcul
|
||||
traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE). La température initialement à 20°C évolue linéairement en fonction du temps pour atteindre 5°C à la fin du calcul
|
||||
|
||||
Loi de comportement :
|
||||
- loi des mélanges en contrainte : 2 lois élastiques ISOELAS
|
||||
|
|
|
|||
|
|
@ -66,8 +66,8 @@ N_arriere UY
|
|||
N_bas UZ
|
||||
N_gauche 'UX=0.1'
|
||||
|
||||
#temperature initialement a 20 passe a 0 (delta = -20)
|
||||
N_tout 'TEMP= -20'
|
||||
#temperature initialement a 20 passe a 5 (delta = -15)
|
||||
N_tout 'TEMP= -15'
|
||||
|
||||
|
||||
initialisation
|
||||
|
|
|
|||
|
|
@ -1,6 +1,6 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.770 #
|
||||
# Visualisation elements finis : Herezh++ V6.774 #
|
||||
# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
|
@ -102,4 +102,4 @@
|
|||
# puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs
|
||||
# etc. pour tous les points d'integration - element
|
||||
|
||||
1.000000000000e+00 1.100000000000e+00 9.426378600697e-01 9.426378600697e-01 1.100000000000e+00 9.426378600697e-01 9.426378600697e-01 0.000000000000e+00 1.000000000000e+00 2.324573519457e-01 7.434350411647e-01 7.434350411647e-01 1.067768595041e-01 -4.270431673240e-02 -4.270431673240e-02 2.135904451183e+02 2.324573519457e-01 7.434350411647e-01 7.434350411647e-01 1.221778435951e-01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 2.135904451183e+02 2.434002422697e-13 4.590763753342e-02 -1.163496563011e-15 6.080083797233e-29 4.590763753326e-02 1.067952225591e+02 1.217001211349e-13 2.295381876671e-02 -5.817482815057e-16 3.040041898616e-29 2.295381876663e-02 2.135904451183e+02 2.434002422697e-13 4.590763753342e-02 -1.163496563011e-15 6.080083797233e-29 4.590763753326e-02
|
||||
1.000000000000e+00 1.100000000000e+00 9.485545136316e-01 9.485545136316e-01 1.100000000000e+00 9.485545136316e-01 9.485545136316e-01 5.000000000000e+00 1.000000000000e+00 2.324573519457e-01 7.481013587089e-01 7.481013587089e-01 1.017768595041e-01 -4.070642038261e-02 -4.070642038261e-02 1.781311212274e+02 2.324573519457e-01 7.481013587089e-01 7.481013587089e-01 1.237164039830e-01 6.361825758123e+00 2.503888242049e-15 9.650488932380e-04 -6.291924042228e-15 1.771277769853e-30 9.650488932333e-04 1.145128636462e+02 4.506998835688e-14 1.737088007828e-02 -1.132546327601e-13 3.188338504334e-29 1.737088007820e-02 1.017892121300e+02 4.006221187278e-14 1.544078229181e-02 -1.006707846756e-13 2.833967394567e-29 1.544078229173e-02 2.035784242599e+02 8.012442374556e-14 3.088156458362e-02 -2.013415693513e-13 5.667934789134e-29 3.088156458347e-02
|
||||
|
|
|
|||
|
|
@ -1,46 +0,0 @@
|
|||
------------------------------------------------------
|
||||
Auteur
|
||||
------------------------------------------------------
|
||||
Julien Troufflard (troufflard@univ-ubs.fr)
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------
|
||||
RESTART
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------
|
||||
Vérification du RESTART pour un élément 1D
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------
|
||||
Calcul simple de traction sur une barre (POUT BIE1) encastrée-libre afin de vérifier le fonctionnement d un RESTART. Le calcul se fait en deux parties. Jusqu'au temps 1, la barre est soumise à une force ponctuelle. Ensuite, le deuxième calcul (RESTART) soumet la barre à la meme force jusqu au temps 2. Si la relecture des informations est faite correctement lors du RESTART, les contraintes, déformations et déplacements doivent donc rester inchangés à l'issue de ce deuxième calcul.
|
||||
|
||||
\latex_brut
|
||||
Test implanté suite à la modif Herezh en réponse à la demande \#119 sur le site d Herezh :
|
||||
|
||||
\url{http://kleger.univ-ubs.fr/Herezh/issues/119} (voir ce lien pour plus d infos)
|
||||
|
||||
\fin_latex_brut
|
||||
|
||||
|
||||
---------------------------------------------------
|
||||
Informations sur les fichiers facultatifs
|
||||
---------------------------------------------------
|
||||
Le script test_RESTART_1D.posttrait1 permet de changer le TEMPSFIN pour le deuxième calcul, et de sauvegarder le .maple du premier calcul dans le but de le comparer au .maple du deuxième calcul. En fait, le .maple du premier calcul est sauvegardé en tant que .maple.ref2, de manière à provoquer une comparaison de .maple classique à l issue du deuxième calcul. Le script test_RESTART_1D.posttrait1 remplace le temps 1 par un temps de 2 dans ce .maple.ref2, de manière à rendre possible la comparaison au .maple du deuxième calcul.
|
||||
|
||||
|
||||
-------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
-------------------------------------------------------------
|
||||
au point d intégration 1 de l élément 1 :
|
||||
- contrainte : SIG11
|
||||
- déformation : EPS11
|
||||
|
||||
au noeud 2 (lieu d application de la force ponctuelle) :
|
||||
- déplacement : X1
|
||||
|
|
@ -1,19 +0,0 @@
|
|||
nom_maillage barre_1elt
|
||||
|
||||
noeuds
|
||||
2 NOEUDS
|
||||
|
||||
1 0 0 0
|
||||
2 1 0 0
|
||||
|
||||
N1 1
|
||||
N2 2
|
||||
|
||||
elements
|
||||
1 ELEMENTS
|
||||
|
||||
1 POUT BIE1 1 2
|
||||
|
||||
A_to 1 1
|
||||
|
||||
E_to 1
|
||||
|
|
@ -1,114 +0,0 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.709 #
|
||||
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
||||
debut_fichier_commande_visu # >>>>>> le mot cle: <debut_fichier_commande_visu>
|
||||
# permet au programme de se positionner au debut du fichier, il est indispensable
|
||||
|
||||
|
||||
# =================================================================================
|
||||
# || ***** demande d une visualisation maple: ***** ||
|
||||
# =================================================================================
|
||||
# un mot cle de debut (debut_visualisation_maple)
|
||||
# un mot cle de fin ( fin_visualisation_maple)
|
||||
# la seule presence de ces deux mots cle suffit a activer la visualisation maple
|
||||
debut_visualisation_maple
|
||||
|
||||
# ----------------------------- definition de la liste des increments a balayer: ----------------
|
||||
debut_list_increment # un mot cle de debut de liste
|
||||
actif 1 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# une liste d entier separee par des blancs, ou le mot cle (tous_les_increments)
|
||||
# un mot cle de fin de liste ( fin_list_increment)
|
||||
dernier_increment fin_list_increment
|
||||
|
||||
# ----------------------------- choix des maillages a visualiser: ----------------
|
||||
# la liste est facultative, par defaut la visualisation concerne le premier maillage
|
||||
debut_choix_maillage # un mot cle de debut,
|
||||
actif 0 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# une liste d entiers , puis <fin_choix_maillage>, sur une meme ligne
|
||||
1 fin_choix_maillage
|
||||
|
||||
# ----------------------------- definition des grandeurs a visualiser (maple): ----------------
|
||||
debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple),
|
||||
actif 1 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin
|
||||
debut_grandeurs_globales fin_grandeurs_globales
|
||||
# ensuite pour chaque maillage:,
|
||||
# le numero du maillage <un entier>,
|
||||
# les infos pour la visualisation eventuelle des torseurs de reaction,
|
||||
# les infos pour la visualisation eventuelle aux noeud,
|
||||
# - tout d abord les ddl principaux (position, deplacements, temperature etc.)
|
||||
# - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses
|
||||
# les infos pour la visualisation eventuelle aux elements,
|
||||
# - tout d abord les grandeurs generiques (deformations, contraintes etc.)
|
||||
# - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def
|
||||
# enfin un mot cle de fin ( fin_grandeurs_maple)
|
||||
1 # le numero de maillage
|
||||
# les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction
|
||||
# une liste de nom <chaine de caracteres >, puis <fin_list_torseur_reaction>
|
||||
debut_list_torseur_reaction fin_list_torseur_reaction
|
||||
debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds
|
||||
# debut de la liste de noeuds, puis une liste de numero de noeud <entier>, puis <fin_list_noeud>
|
||||
deb_list_noeud 2 fin_list_noeud
|
||||
# debut de la liste des ref de noeuds, puis une liste de nom <chaine de caracteres >, puis <fin_list_ref_noeud>
|
||||
deb_list_ref_noeud fin_list_ref_noeud
|
||||
# debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_noeud>
|
||||
deb_list_ddl_noeud X1 fin_list_ddl_noeud
|
||||
type_sortie_ddl_retenue= 0
|
||||
# debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_etendu_noeud>
|
||||
deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds
|
||||
# debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis <fin_list_GrandParticuliere_noeud>
|
||||
deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds
|
||||
fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds
|
||||
debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements
|
||||
# debut de la liste des elements et points d integration, une liste de (un element, un numero de pt d integ), puis <fin_list_NbElement_NbPtInteg>
|
||||
deb_list_NbElement_NbPtInteg fin_list_NbElement_NbPtInteg
|
||||
# debut de la liste des ref d elements, puis une liste de: nom <chaine de caracteres > + numero d integ, puis <fin_list_ref_element>
|
||||
deb_list_ref_element E_tout 1 fin_list_ref_element
|
||||
# debut de la liste des ref de ptinteg d elements, puis une liste de: nom <chaine de caracteres > puis <fin_list_ref_ptinteg_element>
|
||||
deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element
|
||||
# debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis <fin_list_ddl_element>
|
||||
deb_list_ddl_element SIG11 EPS11 fin_list_ddl_element
|
||||
# debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis <fin_list_GrandParticuliere_element>
|
||||
deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
|
||||
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
|
||||
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
|
||||
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
|
||||
# informations particuliere dans le cas ou il y a une animation
|
||||
# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c est au elements)
|
||||
# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent
|
||||
# si oui ou non x1 et x2 represente l accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme.
|
||||
debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier
|
||||
# un parametre de pilotage du style de sortie
|
||||
parametre_style_de_sortie 1
|
||||
# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc
|
||||
# (tangent pour les coques, suivant la fibre moyenne pour les element 1D )
|
||||
tenseur_en_absolue_ 0
|
||||
fin_grandeurs_maple # fin des grandeurs a visualiser au format maple
|
||||
|
||||
# ----------------------------- definition des parametres d animation: ----------------
|
||||
debut_animation # un mot cle de debut de liste (debut_animation)
|
||||
actif 0 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# des parametres avec des valeurs: (sur une meme ligne)
|
||||
cycleInterval 8 # cycleInterval <un reel> (indique le temps en seconde du cycle de l animation)
|
||||
fin_animation # un mot cle de fin
|
||||
|
||||
fin_visualisation_maple
|
||||
# =================================================================================
|
||||
# || fin de la visualisation maple ||
|
||||
# =================================================================================
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
fin_fichier_commande_visu # <<<<<< le mot cle <fin_fichier_commande_visu> permet
|
||||
# l arret de la lecture des commandes, apres ce mot cle, aucune commande n est lu, de plus
|
||||
# sans le mot cle de fin de fichier, le fichier n est pas valide
|
||||
|
||||
###############################################################################################
|
||||
|
|
@ -1,114 +0,0 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.709 #
|
||||
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
||||
debut_fichier_commande_visu # >>>>>> le mot cle: <debut_fichier_commande_visu>
|
||||
# permet au programme de se positionner au debut du fichier, il est indispensable
|
||||
|
||||
|
||||
# =================================================================================
|
||||
# || ***** demande d une visualisation maple: ***** ||
|
||||
# =================================================================================
|
||||
# un mot cle de debut (debut_visualisation_maple)
|
||||
# un mot cle de fin ( fin_visualisation_maple)
|
||||
# la seule presence de ces deux mots cle suffit a activer la visualisation maple
|
||||
debut_visualisation_maple
|
||||
|
||||
# ----------------------------- definition de la liste des increments a balayer: ----------------
|
||||
debut_list_increment # un mot cle de debut de liste
|
||||
actif 1 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# une liste d entier separee par des blancs, ou le mot cle (tous_les_increments)
|
||||
# un mot cle de fin de liste ( fin_list_increment)
|
||||
dernier_increment fin_list_increment
|
||||
|
||||
# ----------------------------- choix des maillages a visualiser: ----------------
|
||||
# la liste est facultative, par defaut la visualisation concerne le premier maillage
|
||||
debut_choix_maillage # un mot cle de debut,
|
||||
actif 0 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# une liste d entiers , puis <fin_choix_maillage>, sur une meme ligne
|
||||
1 fin_choix_maillage
|
||||
|
||||
# ----------------------------- definition des grandeurs a visualiser (maple): ----------------
|
||||
debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple),
|
||||
actif 1 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin
|
||||
debut_grandeurs_globales fin_grandeurs_globales
|
||||
# ensuite pour chaque maillage:,
|
||||
# le numero du maillage <un entier>,
|
||||
# les infos pour la visualisation eventuelle des torseurs de reaction,
|
||||
# les infos pour la visualisation eventuelle aux noeud,
|
||||
# - tout d abord les ddl principaux (position, deplacements, temperature etc.)
|
||||
# - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses
|
||||
# les infos pour la visualisation eventuelle aux elements,
|
||||
# - tout d abord les grandeurs generiques (deformations, contraintes etc.)
|
||||
# - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def
|
||||
# enfin un mot cle de fin ( fin_grandeurs_maple)
|
||||
1 # le numero de maillage
|
||||
# les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction
|
||||
# une liste de nom <chaine de caracteres >, puis <fin_list_torseur_reaction>
|
||||
debut_list_torseur_reaction fin_list_torseur_reaction
|
||||
debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds
|
||||
# debut de la liste de noeuds, puis une liste de numero de noeud <entier>, puis <fin_list_noeud>
|
||||
deb_list_noeud 2 fin_list_noeud
|
||||
# debut de la liste des ref de noeuds, puis une liste de nom <chaine de caracteres >, puis <fin_list_ref_noeud>
|
||||
deb_list_ref_noeud fin_list_ref_noeud
|
||||
# debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_noeud>
|
||||
deb_list_ddl_noeud X1 fin_list_ddl_noeud
|
||||
type_sortie_ddl_retenue= 0
|
||||
# debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_etendu_noeud>
|
||||
deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds
|
||||
# debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis <fin_list_GrandParticuliere_noeud>
|
||||
deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds
|
||||
fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds
|
||||
debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements
|
||||
# debut de la liste des elements et points d integration, une liste de (un element, un numero de pt d integ), puis <fin_list_NbElement_NbPtInteg>
|
||||
deb_list_NbElement_NbPtInteg fin_list_NbElement_NbPtInteg
|
||||
# debut de la liste des ref d elements, puis une liste de: nom <chaine de caracteres > + numero d integ, puis <fin_list_ref_element>
|
||||
deb_list_ref_element E_tout 1 fin_list_ref_element
|
||||
# debut de la liste des ref de ptinteg d elements, puis une liste de: nom <chaine de caracteres > puis <fin_list_ref_ptinteg_element>
|
||||
deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element
|
||||
# debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis <fin_list_ddl_element>
|
||||
deb_list_ddl_element SIG11 EPS11 fin_list_ddl_element
|
||||
# debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis <fin_list_GrandParticuliere_element>
|
||||
deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
|
||||
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
|
||||
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
|
||||
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
|
||||
# informations particuliere dans le cas ou il y a une animation
|
||||
# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c est au elements)
|
||||
# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent
|
||||
# si oui ou non x1 et x2 represente l accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme.
|
||||
debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier
|
||||
# un parametre de pilotage du style de sortie
|
||||
parametre_style_de_sortie 1
|
||||
# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc
|
||||
# (tangent pour les coques, suivant la fibre moyenne pour les element 1D )
|
||||
tenseur_en_absolue_ 0
|
||||
fin_grandeurs_maple # fin des grandeurs a visualiser au format maple
|
||||
|
||||
# ----------------------------- definition des parametres d animation: ----------------
|
||||
debut_animation # un mot cle de debut de liste (debut_animation)
|
||||
actif 0 # <0 ou 1> indique si l ordre est actif ou non
|
||||
# des parametres avec des valeurs: (sur une meme ligne)
|
||||
cycleInterval 8 # cycleInterval <un reel> (indique le temps en seconde du cycle de l animation)
|
||||
fin_animation # un mot cle de fin
|
||||
|
||||
fin_visualisation_maple
|
||||
# =================================================================================
|
||||
# || fin de la visualisation maple ||
|
||||
# =================================================================================
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
fin_fichier_commande_visu # <<<<<< le mot cle <fin_fichier_commande_visu> permet
|
||||
# l arret de la lecture des commandes, apres ce mot cle, aucune commande n est lu, de plus
|
||||
# sans le mot cle de fin de fichier, le fichier n est pas valide
|
||||
|
||||
###############################################################################################
|
||||
|
|
@ -1,71 +0,0 @@
|
|||
dimension 3
|
||||
|
||||
niveau_commentaire 3
|
||||
|
||||
TYPE_DE_CALCUL
|
||||
non_dynamique avec plus lectureCommandesVisu
|
||||
|
||||
|
||||
< barre_1elt.her
|
||||
|
||||
|
||||
les_courbes_1D
|
||||
rampe COURBEPOLYLINEAIRE_1_D
|
||||
Debut_des_coordonnees_des_points
|
||||
Coordonnee dim= 2 0. 0.
|
||||
Coordonnee dim= 2 1. 1.
|
||||
Coordonnee dim= 2 2. 1.
|
||||
Fin_des_coordonnees_des_points
|
||||
|
||||
|
||||
choix_materiaux
|
||||
E_tout MAT
|
||||
|
||||
|
||||
materiaux
|
||||
MAT ISOELAS1D
|
||||
1000. 0.3
|
||||
|
||||
|
||||
masse_volumique
|
||||
E_tout 1.
|
||||
|
||||
|
||||
sections
|
||||
E_tout 0.1
|
||||
|
||||
|
||||
charges
|
||||
N2 PONCTUELLE 10 0 0 COURBE_CHARGE: rampe ECHELLE: 1
|
||||
|
||||
|
||||
blocages
|
||||
N1 UX UY UZ
|
||||
N2 UY UZ
|
||||
|
||||
|
||||
controle
|
||||
DELTAtMINI 1.
|
||||
DELTAtMAXI 1.
|
||||
TEMPSFIN 1.
|
||||
DELTAt 1.
|
||||
ITERATIONS 15
|
||||
PRECISION 1e-3
|
||||
SAUVEGARDE 1
|
||||
MAXINCRE 9999999
|
||||
NORME Residu/Reaction_et_VarRes
|
||||
|
||||
|
||||
para_pilotage_equi_global
|
||||
|
||||
|
||||
para_syteme_lineaire
|
||||
|
||||
|
||||
para_affichage
|
||||
|
||||
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
|
||||
_fin_point_info_
|
||||
|
|
@ -1,98 +0,0 @@
|
|||
#!/usr/bin/env perl
|
||||
use strict;
|
||||
use warnings;
|
||||
use English;
|
||||
use File::Basename;
|
||||
use Cwd;
|
||||
my $NOM_PROG = basename $PROGRAM_NAME;
|
||||
#pattern d un reel pour les regex
|
||||
my $format_reel = '[+-]?[\.]?\d+[\.]?\d*(?:[eE][+-]?\d*)?';
|
||||
|
||||
#
|
||||
#
|
||||
# les buts de ce script sont :
|
||||
# 1- changer le TEMPSFIN (passer de 1 a 2)
|
||||
# 2- conserver le .maple du calcul en vue d une comparaison future avec le .maple du calcul .CVisu2
|
||||
# (ce .maple sera nomme test_RESTART_1D_princ.maple.ref2, c est-a-dire que l on va faire en sorte
|
||||
# qu il devienne un .maple de reference pour le futur deuxieme calcul .CVisu2)
|
||||
# 3- changement du temps dans le .maple test_RESTART_1D_princ.maple.ref2
|
||||
# (comme ce .maple provient du calcul .CVisu1 pour un TEMPSFIN egal a 1, on modifie arbitrairement le
|
||||
# temps ecrit dans le .maple pour pouvoir le comparer proprement au calcul .CVisu2 pour un TEMPSFIN egal
|
||||
# a 2. Sinon, il y a aura forcement ECHEC de la comparaison de maple pour le calcul .CVisu2 a cause d une
|
||||
# difference de grandeur en colonne 1)
|
||||
#
|
||||
#
|
||||
|
||||
#chemin absolu vers l executable Herezh
|
||||
my $hzchemin = shift(@ARGV);
|
||||
|
||||
#nom du fichier .info
|
||||
my $finfo = shift(@ARGV);
|
||||
|
||||
|
||||
#------------------------
|
||||
#changement TEMPSFIN
|
||||
#------------------------
|
||||
#ouverture .info d origine
|
||||
open(FIC1, "<$finfo") or do {
|
||||
print "Erreur (prog:$NOM_PROG) : probleme en ouverture du fichier $finfo ($!)...\n\n";
|
||||
print "resultat posttrait : ECHEC\n";
|
||||
};
|
||||
open(FIC2, ">$finfo.tmp") or do {
|
||||
print "Erreur (prog:$NOM_PROG) : probleme en ouverture du fichier $finfo.tmp ($!)...\n\n";
|
||||
print "resultat posttrait : ECHEC\n";
|
||||
};
|
||||
while(<FIC1>) {
|
||||
|
||||
#modif TEMPSFIN
|
||||
$_ = 'TEMPSFIN 2.' if(/^\s*TEMPSFIN /);
|
||||
|
||||
#ecriture dans le fichier temporaire
|
||||
print FIC2;
|
||||
}
|
||||
close(FIC2);
|
||||
close(FIC1);
|
||||
|
||||
#remplacement du .info d origine par le .info temporaire
|
||||
system("mv -f $finfo.tmp $finfo");
|
||||
|
||||
|
||||
#------------------------
|
||||
#sauvegarde du .maple en .maple.ref2
|
||||
#------------------------
|
||||
(my $fmaple = $finfo) =~ s/.info$/_princ.maple/;
|
||||
#verif existence du fichier .maple
|
||||
(-e $fmaple) or do {
|
||||
print "Erreur (prog:$NOM_PROG) : fichier maple $fmaple introuvable...\n\n";
|
||||
print "resultat posttrait : ECHEC\n";
|
||||
};
|
||||
#recopie du .maple sous le nom sauvegarde_maple_CVisu1.maple
|
||||
system("cp $fmaple $fmaple.ref2");
|
||||
|
||||
#verif existence fichier sauvegarde_maple_CVisu1.maple
|
||||
(-e "$fmaple.ref2") or do {
|
||||
print "Erreur (prog:$NOM_PROG) : fichier $fmaple.ref2 n a pas ete cree...\n\n";
|
||||
print "resultat posttrait : ECHEC\n";
|
||||
};
|
||||
|
||||
#modif temps (on change le temps de 1 en 2. Pour info, le temps est en colonne 1 mais egalement en colonne 6 pour ce test)
|
||||
open(FIC, "<$fmaple.ref2");
|
||||
open(FIC2, ">$fmaple.ref2.tmp");
|
||||
while(<FIC>) {
|
||||
if(/^\s*$format_reel/o) {
|
||||
@_ = split;
|
||||
$_[0] = 2.;
|
||||
$_[5] = 2.;
|
||||
print FIC2 "@_\n";
|
||||
}
|
||||
else {
|
||||
print FIC2;
|
||||
}
|
||||
}
|
||||
close(FIC2);
|
||||
close(FIC);
|
||||
system("mv -f $fmaple.ref2.tmp $fmaple.ref2");
|
||||
|
||||
|
||||
#reussite du script
|
||||
print "resultat posttrait : OK\n";
|
||||
|
|
@ -1,100 +0,0 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.770 #
|
||||
# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
||||
# entete des donnees : informations gererales: on trouve successivement:
|
||||
# >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs
|
||||
# precedes du numero de colonne entre crochet
|
||||
# >> le nombre de maillages m, et dimension de l'espace de travail
|
||||
# puis pour chaque maillage,
|
||||
# >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit
|
||||
# correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions
|
||||
# des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment
|
||||
# ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds
|
||||
# >> le nombre de ref de noeuds, nombre total de grandeurs associees
|
||||
# puis pour chaque maillage
|
||||
# le numero de maillage puis pour chaque reference de noeuds
|
||||
# le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites
|
||||
# puis entre crochet la position suivi de la signification de chaque grandeur
|
||||
#
|
||||
# ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti
|
||||
# >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees
|
||||
# puis pour chaque maillage
|
||||
# le numero de maillage puis pour chaque reference d'element et de pti
|
||||
# le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites
|
||||
# puis entre crochet la position suivi de la signification de chaque grandeur
|
||||
#
|
||||
# puis pour chaque maillage
|
||||
# >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie ,
|
||||
# puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees
|
||||
# et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere)
|
||||
# precedes du numero de colonne correspondant entre crochet
|
||||
# puis pour chaque maillage
|
||||
# >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie ,
|
||||
# les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant
|
||||
# a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques,
|
||||
# on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl
|
||||
# chacun precedes du numero de colonne entre crochet
|
||||
# puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres),
|
||||
# suivi entre crochet, de la plage des numeros de colonnes, correspondant
|
||||
# chacun sur une ligne differentes
|
||||
# ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D:
|
||||
# tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2)
|
||||
# en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3)
|
||||
# non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3)
|
||||
# ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements,
|
||||
# ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique
|
||||
# ** ( cf: parametre_style_de_sortie = 0)
|
||||
|
||||
#====================================================================
|
||||
#|| recapitulatif des differentes grandeurs par colonne ||
|
||||
#====================================================================
|
||||
#---------------------------------- grandeur globales ------------------------------------
|
||||
#0 (nombre de grandeurs globales)
|
||||
#---------------------------------- maillage et dimension --------------------------------
|
||||
#1 3 (nombre de maillages et dimension)
|
||||
#---------------------------------- torseurs de reactions --------------------------------
|
||||
#0 0 (nombre de torseurs et nombre total de grandeurs associees)
|
||||
#
|
||||
#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref ---------------
|
||||
#0 0 (nombre de ref de noeud et nombre total de grandeurs associees)
|
||||
#
|
||||
#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref ---------------
|
||||
#0 0 (nombre de ref d'element et nombre total de grandeurs associees)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#1 1 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
# noeud_2 [2]X [3]Y [4]Z [5]X1 ;
|
||||
#
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 2 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# reference E_tout : element_1 pt_integ_1: [7]X [8]Y [9]Z [10] SIG11 [11] EPS11
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| fin du recapitulatif des differentes grandeurs ||
|
||||
#====================================================================
|
||||
|
||||
# ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant
|
||||
# a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun
|
||||
# correspondant a un maillage. On trouve pour chaque enregistrement successivement :
|
||||
# s'il y a des grandeurs globales: le temps puis les grandeurs globales,
|
||||
# puis s'il y a des torseurs de reaction :
|
||||
# de nouveau le temps, les composantes de la resultante puis les composantes du moments
|
||||
# donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels
|
||||
# puis s'il y a des grandeurs aux noeuds: de nouveau le temps
|
||||
# les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud
|
||||
# puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds
|
||||
# puis s'il y a des grandeur aux elements:
|
||||
# le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques)
|
||||
# suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration
|
||||
# correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes
|
||||
# puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs
|
||||
# etc. pour tous les points d'integration - element
|
||||
|
||||
1.000000000000e+00 1.129446618973e+00 0.000000000000e+00 0.000000000000e+00 1.129446618973e+00 1.000000000000e+00 5.647233094863e-01 0.000000000000e+00 0.000000000000e+00 1.080428555928e+02 1.080428555928e-01
|
||||
|
|
@ -1,64 +0,0 @@
|
|||
------------------------------------------------------
|
||||
Auteur
|
||||
------------------------------------------------------
|
||||
Julien Troufflard (troufflard@univ-ubs.fr)
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------
|
||||
contrainte_individuelle_a_chaque_loi_a_t
|
||||
LOI_ADDITIVE_EN_SIGMA
|
||||
ISOELAS
|
||||
TEMP
|
||||
temperature
|
||||
dilatation thermique
|
||||
dilatation_thermique
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------
|
||||
test de la combinaison "dilatation_thermique + LOI_ADDITIVE"
|
||||
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
**
|
||||
** ----- remarques a supprimer une fois le bug resolu -----
|
||||
**
|
||||
**
|
||||
** MEME PROBLEME QUE LE TEST EN ATTENTE : Test_R_sortie_contrainte_individuelle_a_chaque_loi_a_t_1/
|
||||
** sera OK une fois resolu le probleme sur la sortie contrainte_individuelle_a_chaque_loi_a_t
|
||||
**
|
||||
** le .maple.ref1 ici present est celui obtenu avec la version 6.768
|
||||
**
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
|
||||
------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------
|
||||
traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE). La température initialement à 20°C évolue linéairement en fonction du temps pour atteindre 0°C à la fin du calcul
|
||||
|
||||
Loi de comportement :
|
||||
- loi additive en contrainte : 2 lois élastiques ISOELAS
|
||||
- dilatation thermique
|
||||
|
||||
|
||||
-------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
-------------------------------------------------------------
|
||||
pour le point d integration 1 de l element 1 :
|
||||
- contrainte totale : SIG11
|
||||
- contribution loi 1 : toutes les composantes
|
||||
- contribution loi 2 : toutes les composantes
|
||||
- déformations : EPS11, EPS22, EPS33
|
||||
- volume du point d intégration : VOLUME_PTI
|
||||
|
||||
au noeud 8 :
|
||||
- positions : X1, X2, X3
|
||||
- température : TEMP
|
||||
|
||||
remarque : seuls les positions permettent de rendre compte de la dilatation thermique
|
||||
(car les deformations EPS11, etc... ne representent que la partie mecanique
|
||||
du tenseur des deformations)
|
||||
|
|
@ -1,116 +0,0 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.687 #
|
||||
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
||||
debut_fichier_commande_visu # >>>>>> le mot cle: <debut_fichier_commande_visu>
|
||||
# permet au programme de se positionner au debut du fichier, il est indispensable
|
||||
|
||||
|
||||
# =================================================================================
|
||||
# || ***** demande d'une visualisation maple: ***** ||
|
||||
# =================================================================================
|
||||
# un mot cle de debut (debut_visualisation_maple)
|
||||
# un mot cle de fin ( fin_visualisation_maple)
|
||||
# la seule presence de ces deux mots cle suffit a activer la visualisation maple
|
||||
debut_visualisation_maple
|
||||
|
||||
# ----------------------------- definition de la liste des increments a balayer: ----------------
|
||||
debut_list_increment # un mot cle de debut de liste
|
||||
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# une liste d'entier separee par des blancs, ou le mot cle (tous_les_increments)
|
||||
# un mot cle de fin de liste ( fin_list_increment)
|
||||
dernier_increment fin_list_increment
|
||||
|
||||
# ----------------------------- choix des maillages a visualiser: ----------------
|
||||
# la liste est facultative, par defaut la visualisation concerne le premier maillage
|
||||
debut_choix_maillage # un mot cle de debut,
|
||||
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# une liste d'entiers , puis <fin_choix_maillage>, sur une meme ligne
|
||||
1 fin_choix_maillage
|
||||
|
||||
# ----------------------------- definition des grandeurs a visualiser (maple): ----------------
|
||||
debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple),
|
||||
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin
|
||||
debut_grandeurs_globales fin_grandeurs_globales
|
||||
# ensuite pour chaque maillage:,
|
||||
# le numero du maillage <un entier>,
|
||||
# les infos pour la visualisation eventuelle des torseurs de reaction,
|
||||
# les infos pour la visualisation eventuelle aux noeud,
|
||||
# - tout d'abord les ddl principaux (position, deplacements, temperature etc.)
|
||||
# - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses
|
||||
# les infos pour la visualisation eventuelle aux elements,
|
||||
# - tout d'abord les grandeurs generiques (deformations, contraintes etc.)
|
||||
# - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def
|
||||
# enfin un mot cle de fin ( fin_grandeurs_maple)
|
||||
1 # le numero de maillage
|
||||
# les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction
|
||||
# une liste de nom <chaine de caracteres >, puis <fin_list_torseur_reaction>
|
||||
debut_list_torseur_reaction fin_list_torseur_reaction
|
||||
debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds
|
||||
# debut de la liste de noeuds, puis une liste de numero de noeud <entier>, puis <fin_list_noeud>
|
||||
deb_list_noeud 8 fin_list_noeud
|
||||
# debut de la liste des ref de noeuds, puis une liste de nom <chaine de caracteres >, puis <fin_list_ref_noeud>
|
||||
deb_list_ref_noeud fin_list_ref_noeud
|
||||
# debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_noeud>
|
||||
deb_list_ddl_noeud X1 X2 X3 TEMP fin_list_ddl_noeud
|
||||
type_sortie_ddl_retenue= 0
|
||||
# debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_etendu_noeud>
|
||||
deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds
|
||||
# debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis <fin_list_GrandParticuliere_noeud>
|
||||
deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds
|
||||
fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds
|
||||
debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements
|
||||
# debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis <fin_list_NbElement_NbPtInteg>
|
||||
deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg
|
||||
# debut de la liste des ref d'elements, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_element>
|
||||
deb_list_ref_element fin_list_ref_element
|
||||
# debut de la liste des ref de ptinteg d'elements, puis une liste de: nom <chaine de caracteres > puis <fin_list_ref_ptinteg_element>
|
||||
deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element
|
||||
# debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis <fin_list_ddl_element>
|
||||
deb_list_ddl_element EPS11 EPS22 EPS33 SIG11 fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements
|
||||
# debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis <fin_list_GrandParticuliere_element>
|
||||
deb_list_GrandParticuliere_element VOLUME_PTI contrainte_individuelle_a_chaque_loi_a_t fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
|
||||
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
|
||||
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
|
||||
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
|
||||
# informations particuliere dans le cas ou il y a une animation
|
||||
# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements)
|
||||
# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent
|
||||
# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme.
|
||||
debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier
|
||||
# un parametre de pilotage du style de sortie
|
||||
parametre_style_de_sortie 1
|
||||
# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc
|
||||
# (tangent pour les coques, suivant la fibre moyenne pour les element 1D )
|
||||
tenseur_en_absolue_ 0
|
||||
fin_grandeurs_maple # fin des grandeurs a visualiser au format maple
|
||||
|
||||
# ----------------------------- definition des parametres d'animation: ----------------
|
||||
debut_animation # un mot cle de debut de liste (debut_animation)
|
||||
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# des parametres avec des valeurs: (sur une meme ligne)
|
||||
cycleInterval 8 # cycleInterval <un reel> (indique le temps en seconde du cycle de l'animation)
|
||||
fin_animation # un mot cle de fin
|
||||
|
||||
fin_visualisation_maple
|
||||
# =================================================================================
|
||||
# || fin de la visualisation maple ||
|
||||
# =================================================================================
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
fin_fichier_commande_visu # <<<<<< le mot cle <fin_fichier_commande_visu> permet
|
||||
# l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus
|
||||
# sans le mot cle de fin de fichier, le fichier n'est pas valide
|
||||
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
|
@ -1,91 +0,0 @@
|
|||
dimension 3
|
||||
|
||||
|
||||
niveau_commentaire 3
|
||||
|
||||
|
||||
TYPE_DE_CALCUL
|
||||
|
||||
non_dynamique avec plus lectureCommandesVisu
|
||||
|
||||
|
||||
< cube_1elt_npti8.her
|
||||
< cube_1elt_npti8.lis
|
||||
|
||||
|
||||
les_courbes_1D
|
||||
|
||||
|
||||
choix_materiaux
|
||||
|
||||
E_tout MAT
|
||||
E_tout loi_DILATATION
|
||||
|
||||
materiaux
|
||||
|
||||
MAT LOI_ADDITIVE_EN_SIGMA
|
||||
ISOELAS
|
||||
1000. 0.4
|
||||
|
||||
ISOELAS
|
||||
2000. 0.4
|
||||
fin_liste_lois_elementaires
|
||||
|
||||
|
||||
loi_DILATATION LOI_ISO_THERMO
|
||||
alphaT= 1.e-3 lambda= 1.e-10 cp= 1.e-10
|
||||
fin_thermique_isotrope
|
||||
|
||||
|
||||
|
||||
masse_volumique
|
||||
E_tout 1.
|
||||
|
||||
dilatation_thermique
|
||||
E_tout 1.
|
||||
|
||||
|
||||
charges
|
||||
|
||||
|
||||
blocages
|
||||
|
||||
#traction uniaxiale suivant X par deplacement impose
|
||||
N_droit UX
|
||||
N_arriere UY
|
||||
N_bas UZ
|
||||
N_gauche 'UX=0.1'
|
||||
|
||||
#temperature initialement a 20 passe a 0 (delta = -20)
|
||||
N_tout 'TEMP= -20'
|
||||
|
||||
|
||||
initialisation
|
||||
N_tout 'TEMP= 20'
|
||||
|
||||
controle
|
||||
|
||||
DELTAtMINI 1.e-7
|
||||
DELTAtMAXI 0.1
|
||||
TEMPSFIN 1
|
||||
DELTAt 0.1
|
||||
ITERATIONS 15
|
||||
PRECISION 1e-3
|
||||
SAUVEGARDE DERNIER_CALCUL
|
||||
MAXINCRE 99999999
|
||||
NORME Residu/Reaction_et_VarRes
|
||||
|
||||
|
||||
para_pilotage_equi_global
|
||||
|
||||
|
||||
para_syteme_lineaire
|
||||
|
||||
|
||||
para_affichage
|
||||
|
||||
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
|
||||
_fin_point_info_
|
||||
|
|
@ -1,18 +0,0 @@
|
|||
nom_maillage cube_1elt_npti8
|
||||
|
||||
noeuds
|
||||
8 NOEUDS
|
||||
|
||||
1 0 0 0
|
||||
2 0 0 1
|
||||
3 0 1 0
|
||||
4 0 1 1
|
||||
5 1 0 0
|
||||
6 1 0 1
|
||||
7 1 1 0
|
||||
8 1 1 1
|
||||
|
||||
elements
|
||||
1 ELEMENTS
|
||||
|
||||
1 HEXAEDRE LINEAIRE 5 7 3 1 6 8 4 2
|
||||
|
|
@ -1,158 +0,0 @@
|
|||
# CORRESPONDANCE label => position
|
||||
#
|
||||
# gauche => X=1
|
||||
#
|
||||
# droit => X=0
|
||||
#
|
||||
# avant => Y=1
|
||||
#
|
||||
# arriere => Y=0
|
||||
#
|
||||
# haut => Z=1
|
||||
#
|
||||
# bas => Z=0
|
||||
#
|
||||
|
||||
N_avant 3 4 7 8
|
||||
|
||||
N_arriere 1 2 5 6
|
||||
|
||||
N_haut 2 4 6 8
|
||||
|
||||
N_bas 1 3 5 7
|
||||
|
||||
N_droit 1 2 3 4
|
||||
|
||||
N_gauche 5 6 7 8
|
||||
|
||||
N_haut_droit 2 4
|
||||
|
||||
N_haut_gauche 6 8
|
||||
|
||||
N_haut_arriere 2 6
|
||||
|
||||
N_haut_avant 4 8
|
||||
|
||||
N_bas_droit 1 3
|
||||
|
||||
N_bas_gauche 5 7
|
||||
|
||||
N_bas_arriere 1 5
|
||||
|
||||
N_bas_avant 3 7
|
||||
|
||||
N_arriere_droit 1 2
|
||||
|
||||
N_arriere_gauche 5 6
|
||||
|
||||
N_avant_droit 3 4
|
||||
|
||||
N_avant_gauche 7 8
|
||||
|
||||
N_haut_arriere_droit 2
|
||||
|
||||
N_haut_arriere_gauche 6
|
||||
|
||||
N_haut_avant_gauche 8
|
||||
|
||||
N_haut_avant_droit 4
|
||||
|
||||
N_bas_arriere_droit 1
|
||||
|
||||
N_bas_arriere_gauche 5
|
||||
|
||||
N_bas_avant_gauche 7
|
||||
|
||||
N_bas_avant_droit 3
|
||||
|
||||
N_to 1 2 3 4 5 6 7 8
|
||||
|
||||
E_avant 1
|
||||
|
||||
E_derriere 1
|
||||
|
||||
E_haut 1
|
||||
|
||||
E_bas 1
|
||||
|
||||
E_droit 1
|
||||
|
||||
E_gauche 1
|
||||
|
||||
E_haut_droit 1
|
||||
|
||||
E_haut_gauche 1
|
||||
|
||||
E_haut_arriere 1
|
||||
|
||||
E_haut_avant 1
|
||||
|
||||
E_bas_droit 1
|
||||
|
||||
E_bas_gauche 1
|
||||
|
||||
E_bas_arriere 1
|
||||
|
||||
E_bas_avant 1
|
||||
|
||||
E_arriere_droit 1
|
||||
|
||||
E_arriere_gauche 1
|
||||
|
||||
E_avant_droit 1
|
||||
|
||||
E_avant_gauche 1
|
||||
|
||||
E_haut_arriere_droit 1
|
||||
|
||||
E_haut_arriere_gauche 1
|
||||
|
||||
E_haut_avant_gauche 1
|
||||
|
||||
E_haut_avant_droit 1
|
||||
|
||||
E_bas_arriere_droit 1
|
||||
|
||||
E_bas_arriere_gauche 1
|
||||
|
||||
E_bas_avant_gauche 1
|
||||
|
||||
E_bas_avant_droit 1
|
||||
|
||||
E_to 1
|
||||
|
||||
F_avant 1 5
|
||||
|
||||
F_arriere 1 2
|
||||
|
||||
F_haut 1 4
|
||||
|
||||
F_bas 1 1
|
||||
|
||||
F_droit 1 6
|
||||
|
||||
F_gauche 1 3
|
||||
|
||||
A_haut_droit 1 11
|
||||
|
||||
A_haut_gauche 1 9
|
||||
|
||||
A_haut_arriere 1 12
|
||||
|
||||
A_haut_avant 1 10
|
||||
|
||||
A_bas_droit 1 3
|
||||
|
||||
A_bas_gauche 1 1
|
||||
|
||||
A_bas_arriere 1 4
|
||||
|
||||
A_bas_avant 1 2
|
||||
|
||||
A_arriere_droit 1 8
|
||||
|
||||
A_arriere_gauche 1 5
|
||||
|
||||
A_avant_droit 1 7
|
||||
|
||||
A_avant_gauche 1 6
|
||||
|
|
@ -1,54 +0,0 @@
|
|||
------------------------------------------------------
|
||||
Auteur
|
||||
------------------------------------------------------
|
||||
Julien Troufflard (troufflard@univ-ubs.fr)
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------
|
||||
Traction biaxiale
|
||||
BLOQUAGE_RELATIF_
|
||||
TEMPS_MINI
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------
|
||||
Vérification du fonctionnement de la combinaison BLOQUAGE_RELATIF_ + TEMPS_MINI dans une condition de déplacement imposé
|
||||
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
**
|
||||
** ----- remarques a supprimer une fois le bug resolu -----
|
||||
**
|
||||
** ce calcul tourne avec la version 6.726
|
||||
** il ne tourne pas avec la version 6.727
|
||||
**
|
||||
** selon redmine, la version 6.727 correspond aux modifs suivantes :
|
||||
** <<< introduction des chargements en champs de valeurs >>>
|
||||
**
|
||||
**
|
||||
** le .maple.ref1 ici present est celui obtenu avec la version 6.726
|
||||
**
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
** @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @
|
||||
|
||||
|
||||
------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------
|
||||
Calcul de traction biaxiale sur un cube unitaire.
|
||||
|
||||
La gestion du chargement se fait via des courbes de charge :
|
||||
1) Entre le temps 0 et 1, c est de la traction uniaxiale classique selon X car la condition sur UY de la face Y=1 est inactive jusqu au temps 1 inclus (TEMPS_MINI= 1).
|
||||
|
||||
2) Entre le temps 1 et 2, la courbe dpl_UX permet de figer la face X=1 dans sa position actuelle (X=1.02). La condition sur la face Y=1 devient active, c'est-à-dire un déplacement UY=0.02 (piloté par la courbe dpl_UY_relatif) depuis la position atteinte au temps 1 (BLOQUAGE_RELATIF_).
|
||||
|
||||
|
||||
-------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
-------------------------------------------------------------
|
||||
les positions finales du noeud 8 : X1, X2, X3
|
||||
|
||||
contraintes au point d intégration 1 : SIG11, SIG22, SIG33
|
||||
|
|
@ -1,18 +0,0 @@
|
|||
nom_maillage cube_1elt_npti8
|
||||
|
||||
noeuds
|
||||
8 NOEUDS
|
||||
|
||||
1 0 0 0
|
||||
2 0 0 1
|
||||
3 0 1 0
|
||||
4 0 1 1
|
||||
5 1 0 0
|
||||
6 1 0 1
|
||||
7 1 1 0
|
||||
8 1 1 1
|
||||
|
||||
elements
|
||||
1 ELEMENTS
|
||||
|
||||
1 HEXAEDRE LINEAIRE 5 7 3 1 6 8 4 2
|
||||
|
|
@ -1,158 +0,0 @@
|
|||
# CORRESPONDANCE label => position
|
||||
#
|
||||
# gauche => X=1
|
||||
#
|
||||
# droit => X=0
|
||||
#
|
||||
# avant => Y=1
|
||||
#
|
||||
# arriere => Y=0
|
||||
#
|
||||
# haut => Z=1
|
||||
#
|
||||
# bas => Z=0
|
||||
#
|
||||
|
||||
N_avant 3 4 7 8
|
||||
|
||||
N_arriere 1 2 5 6
|
||||
|
||||
N_haut 2 4 6 8
|
||||
|
||||
N_bas 1 3 5 7
|
||||
|
||||
N_droit 1 2 3 4
|
||||
|
||||
N_gauche 5 6 7 8
|
||||
|
||||
N_haut_droit 2 4
|
||||
|
||||
N_haut_gauche 6 8
|
||||
|
||||
N_haut_arriere 2 6
|
||||
|
||||
N_haut_avant 4 8
|
||||
|
||||
N_bas_droit 1 3
|
||||
|
||||
N_bas_gauche 5 7
|
||||
|
||||
N_bas_arriere 1 5
|
||||
|
||||
N_bas_avant 3 7
|
||||
|
||||
N_arriere_droit 1 2
|
||||
|
||||
N_arriere_gauche 5 6
|
||||
|
||||
N_avant_droit 3 4
|
||||
|
||||
N_avant_gauche 7 8
|
||||
|
||||
N_haut_arriere_droit 2
|
||||
|
||||
N_haut_arriere_gauche 6
|
||||
|
||||
N_haut_avant_gauche 8
|
||||
|
||||
N_haut_avant_droit 4
|
||||
|
||||
N_bas_arriere_droit 1
|
||||
|
||||
N_bas_arriere_gauche 5
|
||||
|
||||
N_bas_avant_gauche 7
|
||||
|
||||
N_bas_avant_droit 3
|
||||
|
||||
N_to 1 2 3 4 5 6 7 8
|
||||
|
||||
E_avant 1
|
||||
|
||||
E_derriere 1
|
||||
|
||||
E_haut 1
|
||||
|
||||
E_bas 1
|
||||
|
||||
E_droit 1
|
||||
|
||||
E_gauche 1
|
||||
|
||||
E_haut_droit 1
|
||||
|
||||
E_haut_gauche 1
|
||||
|
||||
E_haut_arriere 1
|
||||
|
||||
E_haut_avant 1
|
||||
|
||||
E_bas_droit 1
|
||||
|
||||
E_bas_gauche 1
|
||||
|
||||
E_bas_arriere 1
|
||||
|
||||
E_bas_avant 1
|
||||
|
||||
E_arriere_droit 1
|
||||
|
||||
E_arriere_gauche 1
|
||||
|
||||
E_avant_droit 1
|
||||
|
||||
E_avant_gauche 1
|
||||
|
||||
E_haut_arriere_droit 1
|
||||
|
||||
E_haut_arriere_gauche 1
|
||||
|
||||
E_haut_avant_gauche 1
|
||||
|
||||
E_haut_avant_droit 1
|
||||
|
||||
E_bas_arriere_droit 1
|
||||
|
||||
E_bas_arriere_gauche 1
|
||||
|
||||
E_bas_avant_gauche 1
|
||||
|
||||
E_bas_avant_droit 1
|
||||
|
||||
E_to 1
|
||||
|
||||
F_avant 1 5
|
||||
|
||||
F_arriere 1 2
|
||||
|
||||
F_haut 1 4
|
||||
|
||||
F_bas 1 1
|
||||
|
||||
F_droit 1 6
|
||||
|
||||
F_gauche 1 3
|
||||
|
||||
A_haut_droit 1 11
|
||||
|
||||
A_haut_gauche 1 9
|
||||
|
||||
A_haut_arriere 1 12
|
||||
|
||||
A_haut_avant 1 10
|
||||
|
||||
A_bas_droit 1 3
|
||||
|
||||
A_bas_gauche 1 1
|
||||
|
||||
A_bas_arriere 1 4
|
||||
|
||||
A_bas_avant 1 2
|
||||
|
||||
A_arriere_droit 1 8
|
||||
|
||||
A_arriere_gauche 1 5
|
||||
|
||||
A_avant_droit 1 7
|
||||
|
||||
A_avant_gauche 1 6
|
||||
|
|
@ -1,116 +0,0 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.687 #
|
||||
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
||||
debut_fichier_commande_visu # >>>>>> le mot cle: <debut_fichier_commande_visu>
|
||||
# permet au programme de se positionner au debut du fichier, il est indispensable
|
||||
|
||||
|
||||
# =================================================================================
|
||||
# || ***** demande d'une visualisation maple: ***** ||
|
||||
# =================================================================================
|
||||
# un mot cle de debut (debut_visualisation_maple)
|
||||
# un mot cle de fin ( fin_visualisation_maple)
|
||||
# la seule presence de ces deux mots cle suffit a activer la visualisation maple
|
||||
debut_visualisation_maple
|
||||
|
||||
# ----------------------------- definition de la liste des increments a balayer: ----------------
|
||||
debut_list_increment # un mot cle de debut de liste
|
||||
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# une liste d'entier separee par des blancs, ou le mot cle (tous_les_increments)
|
||||
# un mot cle de fin de liste ( fin_list_increment)
|
||||
dernier_increment fin_list_increment
|
||||
|
||||
# ----------------------------- choix des maillages a visualiser: ----------------
|
||||
# la liste est facultative, par defaut la visualisation concerne le premier maillage
|
||||
debut_choix_maillage # un mot cle de debut,
|
||||
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# une liste d'entiers , puis <fin_choix_maillage>, sur une meme ligne
|
||||
1 fin_choix_maillage
|
||||
|
||||
# ----------------------------- definition des grandeurs a visualiser (maple): ----------------
|
||||
debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple),
|
||||
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin
|
||||
debut_grandeurs_globales fin_grandeurs_globales
|
||||
# ensuite pour chaque maillage:,
|
||||
# le numero du maillage <un entier>,
|
||||
# les infos pour la visualisation eventuelle des torseurs de reaction,
|
||||
# les infos pour la visualisation eventuelle aux noeud,
|
||||
# - tout d'abord les ddl principaux (position, deplacements, temperature etc.)
|
||||
# - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses
|
||||
# les infos pour la visualisation eventuelle aux elements,
|
||||
# - tout d'abord les grandeurs generiques (deformations, contraintes etc.)
|
||||
# - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def
|
||||
# enfin un mot cle de fin ( fin_grandeurs_maple)
|
||||
1 # le numero de maillage
|
||||
# les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction
|
||||
# une liste de nom <chaine de caracteres >, puis <fin_list_torseur_reaction>
|
||||
debut_list_torseur_reaction fin_list_torseur_reaction
|
||||
debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds
|
||||
# debut de la liste de noeuds, puis une liste de numero de noeud <entier>, puis <fin_list_noeud>
|
||||
deb_list_noeud fin_list_noeud
|
||||
# debut de la liste des ref de noeuds, puis une liste de nom <chaine de caracteres >, puis <fin_list_ref_noeud>
|
||||
deb_list_ref_noeud N_haut_avant_gauche fin_list_ref_noeud
|
||||
# debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_noeud>
|
||||
deb_list_ddl_noeud X1 X2 X3 fin_list_ddl_noeud
|
||||
type_sortie_ddl_retenue= 0
|
||||
# debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_etendu_noeud>
|
||||
deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds
|
||||
# debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis <fin_list_GrandParticuliere_noeud>
|
||||
deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds
|
||||
fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds
|
||||
debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements
|
||||
# debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis <fin_list_NbElement_NbPtInteg>
|
||||
deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg
|
||||
# debut de la liste des ref d'elements, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_element>
|
||||
deb_list_ref_element fin_list_ref_element
|
||||
# debut de la liste des ref de ptinteg d'elements, puis une liste de: nom <chaine de caracteres > puis <fin_list_ref_ptinteg_element>
|
||||
deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element
|
||||
# debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis <fin_list_ddl_element>
|
||||
deb_list_ddl_element SIG11 SIG22 SIG33 fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements
|
||||
# debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis <fin_list_GrandParticuliere_element>
|
||||
deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
|
||||
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
|
||||
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
|
||||
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
|
||||
# informations particuliere dans le cas ou il y a une animation
|
||||
# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements)
|
||||
# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent
|
||||
# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme.
|
||||
debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier
|
||||
# un parametre de pilotage du style de sortie
|
||||
parametre_style_de_sortie 1
|
||||
# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc
|
||||
# (tangent pour les coques, suivant la fibre moyenne pour les element 1D )
|
||||
tenseur_en_absolue_ 0
|
||||
fin_grandeurs_maple # fin des grandeurs a visualiser au format maple
|
||||
|
||||
# ----------------------------- definition des parametres d'animation: ----------------
|
||||
debut_animation # un mot cle de debut de liste (debut_animation)
|
||||
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
|
||||
# des parametres avec des valeurs: (sur une meme ligne)
|
||||
cycleInterval 8 # cycleInterval <un reel> (indique le temps en seconde du cycle de l'animation)
|
||||
fin_animation # un mot cle de fin
|
||||
|
||||
fin_visualisation_maple
|
||||
# =================================================================================
|
||||
# || fin de la visualisation maple ||
|
||||
# =================================================================================
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
fin_fichier_commande_visu # <<<<<< le mot cle <fin_fichier_commande_visu> permet
|
||||
# l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus
|
||||
# sans le mot cle de fin de fichier, le fichier n'est pas valide
|
||||
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
|
@ -1,83 +0,0 @@
|
|||
dimension 3
|
||||
|
||||
|
||||
niveau_commentaire 3
|
||||
|
||||
|
||||
TYPE_DE_CALCUL
|
||||
non_dynamique avec plus lectureCommandesVisu
|
||||
|
||||
|
||||
< cube_1elt_npti8.her
|
||||
< cube_1elt_npti8.lis
|
||||
|
||||
|
||||
les_courbes_1D
|
||||
|
||||
dpl_UX COURBEPOLYLINEAIRE_1_D
|
||||
Debut_des_coordonnees_des_points
|
||||
Coordonnee dim= 2 0. 0.
|
||||
Coordonnee dim= 2 1. 0.02
|
||||
Coordonnee dim= 2 2. 0.02
|
||||
Fin_des_coordonnees_des_points
|
||||
|
||||
dpl_UY_relatif COURBEPOLYLINEAIRE_1_D
|
||||
Debut_des_coordonnees_des_points
|
||||
Coordonnee dim= 2 0. 0.
|
||||
Coordonnee dim= 2 1. 0.
|
||||
Coordonnee dim= 2 2. 0.02
|
||||
Fin_des_coordonnees_des_points
|
||||
|
||||
|
||||
|
||||
choix_materiaux
|
||||
E_tout MAT
|
||||
|
||||
|
||||
materiaux
|
||||
MAT ISOELAS
|
||||
210000. 0.3
|
||||
|
||||
|
||||
masse_volumique
|
||||
E_tout 1.
|
||||
|
||||
|
||||
charges
|
||||
|
||||
|
||||
blocages
|
||||
|
||||
N_droit UX
|
||||
N_arriere UY
|
||||
N_bas UZ
|
||||
N_gauche 'UX= COURBE_CHARGE: dpl_UX ECHELLE: 1.'
|
||||
N_avant 'UY= COURBE_CHARGE: dpl_UY_relatif ECHELLE: 1.' TEMPS_MINI= 1. BLOQUAGE_RELATIF_
|
||||
|
||||
|
||||
|
||||
controle
|
||||
|
||||
DELTAtMINI 1.e-7
|
||||
DELTAtMAXI 0.1
|
||||
TEMPSFIN 2.
|
||||
DELTAt 0.1
|
||||
ITERATIONS 15
|
||||
PRECISION 1e-3
|
||||
SAUVEGARDE DERNIER_CALCUL
|
||||
MAXINCRE 9999999
|
||||
NORME Residu/Reaction_et_VarRes
|
||||
|
||||
|
||||
para_pilotage_equi_global
|
||||
|
||||
para_syteme_lineaire
|
||||
|
||||
para_affichage
|
||||
|
||||
|
||||
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
|
||||
_fin_point_info_
|
||||
|
|
@ -1,99 +0,0 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.726 #
|
||||
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
||||
# entete des donnees : informations gererales: on trouve successivement:
|
||||
# >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs
|
||||
# precedes du numero de colonne entre crochet
|
||||
# >> le nombre de maillages m, et dimension de l'espace de travail
|
||||
# puis pour chaque maillage,
|
||||
# >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit
|
||||
# correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions
|
||||
# des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment
|
||||
# puis pour chaque maillage
|
||||
# >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie ,
|
||||
# puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees
|
||||
# et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere)
|
||||
# precedes du numero de colonne correspondant entre crochet
|
||||
# puis pour chaque maillage
|
||||
# >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie ,
|
||||
# les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant
|
||||
# a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques,
|
||||
# on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl
|
||||
# chacun precedes du numero de colonne entre crochet
|
||||
# puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres),
|
||||
# suivi entre crochet, de la plage des numeros de colonnes, correspondant
|
||||
# chacun sur une ligne differentes
|
||||
# ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D:
|
||||
# tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2)
|
||||
# en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3)
|
||||
# non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3)
|
||||
# ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements,
|
||||
# ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique
|
||||
# ** ( cf: parametre_style_de_sortie = 0)
|
||||
|
||||
#====================================================================
|
||||
#|| recapitulatif des differentes grandeurs par colonne ||
|
||||
#====================================================================
|
||||
#---------------------------------- grandeur globales ------------------------------------
|
||||
#0 (nombre de grandeurs globales)
|
||||
#---------------------------------- maillage et dimension --------------------------------
|
||||
#1 3 (nombre de maillages et dimension)
|
||||
#---------------------------------- torseurs de reactions --------------------------------
|
||||
#0 0 (nombre de torseurs et nombre total de grandeurs associees)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#1 3 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
# reference N_haut_avant_gauche : noeud_8 [2]X [3]Y [4]Z [5]X1 [6]X2 [7]X3 ;
|
||||
#
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 3 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_1 pt_integ_1: [9]X [10]Y [11]Z [12] SIG11 [13] SIG22 [14] SIG33
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| fin du recapitulatif des differentes grandeurs ||
|
||||
#====================================================================
|
||||
|
||||
# ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant
|
||||
# a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun
|
||||
# correspondant a un maillage. On trouve pour chaque enregistrement successivement :
|
||||
# s'il y a des grandeurs globales: le temps puis les grandeurs globales,
|
||||
# puis s'il y a des torseurs de reaction :
|
||||
# de nouveau le temps, les composantes de la resultante puis les composantes du moments
|
||||
# donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels
|
||||
# puis s'il y a des grandeurs aux noeuds: de nouveau le temps
|
||||
# les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud
|
||||
# puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds
|
||||
# puis s'il y a des grandeur aux elements:
|
||||
# le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques)
|
||||
# suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration
|
||||
# correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes
|
||||
# puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs
|
||||
# etc. pour tous les points d'integration - element
|
||||
|
||||
1.000000000000e-01 1.002000000000e+00 9.994021528470e-01 9.994021528470e-01 1.002000000000e+00 9.994021528470e-01 9.994021528470e-01 1.000000000000e-01 2.117475151360e-01 7.882036274110e-01 7.882036274110e-01 4.186999968413e+02 -7.225796468299e-02 -7.225796473688e-02
|
||||
2.000000000000e-01 1.004000000000e+00 9.988092918389e-01 9.988092918389e-01 1.004000000000e+00 9.988092918389e-01 9.988092918389e-01 2.000000000000e-01 2.121701648668e-01 7.877360526756e-01 7.877360526756e-01 8.349867432157e+02 -5.044148934985e-06 -5.044148849720e-06
|
||||
3.000000000000e-01 1.006000000000e+00 9.982208308086e-01 9.982208308086e-01 1.006000000000e+00 9.982208308086e-01 9.982208308086e-01 3.000000000000e-01 2.125928145976e-01 7.872719480933e-01 7.872719480934e-01 1.248750042103e+03 -3.937182839309e-06 -3.937164706258e-06
|
||||
4.000000000000e-01 1.008000000000e+00 9.976368986747e-01 9.976368986747e-01 1.008000000000e+00 9.976368986747e-01 9.976368986747e-01 4.000000000000e-01 2.130154643284e-01 7.868114153390e-01 7.868114153390e-01 1.660052907544e+03 -4.181044289453e-06 -4.181062195130e-06
|
||||
5.000000000000e-01 1.010000000000e+00 9.970529665408e-01 9.970529665408e-01 1.010000000000e+00 9.970529665408e-01 9.970529665408e-01 5.000000000000e-01 2.134381140592e-01 7.863508825847e-01 7.863508825847e-01 2.067819639942e+03 -1.825287218956e+00 -1.825287218992e+00
|
||||
6.000000000000e-01 1.012000000000e+00 9.964824258645e-01 9.964824258645e-01 1.012000000000e+00 9.964824258645e-01 9.964824258645e-01 6.000000000000e-01 2.138607637900e-01 7.859009113400e-01 7.859009113400e-01 2.475355006163e+03 -3.672612228911e-05 -3.672615832784e-05
|
||||
7.000000000000e-01 1.014000000000e+00 9.959118851882e-01 9.959118851882e-01 1.014000000000e+00 9.959118851882e-01 9.959118851882e-01 7.000000000000e-01 2.142834135209e-01 7.854509400953e-01 7.854509400953e-01 2.879416119725e+03 3.912457109516e-02 3.912457105901e-02
|
||||
8.000000000000e-01 1.016000000000e+00 9.953413445119e-01 9.953413445119e-01 1.016000000000e+00 9.953413445119e-01 9.953413445119e-01 8.000000000000e-01 2.147060632517e-01 7.850009688507e-01 7.850009688507e-01 3.280027871816e+03 -1.697817128008e+00 -1.697817128045e+00
|
||||
9.000000000000e-01 1.018000000000e+00 9.947834816473e-01 9.947834816473e-01 1.018000000000e+00 9.947834816473e-01 9.947834816473e-01 9.000000000000e-01 2.151287129825e-01 7.845609962809e-01 7.845609962809e-01 3.680335473648e+03 -3.314144657907e-05 -3.314146454159e-05
|
||||
1.000000000000e+00 1.020000000000e+00 9.942256187827e-01 9.942256187827e-01 1.020000000000e+00 9.942256187827e-01 9.942256187827e-01 1.000000000000e+00 2.155513627133e-01 7.841210237110e-01 7.841210237110e-01 4.077253376881e+03 -4.098855042264e-02 -4.098855042230e-02
|
||||
1.100000000000e+00 1.020000000000e+00 9.962256187827e-01 9.933721658728e-01 1.020000000000e+00 9.962256187827e-01 9.933721658728e-01 1.100000000000e+00 2.155513627133e-01 7.856983739802e-01 7.834479266224e-01 4.217718765379e+03 4.681740028296e+02 -3.801827027837e-02
|
||||
1.200000000000e+00 1.020000000000e+00 9.982256187827e-01 9.925187129629e-01 1.020000000000e+00 9.982256187827e-01 9.925187129629e-01 1.200000000000e+00 2.155513627133e-01 7.872757242494e-01 7.827748295338e-01 4.356436264079e+03 9.326730854311e+02 -2.145438539390e+00
|
||||
1.300000000000e+00 1.020000000000e+00 1.000225618783e+00 9.916871782055e-01 1.020000000000e+00 1.000225618783e+00 9.916871782055e-01 1.300000000000e+00 2.155513627133e-01 7.888530745186e-01 7.821190187472e-01 4.496139628797e+03 1.396205739923e+03 -2.106395879764e-04
|
||||
1.400000000000e+00 1.020000000000e+00 1.002225618783e+00 9.908556434481e-01 1.020000000000e+00 1.002225618783e+00 9.908556434481e-01 1.400000000000e+00 2.155513627133e-01 7.904304247878e-01 7.814632079605e-01 4.634130600036e+03 1.856089274920e+03 8.615849624653e-02
|
||||
1.500000000000e+00 1.020000000000e+00 1.004225618783e+00 9.900241086907e-01 1.020000000000e+00 1.004225618783e+00 9.900241086907e-01 1.500000000000e+00 2.155513627133e-01 7.920077750570e-01 7.808073971738e-01 4.770419864359e+03 2.312349790221e+03 -1.876809064328e+00
|
||||
1.600000000000e+00 1.020000000000e+00 1.006225618783e+00 9.892127058348e-01 1.020000000000e+00 1.006225618783e+00 9.892127058348e-01 1.600000000000e+00 2.155513627133e-01 7.935851253262e-01 7.801674639172e-01 4.907537790639e+03 2.767532914915e+03 -1.794898118987e-04
|
||||
1.700000000000e+00 1.020000000000e+00 1.008225618783e+00 9.884013029789e-01 1.020000000000e+00 1.008225618783e+00 9.884013029789e-01 1.700000000000e+00 2.155513627133e-01 7.951624755954e-01 7.795275306606e-01 5.042987605215e+03 3.219156984343e+03 -1.248669994093e-01
|
||||
1.800000000000e+00 1.020000000000e+00 1.010225618783e+00 9.875899001230e-01 1.020000000000e+00 1.010225618783e+00 9.875899001230e-01 1.800000000000e+00 2.155513627133e-01 7.967398258646e-01 7.788875974040e-01 5.176779706159e+03 3.667247356173e+03 -2.241569013894e+00
|
||||
1.900000000000e+00 1.020000000000e+00 1.012225618783e+00 9.868000503675e-01 1.020000000000e+00 1.012225618783e+00 9.868000503675e-01 1.900000000000e+00 2.155513627133e-01 7.983171761337e-01 7.782646625417e-01 5.311641903360e+03 4.114546652242e+03 -2.077463709611e-04
|
||||
2.000000000000e+00 1.020000000000e+00 1.014225618783e+00 9.860102006119e-01 1.020000000000e+00 1.014225618783e+00 9.860102006119e-01 2.000000000000e+00 2.155513627133e-01 7.998945264029e-01 7.776417276794e-01 5.444880092637e+03 4.558375434422e+03 2.983076643923e-01
|
||||
Binary file not shown.
Loading…
Reference in a new issue