[Getdp] LinearElasticity2D-how to apply a force field?

[Christophe Geuzaine]

Gianmario Pellegrino wrote:
> Dear Christophe,
> first of all congratulations for you very good job and many thanks for 
> sharing it.
> We're trying to evaluate (my colleague is in copy) the centrifugal 
> stress and deformation of a rotating lamination, but we can't figure out 
> how to express the forcing term as a function of x and y.
> Except for that, the LinearElasticity2d example fits very well.
> The centrifugal force per volume is expressed by:
> 
> f = Vector (k*x , k*y, 0)
> 
> We tried to declare this term as a function, using both $X,$Y and simply 
> X,Y,
> 
> F[Slab] = Vector[k*$X, k*$Y, 0 ]
> 

Right... I almost forgot about this.

Currently, you should use the "current values" $X, $Y and $Z only in 
pre-processing constraints, and use the functions X[], Y[] and Z[] in 
the formulations, i.e.:

F[Slab] = Vector[k*X[], k*Y[], 0]

(We definitely need to clarify the documentation about current values; 
or actually fix the code ;-))

Best,

Christophe

-- 
Christophe Geuzaine
Applied and Computational Mathematics, Caltech
geuzaine at acm.caltech.edu - http://geuz.org