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

Christophe Geuzaine geuzaine at acm.caltech.edu
Tue Oct 26 09:48:53 CEST 2004

```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

```

More information about the getdp mailing list