2
0
Fork 0

modif du test pour prendre en contact la correction du bug sur le contact avec des elements axi

This commit is contained in:
rio 2016-02-25 17:02:28 +00:00
parent d615853c8d
commit fd6cbf097e
4 changed files with 107 additions and 4 deletions

View file

@ -0,0 +1,6 @@
Modif Gérard Rio
----------------
à partir de la version version 6.735, il y a un bug sur le contact qui a été corrigé
ce qui fait que la référence n'est pas bonne.
Donc je laisse l'ancienne pour mémoire mais je change la référence avec un calcul effectué
avec la version V6.739

View file

@ -85,14 +85,17 @@ para_contact ------------------------------------
# > 2 : contact avec une methode de penalisation
CONTACT_TYPE 2
PENALISATION_PENETRATION 0.1
PENALISATION_PENETRATION 0.1#0.1
#TYPE_PENALISATION_PENETRATION (indique la methode pour calculer le facteur beta)
# > 1 : dans ce cas beta=1 et alpha = facteur de penalisation
# > 2 : beta est calcule
# > 3 : meme fonctionnement que 2 mais avec un algorithme supplementaire
# > 5 : meme fonctionnement que 2 mais avec un algorithme supplementaire
TYPE_PENALISATION_PENETRATION 2
TYPE_PENALISATION_PENETRATION 4 #2
PENETRATION_CONTACT_MAXI 0.01
PENETRATION_BORNE_REGULARISATION 0.1
NB_DECOLLEMENT_MAXI 1
DISTANCE_MAXI_AU_PT_PROJETE 2.
TYPE_DE_DECOLLEMENT 1

View file

@ -1,6 +1,6 @@
#fichier au format maple6
###############################################################################################
# Visualisation elements finis : Herezh++ V6.687 #
# Visualisation elements finis : Herezh++ V6.739 #
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
# http://www-lg2m.univ-ubs.fr #
###############################################################################################
@ -91,4 +91,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.699124778880e+01 2.151000000000e+01 0.000000000000e+00 1.699124778880e+01 2.151000000000e+01 0.000000000000e+00 7.836307275110e+03 2.299171571673e+01 1.775534086544e+01 0.000000000000e+00 2.299171571673e+01 1.775534086544e+01 -3.827610525875e+03 -2.733190283370e-01 1.000000000000e+00 2.298875938452e+01 1.574822255247e+01 0.000000000000e+00 2.298875938452e+01 1.574822255247e+01 3.955599632662e+03 -2.857909684829e+00 1.000000000000e+00 2.276827274533e+01 1.681640459547e+01 0.000000000000e+00 6.982193210439e+01 2.881540055102e-02 1.000000000000e+00 2.320903886251e+01 1.687243479857e+01 0.000000000000e+00 1.712718075877e+02 7.068360313144e-04
1.000000000000e+00 1.699227339604e+01 2.151000000000e+01 0.000000000000e+00 1.699227339604e+01 2.151000000000e+01 0.000000000000e+00 7.810506010394e+03 2.299410392537e+01 1.775501775788e+01 0.000000000000e+00 2.299410392537e+01 1.775501775788e+01 -3.879341157508e+03 -2.847089384744e-01 1.000000000000e+00 2.298841524205e+01 1.574825045586e+01 0.000000000000e+00 2.298841524205e+01 1.574825045586e+01 3.905049255510e+03 -2.691105799499e+00 1.000000000000e+00 2.277112712990e+01 1.681615554587e+01 0.000000000000e+00 7.049517889726e+01 2.909324843379e-02 1.000000000000e+00 2.320866829579e+01 1.687248444799e+01 0.000000000000e+00 1.667788013196e+02 6.882934657632e-04

View file

@ -0,0 +1,94 @@
#fichier au format maple6
###############################################################################################
# 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 #
###############################################################################################
# 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 --------------------------------
#2 3 (nombre de maillages et dimension)
#---------------------------------- torseurs de reactions --------------------------------
#0 0 (nombre de torseurs et nombre total de grandeurs associees)
#
#0 0 (nombre de torseurs et nombre total de grandeurs associees)
#
#---------------------------------- grandeurs aux noeuds --------------------------------
#2 4 (nombre de noeuds, nombre total de grandeurs associees)
# maillage_1 :
# noeud_27 [2]X [3]Y [4]Z [5]X1 [6]X2 [7]R_X1 [8]R_X2 ;
# noeud_77 [9]X [10]Y [11]Z [12]X1 [13]X2 [14]R_X1 [15]R_X2 ;
#
#1 4 (nombre de noeuds, nombre total de grandeurs associees)
# maillage_2 :
# noeud_75 [17]X [18]Y [19]Z [20]X1 [21]X2 [22]R_X1 [23]R_X2 ;
#
#---------------------------------- grandeurs aux elements ------------------------------
#1 2 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
# maillage_1 :
# reference G_elt_14_pti_6 : element_14 pt_integ_6: [25]X [26]Y [27]Z [28] contrainte_mises [29] def_duale_mises
#
#
#1 2 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
# maillage_2 :
# reference G_elt_16_pti_8 : element_16 pt_integ_8: [31]X [32]Y [33]Z [34] contrainte_mises [35] def_duale_mises
#
#
#====================================================================
#|| 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.699124778880e+01 2.151000000000e+01 0.000000000000e+00 1.699124778880e+01 2.151000000000e+01 0.000000000000e+00 7.836307275110e+03 2.299171571673e+01 1.775534086544e+01 0.000000000000e+00 2.299171571673e+01 1.775534086544e+01 -3.827610525875e+03 -2.733190283370e-01 1.000000000000e+00 2.298875938452e+01 1.574822255247e+01 0.000000000000e+00 2.298875938452e+01 1.574822255247e+01 3.955599632662e+03 -2.857909684829e+00 1.000000000000e+00 2.276827274533e+01 1.681640459547e+01 0.000000000000e+00 6.982193210439e+01 2.881540055102e-02 1.000000000000e+00 2.320903886251e+01 1.687243479857e+01 0.000000000000e+00 1.712718075877e+02 7.068360313144e-04