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 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. 

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 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
>>
>