2
0
Fork 0

ajout d un test de ref sur la simulation d un cisaillement de membrane avec un amortissement visqueux, une loi elastique avec plis et divers parametres de reglage

This commit is contained in:
rio 2018-06-01 08:38:04 +00:00
parent 31206ee42c
commit 32afb16f1c
6 changed files with 964 additions and 0 deletions

View file

@ -0,0 +1,34 @@
------------------------------------------------------
Auteur
------------------------------------------------------
Gerard Rio (gerard.rio@univ-ubs.fr)
------------------------------------------------------
Mots-cles
------------------------------------------------------
dynamique_relaxation_dynam
parametre_calcul_de_la_viscosite_
LOI_CRITERE
PLISSEMENT_MEMBRANE
permet_affichage_
------------------------------------------------------------
But du test
------------------------------------------------------------
Calcul du cisaillement d'une membrane, avec une loi critère qui prend en compte un comportement élastique linéaire avec plis éventuelles.
L'équilibre global est calculé en relaxation dynamique avec amortissement visqueux.
""
TYPE_DE_CALCUL
non_dynamique
""
------------------------------------------------------------
Description du calcul
------------------------------------------------------------
Le calcul s'effectue ici en 3D sur un maillage rectangulaire, avec des triangles linéaires à un point d'intégration.
La plaque étant maillée avec Stamm, les conditions limites sont: N de l'ouest bloqués, N de l'est bloqués en x et avec un déplacement en y imposé.
L'ensemble des noeuds est bloqué en z.
------------------------------------------------------------
Grandeurs de comparaison
------------------------------------------------------------
On regarde l'énergie de déformation globale, le torseur de réaction sur la ligne Est, et les contraintes et déformations sur : un point approximativement à trois extrémités et un point approximativement central.

View file

