On 03/05/2015 05:01 PM, Ronghai Wu wrote:
在 2015年3月3日星期二 UTC+1下午11:35:05，Robert Cimrman写道：
<snip>
Stress is plotted in this way. A line along x-axis(or y-axis) across the center of hole, plot from the center of hole to boundary, distance(r) is normalized by hole radius(R).
BTW. you can get the stress by using the ev_cauchy_stress term, just like the strain... also check the primer for data probes.
I tried to do this, but got error which I cannot fix (please run "elasticity-periodicBC3.py"). Additionally, I have some questions regarding the strain and stress calculated by pb.evalate(please see "questions about strain and stress.pdf").
You just need to use the probing code from examples/linear_elasticity/its2D_interactive.py, because your script is also interactive. Attached is the fixed script (+ some additional small fixes).
As for the questions, the strain and stress coming from ev_cauchy_*() terms are without the prestrain. You need to subtract that yourself.
r.