
dimension 3
niveau_commentaire 3

TYPE_DE_CALCUL
dynamique_relaxation_dynam avec plus sauveMaillagesEnCours
PARA_TYPE_DE_CALCUL

typeCalRelaxation= 1 lambda= 1 type_calcul_mass= 2 option_recalcul_mass= 0
parametre_calcul_de_la_masse_ casMass_relax= 3
lambdaAvecPonderation_Globale_ f_lambda
parametre_activation_du_contact_ type_activation_contact= 1 # val. par dÃ©faut, voire P. 95)

avec_amortissement_cinetique_
   max_nb_decroit_pourRelaxDyn_   1
   taille_moyenne_glissante_      1
   coef_arret_pourRelaxDyn_       0.
   coef_redemarrage_pourRelaxDyn_ 0.
   max_deltaX_pourRelaxDyn_       0.1
   nb_deb_test_amort_cinetique_  1
   nb_max_dX_OK_pourRelaxDyn_     1 # 10
   nb_deb_testfin_pourRelaxDyn_   1 # 250
fi_parametre_amortissement_cinetique_

ARRET_A_EQUILIBRE_STATIQUE_ 2
#mode_debug_= 1
#------------------------------------------------------------------ 

< maillage.her
#renumerotation_des_noeuds_

# def de mouvements solides
def_mouvement_solide_initiaux_
mouvement_solide_ 
 translation_= 0.0850 0.0000 0.0000
 centre_=      0.0850 0.0000 0.0000
 rotation_=     0 -1.56836999  0 en_radian_
fin_mouvement_solide_


nom_maillage  plan1 # plan de symÃ©trie1 (y<0)
< plan_sym.her
def_mouvement_solide_initiaux_
mouvement_solide_ 
 centre_=      0.  0.  0.
 translation_= -1  -1 -1. # depasse au centre et en bas du ballon
 rotation_=    0.  0. 0.042454 en_gradient_
 translation_= 0 -10.00 0.
fin_mouvement_solide_

nom_maillage  plan2 # plan de symÃ©trie2 (y>0)
< plan_sym.her
def_mouvement_solide_initiaux_
mouvement_solide_ 
 centre_=       0.   0.  0.
 translation_= -1   0  -1.  # depasse en bas du ballon et Ã  l'intÃ©rieur
 rotation_=     0.   0   -0.042454 en_gradient_
 translation_= -0.4244 9.990990 0.
fin_mouvement_solide_

domaine_esclave
1 #maillage 1 defini comme esclave
#------------------------------------------------------------------ 

    les_courbes_1D #------------
# Fonction de pondération pour le loi_melange
f_ponderation COURBEPOLYLINEAIRE_1_D
  Debut_des_coordonnees_des_points
    Coordonnee dim= 2  0.0  1
    Coordonnee dim= 2  0.5  1
    Coordonnee dim= 2  1.0  0
    Coordonnee dim= 2  2.0  0
    Coordonnee dim= 2  3.0  0
  Fin_des_coordonnees_des_points  
  
# Fonction déplacement noeuds pôles
f_dep COURBEPOLYLINEAIRE_1_D
  Debut_des_coordonnees_des_points
    Coordonnee dim= 2  0  0
    Coordonnee dim= 2  1  0
    Coordonnee dim= 2  2  1
    Coordonnee dim= 2  3  1
  Fin_des_coordonnees_des_points

# Fonction de déplacement des surfaces de contact  
f_depPlans COURBEPOLYLINEAIRE_1_D
  Debut_des_coordonnees_des_points
    Coordonnee dim= 2  0.   0
    Coordonnee dim= 2  0.5  0
    Coordonnee dim= 2  1    1#0#1
    Coordonnee dim= 2  2    1#0#1
  Fin_des_coordonnees_des_points  

  # Fonctions complémentaires pour les lois
# Pas de fonction complémentaire

  
    les_fonctions_nD #------------
# Fonction pour le pilotage du lambda
f_lambda FONCTION_EXPRESSION_LITTERALE_nD
  un_argument= temps_courant  
  fct= ((temps_courant <= 1.) ?  0.554 :  5) #2.000  )
  fin_parametres_fonction_expression_litterale_  

# Fonction pour le pilotage du critere    
f_critere FONCTION_EXPRESSION_LITTERALE_nD
  un_argument= temps_courant  
  fct= ((temps_courant < 1) ?  -1e9 :  0  )
  fin_parametres_fonction_expression_litterale_ 
  
  # calcul du volume total
le_volume FONCTION_EXPRESSION_LITTERALE_nD 
   deb_list_var_   X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_   
   fct= -(X1*N_surf_1 + X2*N_surf_2 + X3*N_surf_3)/3 * 76 # A adapter
   fin_parametres_fonction_expression_litterale_  

