On 11/30/2015 08:54 PM, nayan...@gmail.com wrote:
Robert, Thanks a lot for your reply. ps0 = [[0.0, 1.0, 0.0], [ 0.0, 1.0, 0.0]] # <- starting points ps1 = [[5.0, 1.0, 0.0], [ 0.0, 1.0, 1.0]] # <- ending points I know the meaning of it. My question was how to plot the displacement at y=1 from x=0 to x=5.
This is done only by (assuming you want x in [0, 5], y = 1, z = 0):
ps0 = [[0.0, 1.0, 0.0]] ps1 = [[5.0, 1.0, 0.0]]
You do not need the second pair of points.
Does it help?
r.
Sorry for asking too my questions. I am new to SfePy.
Regards , Nayan
On Monday, November 30, 2015 at 4:21:22 PM UTC+1, Robert Cimrman wrote:
On 11/30/2015 03:53 PM, naya...@gmail.com <javascript:> wrote:
Robert, I changed " stiffness_from_youngpoisson(2, young, poisson)
by (it's 3D problem)
stiffness_from_youngpoisson(3, young, poisson) "
Now it's fixed. (Output is attached, Result).
But I tried to plot the displacement, strain and stress only u1 and u2 are on the plot, as per my understanding u1 and u2 are displacements in x and y directions. How do I plot u3 ( displacement in the z direction). I referred
on line 43 of Cuboid_graph.py, change
for ic in range(vals.shape[1] - 1):to
for ic in range(vals.shape[1]):its2D_interactive and made changes according to my problem definition. Also at the top surface I want to plot the displacement from x=0 to x=5 but in the plot it is showing only from x=0 to x=1.
The lines
ps0 = [[0.0, 1.0, 0.0], [ 0.0, 1.0, 0.0]] # <- starting points ps1 = [[5.0, 1.0, 0.0], [ 0.0, 1.0, 1.0]] # <- ending pointsmean, that the first probe goes from [0.0, 1.0, 0.0] to [5.0, 1.0, 0.0], and the second from [ 0.0, 1.0, 0.0] to [ 0.0, 1.0, 1.0] - it corresponds to the xlabels in the images. Try fixing that.
r.
Could you please also help me to fix this aswell?
Thanks in advance, Nayan
On Monday, November 30, 2015 at 3:00:23 PM UTC+1, Robert Cimrman wrote:
Hi Nayan,
when strange things happen, it is usually worth to check your region definitions.
Use:
./simple.py Cuboid_1.py --save-regions-as-groups
and fix the bottom region:
'Bottom' : ('vertices in (y<=-0.9999)', 'facet'),Also, in Cuboid_2.py, replace
stiffness_from_youngpoisson(2, young, poisson)
by (it's 3D problem)
stiffness_from_youngpoisson(3, young, poisson)
Note also, that having the code in the module namespace causes it to be called also before running the simulation (on input file import).
Let me know if that helps.
r.
On 11/28/2015 09:07 PM, naya...@gmail.com <javascript:> wrote:
Hi Robert, I calculated displacement, strain and stress for the block under
loading. I am not convinced with the answer. Displacement and strain are maximum at the corners, but stress is maximum in the middle. I think my problem description file is completely messed up.Is it?
Also when I try plot displacement, strain and stress components, one of the components is missing in the graph.
It would be a great help if you could help me. I have attached allarea the
files.
Regards, Nayan
On Friday, November 27, 2015 at 4:55:00 PM UTC+1, Robert Cimrman wrote:
On 11/27/2015 03:30 PM, naya...@gmail.com <javascript:> wrote: > I am using the regular package now and I got the output. I am
confused
with > the units. What's the unit of stress here? Is it N/M2?
SfePy always uses a consistent unit set, so it is up to you to interpret the units. See also [1].
r.
[1] http://www.continuummechanics.org/cm/techforms/SystemsOfUnits.html