#######################################################################
#                                                                     #
#       |    |  ====  ===   ====  ====  |    |      |        |        #
#       |    |  |     |  |  |        /  |    |      |        |        #
#       |====|  |===  ===   |===    /   |====|   -------  -------     #
#       |    |  |     | \   |      /    |    |      |        |        #
#       |    |  ====  |  \  ====  ====  |    |      |        |        #
#                                                                     #
#######################################################################
#  Copyright (c) 1997-2013 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  7.037
#######################################################################
#                       fichier de commande                           # 
#     --- commentaire sur l'objet du travail ???

#----------------------------------------------------
# definition de la dimension de l'espace de travail |
#----------------------------------------------------
 dimension   3

#---------------------------------------------------------------
# definition facultative du niveau d'impression (entre 0 et 10)|
#---------------------------------------------------------------
 niveau_commentaire   1

#-------------------------------
# definition du type de calcul |
#-------------------------------
TYPE_DE_CALCUL
utilitaires avec plus renumerotation_des_noeuds plus sauveMaillagesEnCours

#--------------------------------------
#| definition du (ou des) maillage(s) |
#--------------------------------------

 #-- def maillage 
< Joint_1x1.her
< Chemise_1x1.her
< Tige_1x1.her

domaine_esclave
1

les_courbes_1D ------
< ../1_Donnees/courbe_1D.hz_courbe

les_fonctions_nD #------------
< ../1_Donnees/fct_nD_ponderation_CD.hz_courbe

#======================================
#| definition des lois de comportement|
#--------------------------------------
 
            choix_materiaux ----------
#---------------------------------------
# nom de maillage |Elements | Materiau |
#---------------------------------------
nom_mail= Joint_1x1         E_tout	graphite
nom_mail= Chemise_1x1       E_tout	acier
nom_mail= Tige_1x1          E_tout	acier

		materiaux #---------

#< graphite.hz_loi	

#####
graphite LOI_ADDITIVE_EN_SIGMA

avec_fonction_de_ponderation_

#---------------------------------------
# loi 1: hyperbulk
#---------------------------------------
les_grandeur_ponderation= def_duale_mises_maxi  AuPti_  fin_grandeur_ponderation_
deb_fonct_ponder=   courbe_constante    fin_fonct_ponder_
avec_ponder_grandeur_locale_
les_grandeurs_=  V_vol  fin_grandeurs_
deb_fonct_=   courbe_bulk   fin_fonct_
# .......  loi de comportement 3D hyperelastique isotrope ISOHYPERBULK_GENE
ISOHYPERBULK3
#-------------------------"
#     K   |  fonction     |
#-------------------------"
550 K_V_dependant_ CPL1D DdlP 0. 1.0   1.0 1.0 FdlP sortie_post_ 1


#---------------------------------------
# loi 2: Orgeas sans partie spherique
#---------------------------------------
les_grandeur_ponderation= def_duale_mises_maxi  AuPti_  fin_grandeur_ponderation_
deb_fonct_ponder=   courbe_constante    fin_fonct_ponder_
avec_ponder_grandeur_locale_
les_grandeurs_=  V_vol  fin_grandeurs_
deb_fonct_=   courbe_orgeas   fin_fonct_
ISOHYPER3DORGEAS1
#K Qs mu1 mu2 mu3 alpha1 alpha2 Qe
0 10.0 0.1 1. 11.5 0.1 0.1 0.04 sortie_post_ 1


#---------------------------------------
# loi 3: Hysteresis deviatorique
#---------------------------------------
#         definition de la premiere loi
#-------- debut cas d une premiere loi hysteresis_3D --------
les_grandeur_ponderation= def_duale_mises_maxi  AuPti_  fin_grandeur_ponderation_
deb_fonct_ponder=   courbe_constante    fin_fonct_ponder_
avec_ponder_grandeur_locale_
    les_grandeurs_=  V_vol  fin_grandeurs_
    deb_fonct_=   courbe_hysteresis_3D   fin_fonct_

HYSTERESIS_3D
#np= 5.75       mu= 1000    Qzero= 3
np= 1.0       mu= 2800    Qzero= 1.75   avec_parametres_de_reglage_
# tous les parametres qui suivent doivent etre sur une seule ligne contrairement a l exemple qui suit
type_de_resolution_  5  type_calcul_comportement_tangent_ 2
cas_kutta_  5  erreurAbsolue_  1.e-4  erreurRelative_ 1.e-3#5
#cas_kutta_  5  erreurAbsolue_  1.e-4  erreurRelative_ 1.e-3#5
nbMaxiAppel_ 3000#3000#1600

nb_iteration_maxi_ 30
nb_dichotomie_maxi_ 20
tolerance_residu_     1.e-4
tolerance_residu_rel_ 1.e-3
nbMaxiAppel_  3600

