|
|
|
#######################################################################
|
|
# #
|
|
# | | ==== === ==== ==== | | | | #
|
|
# | | | | | | / | | | | #
|
|
# |====| |=== === |=== / |====| ------- ------- #
|
|
# | | | | \ | / | | | | #
|
|
# | | ==== | \ ==== ==== | | | | #
|
|
# #
|
|
#######################################################################
|
|
# Copyright (c) 1997-2014 Gerard Rio, gerard.rio@univ-ubs.fr #
|
|
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
|
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
|
# #
|
|
#######################################################################
|
|
# (version avec le plus de verifications pendant le calcul et les I/O )
|
|
version 6.674
|
|
-- initialisation de l'entree des donnees racine du nom de fichier principal = Hourglass
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques complet
|
|
debut de la lecture des maillages
|
|
lecture du premier maillage
|
|
debut de lecture de maillage
|
|
lecture du nom de maillage
|
|
lecture des noeuds
|
|
lecture du nombre de noeud
|
|
33
|
|
lecture des references de noeuds
|
|
lecture des elements du maillage
|
|
lecture du nombre d'element
|
|
20
|
|
lecture des references apres les elements
|
|
fin lecture de maillage
|
|
lecture affinage maillage
|
|
fin lecture affinage maillage
|
|
fin de la lecture des maillages
|
|
debut de la lecture des lois de comportement
|
|
lecture des references de materiaux
|
|
lecture des materiaux
|
|
fin de la lecture des lois de comportement
|
|
debut de la lecture des divers stockages 1
|
|
lecture epaisseurs lecture masse_volumique lecture hourglass_gestion_ fin de la lecture des divers stockages 1
|
|
debut de la lecture du chargement
|
|
lecture PONCTUELLE fin de la lecture du chargement
|
|
debut de la lecture des conditions limites ddl bloque
|
|
lecture blocages fin de la lecture des conditions limites ddl bloque
|
|
debut de la lecture du type d'application du chargement
|
|
|
|
warning : pas de mot cle typecharge, utilisation des valeurs par defaut
|
|
******** Type de Chargement ********
|
|
type de chargement : TYPE1, parametre : 1
|
|
fin de la lecture du du type d'application du chargement
|
|
debut de la lecture des parametres
|
|
|
|
fin de la lecture des parametres
|
|
debut de la lecture de la sortie automatique de resultats
|
|
fin de la lecture de la sortie automatique de resultats
|
|
|
|
======== fin de la lecture, Projet::Lecture( ...===========
|
|
|
|
|
|
raideur: largeur de bande initiale pour les elements (en ddl) 19
|
|
augmentation de la largeur de bande pour les conditions lineaires (en ddl) 0
|
|
largeur de bande finale (en ddl) 19
|
|
stockage symetrique d'ou utilisation de la demi largeur= 10
|
|
stockage raideurBANDE_SYMETRIQUE taille (entiers) :660
|
|
|
|
*** Attention, le nombre total de grandeurs generatrices aux point d'integrationx dans le probleme: 60 est inferieur a 1.2 * le nombre total de ddl du probleme !! c-a-d: 79
|
|
en consequence il est possible que la matrice de raideur soit naturellement singuliere
|
|
|
|
|
|
void Algori::VerifSingulariteRaideurMeca(...
|
|
======================================================================
|
|
INCREMENT DE CHARGE : 1 intensite 1 t= 1 dt= 1
|
|
======================================================================
|
|
MatBand resolution PMIN= 1 PMAX= 7.05155e+06
|
|
MatBand resolution PMIN= 2.11486e-07 PMAX= 1
|
|
|
|
ITERATION : 1
|
|
(assemblage 1) reaction maxi : 0.834725, du ddl X1, du noeud/maillage 31/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 1.00118
|
|
max des reactions = 0.834725
|
|
max du residu total = 2.00118
|
|
**** Norme Residu/Reaction --> 2.39741
|
|
|
|
MatBand resolution PMIN= 1 PMAX= 7.06036e+06
|
|
MatBand resolution PMIN= 2.1124e-07 PMAX= 1
|
|
|
|
(assemblage 1) delta maxi: -1.05217e-05, ddl X2, noeud/mail 18/1 , var total:5.26364e-06, sur:20
|
|
|
|
ITERATION : 2
|
|
(assemblage 1) reaction maxi : 0.833637, du ddl X1, du noeud/maillage 1/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 1
|
|
max des reactions = 0.833637
|
|
max du residu total = 1.60446e-06
|
|
**** Norme Residu/Reaction --> 1.92465e-06
|
|
|
|
erreur en ouverture pour la lecture du fichier Hourglass.CVisu
|
|
UtilLecture::Ouverture_CommandeVisu()
|
|
|
|
**** warning: erreur en ecriture pour la visualisation au fil du calcul
|
|
Algori::VisuAuFilDuCalcul(..
|
|
|
|
>>>> temps fin atteind <<<<<
|
|
=========== choix du module de visualisation interactive ============
|
|
sauvegarde des commandes de visualisation ? (rep 1)
|
|
visualisation automatique ? (rep 2)
|
|
visualisation au format vrml ? (rep 3)
|
|
visualisation par fichier de points, format maple ? (rep 4)
|
|
visualisation au format geomview ? (rep 5)
|
|
visualisation au format Gid ? (rep 6)
|
|
changement de fichier de commande .CVisu ? (rep 7)
|
|
visualisation au format Gmsh ? (rep 8)
|
|
fin (rep 0 ou f)
|
|
reponse ?
|
|
======== module de visualisation par fichiers de points au format maple ========
|
|
|
|
=== choix des increments utilises pour l'initialisation de la visualisation ====
|
|
option par defaut : tous les increments (rep 1)
|
|
choix d'un nombre plus petit d'increment (rep 2)
|
|
reponse ? lecture des elements de contact : 0 elements lus
|
|
----------------------------------------------------------
|
|
|
|
.........................choix numeros d'increment: cni
|
|
............choix du ou des maillages a visualiser: cmv
|
|
...................................choix grandeurs: cg
|
|
...................................animation_maple: ani
|
|
.....................................visualisation: visu
|
|
.arret visualisation interactive pour format maple: f
|
|
reponse ?
|
|
grandeurs globales -> rep : glo
|
|
torseurs de reactions -> rep : tre
|
|
ddl aux noeuds -> rep : noe ?
|
|
ddl etendu aux noeuds -> rep : net ?
|
|
grandeur particuliere aux noeuds -> rep : nop ?
|
|
grandeurs generique aux elements -> rep : ele
|
|
grandeurs particulieres aux elements -> rep : elp
|
|
grandeurs tensorielles aux elements -> rep : elt
|
|
style de sortie -> rep : sty
|
|
pour accepter la valeur par defaut -> rep : o
|
|
pour arreter les questions -> rep : fin (ou f)
|
|
reponse ?
|
|
liste des grandeurs globales disponibles en sortie:
|
|
>>> energie_bilan energie_bulk_viscosity energie_cinetique energie_elastique energie_externe energie_frot_elast energie_frot_plast energie_frot_visq energie_hourglass energie_interne energie_penalisation energie_plastique energie_visco_numerique energie_visqueuse puissance_acceleration puissance_bilan puissance_bulk_viscosity puissance_externe puissance_interne quantite_mouvement volume_total_matiere
|
|
|
|
grandeurs que vous voulez visualiser (rep grandeurs?)
|
|
toutes les grandeurs (rep par defaut) (rep : to)
|
|
effacer la liste actuelle (rep : ef)
|
|
(pour terminer tapez : fin (ou f) )
|
|
grandeur ?
|
|
liste des grandeurs globales disponibles en sortie:
|
|
>>> energie_bilan energie_bulk_viscosity energie_cinetique energie_elastique energie_externe energie_frot_elast energie_frot_plast energie_frot_visq energie_hourglass energie_interne energie_penalisation energie_plastique energie_visco_numerique energie_visqueuse puissance_acceleration puissance_bilan puissance_bulk_viscosity puissance_externe puissance_interne quantite_mouvement volume_total_matiere
|
|
liste des grandeurs actuellement a sortir:
|
|
<<<: energie_hourglass
|
|
|
|
grandeurs que vous voulez visualiser (rep grandeurs?)
|
|
toutes les grandeurs (rep par defaut) (rep : to)
|
|
effacer la liste actuelle (rep : ef)
|
|
(pour terminer tapez : fin (ou f) )
|
|
grandeur ?
|
|
liste des grandeurs globales disponibles en sortie:
|
|
>>> energie_bilan energie_bulk_viscosity energie_cinetique energie_elastique energie_externe energie_frot_elast energie_frot_plast energie_frot_visq energie_hourglass energie_interne energie_penalisation energie_plastique energie_visco_numerique energie_visqueuse puissance_acceleration puissance_bilan puissance_bulk_viscosity puissance_externe puissance_interne quantite_mouvement volume_total_matiere
|
|
liste des grandeurs actuellement a sortir:
|
|
<<<: energie_hourglass energie_elastique
|
|
|
|
grandeurs que vous voulez visualiser (rep grandeurs?)
|
|
toutes les grandeurs (rep par defaut) (rep : to)
|
|
effacer la liste actuelle (rep : ef)
|
|
(pour terminer tapez : fin (ou f) )
|
|
grandeur ?
|
|
liste des grandeurs enregistrees : energie_hourglass energie_elastique
|
|
|
|
grandeurs globales -> rep : glo
|
|
torseurs de reactions -> rep : tre
|
|
ddl aux noeuds -> rep : noe ?
|
|
ddl etendu aux noeuds -> rep : net ?
|
|
grandeur particuliere aux noeuds -> rep : nop ?
|
|
grandeurs generique aux elements -> rep : ele
|
|
grandeurs particulieres aux elements -> rep : elp
|
|
grandeurs tensorielles aux elements -> rep : elt
|
|
style de sortie -> rep : sty
|
|
pour accepter la valeur par defaut -> rep : o
|
|
pour arreter les questions -> rep : fin (ou f)
|
|
reponse ? ----------------------------------------------------------
|
|
|
|
.........................choix numeros d'increment: cni
|
|
............choix du ou des maillages a visualiser: cmv
|
|
...................................choix grandeurs: cg
|
|
...................................animation_maple: ani
|
|
.....................................visualisation: visu
|
|
.arret visualisation interactive pour format maple: f
|
|
reponse ?
|
|
======== fin module de visualisation par fichiers de points au format maple ========
|
|
|
|
=========== choix du module de visualisation interactive ============
|
|
sauvegarde des commandes de visualisation ? (rep 1)
|
|
visualisation automatique ? (rep 2)
|
|
visualisation au format vrml ? (rep 3)
|
|
visualisation par fichier de points, format maple ? (rep 4)
|
|
visualisation au format geomview ? (rep 5)
|
|
visualisation au format Gid ? (rep 6)
|
|
changement de fichier de commande .CVisu ? (rep 7)
|
|
visualisation au format Gmsh ? (rep 8)
|
|
fin (rep 0 ou f)
|
|
reponse ?
|
|
=========== choix du module de visualisation interactive ============
|
|
sauvegarde des commandes de visualisation ? (rep 1)
|
|
visualisation automatique ? (rep 2)
|
|
visualisation au format vrml ? (rep 3)
|
|
visualisation par fichier de points, format maple ? (rep 4)
|
|
visualisation au format geomview ? (rep 5)
|
|
visualisation au format Gid ? (rep 6)
|
|
changement de fichier de commande .CVisu ? (rep 7)
|
|
visualisation au format Gmsh ? (rep 8)
|
|
fin (rep 0 ou f)
|
|
reponse ?
|
|
=============================================================
|
|
| fin HEREZH++ |
|
|
=============================================================
|
|
|
|
#######################################################################
|
|
# #
|
|
# | | ==== === ==== ==== | | | | #
|
|
# | | | | | | / | | | | #
|
|
# |====| |=== === |=== / |====| ------- ------- #
|
|
# | | | | \ | / | | | | #
|
|
# | | ==== | \ ==== ==== | | | | #
|
|
# #
|
|
#######################################################################
|
|
# Copyright (c) 1997-2014 Gerard Rio, gerard.rio@univ-ubs.fr #
|
|
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
|
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
|
# #
|
|
#######################################################################
|
|
# (version avec le plus de verifications pendant le calcul et les I/O )
|
|
version 6.674
|
|
-- initialisation de l'entree des donnees racine du nom de fichier principal = Hourglass
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques complet
|
|
debut de la lecture des maillages
|
|
lecture du premier maillage
|
|
debut de lecture de maillage
|
|
lecture du nom de maillage
|
|
lecture des noeuds
|
|
lecture du nombre de noeud
|
|
33
|
|
lecture des references de noeuds
|
|
lecture des elements du maillage
|
|
lecture du nombre d'element
|
|
20
|
|
lecture des references apres les elements
|
|
fin lecture de maillage
|
|
lecture affinage maillage
|
|
fin lecture affinage maillage
|
|
fin de la lecture des maillages
|
|
debut de la lecture des lois de comportement
|
|
lecture des references de materiaux
|
|
lecture des materiaux
|
|
fin de la lecture des lois de comportement
|
|
debut de la lecture des divers stockages 1
|
|
lecture epaisseurs lecture masse_volumique lecture hourglass_gestion_ fin de la lecture des divers stockages 1
|
|
debut de la lecture du chargement
|
|
lecture PONCTUELLE fin de la lecture du chargement
|
|
debut de la lecture des conditions limites ddl bloque
|
|
lecture blocages fin de la lecture des conditions limites ddl bloque
|
|
debut de la lecture du type d'application du chargement
|
|
|
|
warning : pas de mot cle typecharge, utilisation des valeurs par defaut
|
|
******** Type de Chargement ********
|
|
type de chargement : TYPE1, parametre : 1
|
|
fin de la lecture du du type d'application du chargement
|
|
debut de la lecture des parametres
|
|
|
|
fin de la lecture des parametres
|
|
debut de la lecture de la sortie automatique de resultats
|
|
fin de la lecture de la sortie automatique de resultats
|
|
|
|
======== fin de la lecture, Projet::Lecture( ...===========
|
|
|
|
|
|
raideur: largeur de bande initiale pour les elements (en ddl) 19
|
|
augmentation de la largeur de bande pour les conditions lineaires (en ddl) 0
|
|
largeur de bande finale (en ddl) 19
|
|
stockage symetrique d'ou utilisation de la demi largeur= 10
|
|
stockage raideurBANDE_SYMETRIQUE taille (entiers) :660
|
|
|
|
*** Attention, le nombre total de grandeurs generatrices aux point d'integrationx dans le probleme: 60 est inferieur a 1.2 * le nombre total de ddl du probleme !! c-a-d: 79
|
|
en consequence il est possible que la matrice de raideur soit naturellement singuliere
|
|
|
|
|
|
void Algori::VerifSingulariteRaideurMeca(...
|
|
======================================================================
|
|
INCREMENT DE CHARGE : 1 intensite 1 t= 1 dt= 1
|
|
======================================================================
|
|
MatBand resolution PMIN= 1 PMAX= 7.05155e+06
|
|
MatBand resolution PMIN= 2.11486e-07 PMAX= 1
|
|
|
|
ITERATION : 1
|
|
(assemblage 1) reaction maxi : 0.834725, du ddl X1, du noeud/maillage 31/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 1.00118
|
|
max des reactions = 0.834725
|
|
max du residu total = 2.00118
|
|
**** Norme Residu/Reaction --> 2.39741
|
|
|
|
MatBand resolution PMIN= 1 PMAX= 7.06036e+06
|
|
MatBand resolution PMIN= 2.1124e-07 PMAX= 1
|
|
|
|
(assemblage 1) delta maxi: -1.05217e-05, ddl X2, noeud/mail 18/1 , var total:5.26364e-06, sur:20
|
|
|
|
ITERATION : 2
|
|
(assemblage 1) reaction maxi : 0.833637, du ddl X1, du noeud/maillage 1/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 1
|
|
max des reactions = 0.833637
|
|
max du residu total = 1.60446e-06
|
|
**** Norme Residu/Reaction --> 1.92465e-06
|
|
|
|
======== lecture du fichier de commande de visualisation ========
|
|
|
|
======== fin de la lecture du fichier de commande de visualisation ========
|
|
|
|
>>>> temps fin atteind <<<<<
|
|
=============================================================
|
|
| fin HEREZH++ |
|
|
=============================================================
|
|
|
|
#######################################################################
|
|
# #
|
|
# | | ==== === ==== ==== | | | | #
|
|
# | | | | | | / | | | | #
|
|
# |====| |=== === |=== / |====| ------- ------- #
|
|
# | | | | \ | / | | | | #
|
|
# | | ==== | \ ==== ==== | | | | #
|
|
# #
|
|
#######################################################################
|
|
# Copyright (c) 1997-2014 Gerard Rio, gerard.rio@univ-ubs.fr #
|
|
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
|
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
|
# #
|
|
#######################################################################
|
|
# (version avec le plus de verifications pendant le calcul et les I/O )
|
|
version 6.674
|
|
-- initialisation de l'entree des donnees racine du nom de fichier principal = Hourglass
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques complet
|
|
debut de la lecture des maillages
|
|
lecture du premier maillage
|
|
debut de lecture de maillage
|
|
lecture du nom de maillage
|
|
lecture des noeuds
|
|
lecture du nombre de noeud
|
|
33
|
|
lecture des references de noeuds
|
|
lecture des elements du maillage
|
|
lecture du nombre d'element
|
|
20
|
|
lecture des references apres les elements
|
|
fin lecture de maillage
|
|
lecture affinage maillage
|
|
fin lecture affinage maillage
|
|
fin de la lecture des maillages
|
|
debut de la lecture des lois de comportement
|
|
lecture des references de materiaux
|
|
lecture des materiaux
|
|
fin de la lecture des lois de comportement
|
|
debut de la lecture des divers stockages 1
|
|
lecture epaisseurs lecture masse_volumique lecture hourglass_gestion_ fin de la lecture des divers stockages 1
|
|
debut de la lecture du chargement
|
|
lecture PONCTUELLE fin de la lecture du chargement
|
|
debut de la lecture des conditions limites ddl bloque
|
|
lecture blocages fin de la lecture des conditions limites ddl bloque
|
|
debut de la lecture du type d'application du chargement
|
|
|
|
warning : pas de mot cle typecharge, utilisation des valeurs par defaut
|
|
******** Type de Chargement ********
|
|
type de chargement : TYPE1, parametre : 1
|
|
fin de la lecture du du type d'application du chargement
|
|
debut de la lecture des parametres
|
|
|
|
fin de la lecture des parametres
|
|
debut de la lecture de la sortie automatique de resultats
|
|
fin de la lecture de la sortie automatique de resultats
|
|
|
|
======== fin de la lecture, Projet::Lecture( ...===========
|
|
|
|
|
|
raideur: largeur de bande initiale pour les elements (en ddl) 19
|
|
augmentation de la largeur de bande pour les conditions lineaires (en ddl) 0
|
|
largeur de bande finale (en ddl) 19
|
|
stockage symetrique d'ou utilisation de la demi largeur= 10
|
|
stockage raideurBANDE_SYMETRIQUE taille (entiers) :660
|
|
|
|
*** Attention, le nombre total de grandeurs generatrices aux point d'integrationx dans le probleme: 60 est inferieur a 1.2 * le nombre total de ddl du probleme !! c-a-d: 79
|
|
en consequence il est possible que la matrice de raideur soit naturellement singuliere
|
|
|
|
======================================================================
|
|
INCREMENT DE CHARGE : 1 intensite 1 t= 1 dt= 1
|
|
======================================================================
|
|
|
|
MatBand resolution PMIN= 1 PMAX= 6.30492e+06
|
|
MatBand resolution PMIN= 1.91358e-07 PMAX= 1
|
|
|
|
ITERATION : 1
|
|
(assemblage 1) reaction maxi : 0.671454, du ddl X1, du noeud/maillage 1/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 1.00132
|
|
max des reactions = 0.671454
|
|
max du residu total = 2.00132
|
|
**** Norme Residu/Reaction --> 2.98058
|
|
|
|
MatBand resolution PMIN= 1 PMAX= 6.31373e+06
|
|
MatBand resolution PMIN= 1.91092e-07 PMAX= 1
|
|
|
|
(assemblage 1) delta maxi: -1.26471e-05, ddl X2, noeud/mail 18/1 , var total:6.32617e-06, sur:20
|
|
autres ddl R_X2 0 var 0;
|
|
|
|
ITERATION : 2
|
|
(assemblage 1) reaction maxi : 0.670495, du ddl X1, du noeud/maillage 1/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 0.999999
|
|
max des reactions = 0.670495
|
|
max du residu total = 6.49379e-07
|
|
**** Norme Residu/Reaction --> 9.68507e-07
|
|
|
|
======== lecture du fichier de commande de visualisation ========
|
|
|
|
======== fin de la lecture du fichier de commande de visualisation ========
|
|
|
|
>>>> temps fin atteind <<<<<
|
|
=========== choix du module de visualisation interactive ============
|
|
sauvegarde des commandes de visualisation ? (rep 1)
|
|
visualisation automatique ? (rep 2)
|
|
visualisation au format vrml ? (rep 3)
|
|
visualisation par fichier de points, format maple ? (rep 4)
|
|
visualisation au format geomview ? (rep 5)
|
|
visualisation au format Gid ? (rep 6)
|
|
changement de fichier de commande .CVisu ? (rep 7)
|
|
visualisation au format Gmsh ? (rep 8)
|
|
fin (rep 0 ou f)
|
|
reponse ?
|
|
======== module de visualisation format Gmsh ========
|
|
|
|
=== choix des increments utilises pour l'initialisation de la visualisation ====
|
|
option par defaut : tous les increments (rep 1)
|
|
choix d'un nombre plus petit d'increment (rep 2)
|
|
reponse ? lecture des elements de contact : 0 elements lus
|
|
----------------------------------------------------------
|
|
|
|
.................................maillage initiale: mi
|
|
........................................isovaleurs: iso
|
|
..........................................deformee: de
|
|
.........................choix numeros d'increment: cni
|
|
............choix du ou des maillages a visualiser: cmv
|
|
....................................visualisation : visu
|
|
.............arret de la visualisation interactive: f
|
|
reponse ?
|
|
----> preparation de la visualisation des coordonnees initiales
|
|
parametre par defaut ? : resultat brut du calcul elements finis ,
|
|
pas d'homotheties sur les coordonnees initiales ,
|
|
sortie des references dans un seul fichier,
|
|
---> (rep 'o') pour accepter ces parametres sinon autre
|
|
reponse ? ----------------------------------------------------------
|
|
|
|
.................................maillage initiale: mi
|
|
........................................isovaleurs: iso
|
|
..........................................deformee: de
|
|
.........................choix numeros d'increment: cni
|
|
............choix du ou des maillages a visualiser: cmv
|
|
....................................visualisation : visu
|
|
.............arret de la visualisation interactive: f
|
|
reponse ?
|
|
----> preparation de la visualisation des deformees
|
|
parametre par defaut ? : pas de limites d'alerte sur les deplacements , (rep 'o') pour accepter ces parametres sinon autre
|
|
reponse ? ----------------------------------------------------------
|
|
|
|
.................................maillage initiale: mi
|
|
........................................isovaleurs: iso
|
|
..........................................deformee: de
|
|
.........................choix numeros d'increment: cni
|
|
............choix du ou des maillages a visualiser: cmv
|
|
....................................visualisation : visu
|
|
.............arret de la visualisation interactive: f
|
|
reponse ?
|
|
---- isovaleurs ----
|
|
(0 ou f ou fin) fin modif
|
|
(1) ou (de) parametres par defaut:
|
|
(2) ou (to) toutes les isovaleurs
|
|
(3) parametres generaux pour la sortie
|
|
(4) ddl naturellement defini aux noeuds
|
|
(5) grandeurs scalaires venant des pts integ
|
|
(6) grandeurs tensoriel venant des pts integ
|
|
(7) grandeurs particulieres venant des pts integ
|
|
(8) choix de l'ancien format gmsh (faux par defaut)
|
|
(9) ddl etendu aux noeuds
|
|
(10) grandeurs evoluees aux noeuds
|
|
|
|
reponse ?
|
|
choix d'isovaleur venant des pts d'integ a visualiser (de type grandeurs particulieres)
|
|
(0 ou f ou fin) fin choix grandeurs particulieres
|
|
listes de type de grandeur particuliere a visualiser
|
|
Maillage nb: 1 liste des types de grandeur particuliere disponibles
|
|
MODULE_COMPRESSIBILITE_TOTAL
|
|
MODULE_CISAILLEMENT_TOTAL
|
|
VOLUME_ELEMENT
|
|
ENERGIE_HOURGLASS
|
|
PUISSANCE_BULK
|
|
ENERGIE_BULK
|
|
EPAISSEUR_MOY_INITIALE
|
|
EPAISSEUR_MOY_FINALE
|
|
donner la grandeur particuliere a visualiser
|
|
(to) tous les grandeurs particulieres
|
|
ou une liste de grandeurs particulieres
|
|
(0 ou f ou fin) fin choix de grandeurs particulieres
|
|
reponse ?
|
|
(0 ou f ou fin) fin choix grandeurs particulieres
|
|
listes de type de grandeur particuliere a visualiser
|
|
Maillage nb: 1
|
|
liste des types de grandeurs enregistres: ENERGIE_HOURGLASS
|
|
|
|
Maillage nb: 1 liste des types de grandeur particuliere disponibles
|
|
MODULE_COMPRESSIBILITE_TOTAL
|
|
MODULE_CISAILLEMENT_TOTAL
|
|
VOLUME_ELEMENT
|
|
ENERGIE_HOURGLASS
|
|
PUISSANCE_BULK
|
|
ENERGIE_BULK
|
|
EPAISSEUR_MOY_INITIALE
|
|
EPAISSEUR_MOY_FINALE
|
|
donner la grandeur particuliere a visualiser
|
|
(to) tous les grandeurs particulieres
|
|
ou une liste de grandeurs particulieres
|
|
(0 ou f ou fin) fin choix de grandeurs particulieres
|
|
reponse ?
|
|
(0 ou f ou fin) fin modif
|
|
(1) ou (de) parametres par defaut:
|
|
(2) ou (to) toutes les isovaleurs
|
|
(3) parametres generaux pour la sortie
|
|
(4) ddl naturellement defini aux noeuds
|
|
(5) grandeurs scalaires venant des pts integ
|
|
(6) grandeurs tensoriel venant des pts integ
|
|
(7) grandeurs particulieres venant des pts integ
|
|
(8) choix de l'ancien format gmsh (faux par defaut)
|
|
(9) ddl etendu aux noeuds
|
|
(10) grandeurs evoluees aux noeuds
|
|
|
|
reponse ? ----------------------------------------------------------
|
|
|
|
.................................maillage initiale: mi
|
|
........................................isovaleurs: iso
|
|
..........................................deformee: de
|
|
.........................choix numeros d'increment: cni
|
|
............choix du ou des maillages a visualiser: cmv
|
|
....................................visualisation : visu
|
|
.............arret de la visualisation interactive: f
|
|
reponse ?
|
|
======== fin du module de visualisation format Gmsh ========
|
|
|
|
=========== choix du module de visualisation interactive ============
|
|
sauvegarde des commandes de visualisation ? (rep 1)
|
|
visualisation automatique ? (rep 2)
|
|
visualisation au format vrml ? (rep 3)
|
|
visualisation par fichier de points, format maple ? (rep 4)
|
|
visualisation au format geomview ? (rep 5)
|
|
visualisation au format Gid ? (rep 6)
|
|
changement de fichier de commande .CVisu ? (rep 7)
|
|
visualisation au format Gmsh ? (rep 8)
|
|
fin (rep 0 ou f)
|
|
reponse ?
|
|
=========== choix du module de visualisation interactive ============
|
|
sauvegarde des commandes de visualisation ? (rep 1)
|
|
visualisation automatique ? (rep 2)
|
|
visualisation au format vrml ? (rep 3)
|
|
visualisation par fichier de points, format maple ? (rep 4)
|
|
visualisation au format geomview ? (rep 5)
|
|
visualisation au format Gid ? (rep 6)
|
|
changement de fichier de commande .CVisu ? (rep 7)
|
|
visualisation au format Gmsh ? (rep 8)
|
|
fin (rep 0 ou f)
|
|
reponse ?
|
|
=============================================================
|
|
| fin HEREZH++ |
|
|
=============================================================
|
|
|
|
#######################################################################
|
|
# #
|
|
# | | ==== === ==== ==== | | | | #
|
|
# | | | | | | / | | | | #
|
|
# |====| |=== === |=== / |====| ------- ------- #
|
|
# | | | | \ | / | | | | #
|
|
# | | ==== | \ ==== ==== | | | | #
|
|
# #
|
|
#######################################################################
|
|
# Copyright (c) 1997-2014 Gerard Rio, gerard.rio@univ-ubs.fr #
|
|
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
|
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
|
# #
|
|
#######################################################################
|
|
# (version avec le plus de verifications pendant le calcul et les I/O )
|
|
version 6.674
|
|
-- initialisation de l'entree des donnees racine du nom de fichier principal = Hourglass
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques
|
|
|
|
ATTENTION -> dimension 2, pas de definition d'elements hexaedriques quadratiques complet
|
|
debut de la lecture des maillages
|
|
lecture du premier maillage
|
|
debut de lecture de maillage
|
|
lecture du nom de maillage
|
|
lecture des noeuds
|
|
lecture du nombre de noeud
|
|
33
|
|
lecture des references de noeuds
|
|
lecture des elements du maillage
|
|
lecture du nombre d'element
|
|
20
|
|
lecture des references apres les elements
|
|
fin lecture de maillage
|
|
lecture affinage maillage
|
|
fin lecture affinage maillage
|
|
fin de la lecture des maillages
|
|
debut de la lecture des lois de comportement
|
|
lecture des references de materiaux
|
|
lecture des materiaux
|
|
fin de la lecture des lois de comportement
|
|
debut de la lecture des divers stockages 1
|
|
lecture epaisseurs lecture masse_volumique lecture hourglass_gestion_ fin de la lecture des divers stockages 1
|
|
debut de la lecture du chargement
|
|
lecture PONCTUELLE fin de la lecture du chargement
|
|
debut de la lecture des conditions limites ddl bloque
|
|
lecture blocages fin de la lecture des conditions limites ddl bloque
|
|
debut de la lecture du type d'application du chargement
|
|
|
|
warning : pas de mot cle typecharge, utilisation des valeurs par defaut
|
|
******** Type de Chargement ********
|
|
type de chargement : TYPE1, parametre : 1
|
|
fin de la lecture du du type d'application du chargement
|
|
debut de la lecture des parametres
|
|
|
|
fin de la lecture des parametres
|
|
debut de la lecture de la sortie automatique de resultats
|
|
fin de la lecture de la sortie automatique de resultats
|
|
|
|
======== fin de la lecture, Projet::Lecture( ...===========
|
|
|
|
|
|
raideur: largeur de bande initiale pour les elements (en ddl) 19
|
|
augmentation de la largeur de bande pour les conditions lineaires (en ddl) 0
|
|
largeur de bande finale (en ddl) 19
|
|
stockage symetrique d'ou utilisation de la demi largeur= 10
|
|
stockage raideurBANDE_SYMETRIQUE taille (entiers) :660
|
|
|
|
*** Attention, le nombre total de grandeurs generatrices aux point d'integrationx dans le probleme: 60 est inferieur a 1.2 * le nombre total de ddl du probleme !! c-a-d: 79
|
|
en consequence il est possible que la matrice de raideur soit naturellement singuliere
|
|
|
|
|
|
void Algori::VerifSingulariteRaideurMeca(...
|
|
======================================================================
|
|
INCREMENT DE CHARGE : 1 intensite 1 t= 1 dt= 1
|
|
======================================================================
|
|
|
|
MatBand resolution PMIN= 1 PMAX= 6.30492e+06
|
|
MatBand resolution PMIN= 1.91358e-07 PMAX= 1
|
|
|
|
ITERATION : 1
|
|
(assemblage 1) reaction maxi : 0.671454, du ddl X1, du noeud/maillage 1/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 1.00132
|
|
max des reactions = 0.671454
|
|
max du residu total = 2.00132
|
|
**** Norme Residu/Reaction --> 2.98058
|
|
|
|
MatBand resolution PMIN= 1 PMAX= 6.31373e+06
|
|
MatBand resolution PMIN= 1.91092e-07 PMAX= 1
|
|
|
|
(assemblage 1) delta maxi: -1.26471e-05, ddl X2, noeud/mail 18/1 , var total:6.32617e-06, sur:20
|
|
autres ddl R_X2 0 var 0;
|
|
|
|
ITERATION : 2
|
|
(assemblage 1) reaction maxi : 0.670495, du ddl X1, du noeud/maillage 1/1
|
|
|
|
max puissance exterieure = 1
|
|
max puissance interieurs = 0.999999
|
|
max des reactions = 0.670495
|
|
max du residu total = 6.49379e-07
|
|
**** Norme Residu/Reaction --> 9.68507e-07
|
|
|
|
======== lecture du fichier de commande de visualisation ========
|
|
|
|
======== fin de la lecture du fichier de commande de visualisation ========
|
|
|
|
ouverture du fichier Hourglass_Gmsh/Hourglass_ENERG_HOURG_Gmsh.pos
|
|
|
|
ouverture du fichier Hourglass_Gmsh/Hourglass_deplace_Gmsh.pos
|
|
|
|
>>>> temps fin atteind <<<<<
|
|
=============================================================
|
|
| fin HEREZH++ |
|
|
=============================================================
|