@ -0,0 +1,221 @@
###############################################################################################
# Fichier de commande pour la visualisation elements finis #
# Herezh++ V6.844 #
# Copyright (c) 1997-2018, 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)
tous_les_increments_moins_zero 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 energie_elastique 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 N_E fin_list_torseur_reaction
# les sommes, moyennes etc. sur ref de noeuds: un mot cle de debut: debut_list_SM_sur_refN
# une liste de nom de ref <chaine de caracteres >, puis <fin_list__SM_sur_refN>
debut_list_SM_sur_refN fin_list__SM_sur_refN
# debut de la liste des ddl a considerer <deb_list_ddl_SM_sur_refN>, (une liste de ddl), puis <fin_list_ddl_SM_sur_refN>
deb_list_ddl_SM_sur_refN fin_list_ddl_SM_sur_refN
# debut de la liste des ddl etendus a considerer, (une liste de ddl), puis <fin_list_ddl_etendu__SM_sur_refN>
deb_list_ddl_etendu__SM_sur_refN fin_list_ddl_etendu__SM_sur_refN # fin des ddl etendu pour _SM_sur_refNs
# debut de la liste des grandeurs particulieres a considerer, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refN>
deb_list_GrandParticuliere__SM_sur_refN fin_list_GrandParticuliere__SM_sur_refN # fin des grandeurs particulieres aux _SM_sur_refNs
# les SM sur ref d'elements: un mot cle de debut: debut_list_SM_sur_refE
# une liste de nom de ref d'elements <chaine de caracteres >, le pti associe
# puis <fin_list__SM_sur_refE>
debut_list_SM_sur_refE fin_list__SM_sur_refE
# les SMs sur ref de pti: un mot cle de debut: debut_list_SM_sur_refpti
# une liste de nom de ref <chaine de caracteres >, puis <fin_list__SM_sur_refpti>
debut_list_SM_sur_refpti fin_list__SM_sur_refpti
# debut de la liste des ddl a considerer , (une liste de ddl), puis <fin_list_ddl__SM_sur_refEpti>
deb_list_ddl__SM_sur_refEpti fin_list_ddl__SM_sur_refEpti # fin de la liste de ddl a considerer pour les elements
# debut de la liste des grandeurs particulieres a considerer pour les _SM_sur_refEptis, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refEpti>
deb_list_GrandParticuliere__SM_sur_refEpti fin_list_GrandParticuliere__SM_sur_refEpti # fin des grandeurs particulieres aux _SM_sur_refEptis
# tableau de grandeurs evoluees aux _SM_sur_refEptis a visualiser, un par maillage
deb_list_GrandEvoluee__SM_sur_refEpti fin_list_GrandEvoluee__SM_sur_refEpti
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 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 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 39 1 110 1 199 1 4 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 SIG12 EPS11 EPS22 EPS12 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 ||
# =================================================================================
# =================================================================================
# || ***** demande d'une visualisation Gmsh: ***** ||
# =================================================================================
# un mot cle de debut (debut_visualisation_Gmsh)
# un mot cle de fin ( fin_visualisation_Gmsh) apres tous les ordres particuliers
# la seule presence du premier mots cle suffit a activer la visualisation Gmsh
# la presence du second permet une meilleur lisibilite du fichier, mais n'est pas indispensable
debut_visualisation_Gmsh
# ----------------------------- definition des parametres du maillage initial: ----------------
debut_maillage_initial # un mot cle de debut de liste
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
pseudo-homothetie_sur_les_maillages_ 0 # 0 = aucune homothetie, 1 = il y en a
# --- def eventuelle de la pseudo-homothetie: une par maillage,
# pseudo-homothetie pour le maillage : 1
# maillage_ 1
# mot cle: maillage_ puis le numero du maillage,
# pseudo-homothetie_ 0 # 0 = non active, 1 = active
# ensuite si c'est active, on trouve :
# mot cle: centre_homothetie_ puis les coordonnees du centre d'homothetie
# puis mot cle: fact_mult_ puis les coordonnees donnant les coefs multiplicatifs selon les axes.
# centre_homothetie_
# fact_mult_
visualisation_references_sur_les_maillages_ 1 # 0 = pas de visualisation des reference, 1 = sortie des ref dans fichier unique 2= sortie des ref dans plusieurs fichiers
fin_maillage_initial # le mot cle de fin
# ----------------------------- definition des parametres pour les isovaleurs : ----------------
debut_isovaleur_Gmsh # mot cle de debut des parametres pour les isovaleurs
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
ancien_format_ 0 # 1 = ancien format, 0 = nouveau format
# 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_ 1
1 # le numero de maillage
# tableau des ddl aux noeuds a visualiser, un par maillage
debut_tableau_ddl_aux_noeuds fin_tableau_ddl_aux_noeuds
# tableau des choix_var aux noeuds a visualiser, un par maillage
# choix_var (=1 ou 0) indique si oui ou non il s'agit de la variation
debut_tableau_choix_var_ddl_aux_noeuds fin_tableau_choix_var_ddl_aux_noeuds
# tableau des ddl_etendu aux noeuds a visualiser, un par maillage
debut_tableau_ddl_etendu_aux_noeuds fin_tableau_ddl_etendu_aux_noeuds
# tableau de grandeurs evoluees aux noeuds a visualiser, un par maillage
deb_list_GrandEvoluee_noeud fin_list_GrandEvoluee_noeud
# tableau de ddl aux elements a visualiser, un par maillage
debut_tableau_ddl_aux_elements SIG11 SIG22 SIG12 EPS11 EPS22 EPS12 Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII fin_tableau_ddl_aux_elements
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
deb_list_GrandEvoluee_element CONTRAINTE_COURANTE DEFORMATION_COURANTE fin_list_GrandEvoluee_element
# tableau de grandeurs particulieres aux elements a visualiser, un par maillage
deb_list_GrandParticuliere_element DEF_EPAISSEUR DEF_LARGEUR DIRECTION_PLI fin_list_GrandParticuliere_element
fin_isovaleur_Gmsh # mot cle de fin des parametres pour les isovaleurs
# ----------------------------- definition des parametres de deformee: ----------------
debut_deformee # un mot cle de debut de liste
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
# definition des alertes: deb_list_alerte
# un mot clef de debut
# puis deux nombres: un mini et un maxi, et un nom
# un mot clef de fin: fin_list_alerte
deb_list_alerte
fin_list_alerte
avec_vitesse= 0 avec_acceleration= 0 # 1 ou 0 pour la sortie si disponible
fin_deformee # un mot cle de fin
# ----------------------------- 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)
tous_les_increments 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
fin_visualisation_Gmsh
# =================================================================================
# || fin de la visualisation Gmsh ||
# =================================================================================
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
###############################################################################################

View file

