Mailman 3 April 2016 - SfePy

Material Parameters of Individual Nodes
by derj...＠web.de April 12, 2022

April 12, 2022

Hello! It's a great tool you built there! I'd like to use it for some kind of topology optimisation. Is there an "idiomatic" way to set the material parameters of individual nodes? Thank you very much! JonnyB

3 5

Re: Torque
by Robert Cimrman March 9, 2018

March 9, 2018

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

new interactive example with time stepping
by Robert Cimrman March 13, 2017

March 13, 2017

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

2 8

Need help to start with sfepy
by Nicolas DELAN May 13, 2016

May 13, 2016

Hi everyone, I am a new user of sfepy, i need it to compute the diffusions equations on a turbine blade. But i am a little lost with the boundary conditions. I have tried to impose some dirichlet boundary conditions on some points of the edges of the airfoil but it keeps crashing with exceptions. For example in my .mesh file i have this at the line where my point is defined 0.1 0.5 25 (x) (y) So I have written this in the region definition regions = { 'Omega' : 'all', 'point25' : ('vertices of group 25', 'edge'), } I have tried before with : refions = { 'Omega' : 'all', 'extrados' : ('vertices in (y>0.1) ', 'facet'), 'intrados' : ('vertices in (y<0)', 'facet'), } and it was giving resuls (it was only for test purpose) I don't understand very well how to define the boundary conditions (and to be honest the hole configuration file despite the documentation) So if somebody could help me, i would be grateful ! cordially, Nicolas DELAN

2 6

postprocessing
by Nina Halabi May 5, 2016

May 5, 2016

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

4 5

Region definition
by Dimitri Balasoiu April 25, 2016

April 25, 2016

Hello everyone, I have a question about sfepy's definition of regions. I want to define a region Omega_1, wich is a subregion of Omega. So I type : regions = { 'Omega' : 'all', 'Omega_1' : ('vertices of group 1', 'cell', Omega),} , where group 1 is a group of vertices. Question : wich field element P1 will be created over the region Omega_1 ? More precisely, I have a triangle (A,B,C), with node A in group 1, but the other ones outside of group 1. Will the finite element wich is worth 1 at A, and 0 at B and C be in my field element basis ? In wich file are the basis functions for the fe fields defined ? Dimitri

2 3

Re: [sfepy-devel] install sfepy vcvarshell error
by Robert Cimrman April 20, 2016

April 20, 2016

Hi Mark, yes this error is annoying, and I am not aware of a solution that actually works for sfepy (failed myself to install sfepy in a W10 virtual box). I assume you tried the solutions proposed for example in [1, 2, 3]? I do not think that I have tried everything, so in case there is something you did not try, could you give it a shot? r. [1] https://support.enthought.com/hc/en-us/articles/204469210-Windows-Unable-to… [2] https://stackoverflow.com/questions/24380442/getting-error-unable-to-find-v… [3] https://github.com/cython/cython/wiki/CythonExtensionsOnWindows On 04/07/2016 05:36 PM, 'Mark Shart' via sfepy-devel wrote: > Hi everyone, > > I try to install sfepy for window10 following the instruction with > python(xy). > I can accomplish everything till I need to compile sfepy. The I get the > error:Unable to find vcvarsall.bat. > I installed the MS C++ Compiler for Python 2.7 but its not working. I also > tried other C++ compiler... but still this error drives me crazy. > How can I fix this? Please help :\ >

2 4

issues #317, #318
by Robert Cimrman April 14, 2016

April 14, 2016

Hi, any windows users around that could help with the issues #317, #318? r. https://github.com/sfepy/sfepy/issues/317 https://github.com/sfepy/sfepy/issues/318

1 0

Install error Windows10
by Mark Shart April 7, 2016

April 7, 2016

