Mailman 3 September 2015 - SfePy

Material Parameters of Individual Nodes
by derj...＠web.de April 12, 2022

April 12, 2022

Hello! It's a great tool you built there! I'd like to use it for some kind of topology optimisation. Is there an "idiomatic" way to set the material parameters of individual nodes? Thank you very much! JonnyB

3 5

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

new interactive example with time stepping
by Robert Cimrman March 13, 2017

March 13, 2017

Hi, I have just updated the time stepping solvers in sfepy for interactive use, as demonstrated in the new example [1]. For basic use, ignore the probing code - the time stepper can be used as simply as: tss = SimpleTimeSteppingSolver({'t0' : 0.0, 't1' : 100.0, 'n_step' : 11}, problem=problem) tss.init_time() for step, time, state in tss(): pass r. [1] http://sfepy.org/doc-devel/examples/diffusion/time_poisson_interactive.html

2 8

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

ANN: SfePy 2015.3
by Robert Cimrman Sept. 23, 2015

Sept. 23, 2015

I am pleased to announce release 2015.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 -------------------------- - preliminary support for parallel computing - unified evaluation of basis functions (= isogeometric analysis fields can be evaluated in arbitrary points) - (mostly) fixed finding of reference element coordinates of physical points - several new or improved examples For full release notes see http://docs.sfepy.org/doc/release_notes.html#id1 (rather long and technical). Best regards, Robert Cimrman on behalf of the SfePy development team --- Contributors to this release in alphabetical order: Robert Cimrman Vladimir Lukes

1 0

how to display "cylinder.h5", the result of material_nonlinearity.py
by pyontaku14 Sept. 23, 2015

Sept. 23, 2015

Hi, Could you please tell me how to display "cylinder.h5", the result of material_nonlinearity.py? Usually, the result file is vtk. But the result of material_nonlinearity.py seems only h5. The OS I use is Lubuntu15.04-64bit on VM Ware player. Best regards, Takuo Fujita

2 5

constrained nodal force output
by pyontaku14 Sept. 22, 2015

Sept. 22, 2015

Hi, Is there any way to output constrained nodal force of "examples/linear_elasticity/material_nonlinearity.py"? I would like to confirm material nonlinearity with nodal force. Any help or clarification would be greatly appreciated. Best regards, Takuo Fujita

2 1

coexistence between mayavi2 and paraview
by pyontaku14 Sept. 11, 2015

Sept. 11, 2015

Hi, Is it possible coexistence between mayavi2 and paraview? I installed paraview after installation of sfepy. mayavi2 was deleted automatically. ./postproc.py does not work due to it. Thank you. Takuo Fujita

2 4

navier_stokes2d_iga.py failed
by Sebastian Bascuñan Sept. 10, 2015

Sept. 10, 2015

