HI everyone,

I need your help on a problem of elasticity and more precisely on a problem of integration on a sub-domain after deformation of the whole geometry.

I would like to describe the deformation of a square sub-domain strictly included in a bigger square (on wich I resolve linear elasticity equations). Thus, I would like to calculate the averaging of certain quantities (stress, strain, displacement etc) coming from this sub-region only.

I have thought about using the **Subdomain** class but I donâ€™t really understand how it behaves through the deformation *i.e*

if I locally integrate my quantities with dx(1) (where the intiale subdomain is marked by â€ś1â€ť) following the resolution of the total weak problem, will this result in the correspondence â€śsub-domain â†’ deformed sub-domainâ€ť or simply on the portion corresponding to the undeformed square (which I donâ€™t want)?

Thank you in advance for your answer

AMM

PS: creating a physical surface on GMSH for the sub-square may also be possible but as my geometry is a composite, this will result in a double physical label for a single entity and Iâ€™m not sure GMSH likes it.