ff3d-users
[Top][All Lists]

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

 From: Stephane Del Pino Subject: Re: [ff3d-users] [ Solve an elasticity problem ] Date: Sun, 25 Jan 2004 14:48:53 +0100 User-agent: KMail/1.5.4

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