encore quelques questions...
Christophe Geuzaine
Christophe.Geuzaine at ulg.ac.be
Wed Mar 14 15:07:11 CET 2001
trophime christophe wrote:
>
>
> 1/ Si maintenant je veux coupler
> les calculs thermiques et electromagnetiques (stationnaire et
> transitoire).
> Ca se passe comment?
> Pour la resolution il suffit de suivre l'exemple du &5.9.5 mais je
> vois pas comment passer les pertes joules de la partie électromagnétique
> au calcul thermique.
>
> J'ai introduit le terme source dans l'equation de la chaleur comme :
>
> { Name Hregion_qth_Heat_Axi ; Type Scalar ;
> BasisFunction {
> { Name sr ; NameOfCoef qthsr ; Function BF_Node ;
> Support DomainS_Th ; Entity NodesOf[ All ] ; }
> }
> Constraint {
> { NameOfCoef qthsr ; EntityType NodesOf ;
> NameOfConstraint SourceHeatDensity ; }
> }
> }
>
> et
> Formulation {
> { Name Heat_Axi ; Type FemEquation ;
> Quantity {
> { Name T ; Type Local ; NameOfSpace Hgrad_v_Th ; }
> { Name Qth ; Type Local ; NameOfSpace Hregion_qth_Heat_Axi ; }
> }
> Equation {
> Galerkin { Dt [ rho[] * Cp[] * Dof{d T} , {d T} ] ; In Domain_Th
> ;
> Jacobian VolAxi ; Integration GradGrad ; }
> ...
> Galerkin { [ -Dof{Qth} , {T} ] ; In DomainS_Th ;
> Jacobian VolAxi ; Integration GradGrad ; }
> }
>
> dans le fichier de données associes j'ai defini la contrainte :
>
> { Name SourceHeatDensity ; Type Assign;
> Case {
> { Region Ind ; Value 200. ; }
> { Region Ind1 ; Value 40. ; }
> }
> }
>
> J'ai vu que l'on pouvait choisir un type AssignFromResolution :
>
> { Name SourceHeatDensity ; Type AssignFromResolution;
> Case {
> { Region Ind ; NameOfResolution ??? ; }
> }
> }
>
> Je ne pense pas que ca marche comme ca puisque les resolutions
> sont imbriqués?
>
Non, c'est plus simple. Tu peux simplement faire reference a la grandeur
calculee dans la formulation magnetique. Dans ton exemple (si je
comprends bien), il faudrait ajouter, dans la formulation thermique,
Galerkin { [ - 0.5 * sigma[]* SquNorm[Dt[{a}]] , {T} ] ; In DomainS_Th ;
...
Et tu declares 'a' dans les 'Quantity' de la formulation thermique,
comme appartenant a l'espace fonctionnel approprie.
>
> 2/ Pour le calcul des forces volumiques j \times b, on peut
> mettre en place des methodes du type travaux virtuels?
Je n'ai jamais teste...
--
Christophe Geuzaine
Tel: 32 (0) 4 366 37 10 http://geuz.org
Fax: 32 (0) 4 366 29 10 mailto:Christophe.Geuzaine at ulg.ac.be