ff3d-users
[Top][All Lists]

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

 From: hackervalley Subject: [ff3d-users] [ Solve an elasticity problem ] Date: Sun, 25 Jan 2004 09:56:14 +0100 User-agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.0.0) Gecko/20020623 Debian/1.0.0-0.woody.1

```Hello All,

First thank for ff3d.

Now , I would solve an elasticity problem on a cube.

Here are my data :

cube.pov

box
{
<0,0,0>,
<20,20,20>
pigment { color rgb <1,0,0>}
}
plane
{
<0,1,0>, 0
pigment { color rgb <0,1,0>}
}
plane
{
<0,1,0>, 20
pigment { color rgb <0,0,1>}
}

elasticity.txt

/* Elasticity */

/* Young Module */
double E = 210000;
/* Poisson value */
double poisson = 0.3;
double mu = E/(2*(1+poisson));

vector n = (200,200,200);
vertex a = (0,0,0);
vertex b = (20,20,20);

mesh M = structured(n,a,b);
scene S = pov("cube.pov");
domain O = domain(S,inside(<1,0,0>));

/* st stress d deplacement */
/* PFD Navier */
solve(st,d) in O by M
{
pde(st)
div(st) = 0;
s = E/(1+poisson)*(td+poisson/(1-2*poisson)* Trace(td)*mu);
st = 100 on <0,0,1>;
pde(d)
```
```d = 0 on <0,1,0>;
}

save(medit,"st.mesh",st,M,unix);
save(medit,"d.mesh",d,M,unix);

I don't see how compute Trace(td) and Transpose(d).
I would also notice as save in medit , I have also wrong data for medit

-- Medit,  Release 2.2 (Sept., 2002)
compiled: Tue Jan 21 18:34:49 CET 2003.

## ERR 0003, inmesh, WRONG DATA TYPE.
Input seconds:     0.00

Total running seconds:  0.01
Thank you for using Medit. Bye.

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 .

Thank you.

```