SfePy January 2015

sfepy@python.org

6 participants
9 discussions

Re: Torque
by Robert Cimrman 09 Mar '18

09 Mar '18

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

1D line elements in SfePy
by Nimish 22 Jan '16

22 Jan '16

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

Grabbing the mass matrix
by Sukhbinder Singh 04 Feb '15

04 Feb '15

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

2 7

Scipy was build without umfpack support
by Alicia White 29 Jan '15

29 Jan '15

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.

2 1

Problem installing sfepy with Ananconda
by Alicia White 28 Jan '15

28 Jan '15

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

3 4

Simple Stochastic Simulation
by Sukhbinder Singh 22 Jan '15

22 Jan '15

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.

2 3

pressure boundary condition in elasticity problems
by Constantine 16 Jan '15

16 Jan '15

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.

2 1

Mesh quality evaluation
by Jan Heczko 09 Jan '15

09 Jan '15

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

2 1

new examples
by Robert Cimrman 03 Jan '15

03 Jan '15

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.

1 0

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

SfePy January 2015