Well, this example is somewhat artificial - it just demonstrates, how to implement easily strain or strain rate dependent materials. The particular relation for the elastic coefficient mu, as defined in get_pars(), is:
mu = 1 + norm((strain(t) - strain(t-1)) / dt)
where t is the time (step). I would not search a particular "real" material behind that.
r.
On 09/23/2015 12:10 AM, pyontaku14 wrote:
Hi Robert,
Thank you for your instruction. I think I could obtain residual force as Fx.png. Is this Fx=0.000135 valid for Ux=0.04? I could not understand how "def get_parts" defines material nonlinearity.
Thank you, Takuo
2015年9月9日水曜日 6時02分10秒 UTC+9 pyontaku14:
Hi,
Could you please tell me how to display "cylinder.h5", the result of material_nonlinearity.py? Usually, the result file is vtk. But the result of material_nonlinearity.py seems only h5. The OS I use is Lubuntu15.04-64bit on VM Ware player.
Best regards, Takuo Fujita