#######################################################################
#                                                                     #
#       |    |  ====  ===   ====  ====  |    |      |        |        #
#       |    |  |     |  |  |        /  |    |      |        |        #
#       |====|  |===  ===   |===    /   |====|   -------  -------     #
#       |    |  |     | \   |      /    |    |      |        |        #
#       |    |  ====  |  \  ====  ====  |    |      |        |        #
#                                                                     #
#######################################################################
#  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   3

#-------------------------------
# definition du type de calcul |
#-------------------------------
TYPE_DE_CALCUL
#non_dynamique avec plus visualisation
non_dynamique avec plus uniquement_remontee #plus visualisation

#--------------------------------------
#| definition du (ou des) maillage(s) |
#--------------------------------------
#-- def maillage 
< 0_Maillages/Joint_1x1.her
< 0_Maillages/Joint_1x1.lis
< 0_Maillages/Chemise_1x1.her
< 0_Maillages/Chemise_1x1.lis
< 0_Maillages/Tige_1x1.her
< 0_Maillages/Tige_1x1.lis

    domaine_esclave
1

    les_courbes_1D 
< 1_Donnees/courbe_1D.hz_courbe

    les_fonctions_nD 
< 1_Donnees/fct_nD_ponderation_CD.hz_courbe

  
  fc_com_contact FONCTION_EXPRESSION_LITTERALE_nD
     un_argument= NUM_NOEUD  
     un_argument= NUM_ELEMENT
#     permet_affichage_  5
#    fct= (NUM_NOEUD == 0) ? 5  \ # cas sans numéro de noeud
#          : (NUM_NOEUD == 120) ? 5  \ # cas pour un numéro donné
#            : 0 
#     fct= ((NUM_NOEUD == 123) || (NUM_NOEUD == 117)) ? 6  : 5
     fct= (NUM_NOEUD == 8)  ? 9  : 2#5
  fin_parametres_fonction_expression_litterale_
  
  fc_com_lescontacts FONCTION_EXPRESSION_LITTERALE_nD
     un_argument= NUM_NOEUD  
     un_argument= NUM_ELEMENT
#     permet_affichage_  5
#    fct= (NUM_NOEUD == 0) ? 5  \ # cas sans numéro de noeud
#          : (NUM_NOEUD == 120) ? 5  \ # cas pour un numéro donné
#            : 0 
#     fct= ((NUM_NOEUD == 123) || (NUM_NOEUD == 117)) ? 6  : 5
     fct= (NUM_NOEUD == 8)  ? 9  : 2#5
  fin_parametres_fonction_expression_litterale_
  
    fct_bascule  FONCTION_EXPRESSION_LITTERALE_nD
    un_argument= compteur_iteration_algo_global
    fct= (compteur_iteration_algo_global >0) ? 1 :  0, 2
    fin_parametres_fonction_expression_litterale_



#======================================
#| 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= Chemise_1x1		E_tout	loi_frottement
nom_mail= Tige_1x1			E_tout	acier
#nom_mail= Tige_1x1			E_tout	loi_frottement

    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 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
np= 1.75       mu= 850    Qzero= 16.0   avec_parametres_de_reglage_
#np= 1.75       mu= 85.0    Qzero= 16.0   avec_parametres_de_reglage_
#np= 1.75       mu= 8.50    Qzero= 16.0   avec_parametres_de_reglage_
#np= 1.75       mu= 0.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 --------


#---------------------------------------
# loi 4: 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= 1.0       mu= 2800    Qzero= 1.75
np= 1.0       mu= 2800    Qzero= 1.75   avec_parametres_de_reglage_
#np= 1.0       mu= 280.0    Qzero= 1.75   avec_parametres_de_reglage_
#np= 1.0       mu= 28.00    Qzero= 1.75   avec_parametres_de_reglage_
#np= 1.0       mu= 2.800    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 --------
###


fin_liste_lois_elementaires
		
#------------------------------------------
# Nom Materiau    |      Type loi         |
#------------------------------------------
  acier         ISOELAS
#------------------------------------------
# E    |     Nu
#------------------------------------------
190000      0.3
  
  
# -------------------------------------------------------------------
# |...... loi de comportement de frottement de coulomb ....... |
# -------------------------------------------------------------------
# mu_statique= le coefficient de frottement statique
#loi_frottement LOI_COULOMB
#	mu_statique= 0.3 
#fin_loi_frottement_coulomb_


    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_1x1         N_tout	UZ
nom_mail= Chemise_1x1		N_tout	UZ
nom_mail= Tige_1x1			N_tout	UZ

#nom_mail= Joint_1x1	N_ext   UX
#nom_mail= Joint_1x1	N_int   UX
nom_mail= Joint_1x1     N_bas   UY

nom_mail= Chemise_1x1		N_tout	UX, UY
nom_mail= Tige_1x1			N_tout	UX, UY

nom_mail= Joint_1x1	N_haut	'UY= COURBE_CHARGE: charge_dechargeY ECHELLE: -0.94'



    controle ------------