@ -0,0 +1,159 @@
#######################################################################
# #
# | | ==== === ==== ==== | | | | #
# | | | | | | / | | | | #
# |====| |=== === |=== / |====| ------- ------- #
# | | | | \ | / | | | | #
# | | ==== | \ ==== ==== | | | | #
# #
#######################################################################
# Copyright (c) 1997-2015 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.710
#######################################################################
# fichier de commande #
# --- test de la loi critere pour des plis dans le cas d'un cisaillement
# avec un amortissement visqueux et une loi élastique lineaire
#----------------------------------------------------
# definition de la dimension de l'espace de travail |
#----------------------------------------------------
dimension 3
#---------------------------------------------------------------
# definition facultative du niveau d'impression (entre 0 et 10)|
#---------------------------------------------------------------
niveau_commentaire 3
#-------------------------------
# definition du type de calcul |
#-------------------------------
TYPE_DE_CALCUL
dynamique_relaxation_dynam #avec plus visualisation
< oDR1.algo
#--------------------------------------
#| definition du (ou des) maillage(s) |
#--------------------------------------
# -- def maillage
< tria_2X3_10x10.her
#======================================
#| definition des lois de comportement|
#--------------------------------------
choix_materiaux ----------
#----------------------------------------
# Elements | Nom Materiau |
#----------------------------------------
E_tout plis_elastique
materiaux #----------
#------------------------------------------
# Nom Materiau | Type loi |
#------------------------------------------
plis_elastique LOI_CRITERE
# ....... loi de comportement avec critere ........
TYPE_DE_CRITERE_ PLISSEMENT_MEMBRANE
LOI_CONTRAINTES_PLANES # loi de contrainte plane
NEWTON_LOCAL avec_parametres_de_reglage_
nb_iteration_maxi_ 20
nb_dichotomie_maxi_ 20
tolerance_residu_ 1.e-3
permet_affichage_ 0 #5
fin_parametres_reglage_Algo_Newton_
ISOELAS
2000 0.3
fin_loi_contrainte_plane # --- fin de la loi de contrainte plane
fin_loi_interne # --- fin des lois internes
parametres_controle_pour_contraintes_planes_double_
NEWTON_LOCAL avec_parametres_de_reglage_
nb_iteration_maxi_ 20
nb_dichotomie_maxi_ 20
tolerance_residu_ 0.001 # MPa
tolerance_residu_rel_ 1.e-3
mini_bsurb0_ 0.001
mini_bsurb0_ 0.001
maxi_hsurh0_ 1.0 #4
maxi_bsurb0_ 1.
sortie_post_ 1
permet_affichage_ 0 #4
fin_parametres_reglage_Algo_Newton_
fin_parametres_controle_pour_contraintes_planes_double_
avec_parametres_controle_pour_loi_critere_
permet_affichage_ 0 #5
fin_parametres_controle_pour_loi_critere_
fin_loi_critere # ----- fin de LOI_CRITERE
#--------------------------------------- fin def des lois de comportement ------------------------
# --- divers stockages (1) -------
epaisseurs #--------#
E_tout 1
masse_volumique #--------#
E_tout 1
charges #------------#
blocages #------------#
#-----------------------------------------------------------
# nom du maillage | Ref noeud | Bloquages
#-----------------------------------------------------------
N_O UX UY UZ
N_E ' UY = 0.1 ' UX
N_tout UZ
controle #------------
#---------------------------
# PARAMETRE | VALEUR |
#---------------------------
SAUVEGARDE 1
PRECISION 0.002
DELTAt 0.2#2
MAXINCRE 1#3
ITERATIONS 1000000000
SAUVEGARDE 1
NORME E_cinetique/E_statique_ET_Res/Reac_et_Fext # min(Res,Res/Reaction) #
para_affichage ------------
#--------------------#
# PARAMETRE | VALEUR #
#--------------------#
#FREQUENCE_AFFICHAGE_INCREMENT 100
FREQUENCE_AFFICHAGE_ITERATION 1000
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
# ------------------------------------------------------------------------------------
resultats pas_de_sortie_finale_
COPIE 0
#
_fin_point_info_

View file

