Assistance #6
Forgeage : 3eme année mécatro. Ajout de la partie visqueuse
Description
Bonjour
J'ai un souci avec le fichier .loi du TP sur le forgeage. En effet, sans y toucher le calcul Herezh passe parfaitement mais lorsque je décommente la partie visqueuse (Newton 3D), le calcul ne passe plus et me met le message suivant :
"incoherence, il devrait rester un pointeur de valide sur la liste de pt de ref
Hysteresis3D::Gestion_pointeur_coïncidence(...
non convergence sur une loi de comportement "
Merci d'avance
Mathieu LE NEL & Erwan KERGOURLAY
Fichiers
Mis à jour par Gérard Rio il y a presque 12 ans
Bonjour,
effectivement le fait d'ajouter une viscosité non linéaire peut poser problème. En fait cela provient de ce qu'à l'origine des vitesses (c-a-d pour des vitesses nulles) la dérivée de la contrainte est infinie. Donc cela bloque le système si jamais on atteind des vitesses nulles ou très faible (ce qui est le cas semble-t-il ici à certains moments). Classiquement, on ajoute une vitesse faible mais non nulle pour l'origine de manière à éviter ce problème, mais cela n'est pas présent actuellement dans Herezh++. Je vais donc le rajouter. Pour votre travail, le mieux est alors d'utiliser une viscosité linéaire (sans le coef xn). Pour les grandes vitesses cela risque d'entraîner de très grands efforts ! A vous d'analyser les résultats, et dans le cas où la force visqueuse est trop grande il faudra diminuer la valeur du coefficient K visqueux.
Bon courage,
@+