dimension 3

niveau_commentaire 1#3

TYPE_DE_CALCUL
non_dynamique  #avec plus visualisation
PARA_TYPE_DE_CALCUL
# .........................................................
# / type dĠalgorithme /
#.........................................................
#mode_debug_=  1

< mail_gmsh/joint.her
< mail_gmsh/joint.lis

< mail_gmsh/logement_creux.her
< mail_gmsh/logement_creux.lis

< mail_gmsh/plateau_inferieur.her
< mail_gmsh/plateau_inferieur.lis

< maillage_plateau.her



domaine_esclave
1



les_courbes_1D
  courbe_UY COURBEPOLYLINEAIRE_1_D
    Debut_des_coordonnees_des_points
      Coordonnee dim= 2 0. 0.
      Coordonnee dim= 2 1. 1.
    Fin_des_coordonnees_des_points

les_fonctions_nD

f_com_contact FONCTION_EXPRESSION_LITTERALE_nD

un_argument= NUM_NOEUD un_argument= NUM_MAIL_NOEUD
   fct= ((NUM_NOEUD == 32)&&(NUM_MAIL_NOEUD == 1)) \
           ? 7 : 0 
  fin_parametres_fonction_expression_litterale_


choix_materiaux
nom_mail= joint  E_tout MAT_joint_HYPO
#nom_mail= joint  E_tout loi_frottement
nom_mail= logement_creux  E_tout MAT_outil
nom_mail= plateau_inferieur  E_tout MAT_outil
nom_mail= plateau   E_tout MAT_outil

materiaux
MAT_joint_ISOELAS ISOELAS
  500. 0.45

MAT_joint_HYPO HYPO_ELAS3D
  #conversion : E=500 NU=0.45
  Kc= 5000   mu= 344.827586206897    type_derivee -1
  fin_loi_HYPO_ELAS3D

MAT_outil ISOELAS
  190000. 0.3

loi_frottement LOI_COULOMB
  mu_statique= 0.1
fin_loi_frottement_coulomb_


masse_volumique
nom_mail= joint  E_tout   1.
nom_mail= logement_creux  E_tout   1.
nom_mail= plateau_inferieur  E_tout   1.
nom_mail= plateau  E_tout   1.

zone_contact avec_lissage_de_la_normale
#nom_mail= joint_1 N_contact nom_mail= logement A_contact
#nom_mail= joint_1 N_contact nom_mail= plateau A_contact
nom_mail= joint N_contact nom_mail= logement_creux A_interne
nom_mail= joint N_contact nom_mail= logement_creux A_externe

zone_contact 
nom_mail= joint N_contact nom_mail= plateau A_contact
nom_mail= joint N_contact nom_mail= plateau_inferieur A_contact


charges


blocages

#nom_mail= joint    N_tout UX


nom_mail= logement_creux N_tout UZ
nom_mail= plateau  N_tout UZ
nom_mail= plateau_inferieur  N_tout UZ
nom_mail= joint    N_tout UZ

#nom_mail= joint    N_base UX UY

nom_mail= logement_creux N_LIGNE_axeRevol UX
nom_mail= logement_creux N_LIGNE_encastr_bas UX UY
nom_mail= logement_creux N_LIGNE_encastr_haut UX UY

nom_mail= plateau_inferieur N_LIGNE_axeRevol UX
nom_mail= plateau_inferieur N_LIGNE_encastr_bas UX UY


nom_mail= plateau N_tout UX

nom_mail= plateau N_tout 'UY= COURBE_CHARGE: courbe_UY ECHELLE: -10.'



controle

DELTAtMINI 1.e-6
DELTAtMAXI 0.05
TEMPSFIN 1.
DELTAt 0.01 #0.00565685 #0.001
ITERATIONS 50#15
PRECISION 7.e-3
SAUVEGARDE   1
MAXINCRE 9999999
NORME Residu/Reaction
#RESTART    83
 NORME  min(Res,Res/MaX(Reaction_et_PVExterne))


para_contact
CONTACT_TYPE 2#2
PENALISATION_PENETRATION    1000.
TYPE_PENALISATION_PENETRATION   4
NB_DECOLLEMENT_MAXI         4#1
DISTANCE_MAXI_AU_PT_PROJETE   10.#1.
TYPE_DE_DECOLLEMENT     0
PENETRATION_BORNE_REGULARISATION 1e-4#3
PENETRATION_CONTACT_MAXI 1.e-3#4
OPTIMISATION_NUMEROTATION 1#1
FORCE_CONTACT_NOEUD_MAXI 1e7#7#5#7

PENALISATION_TANGENTIELLE       1.e2
TYPE_PENALISATION_TANGENTIELLE  2
TANGENTIELLE_CONTACT_MAXI        0.1
TANGENTIELLE_BORNE_REGULARISATION 0.01


NIVEAU_COMMENTAIRE_CONTACT  2#FCT_ND_NIVEAU_COMMENTAIRE_CONTACT f_com_contact #0#5#7#4
NIVEAU_COMMENTAIRE_LESCONTACTS  2#5#0#5

para_pilotage_equi_global
#TYPE_DE_PILOTAGE AUCUN_PILOTAGE
#NB_BONNE_CONVERGENCE 1
FACT_DIM_EN_MAUVAISE_CONVERGENCE  1.5 # diminution of the increment if divergence
NB_ITER_POUR_MAUVAISE_CONVERGENCE 20   # number of iteration after what we considere that the convergence is no good
NB_ITER_POUR_BONNE_CONVERGENCE  10    # under this value, the convergence is considered good
NORME_MAXI_INCREMENT  0.2#5 # --- maximum of the norm of the delta increment for a Newton iteration  
#MAXI_VARIATION_DDL_POUR_CONVERGENCE 2.   # if the maximum of variation of dof is up this value => divergence
###MINI_VARIATION_DDL_POUR_CONVERGENCE 1.e-8 # if the maximum of variation of the dof for an iteration is under this value => divergence
#NB_CYCLE_CONTROLE_RESIDU  5   #  if after 4 iterations, the residual balance at node don t decrease, => diversgence 
#PLAGE_CONTROLE_RESIDU     1   # si au bout de 4 iterations, on a toujours pas de diminution du résidue => divergence
#CAS_JACOBIEN_NEGATIF      2  # if a negatif jacobien is detected at a gauss point  => divergence
###VAR_MAXI_JACOBIEN           1.5  # up to 50% of volume  variation, we consider that a divergence appears
INIT_INCRE_AVEC_PAS_PREC  0 # initialisation a the new step with the last direction of displacement

            para_calculs_geometriques #------------
#---------------------------
# PARAMETRE    |  VALEUR   |
#---------------------------
POINT_INTERNE_DELTA_THETAI_MAXI    1.e-4#6#4
POINT_INTERNE_PREC_THETAI_INTERNE  1.e-6#2#6
POINT_INTERNE_NB_BOUCLE_SUR_DELTA_THETAI  16#0
POINT_INTERNE_NB_EXTERNE                  4


para_syteme_lineaire
#OPTIMISATION_POINTEURS_ASSEMBLAGE 1
#TYPE_MATRICE BANDE_SYMETRIQUE_LAPACK
#MATRICE_S_SECONDAIRE_S_   BANDE_SYMETRIQUE   FIN_TYPE_MAT_SECONDAIRES_
#TYPE_RESOLUTION_S_SECONDAIRE_S_  CHOLESKY  FIN_TYPE_RESOLUTION_SECONDAIRES_

para_affichage
FREQUENCE_SORTIE_FIL_DU_CALCUL  1


resultats pas_de_sortie_finale_
COPIE 0


_fin_point_info_
