Does esys-escript mesh a fixed one?

Asked by Aaron Chan

Hello guys,
Recently, I repeated some examples of esys-escript of user guide doc.
One of the example is as a elastic deformation of a cube caused by temperature.

Details are as followed:
Cube with a pre-defined temperature field.
pde is: -\sigma_ij, i = 0
bc: plane x=0, y=0, z=0 stress x y z are constraint to 0 respectively.

I wonder whether the mesh are fixed during the iteration of the pde solver.
When results obtained, can I get the information of the deformed mesh?
As a phd student majoring in geotechnical engineering, I do some simulations on the abaqus or plaxis.
Postprocess of these packages provide the deformed mesh.
That is why I am curious about the deformed mesh of esys-escript.
There is a lot of example of cookbook about the temperature or wave propagation simulations.
For these problems, maybe there is no need to update the mesh each time step, cause the mesh is only a topology tool.
However, for geotechnical problem, deformation of the soil or rocks can be what we focus on.
For instance, if there is rectangle soil region with left, right and bottom fixed.
And there is a point load at the mid-point of the top, can I use esys-escript to model this?
Thank you very much for any reply or advice!

Aaron Chan

Question information

Language:
English Edit question
Status:
Solved
For:
esys-escript Edit question
Assignee:
No assignee Edit question
Solved by:
Lutz Gross
Solved:
Last query:
Last reply:
Revision history for this message
Best Lutz Gross (l-gross) said :
#1

By default the mesh is not deformed but mesh nodes can be moved in a post-processing step:

domain=.....
pde=LinearPDE(domain)

u=pde.getSolution()

newX=domain.getX()+u
domain.setX(newX)

Notice that not all solvers support this (finley does so).

Revision history for this message
Aaron Chan (aaaron) said :
#2

Thank you very much, Lutz.
Does this command domain.setX(newX) imply that,
not only the change of nodes coordinate but the shape funtion.
I mean, once domain reset a new u, shape function and Gauss point coordinate of the mesh changed with u?

Revision history for this message
Lutz Gross (l-gross) said :
#3

Correct.

Revision history for this message
Aaron Chan (aaaron) said :
#4

Thanks!

Revision history for this message
Aaron Chan (aaaron) said :
#5

Thanks Lutz Gross, that solved my question.