207 lines
7.2 KiB
Text
207 lines
7.2 KiB
Text
|
|
|||
|
#######################################################################
|
|||
|
# #
|
|||
|
# | | ==== === ==== ==== | | | | #
|
|||
|
# | | | | | | / | | | | #
|
|||
|
# |====| |=== === |=== / |====| ------- ------- #
|
|||
|
# | | | | \ | / | | | | #
|
|||
|
# | | ==== | \ ==== ==== | | | | #
|
|||
|
# #
|
|||
|
#######################################################################
|
|||
|
# 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_
|
|||
|
|
|||
|
|