@ -0,0 +1,102 @@
#fichier au format maple6
###############################################################################################
# Visualisation elements finis : Herezh++ V6.844 #
# Copyright (c) 1997-2018, 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 ------------------------------------
#1 (nombre de grandeurs globales) [2]energie_elastique
#---------------------------------- maillage et dimension --------------------------------
#1 3 (nombre de maillages et dimension)
#---------------------------------- torseurs de reactions --------------------------------
#1 6 (nombre de torseurs et nombre total de grandeurs associees)
# N_E [4]Rx [5]Ry [6]Rz [7]Mx [8]My [9]Mz ;
#
#-------- 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 --------------------------------
#0 0 (nombre de noeuds, nombre total de grandeurs associees)
#---------------------------------- grandeurs aux elements ------------------------------
#4 6 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
# element_39 pt_integ_1: [11]X [12]Y [13]Z [14] SIG11 [15] SIG22 [16] SIG12 [17] EPS11 [18] EPS22 [19] EPS12
# element_110 pt_integ_1: [20]X [21]Y [22]Z [23] SIG11 [24] SIG22 [25] SIG12 [26] EPS11 [27] EPS22 [28] EPS12
# element_199 pt_integ_1: [29]X [30]Y [31]Z [32] SIG11 [33] SIG22 [34] SIG12 [35] EPS11 [36] EPS22 [37] EPS12
# element_4 pt_integ_1: [38]X [39]Y [40]Z [41] SIG11 [42] SIG22 [43] SIG12 [44] EPS11 [45] EPS22 [46] EPS12
#
#
#====================================================================
#|| 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
2.000000000000e-01 1.004084483619e-01 2.000000000000e-01 -1.319223511620e+01 -1.020405396747e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.711477986779e+00 2.000000000000e-01 2.661416561885e-01 2.799711644625e+00 0.000000000000e+00 2.993635083204e-01 1.038701489699e-02 5.576283009802e-02 1.085659254954e-04 -1.180045845983e-03 2.486591155702e-04 1.133163584008e+00 1.411999454399e+00 0.000000000000e+00 5.531074508322e+00 4.096622612176e+00 4.760118160342e+00 -1.008548440097e-03 -3.047288528760e-03 6.765401995238e-03 1.866041249000e+00 2.817892783088e+00 0.000000000000e+00 5.950282359006e+00 3.016084456303e+00 4.236337348891e+00 1.113213148226e-03 -2.165262740769e-03 4.733399149114e-03 1.334402295306e-01 5.016644995061e-01 0.000000000000e+00 5.537158057183e+00 3.558443761819e+00 4.438881114289e+00 3.434884799205e-04 -1.994368749390e-03 5.244552145503e-03

View file

@ -0,0 +1,11 @@
PARA_TYPE_DE_CALCUL
# ----- oDR
typeCalRelaxation= 2 lambda= 15. type_calcul_mass= 2 option_recalcul_mass= 0.
parametre_calcul_de_la_masse_ casMass_relax= 3
parametre_recalcul_de_la_masse_ fac_epsilon= 100.
parametre_calcul_de_la_viscosite_ type_calcul_visqu_critique= 1 opt_cal_C_critique= 1 f_= 0.9
#mode_debug_= 100
ARRET_A_EQUILIBRE_STATIQUE_ 2

View file