tolerance_coincidence_  1.e-3#4#4#
nb_maxInvCoinSurUnPas_  -20
mini_Qsig_pour_phase_sigma_Oi_tdt_ 0.1 #1 #0 #10
mini_rayon_  1.e-12  min_angle_trajet_neutre_ 1.e-2
possibilite_cosAlphaNegatif_ 0
permet_affichage_ 5
type_de_transport_memorisation_ -1
fin_parametres_reglage_Hysteresis_
#-------- fin cas d une premiere loi hysteresis_3D --------
###

#---------------------------------------
# loi 4: Hysteresis bulk
#---------------------------------------
#         definition de la seconde loi
#-------- debut cas d une premiere loi hysteresis_bulk --------
les_grandeur_ponderation= def_duale_mises_maxi  AuPti_  fin_grandeur_ponderation_
deb_fonct_ponder=   courbe_constante    fin_fonct_ponder_
avec_ponder_grandeur_locale_
les_grandeurs_=  V_vol  fin_grandeurs_
deb_fonct_=   courbe_hysteresis_bulk   fin_fonct_

HYSTERESIS_BULK #hysteresis spherique
np= 1.75       mu= 850    Qzero= 16.0   avec_parametres_de_reglage_
type_de_resolution_  2
cas_kutta_  5
erreurAbsolue_  1.e-4 #1.e-4 #1.e-5
erreurRelative_ 1.e-3 #1.e-3 #1.e-4
#erreurAbsolue_  1.e-8  erreurRelative_ 1.e-8
nbMaxiAppel_  3600
nb_iteration_maxi_  30  nb_dichotomie_maxi_  20
tolerance_residu_ 1.e-3 tolerance_residu_rel_ 1.e-5#5
depassement_Q0_ 1.1 #1.1#1000
sortie_post_ 1
permet_affichage_ 5
tolerance_coincidence_ 1.e-5 #1.e-8 #1.e-10
nb_maxInvCoinSurUnPas_ -20
sortie_post_ 1
permet_affichage_ 5#0
fin_parametres_reglage_Hysteresis_
#-------- debut cas d une premiere loi hysteresis_bulk --------


fin_liste_lois_elementaires
		
#------------------------------------------
# Nom Materiau    |      Type loi         |
#------------------------------------------
  acier         ISOELAS
#------------------------------------------
# E    |     Nu
#------------------------------------------
  190000     0.3

		masse_volumique ------------
nom_mail= Joint_1x1         E_tout	1. 
nom_mail= Chemise_1x1       E_tout	1.
nom_mail= Tige_1x1          E_tout	1.

  zone_contact
#----------------------------------------------------
# ref noeuds | ref face |
#----------------------------------------------------
nom_mail= Joint_1x1 N_ext		nom_mail= Chemise_1x1   A_contact
nom_mail= Joint_1x1 N_int		nom_mail= Tige_1x1      A_contact

		charges ------------
  
            blocages ------------
#-----------------------------------------------------------
# Ref noeud | Bloquages
#-----------------------------------------------------------
#nom_mail= Joint_4x4_quadracompl	N_tout	UZ
#nom_mail= Chemise_quadracompl	N_tout	UZ
#nom_mail= Tige			N_tout	UZ

#nom_mail= Joint_4x4_quadracompl	N_ext		UX
#nom_mail= Joint_4x4_quadracompl	N_int		UX
  
#nom_mail= Joint_4x4_quadracompl	N_bas		UY
  
#nom_mail= Chemise_quadracompl	N_bloq	UX, UY
#nom_mail= Tige			N_tout	UX, UY
  
#nom_mail= Joint_4x4_quadracompl	N_haut	'UY= COURBE_CHARGE: charge_dechargeY ECHELLE: -0.95'

#		condition_limite_lineaire_ ------------
#nom_mail= Joint_4x4_quadracompl	N_tout	enu= UX
#	def_auto_coef_planOuDroite_ centre_noeud_a_t0_ nom_mail= Tige coefficients= 1. 0. 0.4 fin_list_coefficients_
	
#nom_mail= Chemise_quadracompl	N_tout	enu= UZ
#nom_mail= Tige			N_tout	enu= UX, UY, UZ


#nom_mail= mail1 N_haut_avant_gauche enu= UX
#    refs_associe= nom_mail= mail2 N_bas_arriere_droit fin_list_refs_associe_
#    val_condi_lineaire= 0. coefficients= 1. 0. 0.  1. 0. 0. fin_list_coefficients_
  
            controle ------------
#---------------------------
# PARAMETRE    | VALEUR    |
#---------------------------
SAUVEGARDE        1
TEMPSFIN          2.0
DELTAtMINI	  	1.e-6
DELTAtMAXI        0.005
DELTAt            0.0001
ITERATIONS        15
#MAXINCRE  		9999999.
PRECISION         1.e-3
NORME 		Residu/Reaction

            para_affichage ------------
#---------------------------
# PARAMETRE    | VALEUR    |
#---------------------------
FREQUENCE_SORTIE_FIL_DU_CALCUL        1
    
            resultats pas_de_sortie_finale_ ------------
COPIE 0

  _fin_point_info_  
 
  
