Hi Dennis,

could you try [1], if it does what you need? See the updated examples/multi_physics/piezo_elasticity.py example.

r.

[1] https://github.com/rc/sfepy

On 08/12/2015 01:15 AM, Dennis Perchak wrote:

Hi Robert,

Thanks for the response. And for making the code.

Dennis

On Wednesday, August 5, 2015 at 6:11:30 AM UTC-7, Dennis Perchak wrote:

I have been doing some modeling of a piezoelectric system and it works well. One can apply a field and see the resulting deformation. I also get how one can calculate the Cauchy strain and stress using the post-process hook and the terms 'ev_cauchy_strain.2.Omega(u)' and 'ev_cauchy_stress.2.Omega(solid.D, u)'. But since the constitutive law for the stress is the Cauchy term + the piezo-coupling term, how does one evaluate the stress and output it for given displacements and fields?

Dennis