####################################################################### # # # | | ==== === ==== ==== | | | | # # | | | | | | / | | | | # # |====| |=== === |=== / |====| ------- ------- # # | | | | \ | / | | | | # # | | ==== | \ ==== ==== | | | | # # # ####################################################################### # Copyright (c) 1997-2015 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 6.686 ####################################################################### # fichier de commande # # --- écrasement du joint, avec contact entre le joint et la chemise fine #---------------------------------------------------- # definition de la dimension de l'espace de travail | #---------------------------------------------------- dimension 3 #--------------------------------------------------------------- # definition facultative du niveau d'impression (entre 0 et 10)| #--------------------------------------------------------------- niveau_commentaire 0 #------------------------------- # definition du type de calcul | #------------------------------- TYPE_DE_CALCUL non_dynamique #avec plus visualisation #-------------------------------------- #| definition du (ou des) maillage(s) | #-------------------------------------- # -- def maillage joint < Joint17_1x1.her def_mouvement_solide_initiaux_ mouvement_solide_ translation_= 15.015 7.6 0. centre_= 0. 0. 0. rotation_= 0. 0. 0. fin_mouvement_solide_ #----------------------------- # Definition des courbes #----------------------------- les_courbes_1D < courbe_ponderation.txt courbe_polynomiale COURBEPOLYNOMIALE # nom de la courbe puis le type de la courbe # def des coefficients d’un polynome du quatrieme degre a0+a1x+a2x^2+a3x^3+a4x^4 # debut_coef= a0 a1 a2 a3 a4 fin_coef debut_coef= -171.403600607624 1489.547106 -3391.128916 2999.230915 -926.245505244461 fin_coef #w1_4(x) # def du chargement de traction simple_char COURBEPOLYLINEAIRE_1_D Debut_des_coordonnees_des_points Coordonnee dim= 2 0.0 0.00000000 Coordonnee dim= 2 1 0.39 Coordonnee dim= 2 2 0.27 Coordonnee dim= 2 3 0.39 Coordonnee dim= 2 4 0.79 Coordonnee dim= 2 5 0.55 Coordonnee dim= 2 6 0.79 Coordonnee dim= 2 7 1.17 Coordonnee dim= 2 8 0.83 Coordonnee dim= 2 9 1.17 Coordonnee dim= 2 10 1.57 Coordonnee dim= 2 11 1.10 Coordonnee dim= 2 12 1.57 Coordonnee dim= 2 13 1.98 Coordonnee dim= 2 14 1.38 Coordonnee dim= 2 15 1.98 Fin_des_coordonnees_des_points les_Fonctions_nD #------------ # fonction Fonction_expression_litterale_nD # f(P,PR,PRM1) # P : la pression, # PR : la dernière pression de référence # PRM1: la pression de référence avant PR fct_dep_pression FONCTION_EXPRESSION_LITTERALE_nD un_argument= P un_argument= PR un_argument= PRM1 fct= (P > 0.) ? \ # d'abord le cas positif ( ( (PR != 0.) && (P > PR) && (P <= (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ ( 0.4/(PR- (PR+PRM1)/2.) * P + (1. - PR*0.4/(PR- (PR+PRM1)/2.))) \ # première moitiè de la remontée : ( ( (PR != 0.) && (P > (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ ( 0.4/(PRM1 - (PR+PRM1)/2.) * P + (1. - 0.4 * PRM1/(PRM1 - (PR+PRM1)/2.))) \ #deuxième moitié de remontée : 1. \ # tous les autres cas ) \ ) \ : ( ( (PR != 0.) && (P > PR) && (P <= (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ # sinon cas négatif ( -0.4/(PR- (PR+PRM1)/2.) * P + (1. + PR*0.4/(PR- (PR+PRM1)/2.))) \ # première moitiè de la remontée : ( ( (PR != 0.) && (P > (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ ( -0.4/(PRM1 - (PR+PRM1)/2.) * P + (1. + 0.4 * PRM1/(PRM1 - (PR+PRM1)/2.))) \ #deuxième moitié de remontée : 1. \ # tous les autres cas ) \ ) fin_parametres_fonction_expression_litterale_ #====================================== #| definition des lois de comportement| #-------------------------------------- choix_materiaux ---------- #---------------------------------------- # Elements | Nom Materiau | #---------------------------------------- nom_mail= Joint E_tout elastomere materiaux #---------- #------------------------------------------ # Nom Materiau | Type loi | #------------------------------------------ acier ISOELAS # ....... loi de comportement isoelastique 3D ........ # module d'young : coefficient de poisson 210000 0.3 < elastomere_joint.loi #--------------------------------------- fin def des lois de comportement ------------------------ # --- divers stockages (1) ------- masse_volumique #--------# nom_mail= Joint E_tout 1 charges #------------# blocages #------------# #----------------------------------------------------------- # nom du maillage | Ref noeud | Bloquages #----------------------------------------------------------- # le joint nom_mail= Joint N_S UY nom_mail= Joint N_O UX nom_mail= Joint N_N 'UY = COURBE_CHARGE: simple_char ECHELLE: -1', UZ nom_mail= Joint N_tout UZ nom_mail= Joint N_E UX controle #------------ #--------------------------- # PARAMETRE | VALEUR | #--------------------------- SAUVEGARDE 0 DELTAtMAXI 0.05 DELTAtMINI 1.e-8 TEMPSFIN 10 DELTAt 0.01 MAXINCRE 1000000 ITERATIONS 16 PRECISION 5.e-3 para_pilotage_equi_global ------------ #----------------------------------------------- # PARAMETRE | VALEUR | #----------------------------------------------- FACTEUR_DIMINUTION 1.732 #1.732 FACTEUR_AUGMENTATION 1.414 #1.414 NB_BONNE_CONVERGENCE 3 NORME_MAXI_INCREMENT 0.5#0.1 para_syteme_lineaire ------------ #----------------------------------------------- # PARAMETRE | VALEUR | #----------------------------------------------- TYPE_MATRICE BANDE_NON_SYMETRIQUE_LAPACK SYMETRIE_MATRICE 0 para_affichage #------------ #--------------------------- # PARAMETRE | VALEUR | #--------------------------- FREQUENCE_AFFICHAGE_ITERATION 20 FREQUENCE_AFFICHAGE_INCREMENT 20 FREQUENCE_SORTIE_FIL_DU_CALCUL DERNIER_CALCUL # ------------------------------------------------------------------------------------ resultats pas_de_sortie_finale_ COPIE 0 # _fin_point_info_