RDforum.free.fr/Cast3M   Index du Forum
Site Web Cast3M
RDforum.free.fr/Cast3M
Forum d'entraide des utilisateurs de Cast3M
Répondre au sujet
Force ponctuelle en dehors des noeuds
Allo rédo


Inscrit le: 31 Mar 2009
Messages: 6
Répondre en citant
Bonjour,

Je n'arrive pas à trouver de solution pour imposer une force ponctuelle ailleurs que sur un nœud du maillage.

Le mieux que je peux faire c'est une suite d'instructions du type :

Pexc = 0.08 0.07 0.; --> le point ou je veux mettre la force
NoeuExc = SUR3 POIN 'PROC' Pexc; --> le nœud le plus proche
CHRPR = FORC (0. 0. 1.) NoeuExc; --> le chargement dudit nœud

Cette façon de faire est bien sûr inexacte !!!

La solution est pourtant triviale quand on connait les fonctions d'interpolation et qu'on sait dans quel élément se trouve le point Pexc. Je suis étonné que CASTEM ne propose pas une fonction toute faite pour faire ça. Pourtant après de longues recherches, je ne l'ai toujours pas trouvé ?

Connaissez-vous une méthode ?
Voir le profil de l'utilisateurEnvoyer un message privé
Allo rédo


Inscrit le: 31 Mar 2009
Messages: 6
Répondre en citant
J'ai toujours ce problème. Personne pour m'aider ?

N'existe-t-il vraiment pas dans Castem de fonction capable de calculer la projection d'un dirac sur les noeuds via les fonctions d'interpolation ? Ca se fait de tête sur des exemples simples.
Voir le profil de l'utilisateurEnvoyer un message privé
lguenhael
Administrateur

Inscrit le: 25 Nov 2006
Messages: 146
Localisation: Laboratoire Roberval - UTC
Répondre en citant
Bonjour,

On ne peut imposer un effort ponctuel ailleurs que sur un noeud du maillage. S'il n'y a pas de noeud là où vous souhaitez imposer votre effort alors il faut chercher une solution plus en amont, en faisant en sort qu'il y ait justement un noeud à cet endroit.

Cordialement.
Voir le profil de l'utilisateurEnvoyer un message privéEnvoyer un e-mail
Allo rédo


Inscrit le: 31 Mar 2009
Messages: 6
Répondre en citant
Merci pour votre réponse rapide.

C'est bien dommage car la solution est pourtant simple si on se situe au niveau du code (je l'ai déjà programée dans des codes perso) : il suffit de trouver l'élément qui contient le point, trouver les coordonnées locales du point (xi, eta, dzeta) dans le repère de l'élément (cela peut se faire par une méthode itérative simple faisant intervenir la matrice Jacobienne de la transformation) puis exprimer chacune des fonctions de forme pour ce triplet (xi, eta, dzeta) ce qui donne directement la répartition nodale de l'effort.

Réflexion faite, il y a dans Castem des "champs par point (noeuds) "des champs par élément (points de gauss)", il manque les "champs tout court (x,y,z)"
Voir le profil de l'utilisateurEnvoyer un message privé
Force ponctuelle en dehors des noeuds
Vous ne pouvez pas poster de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas voter dans les sondages de ce forum
Vous pouvez joindre des fichiers
Vous pouvez télécharger des fichiers
Toutes les heures sont au format GMT + 2 Heures  
Page 1 sur 1  

  
  
 Répondre au sujet