SfePy April 2017

sfepy@python.org

12 participants
17 discussions

How to properly build a mesh?
by rmk236＠gmail.com 19 Apr '17

19 Apr '17

I am trying to build a custom rectangle mesh, but it seems that I am failing to do a proper one. I wrote a small python script to build the rectangle, using Delaunay triangulation [1]. This results in a file [2] that looks OK when compared with '/meshes/2d/rectangle_tri.mesh'. However, when I try to use this file into the a slightly modified version of '/examples/linear_elasticity/linear_elastic_interactive.py' [3], it raises an orientation error [4]. I really don't undestand what seems to be the issue here. Is it an ordering in the triangular elements? Does Sfepy presupose a specific ordering? [1] http://www.codesend.com/view/b6612f600b38f6f8552dc89ce072c938/ [2] http://www.codesend.com/view/25adf8bbe8e29f151d068b7ba08b4c82/ [3] http://www.codesend.com/view/910a5af4ba2a2b46f8a83572a7280529/ [4] Error below: > sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8] (MyMesh.mesh)... > sfepy: ...done in 0.00 s > sfepy: warning: bad element orientation, trying to correct... > sfepy: warning: bad element orientation, trying to correct... > Traceback (most recent call last): > File "./linear_mod.py", line 98, in <module> > main() > File "./linear_mod.py", line 44, in main > domain = FEDomain('domain', mesh) > File "/home/homeGlobal/lcampos/.local/lib/python2.7/site-packages/sfepy-2017.1-py2.7-linux-x86_64.egg/sfepy/discrete/fem/domain.py", line 61, in __init__ > self.fix_element_orientation() > File "/home/homeGlobal/lcampos/.local/lib/python2.7/site-packages/sfepy-2017.1-py2.7-linux-x86_64.egg/sfepy/discrete/fem/domain.py", line 146, in fix_element_orientation > raise RuntimeError('elements cannot be oriented! (%s)' % key) > RuntimeError: elements cannot be oriented! (2_3)

2 3

new mailing list
by Robert Cimrman 10 Apr '17

10 Apr '17

Hi everyone, I have transitioned our mailing list over to https://mail.python.org/mm3/mailman3/lists/sfepy.python.org/ As you can see, your subscription has been transferred, but you might need to tweak you e-mail delivery settings. Best regards, r.

1 0

migrating to @python.org
by Robert Cimrman 09 Apr '17

09 Apr '17

Hi everyone, as described in [1] I will do the move from Google groups to a new Mailman 3 mailing list on Python.org. I will try to do that on Monday 10.4. - all existing subscriptions will be transferred to the new list, and I will also port the content. Best regards, r. [1] https://github.com/sfepy/sfepy/issues/387

1 0

Calculate diffusive flux from solution in interactive mode
by Scott Calabrese Barton 09 Apr '17

09 Apr '17

I've made a lot of progress learning sfepy this weekend but I'm stuck on post-processing. I'm in interactive mode (using Jupyter) using a toy model that solves the laplace equation on a rectangle with dirichlet conditions on the ends (image below). Now, I'd like to evaluate the flux at designated positions, say on a line or a set of points. I'm able to evaluate the flux field using: flux=pb.evaluate('ev_diffusion_velocity.3.Omega(s.D, c)', copy_materials=False, mode='qp') where 'c' is the solution field variable and s is the material. This produces a set of flux vectors at quadrature points. What I can't seem to do is put this data in a proper field variable so that I can evaluate it with probes, say for example at the boundaries. I've looked at the examples but have found them too complex to reduce to a few lines of code. Can anyone suggest the easiest way to get this data into a field variable? <https://lh3.googleusercontent.com/-IGio8BtAJJE/WOF2SY3wmCI/AAAAAAAAAKI/q9S1…>

2 3

From Solidworks to sfepy
by François Gingras 09 Apr '17

09 Apr '17

Hi all, I am new using Sfepy and I am trying to get a part from solidworks, mesh it (using whatever) and then do some FEA on Sfepy. So right now, I exported the part in STEP and used GMSH to generate the mesh. But the VTK/MESH mesh file contain more dans 1 cell type. So Sfepy pop this error : "meshes with several cell kinds are not supported!" I am trying to figure how to have only one cell kind. Is this seem right? Solidwork -> Step file -> Gmsh -> mesh/vtk file -> sfepy

2 1

Finding versions of VTK and pysparce
by Gary Roach 09 Apr '17

09 Apr '17

Debian linux stretch OS Intel i5750 processor 64bit KDE desktop Directory structure: /root/anaconda3/lib/python3.5/site-packages/sfepy Home Directory structure: /home/gary/Downloads I am still trying to get a clean installation of sfepy. I am still missing vtk and pysparce libraries. I need to get components that are compatible with an amd-64, python3.x system. Any help with finding these packages and help with their installation will be sincerely appreciated. Thank you for your responses Gary R

3 8

migration of mailing list
by Robert Cimrman 05 Apr '17

05 Apr '17

Hi everyone, I consider migrating this mailing list over to mailman on Python.org, see [1]. Let me know if you have any suggestions/recommendations. Best regards, r. [1] https://github.com/sfepy/sfepy/issues/387

1 0

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

SfePy April 2017