search

thegooglecodearchive / sfepy

Automatically exported from code.google.com/p/sfepy
0 stars 0 forks source link

Linear elastic anisotropic materials #116

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago 
We have functions for isotropic materials, can this function be extended to
deal with anisotropy.

Original issue reported on code.google.com by freevryh...@gmail.com on 12 Apr 2010 at 12:20

GoogleCodeExporter commented 9 years ago 
See dw_lin_elastic, which accepts the material parameters in the general tensor 
form
- example use is e.g. in examples/biot/biot_npbc.py ...

Original comment by robert.c...@gmail.com on 12 Apr 2010 at 1:10

GoogleCodeExporter commented 9 years ago 
Please expand - to set the anisotropic material parameters in the 3 dimensions 
would
you need to do something like this:

        out['D'] = nm.tile(stiffness_tensor_lame(dim, lam=1.7, mu=0.3),
                           (coor.shape[0], 1, 1))
        out['D'] = nm.tile(stiffness_tensor_lame(dim, lam=2.7, mu=0.4),
                           (coor.shape[0], 2, 2))
        out['D'] = nm.tile(stiffness_tensor_lame(dim, lam=3.7, mu=0.45),
                           (coor.shape[0], 3, 3))

A

Original comment by freevryh...@gmail.com on 12 Apr 2010 at 2:30

GoogleCodeExporter commented 9 years ago 
You do not need stiffness_tensor_lame() at all - this function just builds the 
tensor
using the two Lame constants.

So to have anisotropic material, set the components of D directly. It is a 6x6 
tensor
in 3D, 3x3 in 2D, with the order of components described in [1] - it is 11 22 
33 12
13 23 in 3D, 11, 22, 12 in 2D. There is also explained the relation of 
isotropic D to
the Lame constants.

As an example of anisotropic material, an orthotropic material can be defined as
described at [2].

[1] http://docs.sfepy.org/doc-devel/src/sfepy/mechanics/matcoefs.html
[2] http://en.wikipedia.org/wiki/Hooke%27s_law#Orthotropic_materials

Original comment by robert.c...@gmail.com on 12 Apr 2010 at 2:40

GoogleCodeExporter commented 9 years ago

Original comment by robert.c...@gmail.com on 13 Apr 2010 at 1:51

GoogleCodeExporter commented 9 years ago 
Migrated to http://github.com/sfepy/sfepy/issues/118

Original comment by robert.c...@gmail.com on 30 Jan 2012 at 10:26