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
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.