Anomalie #325
Problème loi de melange ORGEAS + MAHEO_HYPER et ISOHYPERBULK
100%
Description
Salut Gérard,
J'ai réalisé une simulation en traction uniaxiale sur un treillis avec une loi des mélanges en sigma qui utilise la loi "ISOHYPER3DORGEAS1" pour les poutres sollicitées en compression et une loi "ISOELAS1D" pour les poutres sollicitées en traction.
=> Pas de soucis particulier pour cette simulation
Ensuite, j'ai décidé de changer la loi "ISOELAS1D". Pour cela, j'ai réalisé une simulation sur une seule poutre en compression uniaxiale puis ensuite en traction uniaxiale avec une loi des mélanges en sigma qui utilise la loi "ISOHYPER3DORGEAS1" pour les poutres sollicitées en compression et une loi additive en sigma utilisant la loi "MAHEO_HYPER" et "ISOHYPERBULK3".
=> Pas de soucis particulier, la valeur ma première pente en compression et en traction sont proches.
Enfin, j'ai réutilisé le même treillis que précédemment (cas n°1) et j'ai appliquée ma loi des mélanges en sigma qui utilise la loi "ISOHYPER3DORGEAS1" pour les poutres sollicitées en compression et une loi additive en sigma utilisant la loi "MAHEO_HYPER" et "ISOHYPERBULK3".
=> Ma simulation ne converge pas.
Aurais-tu une idée afin de résoudre ce problème ?
Ci-joint (.zip) tu trouveras les fichiers .her, .info, .CVisu de mes 3 cas. Ainsi que les fichiers comportant mes différentes lois (.txt). Enfin, un fichier .gnu permet de tracer la contrainte en fonction de la déformation pour le cas 1 et 2.
A+
David
Fichiers
Mis à jour par Gérard Rio il y a environ 2 ans
- Statut changé de Nouveau à En cours
oui, je pense qu'il y a un pb dans le choix des paramètres de la loi Mahéo. Voici ce que j'obtiens en débug avec ma version actuelle (sans doute un choua différente de la tienne):
gerardrio@lg2m-2017.local% HZppD_Vn-1 -f Treillis-Orgeas_Maheo-Isohyperbulk
#######################################################################- #
- | | ==== === ==== ==== | | | | #
- | | | | | | / | | | | #
- |====| |=== === |=== / |====| ------- ------- #
- | | | | \ | / | | | | #
- | | ==== | \ ==== ==== | | | | #
- # #######################################################################
- Herezh++ is distributed under GPL 3 license ou ultérieure. #
- Copyright (C) 1997-2022 Université Bretagne Sud (France) #
- AUTHOR : Gérard Rio #
- E-MAIL : gerardrio56@free.fr #
- 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 7.004
-- initialisation de l'entree des donnees racine du nom de fichier principal = Treillis-Orgeas_Maheo-Isohyperbulk======== fin de la lecture, Projet::Lecture( ...===========
======================================================================
INCREMENT DE CHARGE : 1 intensite 0.01 t= 0.01 dt= 0.01
======================================================================
Erreur : (J_r(1)-3.)= -0.667009 est < 0. ! ce n'est pas possible avec les formules actuelles
Maheo_hyper::Calcul_dsigma_deps
passage dans la methode Sortieerreur non identifiee sur l'algo de la resolution de sig22_33(bsurb0,hsurh0)=(0,0)
LoiContraintesPlanesDouble::Calcul_DsigmaHH_tdt (...
Erreur : (J_r(1)-3.)= -0.667009 est < 0. ! ce n'est pas possible avec les formules actuelles
Maheo_hyper::Calcul_dsigma_deps- arret apres une erreur ** *******************************************************
======== erreur (Sortie) detectee ======temps_user:0/00:00:00.10 system:0/00:00:00.02 reel:0/00:00:00.11=============================================================
| fin HEREZH++ |
=============================================================Mis à jour par Gérard Rio il y a presque 2 ans
bonjour David,
est-ce que je peux fermer ce ticket ?
Mis à jour par Gérard Rio il y a presque 2 ans
- Statut changé de En cours à Résolu
- % réalisé changé de 0 à 100