Hi, I installed sfepy 2015.2 on ubuntu 14.04 lts, and ran the tests. The only one that failed was the one mentioned on the subject. This is what I got after running the debug. <<< directory: tests, test files: 1 <<< tests/test_input_navier_stokes2d_iga.py sfepy: left over: ['verbose', '__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] sfepy: left over: ['data_dir', 'verbose', '_filename', '__builtins__', '__package__', '__name__', '__doc__', '__file__'] >>> test instance prepared (1 test(s)) ... solving ../examples/navier_stokes/navier_stokes2d_iga.py... sfepy: creating regions... sfepy: Right sfepy: Bottom sfepy: Top sfepy: Omega sfepy: Left sfepy: Walls sfepy: ...done in 0.00 s >>> <type 'exceptions.ImportError'> Traceback (most recent call last): File "run_tests.py", line 206, in <module> main() File "run_tests.py", line 197, in main run_tests(stats, dirname, [filename]) File "run_tests.py", line 131, in run_tests n_fail, n_total, test_time = run_test(conf_name, options) File "run_tests.py", line 99, in run_test ok, n_fail, n_total = test.run(options.debug) File "/usr/local/src/sfepy-2015.2/sfepy/base/testing.py", line 58, in run ret = test_method() File "/usr/local/src/sfepy-2015.2/tests/tests_basic.py", line 73, in test_input post_process_hook_final=self.post_process_hook_final) File "/usr/local/src/sfepy-2015.2/sfepy/applications/pde_solver_app.py", line 50, in solve_pde app = PDESolverApp(conf, options, output_prefix) File "/usr/local/src/sfepy-2015.2/sfepy/applications/pde_solver_app.py", line 136, in __init__ self.problem = Problem.from_conf(conf, init_equations=is_eqs, **kwargs) File "/usr/local/src/sfepy-2015.2/sfepy/discrete/problem.py", line 96, in from_conf obj.set_fields(conf.fields) File "/usr/local/src/sfepy-2015.2/sfepy/discrete/problem.py", line 344, in set_fields self.fields = fields_from_conf(conf_fields, self.domain.regions) File "/usr/local/src/sfepy-2015.2/sfepy/discrete/common/fields.py", line 69, in fields_from_conf field = Field.from_conf(val, regions) File "/usr/local/src/sfepy-2015.2/sfepy/discrete/common/fields.py", line 153, in from_conf approx_order=approx_order[:2]) File "/usr/local/src/sfepy-2015.2/sfepy/discrete/iga/fields.py", line 72, in __init__ self.nurbs = self.domain.nurbs.elevate(elevate_times) File "/usr/local/src/sfepy-2015.2/sfepy/discrete/iga/domain.py", line 105, in elevate aux = self._to_igakit() File "/usr/local/src/sfepy-2015.2/sfepy/discrete/iga/domain.py", line 71, in _to_igakit import igakit.cad as cad ImportError: No module named igakit.cad Any help or clarification would be greatly appreciated. Regards, Sebastian

2 1

Postproc issues.
by Arthur Campos Sept. 10, 2015

Sept. 10, 2015

Hello, I am new to sfepy and as I was going over the examples I had issues running postproc.py and the error is posted below [arthur@arthurpc sfepy]$ ./postproc.py cylinder.vtk (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window sfepy: point scalars t at [-0.05 -0.02 0. ] sfepy: range: -2.00e+00 2.00e+00 l2 norm range: 1.06e-02 2.00e+00 No handlers could be found for logger "mayavi.core.common" (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window (postproc.py:3360): Gdk-WARNING **: gdkdrawable-x11.c:952 drawable is not a pixmap or window Exception occurred in traits notification handler. Please check the log file for details. Exception occurred in traits notification handler for object: <mayavi.tools.mlab_scene_model.MlabSceneModel object at 0x7effbfbdb3b0>, trait: activated, old value: <undefined>, new value: True Traceback (most recent call last): File "/home/arthur/anaconda/lib/python2.7/site-packages/traits/trait_notifiers.py", line 520, in _dispatch_change_event self.dispatch( handler, *args ) File "/home/arthur/anaconda/lib/python2.7/site-packages/traits/trait_notifiers.py", line 483, in dispatch handler( *args ) File "/home/arthur/sfepy/sfepy/postprocess/viewer.py", line 1217, in _post_init viewer.render_scene(self.scene, viewer.options) File "/home/arthur/sfepy/sfepy/postprocess/viewer.py", line 687, in render_scene self.build_mlab_pipeline(**options) File "/home/arthur/sfepy/sfepy/postprocess/viewer.py", line 527, in build_mlab_pipeline active = mlab.pipeline.set_active_attribute(source) File "/home/arthur/anaconda/lib/python2.7/site-packages/mayavi/tools/pipe_base.py", line 38, in the_function factory = factory_class(*args, **kwargs) File "/home/arthur/anaconda/lib/python2.7/site-packages/mayavi/tools/pipe_base.py", line 163, in __init__ self.set(**traits) File "/home/arthur/anaconda/lib/python2.7/site-packages/mayavi/tools/pipe_base.py", line 180, in set self._anytrait_changed(trait, value) File "/home/arthur/anaconda/lib/python2.7/site-packages/mayavi/tools/pipe_base.py", line 198, in _anytrait_changed setattr(obj, components[-1], value) File "/home/arthur/anaconda/lib/python2.7/site-packages/mayavi/core/trait_defs.py", line 101, in set_value value ) TraitError: The 'point_vectors_name' trait of a SetActiveAttribute instance must be one of [], but a value of '' <type 'str'> was specified. ./src/common/module.cpp(194): assert "module->m_state == State_Initialized" failed in DoCleanUpModules(): not initialized module being cleaned up Any ideas on how I can get it working? Appreciate your help, Arthur.

2 1