Simulation d'une flexion '4' points en axis [Résolu] |
bthomas
|
Bonjour à tous,
Le sujet sur lequel je bosse en ce moment est la modéllisation d'une sorte de flexion 4 points réalisée par 2 tores sur une éprouvette circulaire. Je réalise un modèle axisymétrique du disque pour simplifier les calculs. Les tores sont de simples demi-disque qui sont liés à mon disque-éprouvette par des conditions de contacts. J'impose ma flexion grace à ces conditions. Le premier objectif était de réaliser un modèle fonctionnant sans me soucier des dimensions ni des efforts ou déplacements appliqués. Le calcul a l'air de se faire correctement, mais un problème apparait puisque maintenant que je souhaite jouer sur mes efforts ou déplacements imposés, les résultats ne changent pas! Je pense qu'il doit donc y avoir une erreur soit dans mes conditions de contacts soit dans ma façon d'appliquer mes efforts/déplacements. SI quelqu'un pouvait jeter un coup d'oeil et me dire ce qu'il en pense. Je précise que c'est la première fois que je me préoccupe de contacts sous castem et que je me suis donc largement inspiré de programmes en ma possession. merci d'avance |
||||||||||||||
_________________ Benjamin THOMAS 5 ème année de Mécanique des Structures Institut National des Sciences Appliquées de Rouen Etudiant Master, stagiaire Laboratoire d'Intégrité des Structures et de Normalisation, CEA Saclay |
bthomas
|
Merci de la rapidité de la réponse.
En modifiant la condition ainsi la déformée ne correspond pas du tout à ce que j'attends. cela provient sans doute de ce que tu soulignes. Mais je ne vois pas ce que j'ai fais qui provoque cela. Peux tu me dire comment as-tu noté ce problème? |
||||||||||||
_________________ Benjamin THOMAS 5 ème année de Mécanique des Structures Institut National des Sciences Appliquées de Rouen Etudiant Master, stagiaire Laboratoire d'Intégrité des Structures et de Normalisation, CEA Saclay |
bthomas
|
Je continue à travailler sur ce problème et je constate qu'en apportant la correction que tu propose, hormis la déformée qui ne correspond pas àa ce que je pensais obtenir, le problème d'"indifférence" de mon modèle au niveau de chargement persiste.
SI ce n'est pas le chargement que j'impose qui agit, que ce passe-t-il puisque le calcul se fait quand même! Si quelqu'un a déjà vu ça sur un autre cas ça m'intéresse! |
||||||||||||
_________________ Benjamin THOMAS 5 ème année de Mécanique des Structures Institut National des Sciences Appliquées de Rouen Etudiant Master, stagiaire Laboratoire d'Intégrité des Structures et de Normalisation, CEA Saclay |
lguenhael
Administrateur
|
Le "collage" des surfaces provient du fait que tu inverses l'orientation des lignes de contact alors que celles-ci sont déjà dans l'orientation voulue.
Il faut donc que tu supprimes les lignes : pdisqb = 'INVERSE' pdisqb ; et larch1 = 'INVERSE' larch1 ; Sinon après toutes ces modifications je ne constate pas d’indifférence. |
||||||||||||
|
bthomas
|
Merci de ton aide mais ça ne veut décidément pas changer sur ma machine, alors que j'ai fait les modifications (j'ai supprimé la ligne INVE pdisqb1.. également car dans le cas contraire ça ne converge pas). Le résultat est toujours le même et pour toutes les composantes des contraintes quelque soient les déplacements imposés.
Pourrais tu me mettre le prog que tu as modifié dans un post que je le fasse tourner chez moi ce soir histoire de voir si ça marche? Encore merci du coup de main, je desespere un peu d'arriver à le faire fonctionner! |
||||||||||||
_________________ Benjamin THOMAS 5 ème année de Mécanique des Structures Institut National des Sciences Appliquées de Rouen Etudiant Master, stagiaire Laboratoire d'Intégrité des Structures et de Normalisation, CEA Saclay |
lguenhael
Administrateur
|
Effectivement j’ai été un peut trop vite et je n’avais pas supprimé correctement ces deux lignes.
Cependant je maintient ce que j’ai dit, à savoir que pdisqb = 'INVERSE' pdisqb ; et larch1 = 'INVERSE' larch1 ; sont à supprimer. Mais il ne faut surtout pas supprimer la ligne pdisqb1 = inverse pdisqb1 ; car ça revient au même que de laisser pdisqb = 'INVERSE' pdisqb ; Cependant avec ces modifications (correctement effectuées cette fois) le calcul ne converge effectivement plus. Ceci est du au fait que le contact du tore supérieur est mal défini (je n’ai pas trouver pourquoi exactement) de ce fait le pièce à tendance à vouloir traverser ce tore et puisque plus rien ne la retient suivant z la résolution ne peut plus se faire. J’ai donc (par manque de temps) redéfini ce contact :
Et cette fois ça devrait fonctionner. |
||||||||||||||
|
bthomas
|
Merci encore de ton aide,
Je me doutais qu'il s'agissait d'une mauvaise définition d'un des contacts mais sans voir pourquoi. Je pars donc sur cette piste à partir du code que tu me propose. Je vais travailler là dessus et posterai suivant mes progrès. |
||||||||||||
_________________ Benjamin THOMAS 5 ème année de Mécanique des Structures Institut National des Sciences Appliquées de Rouen Etudiant Master, stagiaire Laboratoire d'Intégrité des Structures et de Normalisation, CEA Saclay |
bthomas
|
Un petit post pour signaler que le problème est réglé.
Comme je n'arrivais pas à trouver exactement d'où cela venait et que mon maillage devait subir de grosses modifications, j'ai recommencé complètement et le calcul a fonctionné du premier coup!! J'ai pourtant suivi exactement la même procédure, mis à part que cette fois je n'avais plus à comprendre comment fonctionnait les contacts linéiques sous castem. Je joins une des versions de mon programme pour ceux qui aurait besoin d'une base pour ce type de contacts. |
||||||||||||||
_________________ Benjamin THOMAS 5 ème année de Mécanique des Structures Institut National des Sciences Appliquées de Rouen Etudiant Master, stagiaire Laboratoire d'Intégrité des Structures et de Normalisation, CEA Saclay |
lguenhael
Administrateur
|
Bonjour et merci pour ce programme.
Je préciserai simplement que l’on peut éventuellement affichée le maillage du contact. Après :
Ajouter :
On peut alors observer tous les noeuds sur lesquels la procédure PASAPAS va vérifier les conditions de contact et ils sont nombreux… Or dans cet exemple, on sait d’avance que les contacts resteront localisés près des disques inférieur et supérieur. Il est donc conseiller de limiter le maillage des contacts aux noeuds dont la distance est inférieur à un certain critère que l’on précise dans IMPOse. Par exemple :
|
||||||||||||||||||
Dernière édition par lguenhael le 11 Mai , 2007 16:43; édité 1 fois |
bthomas
|
Merci de cette précision, je ne connaissais pas cette possibilité de l'opérateur IMPO.
Je comptais réduire mes surfaces de contact du disque manuellement. Les calculs effectués sur les noeuds 'en trop' n'alourdissent pas encore mon calcul mais c'est bien entendu très important dans l'optique de raffiner mon maillage. Sur ce je vous souhaite bon week end à tous! |
||||||||||||
_________________ Benjamin THOMAS 5 ème année de Mécanique des Structures Institut National des Sciences Appliquées de Rouen Etudiant Master, stagiaire Laboratoire d'Intégrité des Structures et de Normalisation, CEA Saclay |
Simulation d'une flexion '4' points en axis [Résolu] |
|
||
Powered by phpBB Forum © phpBB Group
Design by phpBBStyles.com | Styles Database.
Content © RDforum.free.fr/Cast3M
Design by phpBBStyles.com | Styles Database.
Content © RDforum.free.fr/Cast3M