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?
----- Reply message -----
From: "Andre Smit" <freev...(a)gmail.com>
Date: Sat, Dec 18, 2010 05:10
What is the best way to apply a torque load to a model?
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.
I am currrently looking for FEM packages to help me solve a system of
beams and columns, basically a collection of 1D bernoulli/timoshenko
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..
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
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,
> 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
(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 !
FYI: As SciPy 0.12.0 is out and one of the release highlights is "Support for
Python 2 and Python 3 from the same code base (no more 2to3)", we can think
seriously about updating SfePy in this respect as well, cf. .
I'm new to sfepy and having a good experience so far. But I've been having
some issues understanding the d_surface_flux function.. a couple of
questions that might help:
1. In the poisson_neumann example, if I add a post_process function which
evaluates various surface fluxes the results are not what I would expect.
Here is my post_process function:
def post_process(out, pb, state, extend=False):
from sfepy.base.base import Struct
flx = pb.evaluate('d_surface_flux.2.Gamma(coef.sigma, t)',
print "Flux, entire surface:", np.sum(flx), flx.max(), flx.min()
snkf = pb.evaluate('d_surface_flux.2.Gamma_Left(coef.sigma, t)',
print "Sink flux: ", snkf
srcf = pb.evaluate('d_surface_flux.2.Gamma_Right(coef.sigma, t)',
print "Source flux: ", srcf
blkf = pb.evaluate('d_surface_flux.2.Gamma_N(coef.sigma, t)',
print "Block flux: ", blkf
I would expect that "Flux, entire surface" would be approximately zero, and
the sink flux and source flux should be approximately equal but opposite
sign, and block flux to be zero. Since Gamma_N is supposed to have a zero
neumann condition on it my understanding is that flux would only be present
over Gamma_Left and Gamma_Right and that they would be equal but opposite
sign. What am I missing here?
2. If I compute surface_flux using mode='el_avg' and then sum the results,
why is that number not equal to surface_flux computed with mode='eval'?
Whats the difference?
Thanks for your help,
I followed the interactive example in sfepy tutorial and wrote my custom
simulation script for the block domain with body force (t2 term). That
works. But for the point loaded block model (t3 term), I can't make it
work. The script stopped at "assert_(mat.shape[-1]==virtual.dim)" of
terms_point.py with IndexError: tuple index out of range. Searching the
web does not find help information for this.
The reason I write the attached custom simulation script is that I failed
doing the things below:
the "non-custom" script for a point loaded block running with "simple.py"
works and gmsh can mesh different geometries. I would like to finish those
two things together in python like "os.system('gmsh ........')", and
"os.system('simple.py examples/myscript.py") but errors to import config
in simple.py got me there since "path" settings for import in simple.py
seems a little complicated for me. The custom simulation script might be
easier to setup import "path" for me.
Thanks for your help.
For run_tests.py for the latest sfepy in window 8, here is what I do:
1.1 Upgrade python(x,y) from 22.214.171.124 to 126.96.36.199 (packages installed are
listed in the attached file using "yolk -l".)
1.2 Reinstall msysgit:
Since installation of 1.9.4 "NETINSTALL" verion of msysGit failed
(msysGit-netinstall-1.9.4-preview20140611.exe), I installed 1.9.0
"FULLINSTALL" version of msysGit
2.1 Old version of sfepy:
2.1.1 Start msysGit shell with msys.bat
2.1.2 cd to sfepy root directory
2.1.3 python setup.py build_ext --inplace --compiler=mingw32
(I ignored the error info: "msvcr90.dll not found").
2.1.5 This time, all tests are ok except stokes_slip_bc.py which causes
crash. It was ok before with python(x,y)-188.8.131.52 and previous msysgit
2.2 Latest verstion of sfepy:
2.2.1 python setup.py clean
2.2.2 python setup.py build_ext --inplace --compiler=mingw32
2.2.5 No luck again: 79 test file(s) executed in 10.12 s, 63 failure(s)
of 116 test(s)
I am using SfePy and I would like to know if there is a way to get the time
of simulation to have a time dependant boundary condition (if the time of
simulation is stocked in a variable for instance).
I saw that we could manually specify when to change the boundary conditions
at a specific time, I think i's possible to do something like that, but it
would be really more convenient for me directly to get the time variable if
it exists, but I couldn't find it on the examples.
Can I find this variable somewhere?
Thank you very much!
I tried to regenerate latex version of original sfepy documentation
and got encoding error: "ascii" codec can't encode characters.....at
location of src/sfepy/base/parse_conf (see the attached snapshots
below). I checked the files in sfepy/sfepy/base
and sfepy/doc/src/sfepy/base and could not find any non-ascii characters
OS: window 8, make: msysgit/bin/make, sphinx v1.2.2
If it is possible to generate sfepy doc on my machine, it is good for
me to add my sfepy-learning notes into it for myself.
I have attempted to install SfePy and was able to use Git Bash in Windows 7
Pro 32-bit and followed the instructions listed under "
Steps to get a working *SfePy* on Windows using Python(x,y)
Everything seemed nominal but upon running run_tests.py the program replied
with "79 test file(s) executed in 31.40 s, 63 failure(s) of 116 test(s)".
For one of the tests, test_volume.py, I received the output, "test
instance creation failed".
I attempted to run test_volume.py in debug mode as suggested on the webpage
above, but I receive a number of "Import by filename is not supported"
I am not sure if I have explained what I am experiencing enough, but on the
basis of what i have provided, is there a way to solve this issue?
Additionally, is more information necessary to assist me in
troubleshooting this problem?