> Hello, > > If I want to use other forms of material property matrix instead of > isotropic linear elasticity matrix, how can I do it? Is there any methods > in getfem++ to deal with it? Say, if I want to use anisotropic material > property matrix. > > Thank you very much.
The following function in getfem_assembling.h asm_stiffness_matrix_for_vector_elliptic (MAT &M, const mesh_im &mim, const mesh_fem &mf, const mesh_fem &mf_data, const VECT &A, const mesh_region &rg = mesh_region::all_convexes()) do assemble an arbitrary vector elliptic problem with a given order four tensor A. this is very general and thus not optimal if the elastic law can be described in a more simple manner than the whole order four elastic tensor (as for orthotropic laws for instance). The corresponding model brick has not been built, but is very simple to do adapting the one for isotropic elasticity. Yves. ------------------------------------------------------------------------- Yves Renard ([EMAIL PROTECTED]) tel : (33) 04.72.43.80.11 Pole de Mathematiques, INSA de Lyon fax : (33) 04.72.43.85.29 Institut Camille Jordan - CNRS UMR 5208 20, rue Albert Einstein 69621 Villeurbanne Cedex, FRANCE http://math.univ-lyon1.fr/~renard ------------------------------------------------------------------------- _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
