Hello SfePy-ers!<div><br></div><div>I'm starting to work with the tool. I found it great, very nice to use.</div><div><br></div><div>I have a problem that I can not wrap my head around:</div><div><br></div><div>I have a cylinder close to a sphere, both with different elasticity properties (the cylinder is rigid and the sphere is more elastic).</div><div><br></div><div>I'm trying to rotate the cylinder about the centre of its bottom cap and thus hit the sphere deforming it.</div><div><br></div><div>To do this I put the sphere and the cylinder close, I put them in a bounding box representing the media in wich both are inscribed (air, water, etc), I calculate a mesh for this domain with tetgen &nbsp;having a different material for each object (cylinder, baloon, media)</div><div><br></div><div>and then use a function in the ebcs conditions to move the rod. Similar to the large deformations hyperelastic example (rotate_yz), where the rotation angle increases with time.</div><div><br></div><div>My problem is that if my cylinder rotates more than ~.1 degrees, I get a "warp violation" error.</div><div><br></div><div>How could I solve this? By making the elements of my mesh smaller? Taking smaller steps, increasing the order of the balance integral holding the elasticity properties?</div><div><br></div><div>Any help will be great</div><div><br></div><div>Thanks</div><div>Demian</div><div><div>--</div><div>Demian Wassermann, PhD</div><div>dem...@bwh.harvard.edu</div><div>LMI / PNL / SPL Labs</div><div>Harvard Medical School</div><div>Brigham and Women's Hospital</div><div>1249 Boylston, Boston, MA, USA</div><div>http://lmi.bwh.harvard.edu/~wassermann/</div></div>