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

#-------------------------------
# definition du type de calcul |
#-------------------------------
TYPE_DE_CALCUL
##non_dynamique avec plus visualisation
non_dynamique avec plus uniquement_remontee #plus visualisation
PARA_TYPE_DE_CALCUL
#mode_debug_= 1

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

 #-- def maillage 
< Joint_4x4_quadracompl.her
< Joint_4x4_quadracompl.lis
#< Chemise_quadracompl.her
#< Chemise_quadracompl.lis
< Chemise.her
< Chemise.lis
#< Tige.her
#< Tige.lis

domaine_esclave
1

les_courbes_1D ------
< courbe_1D.hz_courbe

les_fonctions_nD #------------
< fct_nD_ponderation_charge.hz_courbe
< fct_nD_ponderation_decharge.hz_courbe

# Choix des pondération de l'hyperélastique sphérique selon charge ou décharge
hyperelastique_spherique FONC_SCAL_COMBINEES_ND
fct_base= courbe_bulk_charge_2
fct_base= courbe_bulk_decharge
fin_fcts_interne_fonction_combinee_
un_argument= temps_courant
#permet_affichage_  4 
fct= (temps_courant > 1) ? courbe_bulk_decharge : courbe_bulk_charge_2
fin_parametres_fonction_combinee_

# Choix des pondération de l'hyperélastique déviatorique selon charge ou décharge
hyperelastique_deviatorique FONC_SCAL_COMBINEES_ND
fct_base= courbe_orgeas_charge_1
fct_base= courbe_orgeas_decharge
fin_fcts_interne_fonction_combinee_
un_argument= temps_courant
#permet_affichage_  4 
fct= (temps_courant > 1) ? courbe_orgeas_decharge : courbe_orgeas_charge_1
fin_parametres_fonction_combinee_

# Choix des pondération de l'hytérésis sphérique selon charge ou décharge
hysteresis_spherique FONC_SCAL_COMBINEES_ND
fct_base= courbe_hysteresis_bulk_charge_1
fct_base= courbe_hysteresis_bulk_decharge
fin_fcts_interne_fonction_combinee_
un_argument= temps_courant
#permet_affichage_  4 
fct= (temps_courant > 1) ? courbe_hysteresis_bulk_decharge : courbe_hysteresis_bulk_charge_1
fin_parametres_fonction_combinee_

# Choix des pondération de l'hytérésis déviatorique selon charge ou décharge
hysteresis_deviatorique FONC_SCAL_COMBINEES_ND
fct_base= courbe_hysteresis_3D_charge_1
fct_base= courbe_hysteresis_3D_decharge
fin_fcts_interne_fonction_combinee_
un_argument= temps_courant
#permet_affichage_  4 
fct= (temps_courant > 1) ? courbe_hysteresis_3D_decharge : courbe_hysteresis_3D_charge_1
fin_parametres_fonction_combinee_

#======================================
#| definition des lois de comportement|
#--------------------------------------
 
            choix_materiaux ----------
#---------------------------------------
# nom de maillage |Elements | Materiau |
#---------------------------------------
 nom_mail= Joint_4x4_quadracompl	E_tout	 graphite
#nom_mail= Chemise_quadracompl	E_tout	acier
 nom_mail= Chemise			E_tout	plastique_mou #acier
#nom_mail= Chemise		E_tout	loi_frottement
#nom_mail= Tige			E_tout	acier
#nom_mail= Tige			E_tout	loi_frottement

		materiaux #---------

#< graphite.hz_loi	

#####
graphite LOI_ADDITIVE_EN_SIGMA


avec_fonction_de_ponderation_ tangent_ddl_via_eps= 1

#---------------------------------------
# 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_=   hyperelastique_spherique    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_=   hyperelastique_deviatorique   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_=   hysteresis_spherique   fin_fonct_

HYSTERESIS_BULK #hysteresis spherique
#np= 1.75	mu= 850	Qzero= 16.0
np= 0.85       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-3#5
erreurRelative_ 1.e-4#7
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_=   hysteresis_deviatorique   fin_fonct_