# Pilotage du contact de type 4   
f_ct4 FONCTION_EXPRESSION_LITTERALE_nD
un_argument= compteur_iteration_algo_global
#un_argument= compteur_increment_charge_algo_global
#fct= rint(compteur_increment_charge_algo_global/100) ? 1 : 0 , 1
fct= rint(compteur_iteration_algo_global/20) ? 1 : 0 , 4
#fct= (compteur_iteration_algo_global > 100 )? 1 : 0 , 4
fin_parametres_fonction_expression_litterale_
#------------------------------------------------------------------ 

	choix_materiaux #---------
    
nom_mail= ballon E_quad	loi_film
nom_mail= ballon E_biel_bord1  loi_sangleBC
nom_mail= ballon E_biel_bord2  loi_sangleBC
nom_mail= ballon E_biel_centre loi_rien1D

nom_mail= plan1 E_to pour_solide_en_contact #rien3D 
nom_mail= plan2 E_to pour_solide_en_contact #rien3D 

#------------------------------------------------------------------ 

    materiaux #---------
    
#-- Loi rien 1D    
loi_rien1D LOI_RIEN1D    
rien3D LOI_RIEN3D
    
#-- Loi sangle bord courbe
loi_sangleBC ISOELAS1D
 5.923e+09  0.3

#-- Loi film
loi_film  LOI_DES_MELANGES_EN_SIGMA
  les_grandeurs_de_controle_= TEMPS  fin_grandeurs_
  deb_fonct_= f_ponderation fin_fonct_
  
    #-- première loi
    ISOELAS2D_C demarre_a_prop_non_nulle_puis_strictement_decroissante_
        1.883e+09  0.30
    
   #-- seconde loi
    LOI_CRITERE demarre_a_prop_non_nulle_puis_strictement_decroissante_
    TYPE_DE_CRITERE_ PLISSEMENT_MEMBRANE avecNiveauSigmaI_mini_pour_plis_
        les_grandeurs_de_controle_= temps_courant  fin_grandeurs_
        deb_fonct_= f_critere fin_fonct_ 
        
        LOI_CONTRAINTES_PLANES 
            NEWTON_LOCAL  avec_parametres_de_reglage_
                nb_iteration_maxi_  20
                nb_dichotomie_maxi_ 20
                tolerance_residu_ 1e1  # tolerance residu 10 Pa
                tolerance_residu_rel_ 1.e-4
                maxi_hsurh0_  4
                mini_hsurh0_  0.25
                permet_affichage_ 0
            fin_parametres_reglage_Algo_Newton_

        # Loi IsoElastique3D
        ISOELAS
            1.8830e+09  3.0000e-01

        fin_loi_contrainte_plane

    fin_loi_interne

    avec_parametres_controle_pour_loi_critere_
        choix_methode_cal_plis_memb_ 2
    fin_parametres_controle_pour_loi_critere_
	
	fin_loi_critere 
fin_liste_lois_elementaires
#-------------------------------------------------------- 

pour_solide_en_contact ISOELAS
    1.8830e+1  3.0000e-01


	masse_volumique #---------
nom_mail= ballon E_quad  9.130000e+02 # en kg/m^3
nom_mail= ballon E_biel  8.680000e+02 # en kg/m^3

nom_mail= plan1 E_to 1.
nom_mail= plan2 E_to 1.

    epaisseurs #------------
nom_mail= ballon E_quad_film    2.500000e-05


	sections  #---------
nom_mail= ballon E_biel_bord    4.793000e-06
nom_mail= ballon E_biel_centre  1.000000e-06 

    integrale_sur_volume_
nom_mail= ballon E_quad une_fonction_nD_ le_volume

#    zone_contact -------------
    glue_contact 
nom_mail= ballon N_41 nom_mail= plan1 F_avant
nom_mail= ballon N_42 nom_mail= plan2 F_arriere
#------------------------------------------------------------------ 

    charges #---------
 # gradient de pression
nom_mail= ballon F_tout  PHYDRO  0. 0. 1.  0. 0. 0.  0.10395698 ATTRIBUT_  sans_limitation_ 

 # poids propre du a la gravite
nom_mail= ballon E_quad VOLUMIQUE  0. 0. -8.862618e+03
nom_mail= ballon E_biel_bord1  VOLUMIQUE 0. 0. -8425.796645
nom_mail= ballon E_biel_bord2  VOLUMIQUE 0. 0. -8425.796645
nom_mail= ballon E_biel_centre VOLUMIQUE 0. 0. 0.000000


 # poids du pole
