I’m trying to simulate uniaxial compression of a 3D cylindrical specimen. the mesh is imported with markers. Loading is in the -z direction and the deformation is expected in the x and y directions for both the top and bottom surfaces (kind of assuming zero friction). However, I get results showing a bias deformation in y direction. Attached is a part of code where boundary conditions are applied. I believe that I’m making a mistake in bcs. Looking forward to some advise. Thank you in advance.
Your lines of code look Ok to me, although it is impossible to say anything definitive without the context in which they are written.
Of course, with these boundary conditions, a rigid-body mode (dispacement in y) remains unconstrained. It could simply be that that mode is somehow triggered in the iterative solve. An easy fix is to constrain V.sub(1) in a single point.
I ran into more errors with the suggested modification of the code. Specifically,
’ ’ ’
TypeError: int() argument must be a string, a bytes-like object or a real number, not 'function
’ ’ ’
So, I opted to model a quarter of a cylinder and constrain the x and y directions. It works as expected. Thanks.