Re: Torque
by Robert Cimrman
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.
1 year, 1 month
1D line elements in SfePy
by Nimish
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
3 years, 3 months
Grabbing the mass matrix
by Sukhbinder Singh
Hi this might be a silly question, but am not able to find a way to get the mass matrix from the problem in sfepy.
I can get the stiffness matrix K but how to get the Mass matrix M.
Any help is appreciated. Thanks
4 years, 2 months
Scipy was build without umfpack support
by Alicia White
In running the tests for sfepy, I had several fails when I ran the debug as
suggested in your installation I recieved the following message: (Full text
will be at the end of this email)
raise ImportError('Scipy was built without UMFPACK support. '
ImportError: Scipy was built without UMFPACK support. You need to install
the UMFPACK library and header files before building scipy.
I have tried pip installing scikit-umfpack but it told me that the
requirement was already met. I tried updating and it said it was up to
date. I also updated scipy but that did not change the error message. Do
you know of other problems that might be causing this error? I am using
Anaconda on a Windows 8 machine.
Thanks,
Alicia White
C:\Users\awhite40\Documents\GitHub\sfepy [master]> python run_tests.py
--debug tests/test_input_stokes.py
<<< directory: tests, test files: 1
<<< tests\test_input_stokes.py
sfepy: left over: ['verbose', '__builtins__', '_filename', '__file__',
'__package__', 'TestInput', 'input_name', '__name
__', 'output_name', '__doc__']
sfepy: left over: ['verbose', 'w2', 'data_dir', '_filename',
'match_y_line', '__package__', '__doc__', '__builtins__', '
save_format', '__file__', '__name__']
>>> test instance prepared (1 test(s))
... solving ../examples/navier_stokes/stokes.py...
sfepy: reading mesh [line2, tri3, quad4, tetra4, hexa8]
(C:\Users\awhite40\Documents\GitHub\sfepy\meshes\2d\special\chan
nels_symm944t.mesh)...
sfepy: ...done in 0.01 s
sfepy: creating regions...
sfepy: Right
sfepy: EBCGamma1
sfepy: EBCGamma2
sfepy: Bottom
sfepy: Top
sfepy: Y1
sfepy: Y2
sfepy: Left
sfepy: Walls
sfepy: Y1Y2
sfepy: ...done in 0.03 s
sfepy: equation "balance":
sfepy: dw_div_grad.i.Y1Y2( fluid.viscosity, v, u )
- dw_stokes.i.Y1Y2( v, p ) = 0
sfepy: equation "incompressibility":
sfepy: dw_stokes.i.Y1Y2( u, q ) = 0
sfepy: using solvers:
ts: no ts
nls: newton
ls: ls
sfepy: updating variables...
sfepy: ...done
C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\discrete\common\fields.py:196:
FutureWarning: comparison to `None` will r
esult in an elementwise object comparison in the future.
nods.remove(nn)
sfepy: setting up dof connectivities...
sfepy: ...done in 0.00 s
sfepy: matrix shape: (1845L, 1845L)
sfepy: assembling matrix graph...
sfepy: ...done in 0.00 s
sfepy: matrix structural nonzeros: 41800 (1.23e-02% fill)
sfepy: updating materials...
sfepy: fluid
sfepy: ...done in 0.00 s
sfepy: nls: iter: 0, residual: 4.566249e+00 (rel: 1.000000e+00)
>>> <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 "C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\base\testing.py",
line 58, in run
ret = test_method()
File "C:\Users\awhite40\Documents\GitHub\sfepy\tests\tests_basic.py",
line 73, in test_input
post_process_hook_final=self.post_process_hook_final)
File
"C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\applications\pde_solver_app.py",
line 55, in solve_pde
return app(nls_status=nls_status)
File
"C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\applications\application.py",
line 29, in call_basic
return self.call(**kwargs)
File
"C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\applications\pde_solver_app.py",
line 220, in call
nls_status=nls_status)
File
"C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\solvers\ts_solvers.py",
line 29, in __call__
state = problem.solve(state0=state0, nls_status=nls_status)
File
"C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\discrete\problem.py", line
984, in solve
vec = solvers.nls(vec0, status=self.nls_status)
File "C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\solvers\nls.py",
line 396, in __call__
eps_a=eps_a, eps_r=eps_r, mtx=mtx_a)
File "C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\solvers\ls.py", line
45, in _standard_call
**kwargs)
File "C:\Users\awhite40\Documents\GitHub\sfepy\sfepy\solvers\ls.py", line
129, in __call__
return self.sls.spsolve(mtx, rhs)
File
"C:\Users\awhite40\AppData\Local\Continuum\Anaconda\lib\site-packages\scipy\sparse\linalg\dsolve\linsolve.py",
li
ne 127, in spsolve
umf = umfpack.UmfpackContext(family[A.dtype.char])
File
"C:\Users\awhite40\AppData\Local\Continuum\Anaconda\lib\site-packages\scikits\umfpack\umfpack.py",
line 282, in _
_init__
raise ImportError('Scipy was built without UMFPACK support. '
ImportError: Scipy was built without UMFPACK support. You need to install
the UMFPACK library and header files before bu
ilding scipy.
4 years, 2 months
Problem installing sfepy with Ananconda
by Alicia White
Hi,
I am having trouble installing sfepy. I am using a Windows 8 machine with
Anaconda installed.
When I try to install I get a series of messages which I have pasted below.
I think the problem is in the part that I have bolded but I am not sure
how to correct that. Any assistance you could provide would be much
appreciated.
Thanks,
Alicia
Missing optional package "mayavi"; you may get run-time errors
running build_ext
running build_src
build_src
building library "sfepy_common" sources
building library "sfepy_terms" sources
building extension "sfepy.discrete.fem.extmods._fmfield" sources
building extension "sfepy.discrete.fem.extmods.mappings" sources
building extension "sfepy.discrete.fem.extmods.assemble" sources
building extension "sfepy.discrete.fem.extmods.bases" sources
building extension "sfepy.discrete.fem.extmods.cmesh" sources
building extension "sfepy.discrete.fem.extmods.lobatto_bases" sources
building extension "sfepy.discrete.fem.extmods._geommech" sources
building extension "sfepy.discrete.iga.extmods.igac" sources
building extension "sfepy.linalg.extmods.crcm" sources
building extension "sfepy.terms.extmods.terms" sources
building data_files sources
build_src: building npy-pkg config files
Looking for python27.dll
customize Mingw32CCompiler
customize Mingw32CCompiler using build_clib
Looking for python27.dll
customize Mingw32CCompiler
customize Mingw32CCompiler using build_ext
extending extension 'sfepy.discrete.fem.extmods._fmfield' defined_macros
with [('__SDIR__', '\'"sfepy\\discrete\\fem\\ex
tmods"\''), ('SFEPY_PLATFORM', 1)]
extending extension 'sfepy.discrete.fem.extmods.mappings' defined_macros
with [('__SDIR__', '\'"sfepy\\discrete\\fem\\ex
tmods"\''), ('SFEPY_PLATFORM', 1)]
extending extension 'sfepy.discrete.fem.extmods.bases' defined_macros with
[('__SDIR__', '\'"sfepy\\discrete\\fem\\extmo
ds"\''), ('SFEPY_PLATFORM', 1)]
extending extension 'sfepy.discrete.fem.extmods.lobatto_bases'
defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\fe
m\\extmods"\''), ('SFEPY_PLATFORM', 1)]
extending extension 'sfepy.discrete.fem.extmods._geommech' defined_macros
with [('__SDIR__', '\'"sfepy\\discrete\\fem\\e
xtmods"\''), ('SFEPY_PLATFORM', 1)]
extending extension 'sfepy.terms.extmods.terms' defined_macros with
[('__SDIR__', '\'"sfepy\\terms\\extmods"\''), ('SFEP
Y_PLATFORM', 1), ('__SDIR__', '\'"sfepy\\discrete\\fem\\extmods"\'')]
building 'sfepy.discrete.fem.extmods._fmfield' extension
compiling C sources
C compiler: gcc -m64 -g -DNDEBUG -DMS_WIN64 -O2 -Wall -Wstrict-prototypes
compile options: '-D__SDIR__='"sfepy\discrete\fem\extmods"'
-DSFEPY_PLATFORM=1 -D__SDIR__='"sfepy\discrete\fem\extmods"'
-DSFEPY_PLATFORM=1 -D__MSVCRT_VERSION__=0x0900
-Isfepy\discrete\fem\extmods -IC:\Users\awhite40\AppData\Local\Continuum
\Anaconda\lib\site-packages\numpy\core\include
-IC:\Users\awhite40\AppData\Local\Continuum\Anaconda\include -IC:\Users\a
white40\AppData\Local\Continuum\Anaconda\PC -c'
extra options: '-g -O2'
gcc -m64 -g -DNDEBUG -DMS_WIN64 -O2 -Wall -Wstrict-prototypes
-D__SDIR__='"sfepy\discrete\fem\extmods"' -DSFEPY_PLATFORM
=1 -D__SDIR__='"sfepy\discrete\fem\extmods"' -DSFEPY_PLATFORM=1
-D__MSVCRT_VERSION__=0x0900 -Isfepy\discrete\fem\extmods
-IC:\Users\awhite40\AppData\Local\Continuum\Anaconda\lib\site-packages\numpy\core\include
-IC:\Users\awhite40\AppData\L
ocal\Continuum\Anaconda\include
-IC:\Users\awhite40\AppData\Local\Continuum\Anaconda\PC -c
sfepy\discrete\fem\extmods\_f
mfield.c -o
build\temp.win-amd64-2.7\Release\sfepy\discrete\fem\extmods\_fmfield.o -g
-O2
Found executable
C:\Users\awhite40\AppData\Local\Continuum\Anaconda\Scripts\gcc.bat
gcc -m64 -g -shared
build\temp.win-amd64-2.7\Release\sfepy\discrete\fem\extmods\_fmfield.o
-LC:\Users\awhite40\AppData\L
ocal\Continuum\Anaconda\libs -*LC:\Users\awhite40\AppData\Local\Continuum\Anaconda\PCbuild\amd64
-Lbuild\temp.win-amd64-2*
*.7 -lsfepy_common -lpython27 -lmsvcr90 -o
sfepy\discrete\fem\extmods\_fmfield.pyd*
*build\temp.win-amd64-2.7/libsfepy_common.a: could not read symbols:
Archive has no index; run ranlib to add one*
collect2.exe: error: ld returned 1 exit status
error: Command "gcc -m64 -g -shared
build\temp.win-amd64-2.7\Release\sfepy\discrete\fem\extmods\_fmfield.o
-LC:\Users\aw
hite40\AppData\Local\Continuum\Anaconda\libs
-LC:\Users\awhite40\AppData\Local\Continuum\Anaconda\PCbuild\amd64 -Lbuild\
temp.win-amd64-2.7 -lsfepy_common -lpython27 -lmsvcr90 -o
sfepy\discrete\fem\extmods\_fmfield.pyd" failed with exit stat
us 1
4 years, 2 months
Simple Stochastic Simulation
by Sukhbinder Singh
Hi,
Has anyone tried any stochastic simulation with sfepy. I would like to
carry out material uncertainty stochastic simulation with sfepy. Any
pointers are welcome.
Thank you.
4 years, 3 months
pressure boundary condition in elasticity problems
by Constantine
Hi, this is probably a silly question but I don't seem to be able to find
an answer.
I want to calculate stresses within a loaded block, like similar
to http://sfepy.org/doc-devel/primer.html
but instead of force on a node I want to be able to specify pressure on the
whole boundary in Pa.
It seems to be possible to do by specifying "load" in the material
properties but it seems to be the force in N instead of pressure. Is there
some other way or should i calculate the force per node and then use the
calculated value in "load" ?
Thanks in advance.
4 years, 3 months
Mesh quality evaluation
by Jan Heczko
Hi everyone,
is there a tool for mesh quality evaluation in sfepy?
I.e. a function that quantifies wether some elements are too distorted.
JH
4 years, 3 months
new examples
by Robert Cimrman
FYI: Two new examples were added to the git repository:
- modal analysis of a linear elastic block in 2D or 3D:
python examples/standalone/interactive/modal_analysis.py -h
- support for fields with different approximation order was added to IGA, as
demonstrated in:
python simple.py examples/navier_stokes/navier_stokes2d_iga.py
Happy new year!
r.
4 years, 3 months
