How to assign initial value to strain (what fe space should be used)

General
1

Hi all,

In 3d solid mechanics, I want to assign initial values to strain. Should I define HCurl space for the strain? Is it correct to use:

fes = HCurl(mesh, order=3)
strain = GridFunction(fes)
strain.set(initial_values)

Regards
HDA

2

I guess it all depends on the variational formulation you are trying to implement. Could you share it with us?

3

Hi Gonzalo

Sorry for my late reply. I got some problems and couldn’t respond. I am actually trying to solve linear viscoelastic problem. I have attached my variation form. My problem is ∆e" has initial value and should be a tensor of 3x3 (similar to elastic strain). Is that correct to use Hdivdiv space for ∆e"?

Regards
HDA

Screenshot_20250318-000100~2
Screenshot_20250318-000100~21075×287 37.6 KB