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
maillage d'un fil de forme sinusoÏdale de rayon r
mc.renoult


Inscrit le: 03 Mar 2009
Messages: 1
Répondre en citant
Bonjour à tous,

Je cherche à modéliser en 3D un fil magnétique de forme sinusoÏdale de rayon r (bien contrôlé), placé au sein d'une autre structure magnétique régie par l'équation de Laplace.
Pour traiter ce problème sous cast3M, je fais une analogie entre le potentiel scalaire magnétique et la Température.
Pour l'instant il m'est possible de créer un fil sinusoidal au sein de la structure par le code suivant:

*******************
*******************

* option de calcul et choix de l'élément

opti dime 3 echo 1 elem tet4;

* construction du domaine limité à droite et à gauche par les pôles de *Faraday


trans1 = 0. ;

lp1 = 32.512 + trans1 ;
lp2 = 8.5 + trans1 ;
lp3 = 12.64 + trans1 ;
lp4 = 15.952 + trans1 ;
lp5 = 22.99 + trans1 ;
lp13 = 11.6969 + trans1 ;
lp16 = 14.8937 + trans1 ;

p1 = lp1 0. -52.578;
p2 = lp2 0. 0.;
p3 = lp3 0. 16.146;
p4 = lp4 0. 26.082;
p5 = lp5 0. 40.572;
p13 = lp13 0. -7. ;
p16 = lp16 0. -14. ;

O1 = 0. 0. 0.;

p11 = 0. 0. 40.572 ;
p14 = 0. 0. -52.578 ;
p12 = 0. 0. -7. ;
p15 = 0. 0 -14. ;

nelem = 27 ;


d116 = droite (4/3*nelem) p1 p16;
d1613 = droite (nelem/2) p16 p13;
d132 = droite (nelem/2) p13 p2;
d23 = droite (nelem/2) p2 p3;
d34 = droite (nelem/3) p3 p4;
d45 = droite (nelem/2) p4 p5;
d511 = droite (2*nelem/3) p5 p11 ;
d11O1 = droite (4*nelem/3) p11 O1 ;
dO112 = droite (nelem/2) O1 p12;
d1215 = droite (nelem/2) p12 p15 ;
d1514 = droite (nelem) p15 p14 ;

d141 = droite (nelem) p14 p1;
dO12 = droite (nelem/2) O1 p2 ;
d1516 = droite (nelem) p15 p16 ;
d1213 = droite (2*nelem/3) p12 p13 ;

* definition du domaine (surface et bords)
* le volume est engendré par translation de la surface 2D


cdom1 = d23 et d34 et d45 et d511 et d11O1 et dO12 ;
s1 = surf (cdom1) 'PLANE' ;
cdom2 = d116 et d1516 et d1514 et d141 ;
s2 = surf (cdom2) 'PLANE' ;
cdom3 = d1613 et d132 et dO12 et dO112 et d1215 et d1516 ;
s3 = surf (cdom3) 'PLANE' ;
s4 = s1 et s2 et s3 ;
volu1 = s4 volu (2*nelem/3) 'TRAN' (0. 22.5 0.) ;


* Traçage pour vérification

trace cach volu1 ;

* Creation des surfaces où sont appliquées les conditions limites
* par translation du contour cdom4
* Creation du maillage associé

cont1 = d116 et d1613 et d132 et d23 et d34 et d45;
cont2 = d11O1 et dO112 et d1215 et d1514;
scont1 = cont1 tran (2*nelem/3) (0. 22.5 0.) ;
scont2 = cont2 tran (2*nelem/3) (0. 22.5 0.) ;

elim (0.0001) volu1 (scont1) ;
elim (0.0001) volu1 (scont2) ;

* creation du fil pour nelem=27 et dy=I*Lx/18 ici Lx=22.5 mm :

PI = 3.14159 ;
F0 = 2.5 - 7. ;
PFIL1 = 3.5 0. F0 ;
PFIL2 = volu1 POINT 'PROC' PFIL1 ;
FIL1 = PFIL2 ;
I=1 ;
REPETER BAUFIL 18 ;
RESU2 = COS (I*3*180/nelem);
F2 = 2.5*RESU2 - 7. ;
F3 = I*1.25 ;
PFIL1= 3.5 F3 F2 ;
PFIL2 = volu1 poin 'PROC' PFIL1 ;
FIL1 = FIL1 et PFIL2 ;
I = I+1 ;
FIN BAUFIL ;


*******************
ELFIL1 = ELEM VOLU1 APPUYE LARGEMENT FIL1 ;
VOLSFIL1 = DIFF VOLU1 ELFIL1 ;
TRAC CACH (ENVE VOLSFIL1) ;
TRAC (VOLSFIL1 ET (ELFIL1 COUL ROUG)) ;

*******************

dom = VOLSFIL1 et ELFIL1 ;

*******************
*******************

Mais je ne trouve pas dans la liste des opérateurs de Cast3M l'opérateur qui me premettrait à partir de cet ensemble de points formant un fil sinusoîdal de former un fil de rayon r.
J'ai essayé d'utiliser l'opérateur point sphere afin d'approximer mon fil de rayon r par un ensemble de sphères contiguËs de rayon r : ça marche mais d'une part il ne s'agit que d'une approximation et d'autre part je ne contrôle pas trés bien le rayon r de mon fil (et je ne peux pas atteindre des rayons r aussi petits que souhaités)


Et puis dans ces divers cas, pour créer mon fil je m'appuie sur les points du maillage existant déjà dans la structure magnétique dans laquelle est placé le fil (utilisation de l'opérateur point proche). Ce qui fait que le rayon du fil est limité à la résolution spatiale imposée par le choix du maillage de la structure globale.
Or je souhaite créer des fils de petite taille (typiquement r=0.44 mm) et donc pouvoir choisir le nombre d'éléments au sein du fil dans son maillage (donc rendre le maillage au sein du fil indépendant de celui de la structure mère comme je le faisais en 2D où simplement le fil était assimilé à un disque en première approximation (car en fait il s'agit d'une ellipse)).



Ce calcul 3D devrait nous servir pour valider une approximation (et surtout en déterminer les limites) sur le potentiel magnétique obtenu lors du même calcul version 2D déjà effectué.

Si quelqu'un a une idée sur l'opérateur à utiliser pour générer le maillage d'un tel fil de rayon r, je suis preneuse.

merci par avance,
marie-charlotte
Voir le profil de l'utilisateurEnvoyer un message privé
maillage d'un fil de forme sinusoÏdale de rayon r
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