SfePy March 2012

sfepy@python.org

5 participants
8 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

Modeling an adjustable X-ray optic
by Tom Aldcroft 21 May '13

21 May '13

I'm working on modeling a next-generation X-ray mirror for which the shape can be actively controlled by use of many thin piezo-electric actuators mounted on the mirror surface. The mirror is basically a glass conical paraboloid with a 1 meter radius and 200 micron thickness (e.g. http://en.wikipedia.org/wiki/X-ray_optics). Our project is currently using a proprietary FEA package, but the model setup and turnaround time is slow, in part because there is only one part-time engineer who can run it. SfePy looks like a great package and we're hoping that it could be used to automate running a large number of different cases. I've spent some time reading the documentation but I have a few questions that I hope can be answered before going too much further. I want to apologize in advance if some of my wording is imprecise, I have a physics background but this topic is a bit outside my realm... - Is SfePy appropriate for this problem? - If a specify a grid with about 800 x 400 points (azimuthal, axial) and about 10 boundary conditions (corresponding to mount points), what is the rough order of magnitude of time to compute the solution? Is it seconds, minutes, hours, or days? - The linear elastic examples show a problem with a specified displacement. How do I specify an input force? The piezo essentially provides a tensile force along the surface. - Is there a way to specify the problem and solve in cylindrical coordinates? This is the natural coordinate system. - How do I specify 6-DOF constraints which correspond to the mirror mounts? Thanks in advance for any help! Tom Aldcroft

3 6

time counting
by Md. Golam Rashed 06 Apr '12

06 Apr '12

is it possible to measure the total time spent in an analysis?

2 12

31 Mar '12

I want to write a paper themed on, doing simulation using scripting language with legacy code. As a civil engineer (structural), i find the tensile strength of concrete test presented in sfepy primer quite suitable to use in my writing. will this cause any copyright problem??? And secondly, can anybody detail the process of wrapping python around C in sfepy?? just to make sure my understanding is right. Thanks!

5 18

SfePy Documentation > HTML file
by Md. Golam Rashed 20 Mar '12

20 Mar '12

Hello, Can SfePy Documentations (./doc/sfepy.pdf) be made in html format so that it can be easily browse-able offline? I tried to download the http://sfepy.org/doc-devel/index.html site but errors ensued. Thanks.

2 1

benchmarks
by freevryheid 07 Mar '12

07 Mar '12

Claim your bragging rights. OS: Archlinux x64 CPU: Intel(R) Core(TM) i7-3820 CPU @ 3.60GHz RAM: 8Gb 63 test file(s) executed in 91.11 s, 0 failure(s) of 107 test(s)

4 5

quick questions
by David Mashburn 07 Mar '12

07 Mar '12

Hello again, I am performing elasticity simulations with boundary forces. I can tell sfepy is capable of this, and starting with the linear_elasticity.py example has already provided a close approximation of what I need, but I am not clear on a couple of things. My main two related questions right now are: 1. Can sfepy apply forces to boundaries instead of to volumes (effectively a von Neumann boundary condition)? What is the prescribed way to do this? I did notice that there is a surface term similar to the volume term, but I was not sure how to use it... 2. Also, what exactly are the "regions" in sfepy? Are they collections of individual mesh elements? If so, is there any way of specifying borders instead? Also, is there a good way way to visualize regions and/or get at the raw data that makes them up? (I am having trouble getting just the) Thanks, -David BTW, I was not as familiar with quadrature, so I made a small visualizer for the different quadrature tables (I know the formatting will probably get eaten, but you should still get the idea): import numpy as np import matplotlib.pyplot as plt from mayavi import mlab import sfepy.fem quadType='3_8' order=5 qp=sfepy.fem.quadratures.quadrature_tables[quadType][order] c=qp.coors w=qp.weights print c print w plt.cla() mlab.clf() if quadType=='1_2': plt.plot(qp.bounds,[0,0]) for i in range(len(c)): plt.plot(c[i],0,'k.',ms=w[i]*100) elif quadType=='2_3': plt.plot([0,1,0,0],[0,0,1,0]) for i in range(len(c)): plt.plot(c[i][0],c[i][1],'k.',ms=w[i]*100) elif quadType=='2_4': plt.plot([0,1,1,0,0],[0,0,1,1,0]) for i in range(len(c)): plt.plot(c[i][0],c[i][1],'k.',ms=w[i]*100) elif quadType=='3_4': mlab.plot3d([0,0,0], [1,0,0], [0,0,1], color=(0,0,1),tube_radius=0.025) mlab.plot3d([0,1,0,0,1], [0,0,1,0,0], [0,0,0,1,0], color=(0,0,1), tube_radius=0.025) for i in range(len(c)): #plt.plot(c[i][0],c[i][1],'k.',ms=w[i]*100) print [c[i][0]],[c[i][1]],[c[i][2]],w[i] mlab.points3d([c[i][0]],[c[i][1]],[c[i][2]],[w[i]],color=(1,1,1),scale_factor=2) elif quadType=='3_8': pts = [[i,j,k] for i in [0,1] for j in [0,1] for k in [0,1]] lines = [[i,j] for i in pts for j in pts if (pts.index(i)<pts.index(j) and sum(abs(np.array(i)-j))==1)] for l in lines: l=np.transpose(l) mlab.plot3d(l[0], l[1], l[2], color=(0,0,1),tube_radius=0.025) for i in range(len(c)): #plt.plot(c[i][0],c[i][1],'k.',ms=w[i]*100) print [c[i][0]],[c[i][1]],[c[i][2]],w[i] mlab.points3d([c[i][0]],[c[i][1]],[c[i][2]],[w[i]],color=(1,1,1),scale_factor=2)

2 6

[PATCH] Add archlinux installation requirements
by Andre Smit 05 Mar '12

05 Mar '12

--- doc/installation.rst | 14 ++++++++++++++ 1 file changed, 14 insertions(+) diff --git a/doc/installation.rst b/doc/installation.rst index 45db2bb..55c0f4f 100644 --- a/doc/installation.rst +++ b/doc/installation.rst @@ -119,6 +119,20 @@ Debian apt-get install python-tables python-pyparsing python-matplotlib python-scipy +Archlinux +^^^^^^^^^ + +:: + + pacman -S python2-numpy python2-scipy python2-matplotlib ipython2 python2-sympy + yaourt -S python-pytables python2-mayavi + +Instructions +"""""""""""" + +Edit Makefile and change all references from python to python2. +Edit scripts and change shebangs to point to python2. + Ubuntu ^^^^^^ -- 1.7.9.2

2 1

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

SfePy March 2012