[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Problem with dirichlet condition and source term

From: Konstantinos Poulios
Subject: Re: Problem with dirichlet condition and source term
Date: Mon, 2 Aug 2021 17:24:26 +0200

So you apply a displacement of [0,1,0] at the top, and your entire body moves by this amount as rigid body, right? This means that the top BC is applied correctly, now you need to check for mistakes in the other BCs/loads. How is CORNER defined? You cannot define a single point as a region.


On Mon, Aug 2, 2021 at 4:13 PM BLANCKAERT Benjamin <> wrote:

Thank you so much for your fast answer. I have already tried to do this :


But when i visualize the results I have the same displacement everywhere on my mesh. I have well written my constitutive law so normally I should have more little displacement when I go deeper but it is not the case.

/// Blanckaert Benjamin

De : Konstantinos Poulios <>
Envoyé : lundi 2 août 2021 15:38:55
Cc :
Objet : Re: Problem with dirichlet condition and source term
Dear Benjamin,

A source term corresponds to a Neumann BC. So the [0,1,0] vector that you apply is a actually a traction, not a displacement. You need to apply a non-homogeneous Dirichlet condition if you want to enforce a displacement of [0,1,0] at the top surface.


On Mon, Aug 2, 2021 at 3:17 PM BLANCKAERT Benjamin <> wrote:

Dear Getfem Users,

I am doing an internship in which I try to modelize a simple case of soil reinforced with rigid inclusions but I have a problem with the interaction between the source term I apply and the dirichlet conditions. I explain :

I apply Dirichlet Condition in the 4 corner at the bottom of my model in the simplest way : md.add_Dirichlet_condition_with_multipliers(mim,'u',mfu,CORNER) (I use Python module)

I also add a source_term_brick at the top of my model, here an example with an unit displacement:

d = [0,1,0]
md.add_source_term_brick(mim, 'u', 'd',TOP)

When I visualize the results with paraview I don't have the expected results. Indeed the displacement at the TOP is 2.2E-6 (which is very far from 1). 

I tried to identify the problem and I am almost sure that it is a problem with the dirichlet conditions which spread too much. 

Do you have any idea about How to fix it ? 

/// Blanckaert Benjamin

reply via email to

[Prev in Thread] Current Thread [Next in Thread]