On 11/30/2015 10:21 PM, naya...@gmail.com wrote:
> Should I take out ps1?  Did u mean that? I'm really sorry for asking this.

No, use exactly this:

ps0 = [[0.0,  1.0, 0.0]]
ps1 = [[5.0,  1.0, 0.0]]

ps0 is the start point, ps1 is the end point of the line probe.

r.

> Regards,
> Nayan
>
> On Monday, November 30, 2015 at 9:13:33 PM UTC+1, Robert Cimrman wrote:
>>
>> On 11/30/2015 08:54 PM, naya...@gmail.com <javascript:> 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 points
>>>>
>>>> mean, 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
>>>> area
>>>>>>> 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
>> all
>>>> 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
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>
>>>>
>>>
>>
>>
>