Re: [ff3d-users] [ Solve an elasticity problem ]

[Stephane Del Pino]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

On Sunday 25 January 2004 09:56, hackervalley wrote:
Hello.
I will first go with your second problem.
> Here my data .
>
> vector n = (10,10,10);
> vertex a = (0,0,0);
> vertex b = (2,2,2);
> mesh M =structured(n,a,b);
> scene S = pov("void.pov");
> domain O = domain(S);
>
> solve(u) in O by M
> {
> pde(u)
> -div(grad(u)) = 1;
> u =0 on M;
> };
>
> save(medit,"u.mesh",u,M,unix);
>
> I hope you can help me .
Yes. The save instruction may be quite non-intuitive. You have to save the 
mesh *and* the solution (this separation was done to deal with transient 
problems).

so write

save(medit,"u",u,M,unix); // generates the 'u.bb' file
save(medit,"u",M,unix);   // creates 'u.mesh' file

Note that extensions are added automatically and that unix keyword is the 
default value.

Concerning your elasticity problem. You cannot defined it such a way. The 
language do not allow this kind of construction and will not before a long 
time --- this is version 2 of the language and we may wait for version 3 
which is not planed up to now since there are others priority.

If you want to deal with elasticity. You have to use variationnal formula.
Try to mimic the 'testsuite' example in attachement

Regards,
Stephane.
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.4 (GNU/Linux)

iD8DBQFAE8lFbZyhKhZwdc0RArPkAJ97ZbO55+dtveI9ojfcxP2z4NlM2ACgnIma
H4xIHb5MosPM3ZjWG46gjbA=
=mmEv
-----END PGP SIGNATURE-----

elasticity.ff
Description: Text Data