HYSTERESIS_3D
#np= 1.0       mu= 2800    Qzero= 1.75
np= 0.85       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-3#5
erreurRelative_ 1.e-4#7
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
depassement_Q0_ 2.2
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
  
#------------------------------------------
# Nom Materiau    |      Type loi         |
#------------------------------------------
  plastique_mou         ISOELAS
#------------------------------------------
# E    |     Nu
#------------------------------------------
  190     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_4x4_quadracompl 	E_tout	1. 
#nom_mail= Chemise_quadracompl		E_tout	1.
nom_mail= Chemise				E_tout	1.
#nom_mail= Tige				E_tout	1.

  zone_contact
#----------------------------------------------------
# ref noeuds | ref face |
#----------------------------------------------------
#nom_mail= Joint_4x4_quadracompl N_ext		nom_mail= Chemise_quadracompl	A_contact
nom_mail= Joint_4x4_quadracompl N_ext		nom_mail= Chemise			A_contact
#nom_mail= Joint_4x4_quadracompl N_int		nom_mail= Tige			A_contact

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

#nom_mail= Joint_4x4_quadracompl	N_ext		UX # sans contact
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= Chemise			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.94'
  


             controle ------------
#---------------------------
# PARAMETRE    | VALEUR    |
#---------------------------
SAUVEGARDE	1			# defaut : 1
DELTAt	0.002#459560 #0.001 #0.0001		# defaut : 1.
DELTAtMAXI	0.01#6 #0.002			# defaut : 1.
DELTAtMINI	1.e-7			# defaut : 0.
TEMPSFIN	0.9 #0.59810066 #1.49 #2.00			# defaut : 1.
PRECISION	1.e-3			# defaut : 1.e-3
#NORME		Residu/Reaction	# defaut : Residu/Reaction
NORME		Residu/Reaction_et_verification_ddouble
ITERATIONS	20			# defaut : 100
MAXINCRE	999999		# defaut : 400
#RESTART  25

	para_syteme_lineaire ------------
TYPE_MATRICE	BANDE_NON_SYMETRIQUE_LAPACK	# defaut : BANDE_SYMETRIQUE
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.01#0.001		# defaut : +infini
CAS_JACOBIEN_NEGATIF      		2		# defaut : 1

	para_contact ------------
CONTACT_TYPE				0#2# 0 # 4		# defaut : 0
PENALISATION_PENETRATION		0.1		# defaut : 0.1
TYPE_PENALISATION_PENETRATION 	1#8		# defaut : 2
PENALISATION_TANGENTIELLE		1.e6		# defaut : 1.e6
TYPE_PENALISATION_TANGENTIELLE	8		# defaut : 2
TYPE_DE_DECOLLEMENT     		1		# defaut : 0
NB_DECOLLEMENT_MAXI			10		# defaut : 1
PRECISION_BOITE_PRELOCALISATION	1.05		# defaut : 1.05
MINI_EXTRA_BOITE_PRELOCALISATION	0.01	# defaut : 0.001
AJOUT_EXTRA_BOITE_PRELOCALISATION	1.0		# defaut : 0.
PENETRATION_CONTACT_MAXI 		0.00001	# defaut : 0.1
TANGENTIELLE_CONTACT_MAXI		1000.1	# defaut : 0.1
PENETRATION_BORNE_REGULARISATION 	0.01		# 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

	para_affichage ------------
FREQUENCE_SORTIE_FIL_DU_CALCUL        1	

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

 _fin_point_info_  
 ##_suite_point_info_

	controle ------------
	
#---------------------------
# PARAMETRE    | VALEUR    |
#---------------------------
DELTAt            0.001#0001	# defaut : 1.
DELTAtMAXI        0.06#002		# defaut : 1.
TEMPSFIN          1.49		# defaut : 1.

  _fin_point_info_  
 
  
