Greetings!
I am trying to run a structural model with multipoint constraints and
Dirichlet constraints.
I know that the Dirichlet constraints can be added through the respective
DofMap method, but if I wish to add constraints of the form:
\sum_k a_k u_k = 0
or
u_j = 0 (Dirichlet)
using the DofMap::add_constraint_row() method, is it doable?
Currently I am trying to do for a Dirichlet constraint (hoping that the
penalty method would work) to get u_{dof_id} = 0
libMesh::DofConstraintRow row;
row[dof_id] = 1.e12;
dof_map.add_constraint_row(dof_id, row);
and for the multipoint constraint the following to get u_{ref_dof_id} =
u_{dependent_dof_id}
libMesh::DofConstraintRow row;
row[ref_dof_id] = 1.e12;
row[dependent_dof_id] = -1.e12;
dof_map.add_constraint_row(dependent_dof_id, row);
But this is not working.
Any advice would be appreciated.
Thanks,
Manav
------------------------------------------------------------------------------
New Year. New Location. New Benefits. New Data Center in Ashburn, VA.
GigeNET is offering a free month of service with a new server in Ashburn.
Choose from 2 high performing configs, both with 100TB of bandwidth.
Higher redundancy.Lower latency.Increased capacity.Completely compliant.
http://p.sf.net/sfu/gigenet
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