#---------------------------
# PARAMETRE    | VALEUR    |
#---------------------------
SAUVEGARDE      1			                                # defaut : 1
DELTAt          0.0001		                                # defaut : 1.
DELTAtMAXI      0.002			                            # defaut : 1.
DELTAtMINI      1.e-6			                            # defaut : 0.
TEMPSFIN        1.00			                            # defaut : 1.
PRECISION       1.e-3			                            # defaut : 1.e-3
#NORME           Residu/Reactionet_verification_ddouble	    # defaut : Residu/Reaction
NORME           min(Res,Res/MaX(Reaction_et_PVExterne))_verification_ddouble	
ITERATIONS      50#20			                                # defaut : 100
MAXINCRE  		10000#999999		                                # defaut : 400

    para_syteme_lineaire ------------
#TYPE_MATRICE        BANDE_NON_SYMETRIQUE_LAPACK	    # defaut : BANDE_SYMETRIQUE
TYPE_MATRICE     CARREE_LAPACK
SYMETRIE_MATRICE    0					            # defaut : 1
TYPE_RESOLUTION     GAUSS					        # defaut : CHOLESKY

	para_pilotage_equi_global ------------
FACTEUR_DIMINUTION              	1.732		# defaut : 1.732
FACTEUR_AUGMENTATION            	1.414		# defaut : 1.414
NB_BONNE_CONVERGENCE 			    3		    # defaut : 3
NB_ITER_POUR_BONNE_CONVERGENCE  	5		    # defaut : 0.25 ITERATIONS
NB_ITER_POUR_MAUVAISE_CONVERGENCE 	10		    # defaut : 0.5 ITERATIONS
NORME_MAXI_INCREMENT  			    0.1#0.001		# defaut : +infini		<--
CAS_JACOBIEN_NEGATIF      		    2		    # defaut : 1

    para_contact ------------
CONTACT_TYPE				            2
PROP_VALPROPRE_RAIDEUR      0.01
TYPE_PENALISATION_PENETRATION  12

#CONTACT_TYPE				                   4		    # defaut : 0
PENETRATION_CONTACT_MAXI 		        0.0001	    # defaut : 0.1
#FCT_ND_BASCUL_CONTACT_TYPE_4       fct_bascule
PENETRATION_BORNE_REGULARISATION_PLUS  0.

#PENALISATION_PENETRATION		        0.1		    # defaut : 0.1
#TYPE_PENALISATION_PENETRATION 	        2#8		    # defaut : 2

PENALISATION_TANGENTIELLE		        1.e6		# defaut : 1.e6
TYPE_PENALISATION_TANGENTIELLE	       2# 8		    # defaut : 2
TYPE_DE_DECOLLEMENT     		        1		    # defaut : 0
NB_DECOLLEMENT_MAXI			            5		    # defaut : 1
PRECISION_BOITE_PRELOCALISATION	        1.05		# defaut : 1.05
MINI_EXTRA_BOITE_PRELOCALISATION	    0.001	    # defaut : 0.001
AJOUT_EXTRA_BOITE_PRELOCALISATION	    0.		    # defaut : 0.
TANGENTIELLE_CONTACT_MAXI		        0.01 	    # defaut : 0.1
PENETRATION_BORNE_REGULARISATION 	    0.001		# defaut : +infini
TANGENTIELLE_BORNE_REGULARISATION	    100.0		# defaut : +infini
FORCE_CONTACT_NOEUD_MAXI 		        1.e8		# defaut : +infini
FORCE_TANGENTIELLE_NOEUD_MAXI 	        1.e8		# defaut : +infini
DISTANCE_MAXI_AU_PT_PROJETE		        10.0		# defaut : +infini
NB_CHANGE_FRONTIERE_MAX			        4		    # defaut : /
OPTIMISATION_NUMEROTATION 		        1		    # defaut : 0
#NIVEAU_COMMENTAIRE_LESCONTACTS	        5		    # defaut : 0
#NIVEAU_COMMENTAIRE_CONTACT		        5		    # defaut : 0

 NIVEAU_COMMENTAIRE_LESCONTACTS     2#FCT_ND_NIVEAU_COMMENTAIRE_LESCONTACTS  fc_com_lescontacts #6#3
 NIVEAU_COMMENTAIRE_CONTACT         2#FCT_ND_NIVEAU_COMMENTAIRE_CONTACTS  fc_com_contact #6#3# 3


	para_affichage ------------
FREQUENCE_SORTIE_FIL_DU_CALCUL        1	
FREQUENCE_AFFICHAGE_INCREMENT	30 
#FREQUENCE_AFFICHAGE_ITERATION	30


	para_calculs_geometriques ------------
#POINT_INTERNE_DELTA_THETAI_MAXI			1.e-2		# defaut : 1e-4
POINT_INTERNE_PREC_THETAI_INTERNE			1.e-4		# defaut : 1e-6
#POINT_INTERNE_NB_BOUCLE_SUR_DELTA_THETAI   16		    # defaut : 10
#POINT_INTERNE_NB_EXTERNE				    4		    # defaut : 3
    
	resultats pas_de_sortie_finale_------------
COPIE 0
#  _fin_point_info_  


_suite_point_info_

    controle ------------
#---------------------------
# PARAMETRE    | VALEUR    |
#---------------------------
DELTAt            0.0001	# defaut : 1.
DELTAtMAXI        0.001		# defaut : 1.
TEMPSFIN          1.49		# defaut : 1.

  _fin_point_info_  
 
