Hi Takuo,

the dw_lin_elastic_iso term takes the Lame's parameters as arguments (two scalars, see tests/test_elasticity_small_strain.py), if you want to specify the stiffness tensor D, you need to use the dw_lin_elastic term.

r.

On 23. 08. 20 3:59, Takuo Fujita wrote:

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)

---

SfePy mailing list -- sfepy@python.org To unsubscribe send an email to sfepy-leave@python.org https://mail.python.org/mailman3/lists/sfepy.python.org/ Member address: cimrman3@ntc.zcu.cz