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:
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 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.
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 coordinatesout['cauchy_stress'] = Struct(name='output_data', mode='cell', data=stress, dofs=None)
print out['cauchy_stress'].data
any suggestions would be appreciated very much.
Thanks
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 loadedrock 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.
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
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
point 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.
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
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
point 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.
In our FEM textbook for engineering students, we are taught that the nodes
of a element can be numbered both locally and globally starting from
- Based on the local numbering and global numbering, a element stiffness matrix can be assembled into corresponding position in a global stiffness matrix. Because of this reason, most engineers might think that mesh node in FEM have an id which is numbered starting from 1.
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.
Yes, post-processing/displaying the result. Following your suggestion, I
just know how to do that in mayavi.
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.
Thx again, Robert
ouyang
On 07/17/2014 11:59 AM, Ouyang wrote:
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
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
point 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.
In our FEM textbook for engineering students, we are taught that the nodes
of a element can be numbered both locally and globally starting from
- Based on the local numbering and global numbering, a element stiffness matrix can be assembled into corresponding position in a global stiffness matrix. Because of this reason, most engineers might think that mesh node in FEM have an id which is numbered starting from 1.
I see. We just use the Python conventions :) (= from 0 everywhere).
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.
Yes, post-processing/displaying the result. Following your suggestion, I
just know how to do that in mayavi.
Good. If you use that often, you might consider writing a domain-specific plot class (see [1]).
[1] http://sfepy.org/doc-devel/src/sfepy/postprocess/domain_specific.html
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.
Thx again, Robert
Hth, r.
participants (2)
-
Ouyang -
Robert Cimrman