SfePy November 2014

sfepy@python.org

5 participants
7 discussions

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.

3 years, 4 months

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

5 years, 6 months

ANN: SfePy 2014.4

by Robert Cimrman

I am pleased to announce release 2014.4 of SfePy. Description ----------- SfePy (simple finite elements in Python) is a software for solving systems of coupled partial differential equations by the finite element method or by the isogeometric analysis (preliminary support). It is distributed under the new BSD license. Home page: http://sfepy.org Mailing list: http://groups.google.com/group/sfepy-devel Git (source) repository, issue tracker, wiki: http://github.com/sfepy Highlights of this release -------------------------- - preliminary support for 1D problems - data probes using pyVTK library For full release notes see http://docs.sfepy.org/doc/release_notes.html#id1 (rather long and technical). Best regards, Robert Cimrman and Contributors (*) (*) Contributors to this release (alphabetical order): Lubos Kejzlar, Vladimir Lukes

6 years, 8 months

"MemoryError" when using large mesh and how does integral order change prestress shape?

by Ronghai Wu

Dear sfepy users and developers, I am using sfepy for solving an 2D inclusion elastic problem. So, I implemented based on "linear_elasticity/prestress_fibres" example in interactive way. However, the prestress in my problem is not constant, so I manually assign values. I have two questions: 1. why when the integral order is 3, the array shape of prestress.datas[( 'Omega', 'i')][0]['val'] is (nx*ny,4,3,1), so when I assign by prestress_value, I have to make it the same shape. If change integral order to 1, the shape changes to (nx*ny,8,3,1), why? 2. When I use mesh smaller than 70*70, no problem. When 90*90 very slow. When 120*120 "MemoryError" comes out. Thanks a lot. My codes and errors are as follow: field = Field.from_args('displacement', np.float64, 'vector', omega, 1) u = FieldVariable('u', 'unknown', field) v = FieldVariable('v', 'test', field, primary_var_name='u') prestress_value = np.zeros((nx*ny,3,1)) m = Material('m', lam=lambda_, mu=mu) prestress = Material('prestress', val=prestress_value) integral = Integral('i', order=3) t1 = Term.new('dw_lin_elastic_iso(m.lam, m.mu, v, u )', integral, omega, m=m , v=v, u=u) t2 = Term.new('dw_lin_prestress(prestress.val, v )',integral, omega, prestress=prestress, v=v) eq = Equation('balance', t1-t2 ) eqs = Equations([eq]) ls = ScipyDirect({}) nls_status = IndexedStruct() nls = Newton({'i_max' : 1,'eps_a' : 1e-6,'problem' : 'nonlinear'}, lin_solver=ls, status=nls_status) pb = Problem('elasticity', equations=eqs, nls=nls, ls=ls) pb.update_materials() prestress_value = np.zeros((nx*ny,4,3,1)) ####################################################################### this is how I assign values to prestress, eta is a known scalar field. ####################################################################### for i in range(nx*ny): for j in range(4): prestress_value[i][j] = np.array([[stress00_11], [stress00_22], [ stress00_12]]) * eta.value[i] prestress.datas[('Omega', 'i')][0]['val'] = prestress_value pb.update_materials() ######################################################################## following is error: sfepy: updating materials... sfepy: m sfepy: prestress Traceback (most recent call last): File "<stdin>", line 1, in <module> File "/usr/local/lib/python2.7/dist-packages/spyderlib/widgets/externalshell/sitecustomize.py" , line 580, in runfile execfile(filename, namespace) File "/home/ronghaiwu/temporary/test.py", line 304, in <module> main() File "/home/ronghaiwu/temporary/test.py", line 179, in main pb.update_materials() File "/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/problem.py" , line 522, in update_materials problem=self, verbose=verbose) File "/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/equations.py" , line 313, in time_update_materials verbose=verbose) File "/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/materials.py" , line 70, in time_update mat.time_update(ts, equations, mode=mode, problem=problem) File "/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/materials.py" , line 339, in time_update self.update_data(key, ts, equations, term, problem=problem) File "/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/materials.py" , line 251, in update_data **self.extra_args) File "/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/functions.py" , line 34, in __call__ return self.function(*args, **_kwargs) File "/home/ronghaiwu/mypylibs/lib/python2.7/site-packages/sfepy/discrete/functions.py" , line 66, in get_constants out[key] = nm.tile(val, (coors.shape[0], 1, 1)) File "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py", line 860, in tile c = c.reshape(-1, n).repeat(nrep, 0) MemoryError Regards Ronghai

6 years, 8 months

Problem with a hexahedral mesh

by patricia...＠upf.edu