Hi everyone, I'm trying to install sfepy on windows10 with python(xy) as suggested, but get the following error when compiling python setup.py build_ext --inplace --compiler=mingw32 $ python setup.py build_ext --inplace igakit (optional) is missing; you may get run-time errors petsc4py (optional) is missing; you may get run-time errors mpi4py (optional) is missing; you may get run-time errors pymetis (optional) is missing; you may get run-time errors pysparse (optional) is missing; 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.common.extmods._fmfield" sources building extension "sfepy.discrete.common.extmods.mappings" sources building extension "sfepy.discrete.common.extmods.assemble" sources building extension "sfepy.discrete.common.extmods.cmesh" sources building extension "sfepy.discrete.common.extmods.crefcoors" sources building extension "sfepy.discrete.common.extmods._geommech" sources building extension "sfepy.discrete.fem.extmods.bases" sources building extension "sfepy.discrete.fem.extmods.lobatto_bases" 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 customize Mingw32CCompiler customize Mingw32CCompiler using build_clib customize Mingw32CCompiler customize Mingw32CCompiler using build_ext extending extension 'sfepy.discrete.common.extmods._fmfield' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.common.extmods.mappings' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.common.extmods.crefcoors' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.common.extmods._geommech' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.fem.extmods.bases' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.fem.extmods.lobatto_bases' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.terms.extmods.terms' defined_macros with [('__SDIR__', '\'"sfepy\\terms\\extmods"\''), ('SFEPY_PLATFORM', 1), ('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\'')] building 'sfepy.discrete.common.extmods._fmfield' extension compiling C sources C compiler: gcc -O2 -Wall -Wstrict-prototypes compile options: '-D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -DNPY_MINGW_USE_CUSTOM_MSVCR -D__MSVCRT_VERSION__=0x0900 -Isfepy\discrete\common\extmods -IC:\Python27\lib\site-packages\numpy\core\include -IC:\Python27\include -IC:\Python27\PC -c' extra options: '-g -O2' gcc -O2 -Wall -Wstrict-prototypes -D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -DNPY_MINGW_USE_CUSTOM_MSVCR -D__MSVCRT_VERSION__=0x0900 -Isfepy\discrete\common\extmods -IC:\Python27\lib\site-packages\numpy\core\include -IC:\Python27\include -IC:\Python27\PC -c sfepy\discrete\common\extmods\_fmfield.c -o build\temp.win32-2.7\Release\sfepy\discrete\common\extmods\_fmfield.o -g -O2 Found executable C:\MinGW\bin\gcc.exe g++ -shared build\temp.win32-2.7\Release\sfepy\discrete\common\extmods\_fmfield.o -LC:\Python27\libs -LC:\Python27\PCbuild -Lbuild\temp.win32-2.7 -lsfepy_common -lpython27 -lmsvcr90 -o sfepy\discrete\common\extmods\_fmfield.pyd Found executable C:\MinGW\bin\g++.exe C:\Python27\libs/libpython27.a: error adding symbols: File format not recognized collect2.exe: error: ld returned 1 exit status error: Command "g++ -shared build\temp.win32-2.7\Release\sfepy\discrete\common\extmods\_fmfield.o -LC:\Python27\libs -LC:\Python27\PCbuild -Lbuild\temp.win32-2.7 -lsfepy_common -lpython27 -lmsvcr90 -o sfepy\discrete\common\extmods\_fmfield.pyd" failed with exit status 1 Can anyone help me please?

1 0

Re: [sfepy-devel] Obtain the required data from a *.vtk file?
by Robert Cimrman April 4, 2016

April 4, 2016

On 04/02/2016 02:50 PM, Kid Guo wrote: > Thank you for your reply. > Last time I did not have a good description of my problem,but I made some > changes on your script and got some useful results for me. > > Now I have a new question,I have got all points(such as,x=0.1) > and corresponding cells,if you run the modified script,you can see it, > there are 11 cells on the face where x=0.1, OK, so you want cells containing faces that are in x=0.1 plane, right? I have modified the script to find those - there are 16 of them - there was a mistake in your script, as your cells have 4 vertices, not only three. > I want to get the values(p and u) of the four points of each unit instead > of all cells. See the attachment. r. > Any ideas? > > Thank you advance. > > 在 2016年4月1日星期五 UTC+8下午11:07:34，Robert Cimrman写道： >> >> On 04/01/2016 05:08 AM, Kid Guo wrote: >>> I am trying to Obtain the required data from a *.vtk file. >>> For example, >>> 1. I want to output the all points where{x=0.1},such as (0.1,0.1,0.1) >>> 2. I want to obtain the CELLS including points in step1;such as >>> (12,13,56,23) >>> 3. I want to obtain the values(p,u) in the *.vtk file >> >> The attached script should show you how to do what you want. >> >>> < >> https://lh3.googleusercontent.com/-edlklG7xlrM/Vv3mDxDgzDI/AAAAAAAAACo/BajI…> >> >> >> This seems to correspond to an old version of sfepy - the script is w.r.t. >> the >> latest version. >> >>> I am using ipython,but I have no idea. >> >> Yeah, these low-level features are not extensively documented - one has to >> know >> what to look for. >> >> r. >> >>> Here is the attachment >>> >>> >>> What should I do? >>> >>> Thank you advance. >> >> >

1 0