Dear Yves,
Thank you for your help.
Ok so I have to declare a matrix field H and reuse it in the function
you have proposed.
I do not exactly understand the purpose of this matrix and how to
declare it ?
Best regards,
On 07/02/2017 02:34 PM, Yves Renard wrote:
Dear Antoine,
You can use
add_generalized_Dirichlet_condition_with_multipliers(mim, varname,
mult_description, region, dataname, Hname)
instead, see
http://getfem.org/python/cmdref_Model.html?highlight=dirichlet_condition#getfem.Model.add_generalized_Dirichlet_condition_with_multipliers
Best regards,
Yves.
----- Original Message -----
From: "Antoine Mazuyer" <[email protected]>
To: [email protected]
Sent: Sunday, July 2, 2017 1:30:26 PM
Subject: [Getfem-users] Prescribed Displacement in one direction only (free
slip)
Dear all,
I work with 3D linear elastic problems and I have to prescribed
displacements on some boundaries (Dirichlet Conditions)
I use this sample code:
std::vector< bgeot::scalar_type > F( mf_rhs.nb_dof() * 3 );
getfem::interpolation_function( mf_rhs, F,
my_class_function, attribute_id );
getfem::add_Dirichlet_condition_with_multipliers(
model, mim, "u",
mf_u, attribute_id, "Dirichlet" );
to add a Dirichlect condition to my region "attribute_id"
The problem is that my_class_function has an operator() that returns a
bgeot::base_vector with three coordinates, containing the three
prescribed displacements. With this kind of code, I have to prescribed
the displacement in the 3 directions
How can I just precribed the displacement in only one direction, and let
the others free ?
Thank you in advance for your advices !
--
--
*Antoine Mazuyer
<http://www.ring-team.org/component/liad/?view=user&limitstart=0&id=589>
* | *PhD. Candidate* | ASGA
*Team:* RING <http://ring.georessources.univ-lorraine.fr/home> |
GeoModeles | UMR7359 GeoRessources
*Address:* F221 | ENSG | Campus Brabois | 2, rue du Doyen Marcel Roubault
TSA70605 | 54518 Vandoeuvre-les-Nancy Cedex | France
*Email:* [email protected]
*Tel.:* (+33)3 83 59 64 50
