SfePy August 2020

sfepy@python.org

2 participants
4 discussions

Material Parameters of Individual Nodes
by derj...＠web.de 11 Apr '22

11 Apr '22

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

[sfepy]
by Takuo Fujita 02 Sep '20

02 Sep '20

Hi sfepy expert, I tried 1 hexa analysis with the attached file, but encountered the error below. The error seems due to unmatched load shape. Coud you give me an instruction? Best regards, Takuo $python linear_elastic_1hexa.py f Material:f constant_names: set() datas: dict with keys: [] extra_args: dict with keys: [] flags: dict with keys: [] function: ConstantFunction:get_constants_val is_constant: True kind: time-dependent mode: None name: f special_names: set() sfepy: saving regions as groups... sfepy: Omega sfepy: origin sfepy: xyz010 sfepy: xyz001 sfepy: xyz011 sfepy: xyz100 sfepy: x1 sfepy: ...done sfepy: updating variables... sfepy: ...done sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: matrix shape: (71, 71) sfepy: assembling matrix graph... sfepy: ...done in 0.00 s sfepy: matrix structural nonzeros: 5041 (1.00e+00% fill) sfepy: updating variables... sfepy: ...done sfepy: updating materials... sfepy: m sfepy: f sfepy: ...done in 0.00 s sfepy: allowed argument shapes for term "-1.0 * dw_point_load.0.x1(f.val, v)": sfepy: [{'material': '.: N', 'virtual': ('N', None)}] sfepy: actual argument shapes: sfepy: {"material:('f', 'val')": (0, 0, 0, 0), 'virtual:v': (3,)} sfepy: giving up! Traceback (most recent call last): File "linear_elastic_1hexa.py", line 93, in <module> state = pb.solve(status=status) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/discrete/problem.py", line 1414, in solve status=status) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/solvers/ts_solvers.py", line 34, in _standard_ts_call status=status, **kwargs) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/solvers/ts_solvers.py", line 73, in __call__ vec = nls(vec0) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/solvers/nls.py", line 257, in __call__ vec_r = fun(vec_x) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/discrete/evaluate.py", line 73, in eval_residual vec_r = self.problem.equations.eval_residuals(vec) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/discrete/equations.py", line 701, in eval_residuals self.evaluate(mode='weak', dw_mode='vector', asm_obj=out) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/discrete/equations.py", line 634, in evaluate term_mode=term_mode, asm_obj=asm_obj) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/discrete/equations.py", line 853, in evaluate ret_status=True) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/terms/terms.py", line 1436, in evaluate self.check_shapes(*args) File "/home/pyontaku14/anaconda3/lib/python3.7/site-packages/sfepy/terms/terms.py", line 1274, in check_shapes % term_str) ValueError: wrong arguments shapes for "-1.0 * dw_point_load.0.x1(f.val, v)" term! (see above) $

2 1

Error in 5 hexas linear static analysis
by Takuo Fujita 23 Aug '20

23 Aug '20

Hi sfepy expert, I tried 5 hexas linear static analysis and encountered the following error. Could you please tell me what should be corrected? Best regards, Takuo sfepy: allowed argument shapes for term "+1.0 * dw_lin_elastic_iso.2.Omega(solid.D, v, u)": sfepy: [{'material_1': '1, 1', 'material_2': '1, 1', 'virtual': ('D', 'state'), 'state': 'D', 'parameter_1': 'D', 'parameter_2': 'D'}] sfepy: actual argument shapes: sfepy: {"material_2:('solid', 'D')": (6, 6), 'virtual:v': (3,), 'state:u': (3,)} sfepy: giving up! Traceback (most recent call last): File "/home/pyontaku14/sfepy/simple.py", line 182, in <module> main() File "/home/pyontaku14/sfepy/simple.py", line 179, in main app() File "/home/pyontaku14/sfepy/sfepy/applications/application.py", line 29, in call_basic return self.call(**kwargs) File "/home/pyontaku14/sfepy/sfepy/applications/pde_solver_app.py", line 228, in call post_process_hook_final=self.post_process_hook_final) File "/home/pyontaku14/sfepy/sfepy/discrete/problem.py", line 1435, in solve status=status) File "/home/pyontaku14/sfepy/sfepy/solvers/ts_solvers.py", line 34, in _standard_ts_call status=status, **kwargs) File "/home/pyontaku14/sfepy/sfepy/solvers/ts_solvers.py", line 73, in __call__ vec = nls(vec0) File "/home/pyontaku14/sfepy/sfepy/solvers/nls.py", line 257, in __call__ vec_r = fun(vec_x) File "/home/pyontaku14/sfepy/sfepy/discrete/evaluate.py", line 73, in eval_residual vec_r = self.problem.equations.eval_residuals(vec) File "/home/pyontaku14/sfepy/sfepy/discrete/equations.py", line 702, in eval_residuals self.evaluate(mode='weak', dw_mode='vector', asm_obj=out) File "/home/pyontaku14/sfepy/sfepy/discrete/equations.py", line 635, in evaluate term_mode=term_mode, asm_obj=asm_obj) File "/home/pyontaku14/sfepy/sfepy/discrete/equations.py", line 854, in evaluate ret_status=True) File "/home/pyontaku14/sfepy/sfepy/terms/terms.py", line 1436, in evaluate self.check_shapes(*args) File "/home/pyontaku14/sfepy/sfepy/terms/terms.py", line 1274, in check_shapes % term_str) ValueError: wrong arguments shapes for "+1.0 * dw_lin_elastic_iso.2.Omega(solid.D, v, u)" term! (see above)

2 2

Query on sfepy for Finite element analysis
by Aanchal Jasuja 03 Aug '20

03 Aug '20

Hello, We are trying to use sfepy for FEM analysis. However, we have a couple of queries about sfepy. 1. We couldn't find any examples for 1-D line element modelling, although it says so in the description. Can sfepy be used for it or just for plates & shells? 2. According to the manual, the outputs are in the form of stresses & strains. Is there any possible way to directly get results in the form of bending moment, shear force, reactions etc. ? We would really appreciate your support response at the earliest convenience. Thank you Regards Aanchal Jasuja

2 1

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

SfePy August 2020