Hi! I am trying to simulate a hyperelastic 3D cylinder with hexahedral elements. I did the mesh with gmsh and save it as a .mesh file. When I do it with tetrahedral elments everything works perfect, but when I try to simulate the mesh with hexahredal elemnts, SfePy can not read well the mesh and give me this error: sfepy: warning: bad element orientation, trying to correct... sfepy: warning: bad element orientation, trying to correct... RuntimeError: elements cannot be oriented! (0, 3_4) I've attached the .mesh file with hexahedral elements Thanks a lot for the help!

6 years, 8 months

Installation test failed

by Serge-Étienne Parent

Hi, I tried to install SfePy on my openSUSE 13.2 64bit machine. I had a hard time installing PyTables, but I could handle this by installing it using /usr/bin/pip2.7 install tables ...after installing python-pip, hdf5-devel, python-devel and python-numpy-devel. Everything else went smooth and SfePy compiled successfully. But most tests failed and importing in python did not work. *Importing in Python:* Python 2.7.8 (default, Sep 30 2014, 15:34:38) [GCC] on linux2 Type "help", "copyright", "credits" or "license" for more information. >>> from sfepy.base.base import * Traceback (most recent call last): File "<stdin>", line 1, in <module> ImportError: No module named sfepy.base.base *Test:* useraccount@linux:~/Programmes/sfepy> ./run_tests.py <<< directory: tests, test files: 82 tests/test_assembling.py --- test_assemble_matrix: failed! --- test_assemble_matrix_complex: failed! --- test_assemble_vector: failed! --- test_assemble_vector_complex: failed! tests/test_base.py +++ test_resolve_deps: ok +++ test_struct_add: ok +++ test_struct_i_add: ok +++ test_verbose_output: ok tests/test_cmesh.py --- test_cmesh_counts: failed! tests/test_domain.py --- test instance creation failed tests/test_ebcs.py --- test instance creation failed tests/test_elasticity_small_strain.py --- test_get_solution: failed! --- test_linear_terms: failed! tests/test_fem.py --- test instance creation failed tests/test_functions.py --- test instance creation failed tests/test_high_level.py --- test instance creation failed tests/test_homogenization_perfusion.py --- test_solution: failed! tests/test_hyperelastic_tlul.py --- test_solution: failed! tests/test_input_acoustics.py --- test instance creation failed tests/test_input_acoustics3d.py --- test instance creation failed tests/test_input_active_fibres.py --- test instance creation failed tests/test_input_biot.py --- test instance creation failed tests/test_input_biot_npbc.py --- test instance creation failed tests/test_input_biot_npbc_lagrange.py --- test instance creation failed tests/test_input_cube.py --- test instance creation failed tests/test_input_elastic_contact_planes.py --- test instance creation failed tests/test_input_elastic_contact_sphere.py --- test instance creation failed tests/test_input_elastic_shifted_periodic.py --- test instance creation failed tests/test_input_its2D_2.py --- test instance creation failed tests/test_input_laplace_coupling_lcbcs.py --- test instance creation failed tests/test_input_laplace_time_ebcs.py --- test instance creation failed tests/test_input_linear_elastic.py --- test instance creation failed tests/test_input_linear_elastic_damping.py --- test instance creation failed tests/test_input_linear_elastic_mM.py --- test instance creation failed tests/test_input_linear_elastic_probes.py --- test instance creation failed tests/test_input_linear_elastic_tractions.py --- test instance creation failed tests/test_input_linear_elastic_up.py --- test instance creation failed tests/test_input_linear_viscoelastic.py --- test instance creation failed tests/test_input_material_nonlinearity.py --- test instance creation failed tests/test_input_navier_stokes.py --- test instance creation failed tests/test_input_navier_stokes2d.py --- test instance creation failed tests/test_input_perfusion_tl.py --- test instance creation failed tests/test_input_piezo.py --- test instance creation failed tests/test_input_poisson.py --- test instance creation failed tests/test_input_poisson_field_dependent_material.py --- test instance creation failed tests/test_input_poisson_functions.py --- test instance creation failed tests/test_input_poisson_neumann.py --- test instance creation failed tests/test_input_poisson_periodic_boundary_condition.py --- test instance creation failed tests/test_input_prestress_fibres.py --- test instance creation failed tests/test_input_sinbc.py --- test instance creation failed tests/test_input_stabilized_navier_stokes.py --- test instance creation failed tests/test_input_stokes.py --- test instance creation failed tests/test_input_stokes_slip_bc.py --- test instance creation failed tests/test_input_thermo_elasticity.py --- test instance creation failed tests/test_input_thermo_elasticity_ess.py --- test instance creation failed tests/test_input_time_poisson.py --- test instance creation failed tests/test_input_vibro_acoustic3d.py --- test instance creation failed tests/test_io.py +++ test_recursive_dict_hdf5: ok +++ test_sparse_matrix_hdf5: ok tests/test_laplace_unit_disk.py --- test instance creation failed tests/test_laplace_unit_square.py --- test instance creation failed tests/test_lcbc_2d.py --- test_linear_rigid_body_bc: failed! tests/test_lcbc_3d.py --- test_linear_rigid_body_bc: failed! tests/test_lcbcs.py --- test_laplace_shifted_periodic: failed! tests/test_linalg.py --- test_assemble1d: failed! --- test_geometry: failed! --- test_tensors: failed! --- test_unique_rows: failed! tests/test_linear_solvers.py --- test instance creation failed tests/test_linearization.py --- test_linearization: failed! tests/test_matcoefs.py +++ test_conversion_functions: ok +++ test_elastic_constants: ok +++ test_stiffness_tensors: ok tests/test_mesh_generators.py --- test_gen_block_mesh: failed! --- test_gen_cylinder_mesh: failed! --- test_gen_extended_block_mesh: failed! tests/test_mesh_interp.py --- test_interpolation: failed! --- test_interpolation_two_meshes: failed! --- test_invariance: failed! --- test_invariance_qp: failed! tests/test_mesh_smoothing.py --- test_mesh_smoothing: failed! tests/test_meshio.py --- test instance creation failed tests/test_msm_laplace.py --- test instance creation failed tests/test_msm_symbolic.py --- test instance creation failed tests/test_normals.py --- test_normals: failed! tests/test_parsing.py --- test_parse_equations: failed! --- test_parse_regions: failed! tests/test_periodic_bc_2d.py --- test instance creation failed tests/test_periodic_bc_3d.py --- test instance creation failed tests/test_permutations.py --- test_rcm: failed! tests/test_poly_spaces.py --- test instance creation failed tests/test_projections.py --- test instance creation failed tests/test_quadratures.py --- test_quadratures: failed! --- test_weight_consistency: failed! tests/test_regions.py --- test instance creation failed tests/test_semismooth_newton.py --- test_semismooth_newton: failed! tests/test_sparse.py --- test_compose_sparse: failed! tests/test_tensors.py --- test_stress_transform: failed! --- test_tensors: failed! --- test_transform_data: failed! --- test_transform_data4: failed! tests/test_term_call_modes.py --- test instance creation failed tests/test_term_consistency.py --- test instance creation failed tests/test_units.py +++ test_consistent_sets: ok +++ test_units: ok tests/test_volume.py --- test instance creation failed 82 test file(s) executed in 0.46 s, 95 failure(s) of 106 test(s) *Test debug for test_term_consistency.py:* useraccount@linux:~/Programmes/sfepy> python run_tests.py --debug tests/test_term_consistency.py <<< directory: tests, test files: 1 <<< tests/test_term_consistency.py sfepy: left over: ['verbose', '__builtins__', '_integrate', 'test_terms', '__doc__', '__name__', 'get_pars', 'TestCommon', 'dim', 'data_dir', 'nm', '__package__', '_filename', '__file__'] --- test instance creation failed 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 83, in run_test test = conf.funmod.Test.from_conf(conf, options) File "/home/sergeetienne/Programmes/sfepy/tests/test_term_consistency.py", line 94, in from_conf from sfepy.discrete import Problem File "/home/sergeetienne/Programmes/sfepy/sfepy/discrete/__init__.py", line 8, in <module> from variables import Variables, Variable, FieldVariable, create_adof_conns File "/home/sergeetienne/Programmes/sfepy/sfepy/discrete/variables.py", line 12, in <module> import sfepy.linalg as la File "/home/sergeetienne/Programmes/sfepy/sfepy/linalg/__init__.py", line 1, in <module> from extmods.crcm import rcm, permute_in_place ImportError: No module named crcm

6 years, 8 months

Re: [sfepy-devel] How to implement non-constant prestress problem?

by Robert Cimrman

Hi Ronghai, On 11/05/2014 09:42 PM, Ronghai Wu wrote: > Dear Sfepy developers and users, > > I am wondering if it is possible to solve a non-constant prestress problem > with sfepy, as show > > > > where all fields exist in the whole domain(no subdomain is considered). > is a constant stress field, however is a non-constant scale field(having > different values at different nodes, but values are known). So, therefore > prestress is also non-constant. > I know the "linear_elasticity/prestress_fibres" example in Sfepy > documentation. But the prestress in that example is constant(having same > values at different nodes). So, I am wondering if it is possible to > implement non-constant prestress problem. Thanks. It should work out of the box. Look at [1] to see how to define a non-constant material parameter (the prestress here) by using functions. r. [1] http://sfepy.org/doc-devel/users_guide.html#defining-material-parameters

6 years, 8 months

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

SfePy November 2014