From: Justin Perry on
does anyone know how to optimize 3-d structures using this code which can be found at http://www.topopt.dtu.dk/files/matlab.pdf if you are unfamiliar with it.

any help would be much appreciated.

Regards

Justin
From: Justin Perry on
"Justin Perry" <perry_282(a)hotmail.com> wrote in message <hnb64e$a60$1(a)fred.mathworks.com>...
> does anyone know how to optimize 3-d structures using this code which can be found at http://www.topopt.dtu.dk/files/matlab.pdf if you are unfamiliar with it.
>
> any help would be much appreciated.
>
> Regards
>
> Justin

I have managed to get further with my efforts but am now stuck on lines 90 and 91

90 k=[ 1/2-nu/6 1/8+nu/8 -1/4-nu/12 -1/8+3*nu/8 ...
91 -1/4+nu/12 -1/8-nu/8 nu/6 1/8-3*nu/8];

if anyone could possible explain these lines it would be of great help.

Regards

Justin
From: Bruno Luong on
"Justin Perry" <perry_282(a)hotmail.com> wrote in message <hndo1h$fjt$1(a)fred.mathworks.com>...
> "Justin Perry" <perry_282(a)hotmail.com> wrote in message <hnb64e$a60$1(a)fred.mathworks.com>...
> > does anyone know how to optimize 3-d structures using this code which can be found at http://www.topopt.dtu.dk/files/matlab.pdf if you are unfamiliar with it.
> >
> > any help would be much appreciated.
> >
> > Regards
> >
> > Justin
>
> I have managed to get further with my efforts but am now stuck on lines 90 and 91
>
> 90 k=[ 1/2-nu/6 1/8+nu/8 -1/4-nu/12 -1/8+3*nu/8 ...
> 91 -1/4+nu/12 -1/8-nu/8 nu/6 1/8-3*nu/8];
>
> if anyone could possible explain these lines it would be of great help.

It looks like this array is used to compute the elasticity stiffness (workload of the Galerking) matrix from Young and Poisson modulus of the materia (line 92-99). If you are familiar with mechanical FEM, these lines become trivial. If you are not, you need to read the literature on the subject.

Bruno