Mailman 3 April 2012 - SfePy

Re: Torque
by Robert Cimrman 09 Mar '18

09 Mar '18

In examples/large_deformation/hyperelastic.py a rotation by displacements is applied. By using a similar function the vectors defining the force couples could be defined for dw_surface_ltr (IMHO). Does it make sense? r. ----- Reply message ----- From: "Andre Smit" <freev...(a)gmail.com> To: <sfepy...(a)googlegroups.com> Subject: Torque Date: Sat, Dec 18, 2010 05:10 What is the best way to apply a torque load to a model? -- Andre -- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To post to this group, send email to sfepy...(a)googlegroups.com. To unsubscribe from this group, send email to sfepy-devel...(a)googlegroups.com. For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.

4 7

Modeling an adjustable X-ray optic
by Tom Aldcroft 21 May '13

21 May '13

I'm working on modeling a next-generation X-ray mirror for which the shape can be actively controlled by use of many thin piezo-electric actuators mounted on the mirror surface. The mirror is basically a glass conical paraboloid with a 1 meter radius and 200 micron thickness (e.g. http://en.wikipedia.org/wiki/X-ray_optics). Our project is currently using a proprietary FEA package, but the model setup and turnaround time is slow, in part because there is only one part-time engineer who can run it. SfePy looks like a great package and we're hoping that it could be used to automate running a large number of different cases. I've spent some time reading the documentation but I have a few questions that I hope can be answered before going too much further. I want to apologize in advance if some of my wording is imprecise, I have a physics background but this topic is a bit outside my realm... - Is SfePy appropriate for this problem? - If a specify a grid with about 800 x 400 points (azimuthal, axial) and about 10 boundary conditions (corresponding to mount points), what is the rough order of magnitude of time to compute the solution? Is it seconds, minutes, hours, or days? - The linear elastic examples show a problem with a specified displacement. How do I specify an input force? The piezo essentially provides a tensile force along the surface. - Is there a way to specify the problem and solve in cylindrical coordinates? This is the natural coordinate system. - How do I specify 6-DOF constraints which correspond to the mirror mounts? Thanks in advance for any help! Tom Aldcroft

3 6

Re: Fwd:
by Robert Cimrman 27 Apr '12

27 Apr '12

ok, I think I have a suspect: the working mesh has element type 12 (VTK_HEXAHEDRON), while the two that do not work have 11 (VTK_VOXEL). I have updated SfePy to support pixels and voxels - check out the current git master. The sinus.vtk seems to work, about the other file I am not sure, as it takes a long time, but it can be read... r. On 04/27/2012 10:27 AM, Xiaogai Li wrote: > Hi again, > > Here is one more model with the same problem. This geometry is much simpler > and only a few elements. > Thanks. > BR, > XG > > > ---------- Forwarded message ---------- > From: Xiaogai Li<xiaog...(a)gmail.com> > Date: Fri, Apr 27, 2012 at 9:54 AM > Subject: > To: cimr...(a)ntc.zcu.cz > > > Hi Robert, > > Thanks and please find the attached file. > > smooth_voxel_mesh.py > is the python script I used to run the smoothing. > > atlas_csf-label_legacy_notwork.vtk > is the file that with the problem "IndexError: index out of bounds". > > atlas_csf-label_legacy_part_works.vtk > > is the file that is part of the whole image, > and this image works with the python script. > > So I wonder what is the problem in file atlas_csf-label_legacy_notwork.vtk. > > Thanks a lot for your time. > BR > Xiaogai Li > > > > > > > On Fri, Apr 27, 2012 at 9:08 AM, Robert Cimrman<cimr...(a)ntc.zcu.cz> wrote: > >> Hi, >> >> could you send me (off-list) the file? It's a strange error. >> >> r. >> >> >> ----- Reply message ----- >> From: "Xiaogai Li"<xiaog...(a)gmail.com> >> To:<sfepy...(a)googlegroups.com> >> Subject: About the mesh_tools.smooth >> Date: Fri, Apr 27, 2012 00:37 >> >> >> Hello everyone, >> >> I met a strange error when using smooth_mesh. It seems that I could not >> create the object Domain, however, I have verified that the voxel-mesh was >> quite nice. >> >> coordssmooth = smooth_mesh(m, 4, 0.6307, -0.6347, None, True, False) >> File "/usr/local/lib/python2.7/dist-packages/sfepy/mesh/mesh_tools.py", >> line 505, in smooth_mesh >> domain = Domain('mesh', mesh) >> File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line >> 223, in __init__ >> self.setup_facets(verbose=verbose) >> File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line >> 375, in setup_facets >> obj.setup_unique() >> File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/facets.py", line >> 281, in setup_unique >> self.uid[0], self.uid[1:] = 0, ii[:-1] >> IndexError: index out of bounds

