Dear all,
I have managed to solve my problem using a H matrix which in my
case is 3*3 sized and diagonal (I work in 3D).
The diagonal terms are 1 if the corresponding direction is fixed,
and 0 if the corresponding direction is free.
Hope it can be useful for others !
Have a good week-end,
Antoine
On 07/02/2017 03:47 PM, Antoine Mazuyer
wrote:
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" <address@hidden>
To: address@hidden
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 !
--
--
Team: RING | GeoModeles | UMR7359
GeoRessources
Address: F221 | ENSG | Campus Brabois | 2, rue du
Doyen Marcel Roubault
TSA70605 | 54518 Vandoeuvre-les-Nancy Cedex | France
Email: address@hidden
Tel.: (+33)3 83 59 64 50
--
--
Team: RING | GeoModeles | UMR7359
GeoRessources
Address: F221 | ENSG | Campus Brabois | 2, rue du
Doyen Marcel Roubault
TSA70605 | 54518 Vandoeuvre-les-Nancy Cedex | France
Email: address@hidden
Tel.: (+33)3 83 59 64 50
|