From: Xi YUAN <hill_yuan@ho...>  20070306 05:21:17

I am considering add shell elements to libmesh but not quite sure how to do it. The first probelm I encounter is how to build the shape function of such kind of elements. For shell elements, we need know the thickness and fiber vector of all FE nodes. But in libmesh, no such informations. I am considering add those variables into mesh_data and rewrite the function of shape function calculation like Real FE<2,LAGRANGE>::shape(const ElemType type, const Order order, const unsigned int i, const Point& p, const real thickness, real* fibervector) But it looks a little awkward because it changes the basic structure of this soft. Is there any suggestion about it? _________________________________________________________________ Don't just search. Find. Check out the new MSN Search! http://search.msn.com/ 