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)
$

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)

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