On 11/30/2015 03:53 PM, naya...@gmail.com 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
>>>>
>>>
>>
>>
>