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.
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
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
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
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
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
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
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
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.
Hello,
I am new to sfepy, so as I as going over the first example a bumped into
this problem while runing ./postproc.py its2D.vtk
[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 of how can I get this fixed?
Appreciate your help!
Arthur.
