> Dear all,
>
> I have two questions:
>
> 1- Is there any support for vector basis ? The reason I'm asking is I want to
> implement a full elasticity and without vector basis for finite elements it's
> a headache? Is there a way for me to build it using a scalar basis?
Not too much of a headache... Here is a snippet from a linear elasticity
application code. You will find it similar to the way we do Stokes in the
examples:
for (unsigned int qp=0; qp<qrule.n_points(); qp++)
{
const Real
mu_qp = 1,
lambda_qp = 1;
// Matrix
for (unsigned int i=0; i<n_u_dofs; i++)
for (unsigned int j=0; j<n_u_dofs; j++)
for (unsigned int i_dim=0, ioff=0; i_dim<dim; i_dim++, ioff
+= n_u_dofs)
for (unsigned int j_dim=0, joff=0; j_dim<dim; j_dim++, joff
+= n_u_dofs)
Ke(i+ioff,j+joff) +=
JxW[qp]*(
lambda_qp*
dphi[i][qp](i_dim)*
dphi[j][qp](j_dim)
+
mu_qp*
dphi[i][qp](j_dim)*
dphi[j][qp](i_dim)
+
((i_dim == j_dim) ?
mu_qp*(dphi[i][qp]*dphi[j][qp]) : 0.)
);
}
> 2- Is there a support for blocks ? say I want one part of mesh to have ID-1
> and second part to have ID-2?
In libMesh we specify these as different subdomains. You can then add
variables which live only on a given subdomain if need be.
------------------------------------------------------------------------------
EditLive Enterprise is the world's most technically advanced content
authoring tool. Experience the power of Track Changes, Inline Image
Editing and ensure content is compliant with Accessibility Checking.
http://p.sf.net/sfu/ephox-dev2dev
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