1 0

Re: About the mesh_tools.smooth
by Robert Cimrman 27 Apr '12

27 Apr '12

Hi, could you send me (off-list) the file? It's a strange error. r. ----- Reply message ----- From: "Xiaogai Li" <xiaog...(a)gmail.com> To: <sfepy...(a)googlegroups.com> Subject: About the mesh_tools.smooth Date: Fri, Apr 27, 2012 00:37 Hello everyone, I met a strange error when using smooth_mesh. It seems that I could not create the object Domain, however, I have verified that the voxel-mesh was quite nice. coordssmooth = smooth_mesh(m, 4, 0.6307, -0.6347, None, True, False) File "/usr/local/lib/python2.7/dist-packages/sfepy/mesh/mesh_tools.py", line 505, in smooth_mesh domain = Domain('mesh', mesh) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line 223, in __init__ self.setup_facets(verbose=verbose) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line 375, in setup_facets obj.setup_unique() File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/facets.py", line 281, in setup_unique self.uid[0], self.uid[1:] = 0, ii[:-1]IndexError: index out of bounds -- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To post to this group, send email to sfepy...(a)googlegroups.com. To unsubscribe from this group, send email to sfepy-devel...(a)googlegroups.com. For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.

1 0

About the mesh_tools.smooth
by Xiaogai Li 27 Apr '12

27 Apr '12

Hello everyone, I met a strange error when using smooth_mesh. It seems that I could not create the object Domain, however, I have verified that the voxel-mesh was quite nice. coordssmooth = smooth_mesh(m, 4, 0.6307, -0.6347, None, True, False) File "/usr/local/lib/python2.7/dist-packages/sfepy/mesh/mesh_tools.py", line 505, in smooth_mesh domain = Domain('mesh', mesh) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line 223, in __init__ self.setup_facets(verbose=verbose) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line 375, in setup_facets obj.setup_unique() File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/facets.py", line 281, in setup_unique self.uid[0], self.uid[1:] = 0, ii[:-1] IndexError: index out of bounds

1 0

Problem about "sfepy.mesh.mesh_tools.smooth_mesh" module for large datasets
by Xiaogai Li 26 Apr '12

26 Apr '12

Hi everyone, I have a problem about "sfepy.mesh.mesh_tools.smooth_mesh" module, when it comes to large data size, the running time is very long. I have a model with about 300,000 cells, now already taken more than 10 hours, but still not finished yet. So I wonder, if anyone has experience know the some tricks to make it running faster? Thanks for your time. BR, Xiaogai

2 2

Re: isfepy currently doesn't work
by Robert Cimrman 19 Apr '12

19 Apr '12

Hi, I cannot test it right now, but it seems something is wrong with mayavi installation. In regular ipython, is it possible to import mayavi and traits? r. ----- Reply message ----- From: "kassbohm" <kass...(a)gmail.com> To: "sfepy-devel" <sfepy...(a)googlegroups.com> Subject: isfepy currently doesn't work Date: Wed, Apr 18, 2012 17:35 Hi, with the latest git-resource: $ ./isfepy Traceback (most recent call last): File "./isfepy", line 102, in <module> main() File "./isfepy", line 99, in main init_session(ipython, **args) File "/home/kassbohm/SFePy_Install/sfepy/sfepy/interactive/ session.py", line 225, in init_session ip = _init_ipython_session(is_wx, _preexec_source, argv) File "/home/kassbohm/SFePy_Install/sfepy/sfepy/interactive/ session.py", line 131, in _init_ipython_session exec source in aux File "<string>", line 10, in <module> File "/home/kassbohm/SFePy_Install/sfepy/sfepy/postprocess/ __init__.py", line 1, in <module> from viewer import Viewer, ViewerGUI File "/home/kassbohm/SFePy_Install/sfepy/sfepy/postprocess/ viewer.py", line 17, in <module> from traits.api \ ImportError: No module named traits.api What can I do to solve this one? -- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To post to this group, send email to sfepy...(a)googlegroups.com. To unsubscribe from this group, send email to sfepy-devel...(a)googlegroups.com. For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.

1 1

isfepy currently doesn't work
by kassbohm 18 Apr '12

18 Apr '12