nom_mail= ballon N_poleN_ext PONCTUELLE 0. 0. -0.603417
#------------------------------------------------------------------ 

    blocages #---------
    
nom_mail= plan1 N_to mouvement_solide_ 
translation_= 0 10.00 0.
fin_mouvement_solide_ COURBE_CHARGE: f_depPlans 
    
nom_mail= plan2 N_to mouvement_solide_ 
translation_= 0.4244 -9.990990 0.
fin_mouvement_solide_ COURBE_CHARGE: f_depPlans        
    
nom_mail= ballon N_bord1_horsPS  UY TEMPS_MAXI= 1
nom_mail= ballon N_bord2_horsPS  UY TEMPS_MAXI= 1

nom_mail= ballon N_poleS_1 UX, 'UY= COURBE_CHARGE: f_dep ECHELLE:  0.1463892443' UZ
nom_mail= ballon N_poleS_2 UX, 'UY= COURBE_CHARGE: f_dep ECHELLE:  -0.1463892443' UZ


# Encastrement t=3 pour force de réaction
nom_mail= ballon N_tout UX UY UZ TEMPS_MINI= 2

# Encastrement des plans
# nom_mail= plan1 N_tout UX UY UZ
# nom_mail= plan2 N_tout UX UY UZ

    condition_limite_lineaire_  #--------    
# plan de symétrie cyclique 
nom_mail= ballon N_bord2_horsPS enu= X1
TEMPS_MINI= 1 TEMPS_MAXI= 2
def_auto_coef_planOuDroite_ centre_fixe_ 0. 0. 0. \
coefficients= -0.0424412031961483 0.9990989662046815  0 fin_list_coefficients_

nom_mail= ballon N_bord1_horsPS enu= X1
TEMPS_MINI= 1  TEMPS_MAXI= 2
def_auto_coef_planOuDroite_ centre_fixe_ 0. 0. 0. \
coefficients= -0.0424412032 -0.9990989662  0 fin_list_coefficients_


# blocage en x noeuds poleN_bord
nom_mail= ballon  N_poleN_bord enu= UX
TEMPS_MAXI= 2
val_condi_lineaire= 0 coefficients= 1 0 0 fin_list_coefficients_

#-- ajout de masse additionnelle sur les noeuds bloqués
# permet de ne pas avoir de message d'erreur de masse nulle
# mais cela ne change pas le calcul car les noeuds sont totalement bloqués
masse_addi
 nom_mail= plan1 N_tout 1.
 nom_mail= plan2 N_tout 1.

#------------------------------------------------------------------ 
    typecharge #-----------
TYPE3

    controle #------------ 
TEMPSFIN    1
DELTAt 	    0.5
ITERATIONS 100000
MAXINCRE    10
NORME E_cinetique/E_statique_ET_Res/Reac_et_Fext
PRECISION 1.00e-03
SAUVEGARDE  1
RESTART 0

    para_contact ----------
CONTACT_TYPE                   0#2#4 # 2

TYPE_PENALISATION_PENETRATION   8 # 1
PENALISATION_PENETRATION        0.01
PENETRATION_CONTACT_MAXI        0.1 #0.5
FCT_ND_BASCUL_CONTACT_TYPE_4    f_ct4

 NB_DECOLLEMENT_MAXI              10000
 DISTANCE_MAXI_AU_PT_PROJETE      20. #1.
# TYPE_DE_DECOLLEMENT              1
# PENETRATION_BORNE_REGULARISATION 0.001
# FORCE_CONTACT_NOEUD_MAXI         100000

PENALISATION_TANGENTIELLE          0.
TYPE_PENALISATION_TANGENTIELLE     8

NIVEAU_COMMENTAIRE_CONTACT         3#3#7
MINI_EXTRA_BOITE_PRELOCALISATION   0.5
AJOUT_EXTRA_BOITE_PRELOCALISATION  0.1


    para_affichage #------------ 
FREQUENCE_AFFICHAGE_INCREMENT  1
FREQUENCE_AFFICHAGE_ITERATION  100
FREQUENCE_SORTIE_FIL_DU_CALCUL 1

resultats  pas_de_sortie_finale_ 
COPIE 0   

# _fin_point_info_  
 
 _suite_point_info_

    controle #------------ 
TEMPSFIN    2
DELTAt      1 
DELTAtMINI  1 
DELTAtMAXI  1 
PRECISION   5.00e-03   
RESTART 0
    para_contact --------------
CONTACT_TYPE   2

_suite_point_info_

    para_contact --------------
CONTACT_TYPE   0

    controle #------------ 
TEMPSFIN    3
RESTART 0

  _fin_point_info_  