Hi, Robert,
    Thanks a lot for your answer. 

On Saturday, July 12, 2014 3:05:08 PM UTC+8, Robert Cimrman wrote:
Hi Ouyang,

On 07/12/2014 05:25 AM, Ouyang wrote:
> Hi there,
>     I follow the point loaded disk example: its2D and simulate the 3d point
> loaded block as show in the attached figure.
> Here are my questions:
>     1. When I use "vertex <id>", the <id> seems not the node id in mesh and
> looks like an index in a list. My node id is 5 and I have to put 4 there.

Are you talking about the medit mesh format? The vertices in the mesh file are
not numbered at all, so the id is indeed given by their order. When a mesh is
displayed in medit, it labels the vertices starting from 1. In Python, the
indices start from 0.

>     2. Are there any way to customize contour bands? Since the point loaded
> rock block could fails at the center  point (0,0,0), I want to see the
> stress or strain details at the point (0,0,0) and would like to add more
> contour bands in.

Are you talking about post-processing / displaying the results? The postproc.py
script cannot do that, but you can open the mayavi pipeline after the scene is
displayed, and add/edit the visualization modules and filters there as you wish.

>     3. How to output a contour map for stress or strain components,
> \sigma_{ij} or \epsilon_{ij},  on  any cross section in sfepy?
>         The following script can give me stress components but I am still
> looking for their corresponding coordinates
>        out['cauchy_stress'] = Struct(name='output_data', mode='cell',
>   data=stress, dofs=None)
>        print out['cauchy_stress'].data
> any suggestions would be appreciated very much.

This can be done in mayavi or paraview using the appropriate filters (cutting
plane etc).