I am attempting to simulate a peristaltic pump. The model consists of a bar with rounded ends and a tube that goes around a diameter equivalent to the length of the bar. The bar is supposed to be compressing the tube at the contact region but since it has been modeled using actual dimensions, the bar and the tube are overimposed. A surface-to-surface contact has been defined between the rounded end of the bar and the inner wall of the tube. Both tube ends are fixed as well as the bent region (to simulate contact against a wall). The center of the bar is also fixed for the static analysis before turning the pump.
Before attempting any motion or dynamic analysis, I wish to try a static analysis to determine the force, stresses and deformation resulting from the contact only. Afterwards, I'll try to complete a motion simulation (something like this: https://youtu.be/DGBfLxwyS6Y). When I solve with the Sol101 solver, I get some unexpected results. It looks weird like the bar penetrating the tube. I am pretty sure that the issue might be related to the contact definition or some related parameter. I have included the simulation files.
Thank you for your assistance
I am sorry to tell you but linear static analysis (SOL101) is useless in this case, you have a nonlinear problem because large displacement effect are present in the problem, you need to use ADVANCED NONLINEAR ANALYSIS (SOL601). This is an interference-fit problem, then you need to use a negative minimum search distance to account for the interference and get the solver to move any component to its final position.
Thank you, Blas
I appreciate your response. I tried the Sol 601 but it did not run. I suspected something was wrong with the model so I decided to remake the model. Afterwards, at least it started running but now I'm getting a negative jacobian error. I've tried remeshing with smaller element sizes to no avail. I know I am missing something simple. It even tells me which is the element but I do not know how to identify and correct that element.
It has come to my mind to split the bodies to refine the mesh in particular regions but I would like to make sure which elements are causing the error before trying something.
Any thoughts? I included the model files.
First at all you need to prepare the geometry to mesh effciently with HEX elements, avoid to have 3-D solid TETRA or 3-D solid triangular prismatic elements, the best one is the CHEXA elements. If you prepare the geometry (either is the master or idealized model) like this way slicing the solid tube & dividing the surface of the block using a 2D SEED mesh to control the quality of the 3-D solid mesh then you will have a correct mesh to work with the analysis options.
The 2-D Seed mesh using command "2D Mapped Mesh" is the following (make sure to desactivate the option "export mesh to solver"):
And then the 3-D solid Hexaedral mesh is really easy to create using command "3D Swept Mesh", here a CHEXA 20-node elements were create to account for effective contact between parts (later in NXSTRAT command of NX NASTRAN ADVANCED NONLINEAR (SOL601,106) make sure to convert CHEXA 20-nodes to 27-nodes (with extra midside nodes) using ELCV=-1 under NXSTRAT > TRANSLATION):
Also, in the FEM environment you need to define the material properties of the components: please note if the TUBE is a rubber-like material the you need to enter the nonlinear HYPERELASTIC material properties using, for instance, either the MOONEY-RIVLIN, ODGEN or ARRUDA-BOYCE constants; this is a critical steep in the simulation because is difficult to have the correct information available.
Later in the SIM environment you need to define important data using the Modeling Objects Manager like Time Step, SOL601 NXSTRAT Strategy parameters, Surface-to-Surface Contact parameters for the Advanced Nonlinear Pair, the Structural Output Request, etc..
For the Contact setup, here you have two options:
The solution is not easy, the convergence in this type of nonlinear contact problems is critical. In the NXSTRAT card select Automatic Time Stepping (ATS) and select the 3D-Iterative solver, because the model has high-order CHEXA 27-nodes elements it should speed the solution.
To learn more about Advanced Nonlinear I suggest to run the TUTORIAL EXAMPLE under NX 10 HELP > CAE > Tutorials: Advanced Simulation, this way you can make us specific questions, not how to run the full model, OK?.