@ -0,0 +1,437 @@
#############################################################################
# #
# STAMM #
# (maillage automatique de pieces simples) #
#############################################################################
# Copyright © 1997-2013 Gerard Rio, gerard.rio@univ-ubs.fr #
# All rights reserved. http://www-lg2m.univ-ubs.fr/~rio #
# #
#############################################################################
# version 03.29 , compatible avec la version Herezh++6.576
# ----------------------------------------------------------------
# | construction automatique d'un maillage 2D |
# | pour une utilisation par herezh++ |
# ----------------------------------------------------------------
# rectangle de dimension : 2 x 3 ( 10 x 10 elements )
# geometrie rectangulaire, decoupage triangulaire, interpolation lineaire.
nom_maillage tria_2X3_10x10
noeuds ------------
121 NOEUDS
#---------------------------------------------------------------
#|NO DU| X | Y | Z |
#|NOEUD| | | |
#---------------------------------------------------------------
1 0 0 0
2 0 0.3 0
3 0 0.6 0
4 0 0.8999999999999999 0
5 0 1.2 0
6 0 1.5 0
7 0 1.8 0
8 0 2.1 0
9 0 2.4 0
10 0 2.7 0
11 0 3 0
12 0.2 0 0
13 0.2 0.3 0
14 0.2 0.6 0
15 0.2 0.8999999999999999 0
16 0.2 1.2 0
17 0.2 1.5 0
18 0.2 1.8 0
19 0.2 2.1 0
20 0.2 2.4 0
21 0.2 2.7 0
22 0.2 3 0
23 0.4 0 0
24 0.4 0.3 0
25 0.4 0.6 0
26 0.4 0.8999999999999999 0
27 0.4 1.2 0
28 0.4 1.5 0
29 0.4 1.8 0
30 0.4 2.1 0
31 0.4 2.4 0
32 0.4 2.7 0
33 0.4 3 0
34 0.6000000000000001 0 0
35 0.6000000000000001 0.3 0
36 0.6000000000000001 0.6 0
37 0.6000000000000001 0.8999999999999999 0
38 0.6000000000000001 1.2 0
39 0.6000000000000001 1.5 0
40 0.6000000000000001 1.8 0
41 0.6000000000000001 2.1 0
42 0.6000000000000001 2.4 0
43 0.6000000000000001 2.7 0
44 0.6000000000000001 3 0
45 0.8 0 0
46 0.8 0.3 0
47 0.8 0.6 0
48 0.8 0.8999999999999999 0
49 0.8 1.2 0
50 0.8 1.5 0
51 0.8 1.8 0
52 0.8 2.1 0
53 0.8 2.4 0
54 0.8 2.7 0
55 0.8 3 0
56 1 0 0
57 1 0.3 0
58 1 0.6 0
59 1 0.8999999999999999 0
60 1 1.2 0
61 1 1.5 0
62 1 1.8 0
63 1 2.1 0
64 1 2.4 0
65 1 2.7 0
66 1 3 0
67 1.2 0 0
68 1.2 0.3 0
69 1.2 0.6 0
70 1.2 0.8999999999999999 0
71 1.2 1.2 0
72 1.2 1.5 0
73 1.2 1.8 0
74 1.2 2.1 0
75 1.2 2.4 0
76 1.2 2.7 0
77 1.2 3 0
78 1.4 0 0
79 1.4 0.3 0
80 1.4 0.6 0
81 1.4 0.8999999999999999 0
82 1.4 1.2 0
83 1.4 1.5 0
84 1.4 1.8 0
85 1.4 2.1 0
86 1.4 2.4 0
87 1.4 2.7 0
88 1.4 3 0
89 1.6 0 0
90 1.6 0.3 0
91 1.6 0.6 0
92 1.6 0.8999999999999999 0
93 1.6 1.2 0
94 1.6 1.5 0
95 1.6 1.8 0
96 1.6 2.1 0
97 1.6 2.4 0
98 1.6 2.7 0
99 1.6 3 0
100 1.8 0 0
101 1.8 0.3 0
102 1.8 0.6 0
103 1.8 0.8999999999999999 0
104 1.8 1.2 0
105 1.8 1.5 0
106 1.8 1.8 0
107 1.8 2.1 0
108 1.8 2.4 0
109 1.8 2.7 0
110 1.8 3 0
111 2 0 0
112 2 0.3 0
113 2 0.6 0
114 2 0.8999999999999999 0
115 2 1.2 0
116 2 1.5 0
117 2 1.8 0
118 2 2.1 0
119 2 2.4 0
120 2 2.7 0
121 2 3 0
# references des noeuds
N_O 1 2 3 4 5 6 7 8 9 10 11
N_S 1 12 23 34 45 56 67 78 89 100 111
N_E 111 112 113 114 115 116 117 118 119 120 121
N_N 11 22 33 44 55 66 77 88 99 110 121
N_SO 1
N_SE 111
N_NE 121
N_NO 11
N_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96
97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112
113 114 115 116 117 118 119 120 121
# les elements
elements ----------
200 ELEMENTS
#----------------------------------------------------------------------
#| NO | | |
#|ELTS | type element | Noeuds |
#----------------------------------------------------------------------
1 TRIANGLE LINEAIRE 1 12 2
2 TRIANGLE LINEAIRE 12 13 2
3 TRIANGLE LINEAIRE 2 13 3
4 TRIANGLE LINEAIRE 13 14 3
5 TRIANGLE LINEAIRE 3 14 4
6 TRIANGLE LINEAIRE 14 15 4
7 TRIANGLE LINEAIRE 4 15 5
8 TRIANGLE LINEAIRE 15 16 5
9 TRIANGLE LINEAIRE 5 16 6
10 TRIANGLE LINEAIRE 16 17 6
11 TRIANGLE LINEAIRE 6 17 7
12 TRIANGLE LINEAIRE 17 18 7
13 TRIANGLE LINEAIRE 7 18 8
14 TRIANGLE LINEAIRE 18 19 8
15 TRIANGLE LINEAIRE 8 19 9
16 TRIANGLE LINEAIRE 19 20 9
17 TRIANGLE LINEAIRE 9 20 10
18 TRIANGLE LINEAIRE 20 21 10
19 TRIANGLE LINEAIRE 10 21 11
20 TRIANGLE LINEAIRE 21 22 11
21 TRIANGLE LINEAIRE 12 23 13
22 TRIANGLE LINEAIRE 23 24 13
23 TRIANGLE LINEAIRE 13 24 14
24 TRIANGLE LINEAIRE 24 25 14
25 TRIANGLE LINEAIRE 14 25 15
26 TRIANGLE LINEAIRE 25 26 15
27 TRIANGLE LINEAIRE 15 26 16
28 TRIANGLE LINEAIRE 26 27 16
29 TRIANGLE LINEAIRE 16 27 17
30 TRIANGLE LINEAIRE 27 28 17
31 TRIANGLE LINEAIRE 17 28 18
32 TRIANGLE LINEAIRE 28 29 18
33 TRIANGLE LINEAIRE 18 29 19
34 TRIANGLE LINEAIRE 29 30 19
35 TRIANGLE LINEAIRE 19 30 20
36 TRIANGLE LINEAIRE 30 31 20
37 TRIANGLE LINEAIRE 20 31 21
38 TRIANGLE LINEAIRE 31 32 21
39 TRIANGLE LINEAIRE 21 32 22
40 TRIANGLE LINEAIRE 32 33 22
41 TRIANGLE LINEAIRE 23 34 24
42 TRIANGLE LINEAIRE 34 35 24
43 TRIANGLE LINEAIRE 24 35 25
44 TRIANGLE LINEAIRE 35 36 25
45 TRIANGLE LINEAIRE 25 36 26
46 TRIANGLE LINEAIRE 36 37 26
47 TRIANGLE LINEAIRE 26 37 27
48 TRIANGLE LINEAIRE 37 38 27
49 TRIANGLE LINEAIRE 27 38 28
50 TRIANGLE LINEAIRE 38 39 28
51 TRIANGLE LINEAIRE 28 39 29
52 TRIANGLE LINEAIRE 39 40 29
53 TRIANGLE LINEAIRE 29 40 30
54 TRIANGLE LINEAIRE 40 41 30
55 TRIANGLE LINEAIRE 30 41 31
56 TRIANGLE LINEAIRE 41 42 31
57 TRIANGLE LINEAIRE 31 42 32
58 TRIANGLE LINEAIRE 42 43 32
59 TRIANGLE LINEAIRE 32 43 33
60 TRIANGLE LINEAIRE 43 44 33
61 TRIANGLE LINEAIRE 34 45 35
62 TRIANGLE LINEAIRE 45 46 35
63 TRIANGLE LINEAIRE 35 46 36
64 TRIANGLE LINEAIRE 46 47 36
65 TRIANGLE LINEAIRE 36 47 37
66 TRIANGLE LINEAIRE 47 48 37
67 TRIANGLE LINEAIRE 37 48 38
68 TRIANGLE LINEAIRE 48 49 38
69 TRIANGLE LINEAIRE 38 49 39
70 TRIANGLE LINEAIRE 49 50 39
71 TRIANGLE LINEAIRE 39 50 40
72 TRIANGLE LINEAIRE 50 51 40
73 TRIANGLE LINEAIRE 40 51 41
74 TRIANGLE LINEAIRE 51 52 41
75 TRIANGLE LINEAIRE 41 52 42
76 TRIANGLE LINEAIRE 52 53 42
77 TRIANGLE LINEAIRE 42 53 43
78 TRIANGLE LINEAIRE 53 54 43
79 TRIANGLE LINEAIRE 43 54 44
80 TRIANGLE LINEAIRE 54 55 44
81 TRIANGLE LINEAIRE 45 56 46
82 TRIANGLE LINEAIRE 56 57 46
83 TRIANGLE LINEAIRE 46 57 47
84 TRIANGLE LINEAIRE 57 58 47
85 TRIANGLE LINEAIRE 47 58 48
86 TRIANGLE LINEAIRE 58 59 48
87 TRIANGLE LINEAIRE 48 59 49
88 TRIANGLE LINEAIRE 59 60 49
89 TRIANGLE LINEAIRE 49 60 50
90 TRIANGLE LINEAIRE 60 61 50
91 TRIANGLE LINEAIRE 50 61 51
92 TRIANGLE LINEAIRE 61 62 51
93 TRIANGLE LINEAIRE 51 62 52
94 TRIANGLE LINEAIRE 62 63 52
95 TRIANGLE LINEAIRE 52 63 53
96 TRIANGLE LINEAIRE 63 64 53
97 TRIANGLE LINEAIRE 53 64 54
98 TRIANGLE LINEAIRE 64 65 54
99 TRIANGLE LINEAIRE 54 65 55
100 TRIANGLE LINEAIRE 65 66 55
101 TRIANGLE LINEAIRE 56 67 57
102 TRIANGLE LINEAIRE 67 68 57
103 TRIANGLE LINEAIRE 57 68 58
104 TRIANGLE LINEAIRE 68 69 58
105 TRIANGLE LINEAIRE 58 69 59
106 TRIANGLE LINEAIRE 69 70 59
107 TRIANGLE LINEAIRE 59 70 60
108 TRIANGLE LINEAIRE 70 71 60
109 TRIANGLE LINEAIRE 60 71 61
110 TRIANGLE LINEAIRE 71 72 61
111 TRIANGLE LINEAIRE 61 72 62
112 TRIANGLE LINEAIRE 72 73 62
113 TRIANGLE LINEAIRE 62 73 63
114 TRIANGLE LINEAIRE 73 74 63
115 TRIANGLE LINEAIRE 63 74 64
116 TRIANGLE LINEAIRE 74 75 64
117 TRIANGLE LINEAIRE 64 75 65
118 TRIANGLE LINEAIRE 75 76 65
119 TRIANGLE LINEAIRE 65 76 66
120 TRIANGLE LINEAIRE 76 77 66
121 TRIANGLE LINEAIRE 67 78 68
122 TRIANGLE LINEAIRE 78 79 68
123 TRIANGLE LINEAIRE 68 79 69
124 TRIANGLE LINEAIRE 79 80 69
125 TRIANGLE LINEAIRE 69 80 70
126 TRIANGLE LINEAIRE 80 81 70
127 TRIANGLE LINEAIRE 70 81 71
128 TRIANGLE LINEAIRE 81 82 71
129 TRIANGLE LINEAIRE 71 82 72
130 TRIANGLE LINEAIRE 82 83 72
131 TRIANGLE LINEAIRE 72 83 73
132 TRIANGLE LINEAIRE 83 84 73
133 TRIANGLE LINEAIRE 73 84 74
134 TRIANGLE LINEAIRE 84 85 74
135 TRIANGLE LINEAIRE 74 85 75
136 TRIANGLE LINEAIRE 85 86 75
137 TRIANGLE LINEAIRE 75 86 76
138 TRIANGLE LINEAIRE 86 87 76
139 TRIANGLE LINEAIRE 76 87 77
140 TRIANGLE LINEAIRE 87 88 77
141 TRIANGLE LINEAIRE 78 89 79
142 TRIANGLE LINEAIRE 89 90 79
143 TRIANGLE LINEAIRE 79 90 80
144 TRIANGLE LINEAIRE 90 91 80
145 TRIANGLE LINEAIRE 80 91 81
146 TRIANGLE LINEAIRE 91 92 81
147 TRIANGLE LINEAIRE 81 92 82
148 TRIANGLE LINEAIRE 92 93 82
149 TRIANGLE LINEAIRE 82 93 83
150 TRIANGLE LINEAIRE 93 94 83
151 TRIANGLE LINEAIRE 83 94 84
152 TRIANGLE LINEAIRE 94 95 84
153 TRIANGLE LINEAIRE 84 95 85
154 TRIANGLE LINEAIRE 95 96 85
155 TRIANGLE LINEAIRE 85 96 86
156 TRIANGLE LINEAIRE 96 97 86
157 TRIANGLE LINEAIRE 86 97 87
158 TRIANGLE LINEAIRE 97 98 87
159 TRIANGLE LINEAIRE 87 98 88
160 TRIANGLE LINEAIRE 98 99 88
161 TRIANGLE LINEAIRE 89 100 90
162 TRIANGLE LINEAIRE 100 101 90
163 TRIANGLE LINEAIRE 90 101 91
164 TRIANGLE LINEAIRE 101 102 91
165 TRIANGLE LINEAIRE 91 102 92
166 TRIANGLE LINEAIRE 102 103 92
167 TRIANGLE LINEAIRE 92 103 93
168 TRIANGLE LINEAIRE 103 104 93
169 TRIANGLE LINEAIRE 93 104 94
170 TRIANGLE LINEAIRE 104 105 94
171 TRIANGLE LINEAIRE 94 105 95
172 TRIANGLE LINEAIRE 105 106 95
173 TRIANGLE LINEAIRE 95 106 96
174 TRIANGLE LINEAIRE 106 107 96
175 TRIANGLE LINEAIRE 96 107 97
176 TRIANGLE LINEAIRE 107 108 97
177 TRIANGLE LINEAIRE 97 108 98
178 TRIANGLE LINEAIRE 108 109 98
179 TRIANGLE LINEAIRE 98 109 99
180 TRIANGLE LINEAIRE 109 110 99
181 TRIANGLE LINEAIRE 100 111 101
182 TRIANGLE LINEAIRE 111 112 101
183 TRIANGLE LINEAIRE 101 112 102
184 TRIANGLE LINEAIRE 112 113 102
185 TRIANGLE LINEAIRE 102 113 103
186 TRIANGLE LINEAIRE 113 114 103
187 TRIANGLE LINEAIRE 103 114 104
188 TRIANGLE LINEAIRE 114 115 104
189 TRIANGLE LINEAIRE 104 115 105
190 TRIANGLE LINEAIRE 115 116 105
191 TRIANGLE LINEAIRE 105 116 106
192 TRIANGLE LINEAIRE 116 117 106
193 TRIANGLE LINEAIRE 106 117 107
194 TRIANGLE LINEAIRE 117 118 107
195 TRIANGLE LINEAIRE 107 118 108
196 TRIANGLE LINEAIRE 118 119 108
197 TRIANGLE LINEAIRE 108 119 109
198 TRIANGLE LINEAIRE 119 120 109
199 TRIANGLE LINEAIRE 109 120 110
200 TRIANGLE LINEAIRE 120 121 110
# references des elements
A_O 1 3 3 3 5 3 7 3 9 3 11 3 13 3 15 3 17 3 19 3
A_S 1 1 21 1 41 1 61 1 81 1 101 1 121 1 141 1 161 1 181 1
A_E 182 1 184 1 186 1 188 1 190 1 192 1 194 1 196 1 198 1 200 1
A_N 20 2 40 2 60 2 80 2 100 2 120 2 140 2 160 2 180 2 200 2
F_to 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1
12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1
23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1
34 1 35 1 36 1 37 1 38 1 39 1 40 1 41 1 42 1 43 1 44 1
45 1 46 1 47 1 48 1 49 1 50 1 51 1 52 1 53 1 54 1 55 1
56 1 57 1 58 1 59 1 60 1 61 1 62 1 63 1 64 1 65 1 66 1
67 1 68 1 69 1 70 1 71 1 72 1 73 1 74 1 75 1 76 1 77 1
78 1 79 1 80 1 81 1 82 1 83 1 84 1 85 1 86 1 87 1 88 1
89 1 90 1 91 1 92 1 93 1 94 1 95 1 96 1 97 1 98 1 99 1
100 1 101 1 102 1 103 1 104 1 105 1 106 1 107 1 108 1 109 1 110 1
111 1 112 1 113 1 114 1 115 1 116 1 117 1 118 1 119 1 120 1 121 1
122 1 123 1 124 1 125 1 126 1 127 1 128 1 129 1 130 1 131 1 132 1
133 1 134 1 135 1 136 1 137 1 138 1 139 1 140 1 141 1 142 1 143 1
144 1 145 1 146 1 147 1 148 1 149 1 150 1 151 1 152 1 153 1 154 1
155 1 156 1 157 1 158 1 159 1 160 1 161 1 162 1 163 1 164 1 165 1
166 1 167 1 168 1 169 1 170 1 171 1 172 1 173 1 174 1 175 1 176 1
177 1 178 1 179 1 180 1 181 1 182 1 183 1 184 1 185 1 186 1 187 1
188 1 189 1 190 1 191 1 192 1 193 1 194 1 195 1 196 1 197 1 198 1
199 1 200 1
E_O 1 3 5 7 9 11 13 15 17 19
E_S 1 21 41 61 81 101 121 141 161 181
E_E 182 184 186 188 190 192 194 196 198 200
E_N 20 40 60 80 100 120 140 160 180 200
E_SO 1
E_SE 181 182
E_NE 200
E_NO 19 20
F_O 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1
12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1
F_S 1 1 2 1 21 1 22 1 41 1 42 1 61 1 62 1 81 1 82 1 101 1
102 1 121 1 122 1 141 1 142 1 161 1 162 1 181 1 182 1
F_E 91 1 92 1 93 1 94 1 95 1 96 1 97 1 98 1 99 1 100 1 101 1
102 1 103 1 104 1 105 1 106 1 107 1 108 1 109 1 110 1
F_N 19 1 20 1 39 1 40 1 59 1 60 1 79 1 80 1 99 1 100 1 119 1
120 1 139 1 140 1 159 1 160 1 179 1 180 1 199 1 200 1
E_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96
97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112
113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128
129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160
161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176
177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192
193 194 195 196 197 198 199 200