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
Dear Robert
I could overcome my previous problems of not having root privileges, so it
all worked up till the command:
./postproc.py cylinder.vtk
This echoed the below message:
Exception
In
/home/gmh/.local/lib/python2.7/site-packages/mayavi-4.4.4.dev0-py2.7-linux-x86_64.egg/mayavi/sources/vtk_xml_file_reader.py:86
AttributeError: 'AlgorithmOutput' object has no attribute 'point_data' (in
get_all_attributes)
Traceback (most recent call last):
File
"/home/gmh/.local/lib/python2.7/site-packages/traits-4.6.0.dev354-py2.7-linux-x86_64.egg/traits/trait_notifiers.py",
line 340, in __call__
self.handler( *args )
File "/home/gmh/sfepy/sfepy/postprocess/dataset_manager.py", line 177, in
_dataset_changed
self._assign_attribute.input = value
File
"/home/gmh/.local/lib/python2.7/site-packages/traits-4.6.0.dev354-py2.7-linux-x86_64.egg/traits/trait_handlers.py",
line 104, in _read_only
name, class_of( object ) )
TraitError: The 'input' trait of an AssignAttribute instance is 'read only'.
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"
I'm using Python 2.7.8 if that's relevant to the current issue.
Thanks again
Cheers
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
Hello!
I'm trying to solve Poisson equation with Sfepy in a very simple domain,
but I found some problem:
1) I need triangular mesh on a rectangular/square domain
2) I have null von neumann conditions ( and I omitted them in boundary
condition) like in the
discussion https://groups.google.com/forum/#!searchin/sfepy-devel/null$20von$20neumann…
3)I know that this problem have many convergence problem, in one of this
case is the simulation successful? Or I find some error?
I don't know because the simulation is successful but in the postprocessing
the result isn't correct
I attach the problem description file too. If you want I would like some
advice about my problem
Thank you very much
-Francesca
I have made an inline instalation of sfepy-2013.3 on fedora22 , I run the
tests and I have only 3 failed, they maybe related to the parallel option
since I don't have instaled that option. I have python 2.7.10, mayavi
4.4.3, VTK 6.1.0, Qt 4.8.6
Either on the python or ipython window when I try to import:
*from sfepy.postprocess.viewer import Viewe*r
I get the following message,
*cannot set toolkit to wx because it has already been set to qt4*
and of course any postprocessing will fail
Hi Francesca,
On 12/10/2015 11:41 AM, Francesca Camorrino wrote:
> Hello!
> I'm trying to solve Poisson equation with Sfepy in a very simple domain,
> but I found some problem:
> 1) I need triangular mesh on a rectangular/square domain
For this, you can try gmsh. Why do you need triangles, by the way?
> 2) I have null von neumann conditions ( and I omitted them in boundary
> condition) like in the
> discussion https://groups.google.com/forum/#!searchin/sfepy-devel/null$20von$20neumann…
OK
> 3)I know that this problem have many convergence problem, in one of this
> case is the simulation successful? Or I find some error?
You fix the potential in some, so it should work fine.
But your definition for region 'p' selects nothing for the mesh you generate - try
'p' : ('vertices in (x <= 0.99951 )', 'vertex'),
instead.
Always verify your regions using --save-regions-as-groups option:
./simple.py prob_desc_file.py --save-regions-as-groups
- then check the 0_regions.vtk file using postproc.py.
> I don't know because the simulation is successful but in the postprocessing
> the result isn't correct
> I attach the problem description file too. If you want I would like some
> advice about my problem
When I run your code, it does not do anything, because there is no load - maybe
some right-hand side term is missing? Or provide a Dirichlet BC also in the
region 'n' (again, fix the region definition).
r.
PS: I missed your message, so my reply is late.
The example
http://sfepy.org/doc-devel/examples/navier_stokes/navier_stokes.html clearly
illustrated a NS problem with inlet velocity normal to y axis.
Is it possible to set the inlet with arbitrary plane that is not normal to
any axis? e.g. inlet normal velocity = (u1, u2, u3)
I see that in this example the inlet and outlet boundaries are defined by
the cinc functions:
cinc_name = 'cinc_' + op.splitext(op.basename(filename_mesh))[0]cinc = getattr(utils, cinc_name)
functions = {
'cinc0' : (lambda coors, domain=None: cinc(coors, 0),),
'cinc1' : (lambda coors, domain=None: cinc(coors, 1),),}
how would this utility work with my arbitrary inlet/outlet boundaries?