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