I work on geological problems and I use the Finite Element Method. But this question can be applied on other classical mechanical problems.

I work on implicit 3D surfaces (which represent the limits between two geological layers aka two media).

I have to impose displacement on these surfaces (Dirichlet Condition). The condition that I impose on surfaces are transfered on the surrounding nodes. I make the assumption that I can use the " classical" way to apply Dirichlet Condition (penalty method in my case) on these nodes, that are NOT on the the boundaries.

My question is: do you think this assumption is valid? Have you got some references about the subject? As I understand the problem I try to apply Boundary Conditions not on the boundaries...

Thank you in advance.



We try to solve this problem: $$ \left\{ \begin{aligned} \sigma_{ij,j} + F_{i} &= 0 & &\text{in the domain $\Omega$}&\\ u_i&= q_{i}& &\text{on the boundary $\Gamma_{q}$}&\\ \sigma_{ij}n_{j} &= h_{i}& &\text{on the boundary $\Gamma_{h}$}&\\ u_i &= b_{i}& &\text{punctually in the domain $\Omega$}&\\ \end{aligned} \right. $$

We run the FEA on a mesh. We basically try to impose Dirichlet Condition on the nodes which are inside the domain $\Omega$ (fourth line). We constrain the displacement value of $u_i$ to $b_i$ for these nodes.

  • $\begingroup$ Can you explain in mathematical terms what the condition is you want to impose? $\endgroup$ – Wolfgang Bangerth Oct 1 '15 at 11:17
  • $\begingroup$ I have edited the question. $\endgroup$ – AntoineMazuyer Oct 1 '15 at 13:34
  • $\begingroup$ Yes, it's possible to prescribe "boundary conditions" in the domain and that's one of the greatest advantages of FEM. We can do that because in FEM our "domains" are the finite elements. $\endgroup$ – Gustavo Dec 5 '18 at 15:38

The solution of the equation you are looking for is in the space $H^1$ of functions that have one weak derivative, but in 2d and 3d, this does not imply that the solution is in fact continuous. As a consequence, it is not possible to define the value of a solution at individual points, and equally consequentially, it is mathematically not possible to prescribe (boundary) values at individual points.

But I suspect that you're not actually interested in prescribing values for the displacement at individual points, but for all points along a line (in 2d) or surface (in 3d), and this is perfectly valid. (Mathematically speaking, this is so because the trace operator is well defined on $H^1$). At the discrete level, this is then equivalent to prescribing the values on all nodes that lie on this line/surface.

Conceptually, you can think of this as a crack in your domain that is occupied by an infinitely thin device whose displacement you can manipulate and that attached to the two sides of the crack.

  • $\begingroup$ OK thank you very much for these informations ! $\endgroup$ – AntoineMazuyer Oct 5 '15 at 12:15

Well you can model slip between geologic faults (a discontinuity) using constraint equations which can be thought of as loading at the interface instead of the boundary. You simply specify the relative slip and solve the problem. E.g., see the top-left figure in https://bitbucket.org/stali/defmod/wiki/Gallery

If that is indeed what you want then I would suggest that you use Lagrange multipliers instead of the penalty method. Off course you can use any constraint equation you want but it should make physical sense.

  • $\begingroup$ Thank your for suggestions but my initial question is more like "Is FEA, which is based on variational formulation of boundary values problem, still a valid tool to solve a constrained equation with constraints on values inside the domain (which is no longer a boundary value problem)?" EDIT: I take a look to tour code it's pretty interesting $\endgroup$ – AntoineMazuyer Oct 1 '15 at 14:37
  • $\begingroup$ Yes it is. In theory you can think of the fault interface as a boundary (there are coincident nodes). $\endgroup$ – stali Oct 1 '15 at 14:54
  • $\begingroup$ Ok but the surfaces on which I impose the boundary conditions are horizons. (interface between two layers) $\endgroup$ – AntoineMazuyer Oct 1 '15 at 15:10

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.