Mailman 3 September 2014 - SfePy

Re: Torque
by Robert Cimrman March 9, 2018

March 9, 2018

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

1D line elements in SfePy
by Nimish Jan. 22, 2016

Jan. 22, 2016

I am currrently looking for FEM packages to help me solve a system of beams and columns, basically a collection of 1D bernoulli/timoshenko line elements. I started reading SfePy docs and i am getting the idea that doing the above is not really possible here, am i right? Are only 2D area elements permitted in SfePy? Or is there any direct support for solving 1D line elements too.. Cheers Nimish

4 9

SfePy system-wide installation
by G.L. Oct. 27, 2014

Oct. 27, 2014

Hi ! First of all, congratulation for this great piece of software ! I was trying to build and install SfePy system-wide (2014.1 or 2014.2), but ran into troubles. I followed the instructions on the website ("python setup.py build", then "python setup.py install"). There was no special error during the build (some warnings from gcc) or the installation. However, when I try to run the tests, nothing is working : > 79 test file(s) executed in 0.55 s, 92 failure(s) of 103 test(s) The errors are due to python that fails to import SfePy modules : ./run_tests.py --debug tests/test_cmesh.py > .... > .... > .... > File "tests/test_cmesh.py", line 96, in test_cmesh_counts > from sfepy.discrete.fem import Mesh > File "/scratch/packages/sfepy-2014.2/sfepy/discrete/__init__.py", line > 8, in <module> > from variables import Variables, Variable, FieldVariable, > create_adof_conns > File "/scratch/packages/sfepy-2014.2/sfepy/discrete/variables.py", line > 12, in <module> > import sfepy.linalg as la > File "/scratch/packages/sfepy-2014.2/sfepy/linalg/__init__.py", line 1, > in <module> > from extmods.crcm import rcm, permute_in_place > ImportError: No module named crcm I don't understant, because everything is installed in my /usr/local/lib/python2.7/ (I run an ubuntu 14.04) If I do an inplace build, everything works fine (but I really want to install SfePy system-wide...) I really don't know how to fix that... any idea ? Thank you very much ! G.

3 6

Working with several materials simultaneously
by Santiago Puerto Sept. 27, 2014

Sept. 27, 2014

Hello! Im Santiago Puerto, student of engineering from the Universidad Distrital Francisco José de Caldas - Bogotá - Colombia. This is my first Sfepy implementation, i wonder if a FEM mesh can contain multiple material ids (with different physical properties each one). My project is an asphalt mixtures aging simulator, i use X-Ray CAT images to reconstruct an asphalt mixture sample, then run a scaling and segmenting process in order to store the model in three dimensional numpy array. Now im working to transform the digital reconstruction on a FEM mesh but i don't want to lose the material detail (Identifying the material from the mixture, it's position and properties) gained with the X-Ray CAT images processing.

2 1

ANN: SfePy 2014.3
by Robert Cimrman Sept. 25, 2014

Sept. 25, 2014

I am pleased to announce release 2014.3 of SfePy. Description ----------- SfePy (simple finite elements in Python) is a software for solving systems of coupled partial differential equations by the finite element method or by the isogeometric analysis (preliminary support). It is distributed under the new BSD license. Home page: http://sfepy.org Mailing list: http://groups.google.com/group/sfepy-devel Git (source) repository, issue tracker, wiki: http://github.com/sfepy Highlights of this release -------------------------- - isogeometric analysis (IGA) speed-up by C implementation of NURBS basis evaluation - generalized linear combination boundary conditions that work between different fields/variables and support non-homogeneous periodic conditions - non-constant essential boundary conditions given by a function in IGA - reorganized and improved documentation For full release notes see http://docs.sfepy.org/doc/release_notes.html#id1 (rather long and technical). Best regards, Robert Cimrman and Contributors (*) (*) Contributors to this release (alphabetical order): Vladimir Lukes, Matyas Novak, Zhihua Ouyang, Jaroslav Vondrejc

1 0

Thanks!
by Daniel Wheeler Sept. 8, 2014

Sept. 8, 2014

I'd just like to say thanks for the help with the periodic displacement boundary conditions. This was the major stumbling block to our first release of PyMKS (http://openmaterials.github.io/pymks/). Sfepy is now the data generator for the learning algorithm in PyMKS. We did have other solver choices, but we decided on Sfepy because it is both scriptable from Python and can solve linear elasticity (and now the periodic displacement BC). However, as a third party requirement, it needs to be relatively easy to install, which was the most important factor -- Daniel Wheeler

2 1

3D displacement issues with order=3
by Daniel Wheeler Sept. 4, 2014

Sept. 4, 2014

I'm trying to extend the linear combination example to 3D. During the debugging, I noticed that I couldn't solve a simple displacement problem in 3D with "order=3" or "order=2" integrals. The gist contains the code I'm trying to run https://gist.github.com/wd15/8274429e3b4179c3d74f The gist only has displacements with no linear combination or periodic boundary conditions. The image produced seems to be badly distorted without the expected symmetry. Should this work with an "order=3" integral or does it need to be "order=1"? Does this require changing the solver parameters in some way? Thanks. -- Daniel Wheeler

2 3