Projet

Général

Profil

Wiki » Historique » Version 3

Gérard Rio, 22/09/2015 12:16

1 1 Gérard Rio
À partir de la version V 6.719:
2
  - Mise en place de la possibilité d'utiliser plusieurs matrices de raideurs et/ou plusieurs méthodes différentes de résolution, dans le cas d'un comportement transitoire non dynamique: algorithme: "non_dynamique" . 
3
Fonctionnement: Si la résolution du système linéaire avec la matrice normale échoue, la seconde matrice est alors utilisée. Si la seconde échoue, la troisième est utilisée, etc. La doc sera mise à jour dès que possible, mais en attendant voici un exemple de syntaxe qui permet de déclarer une matrice principale de type Lapack symétrique, et une matrice secondaire bande symétrique classique.
4
5
para_syteme_lineaire #----------------------------------------------
6
TYPE_MATRICE BANDE_SYMETRIQUE_LAPACK
7
MATRICE_S_SECONDAIRE_S_  BANDE_SYMETRIQUE   FIN_TYPE_MAT_SECONDAIRES_
8
TYPE_RESOLUTION_S_SECONDAIRE_S_  CHOLESKY FIN_TYPE_RESOLUTION_SECONDAIRES_
9 2 Gérard Rio
10
A voir dans les essais futurs !
11 3 Gérard Rio
12
---------------------------------
13
Version 6.720
14
- mise en place d'une loi d'hystérésis sur un comportement purement volumique. Dans le cas d'un matériau, il est donc nécessaire d'adjoindre à cette loi, une partie déviatorique. En fait la nouvelle loi : Hysteresis_Bulk est prévue pour représenter le comportement de matériaux poreux type mousses - matériaux architecturés. La doc sera mise à jour dès que possible, en attendant, il est possible d'avoir une information via la constitution interactive d'un fichier info et en choisissant comme loi Hysteresis_Bulk, un exemple de loi que j'ai utilisé pour la mise au point suit:
15
16
#------------------------------------------
17
# Nom Materiau    |      Type loi         |
18
#------------------------------------------
19
  hyste_bulk         HYSTERESIS_BULK
20
#-------------------------------------------------------------------
21
# .......  loi_de_comportement d'hysteresis spherique ........      |
22
# para de prager(>=0)  :         mu         : limite de plasticite |
23
#-------------------------------------------------------------------
24
   np= 2       mu= 8000     Qzero= 100 avec_parametres_de_reglage_
25
     type_de_resolution_  2 
26
     cas_kutta_  5  
27
#     erreurAbsolue_  1.e-3  erreurRelative_ 1.e-5
28
     erreurAbsolue_  1.e-8  erreurRelative_ 1.e-8
29
     nbMaxiAppel_  10000
30
     nb_iteration_maxi_  20  nb_dichotomie_maxi_  1  
31
     tolerance_residu_ 1.e-8 tolerance_residu_rel_ 1.e-8#5
32
     depassement_Q0_ 1000
33
     sortie_post_ 1
34
#     permet_affichage_ 8 
35
   fin_parametres_reglage_Hysteresis_
Redmine Appliance - Powered by TurnKey Linux