Hi, with the latest git-resource: $ ./isfepy Traceback (most recent call last): File "./isfepy", line 102, in <module> main() File "./isfepy", line 99, in main init_session(ipython, **args) File "/home/kassbohm/SFePy_Install/sfepy/sfepy/interactive/ session.py", line 225, in init_session ip = _init_ipython_session(is_wx, _preexec_source, argv) File "/home/kassbohm/SFePy_Install/sfepy/sfepy/interactive/ session.py", line 131, in _init_ipython_session exec source in aux File "<string>", line 10, in <module> File "/home/kassbohm/SFePy_Install/sfepy/sfepy/postprocess/ __init__.py", line 1, in <module> from viewer import Viewer, ViewerGUI File "/home/kassbohm/SFePy_Install/sfepy/sfepy/postprocess/ viewer.py", line 17, in <module> from traits.api \ ImportError: No module named traits.api What can I do to solve this one?

1 0

Re: Trying Installation... AttributeError: LooseVersion instance has no attribute 'version'
by Robert Cimrman 16 Apr '12

16 Apr '12

Hi Sven, Does any of these work: from enthought.mayavi import version or from mayavi import version I have mayavi 4.0.0, where the latter works. In 3.* versions, the former way worked, but you may have something older. Note this is the mayavi version, not VTK. Anyway, mayavi is not a hard dependency, so you should be able run sfepy without it, or with older version. r. On 04/16/2012 06:48 AM, KASSBOHM wrote: > Hello Robert, > > Mayavi version - VTK version 5.4.2 > Python 2.7.1+ > wxPython 2.8.11.0 > > Linux thnkpd 2.6.38-11-generic #48-Ubuntu SMP Fri Jul 29 19:02:55 UTC 2011 > x86_64 x86_64 x86_64 GNU/Linux > > Thanks for your help. > Sven > > On Sun, 15 Apr 2012, Robert Cimrman wrote: > >> What is your system and mayavi version? Is mayavi installed? Anyway, thanks >> for creating the >> issue, it looks like a bug in setup. >> >> In the meantime, you can try commenting out the mayavi version check in >> setup.py, around line >> 110. >> >> r. >> >> ----- Reply message ----- >> From: "kassbohm" <kass...(a)gmail.com> >> To: "sfepy-devel" <sfepy...(a)googlegroups.com> >> Subject: Trying Installation... AttributeError: LooseVersion instance has no >> attribute 'version' >> Date: Sun, Apr 15, 2012 14:43 >> >> >> Sorry, I also posted this question here: >> https://github.com/sfepy/sfepy/issues/180 >> >> Getting the following error during installation. How do I fix this? >> >> $ python setup.py buildext --inplace >> Traceback (most recent call last): >> File "setup.py", line 112, in >> version_getter=_mayavi_version) >> File "/home/kassbohm/SFePy_Install/sfepy-2012.1/build_helpers.py", >> line 320, in package_check >> if checker(have_version) < checker(version): >> File "/usr/lib/python2.7/distutils/version.py", line 296, in __cmp_ >> return cmp(self.version, other.version) >> AttributeError: LooseVersion instance has no attribute 'version' >>

2 4

Re: Trying Installation... AttributeError: LooseVersion instance has no attribute 'version'
by Robert Cimrman 16 Apr '12

16 Apr '12

What is your system and mayavi version? Is mayavi installed? Anyway, thanks for creating the issue, it looks like a bug in setup. In the meantime, you can try commenting out the mayavi version check in setup.py, around line 110. r. ----- Reply message ----- From: "kassbohm" <kass...(a)gmail.com> To: "sfepy-devel" <sfepy...(a)googlegroups.com> Subject: Trying Installation... AttributeError: LooseVersion instance has no attribute 'version' Date: Sun, Apr 15, 2012 14:43 Sorry, I also posted this question here: https://github.com/sfepy/sfepy/issues/180 Getting the following error during installation. How do I fix this? $ python setup.py buildext --inplace Traceback (most recent call last): File "setup.py", line 112, in version_getter=_mayavi_version) File "/home/kassbohm/SFePy_Install/sfepy-2012.1/build_helpers.py", line 320, in package_check if checker(have_version) < checker(version): File "/usr/lib/python2.7/distutils/version.py", line 296, in __cmp_ return cmp(self.version, other.version) AttributeError: LooseVersion instance has no attribute 'version' -- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To post to this group, send email to sfepy...(a)googlegroups.com. To unsubscribe from this group, send email to sfepy-devel...(a)googlegroups.com. For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.

2 3