Hi Dennis,
On 08/18/2015 01:45 PM, Dennis Perchak wrote:
Hi Robert,
Sorry for being slow to get back to this but I was busy on another project and then since I was using the version of SfePy that comes with the Enthought Canopy distribution I had to install the version you pointed me to and it took a bit to get everything to compile and work on my Mac. But I got it to work and the updated version with evaluating the piezo stress works. I have a simple test that I found in someone's thesis. If one has a material of large extent in the x and y directions compared to the z directions, then it will behave somewhat like a 1D system. In that case, if one sets phi on the bottom to be zero and some value on the top, then the strain that will arise due to the resulting field (Ez = -(phi_top - phi_bot)/dz can be estimated from the fact that the stress is zero. Thus, the strain eps_zz = (e_zz/C_zz)*E_z where e and C are the piezo and elastic constants. And this is in fact what happens when you do this in SfePy.
Good that it works.
Thanks for adding this to the code.
Another quick question if you don't mind. When using probes, is there a way to output the coordinate values of the probe positions rather than plot the probe values versus (essentially) the index number of the probe?
Not currently, but it will be easy to add - seems like a useful feature to have.
r.