Anomalie #207
masses nulle en relaxation dynamique et loi DEFORMATION PLANE
Description
Bonjour Gérard,
j'ai testé un calcul de relaxation dynamique sur une membrane 2D avec une loi LOI_DEFORMATIONS_PLANES. J'obtiens le message suivant :
--------------------------------
matrice masse principale:
stockage initial matrice masse --> diagonal, nbddl= 363
matrice masse secondaire:
stockage initial matrice masse --> diagonal, nbddl= 363
*** attention : le noeud 1 a une masse nulle en coordonnee 1
*** attention : le noeud 1 a une masse nulle en coordonnee 2
*** attention : le noeud 1 a une masse nulle en coordonnee 3
*** attention : le noeud 2 a une masse nulle en coordonnee 1
*** attention : le noeud 2 a une masse nulle en coordonnee 2
*** attention : le noeud 2 a une masse nulle en coordonnee 3
etc... (idem pour tous les noeuds)
--------------------------------
Ce ne sont que des warning donc le calcul tourne quand même. Mais il produit une erreur de loi de comportement sur la loi d'hysteresis dès la première itération de chaque incrément. Ce qui montre bien que les masses sont réellement nulles.
Je t'ai joint une archive du calcul (un quart de coussin carré sous pression uniforme avec une loi additive encapsulée dans une option LOI_DEFORMATIONS_PLANES).
Comme roue de secours, je me suis dit que ça pourrait être pas mal de pouvoir faire tourner la relaxation dynamique avec une matrice masse physique. L'utilisateur pourrait quand même faire tourner un calcul en attendant la correction de tel ou tel type d'erreur liée au calcul de masse via la raideur.
Je ne sais pas si ça existe déjà. J'ai fouillé la doc mais je n'ai pas trouvé.
Fichiers
Mis à jour par Gérard Rio il y a environ 5 ans
- Fichier HH_20C_3D_DP_test.hz_loi HH_20C_3D_DP_test.hz_loi ajouté
- Fichier coussin_HH_3D_DP_test.info coussin_HH_3D_DP_test.info ajouté
- Statut changé de Nouveau à Résolu
- % réalisé changé de 0 à 100
bonjour Julien,
1) le pb de fonctionnement vient de la forme du potentiel utilisé pour la loi ISOHYPERBULK_GENE
Dans la définition du potentiel on utilise une puissance non entière de (V-1) or cette grandeur peut-être négative ce qui entraine de NAN
J'ai modifié le potentiel (cf. fichiers joins) et cela fonctionne (mais c'est pour un potentiel qui ne répond sans doute pas à ton objectif, donc à voir: cf. remarque dans la doc qui rappel qu'il faut que le potentiel soit toujours calculable notamment pour pouvoir calculer les dérivées première et seconde).
2) concernant les masses nulles (ce qui n'est finalement pas le cas ici), on peut imposer de remplacer les masses nulles par une valeur particulière: cf. doc: "8.3.8 Traitement du cas particulier de masse nulle"). Mais par défaut aucun traitement n'est fait de manière à ce que le calcul diverge et met en évidence un pb éventuel (par exemple de noeud libre ...)