Hi, Robert,
Thanks a lot for your answer.
ouyang
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 pointloaded
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).
r.