# Modeling deformation in Julia

**URL:** <https://discourse.julialang.org/t/modeling-deformation-in-julia/111941>\
**Category:** Performance\
**Tags:** question, package, fem, finite-element\
**Created:** [March 21, 2024, 5:30pm UTC](https://discourse.julialang.org/t/modeling-deformation-in-julia/111941 "2024-03-21T17:30:42Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![marianoarnaiz](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/marianoarnaiz/32/19377_2.png) [@marianoarnaiz](https://discourse.julialang.org/u/marianoarnaiz)\
**Post date:** [March 21, 2024, 5:30pm UTC](https://discourse.julialang.org/t/modeling-deformation-in-julia/111941/1 "2024-03-21T17:30:42Z")

</div>

Hi everyone, a general question here about Finite elements in Julia.

I am trying to do something rather “simple”. Imaging a large body (3D), fully constrained.  
In side there is a fault/fracture/weak zone inside. There is some displacement along this plane, and the medium/body deforms. I want to model this with FEM, it is very usual in my field but I though maybe one of the FEM packages of Julia could handle this.

Something like this:

 ![1-s2.0-S0040195120303681-gr3_lrg](https://global.discourse-cdn.com/julialang/original/3X/5/d/5dd50096c2ff0e9741b531da5106640e7b59033f.jpeg)

Any idea where to start?  
Cheers.  
M.

---

<div class="post-metadata">

**Author:** ![PetrKryslUCSD](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/petrkryslucsd/32/215825_2.png) [@PetrKryslUCSD](https://discourse.julialang.org/u/PetrKryslUCSD)\
**Post date:** [March 21, 2024, 5:37pm UTC](https://discourse.julialang.org/t/modeling-deformation-in-julia/111941/2 "2024-03-21T17:37:56Z")

</div>

Is there an actual fault, so that the left and right part of the domain slide alongside the fault?

---

<div class="post-metadata">

**Author:** ![marianoarnaiz](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/marianoarnaiz/32/19377_2.png) [@marianoarnaiz](https://discourse.julialang.org/u/marianoarnaiz)\
**Post date:** [March 21, 2024, 6:39pm UTC](https://discourse.julialang.org/t/modeling-deformation-in-julia/111941/3 "2024-03-21T18:39:48Z")

</div>

Yes! There might also me asperity on the fault. Although the displacement along the fault would be very small compare to the dimensions of the domain.

---

<div class="post-metadata">

**Author:** ![PetrKryslUCSD](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/petrkryslucsd/32/215825_2.png) [@PetrKryslUCSD](https://discourse.julialang.org/u/PetrKryslUCSD)\
**Post date:** [March 23, 2024, 4:57pm UTC](https://discourse.julialang.org/t/modeling-deformation-in-julia/111941/4 "2024-03-23T16:57:43Z")

</div>

In that case I regret but I have nothing: As far as I know no one has programmed that in Julia yet.

---

<div class="post-metadata">

**Author:** ![marianoarnaiz](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/marianoarnaiz/32/19377_2.png) [@marianoarnaiz](https://discourse.julialang.org/u/marianoarnaiz)\
**Post date:** [March 23, 2024, 11:58pm UTC](https://discourse.julialang.org/t/modeling-deformation-in-julia/111941/5 "2024-03-23T23:58:50Z")

</div>

Now that I think about it, for my case the fault is so small compare the domain, that I know the nodes at each side of the fault will slide, but not the entire domain.  
I found this package with two example that approach my setting:

[https://kailaix.github.io/AdFem.jl/dev/viscoelasticity\_earth/](https://kailaix.github.io/AdFem.jl/dev/viscoelasticity_earth/)

[https://kailaix.github.io/AdFem.jl/dev/earthquake/](https://kailaix.github.io/AdFem.jl/dev/earthquake/)

I will have to toy with them. Thanks for the help