206 lines
7.2 KiB
Text
206 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# 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 20000
|
||
FREQUENCE_AFFICHAGE_INCREMENT 20000
|
||
FREQUENCE_SORTIE_FIL_DU_CALCUL DERNIER_CALCUL
|
||
|
||
# ------------------------------------------------------------------------------------
|
||
resultats pas_de_sortie_finale_
|
||
COPIE 0
|
||
#
|
||
|
||
|
||
_fin_point_info_
|
||
|
||
|