Contact problems from the FEBio test suite. These files are run through FEBio after each compilation to ensure that it is working as intended.
Owner: FEBio Team
contact
Files
Filename
Description
Tags
Download
co01.feb
Two elements with sliding contact are compressed against each other. The nodes on the left face of the left element are constrained in the xdirection. The nodes on the right face of the right element are given a displacement of 0.2 in the xdirection in time t=1.
Not using the augmented Lagrangian gives slightly better convergence, but the contact penetration is greater. This problem converges with fewer equilibrium iterations if the automatic time stepper is turned off.
contact
co02.feb
An element with prescribed displacement of 0.5 units in the zdirection compresses a 5x5x4 element block with sliding contact.
Convergence and contact penetration are improved by setting auto='on' in the contact penalty.
contact
co04.feb
A stiff block with prescribed displacement of 0.14 units in the zdirection from time t=0 to t=1 and then 0.8 units in the xdirection from time t=1 to t=2 compresses then slides across a larger softer block.
Does not converge with auto='on' in the contact penalty. Turning the auto time stepper off leads to better convergence and contact penetration using augmented Lagrangian. Not using augmented Lagrangian creates too much penetration.
contact
co07.feb
A stiff section of a cylinder has a prescribed displacement in the zdirection 0f 0.35 units from time t=0 to t=1 and then a displacement of 3 units in the xdirection from time t=1 to t=4. It presses down then slides along a softer block.
This problem demonstrated the sliding2 symmetric contact with augmented Lagrangian turned off. It did not run well without using fullNewton (setting max_ups to 0). Even though fullNewton converged with fewer iterations than BFGS, it took more than four times as long to do so. However, with these setting, this contact method did better than the other three.
contact
co08.feb
A quarter disk is compressed by a rigid wall with a prescribed displacement of 0.5 units in the ydirection from time t=0 to t=0.5 and then returning to 0 displacement from time t=0.5 to t=1.
contact
co09.feb
Block with prescribed displacement of 2.5 units in zdirection from time t=0 to t=1 sliding and compressing against inclined rigid body plane.
contact
co10.feb
Block with prescribed displacement of 2.5 units in zdirection from time t=0 to t=1 sliding and compressing against inclined rigid body plane.
contact
co11.feb
A Rigid ball with a prescribed displacement of 5.0 units in the zdirection pushes down and slides on a cantilever board.
Convergence and contact penetration are better without using augmented Lagrangian.
contact
co12.feb
A Rigid ball with a prescribed displacement of 5.0 units in the zdirection pushes down and slides on a cantilever board.
Convergence and contact penetration are better without using augmented Lagrangian. This is the same problems as co11 (contact11) using the 'element face' format for the contact surface definition.
contact
co13.feb
Two cylindrical tubes at 90 degrees are compressed against each other. The ends of one tube are constrained in x,y,z and the ends of the other tube are attached to a rigid body which has a prescribed displacement of 6.0 units in the xdirection.
Penetration is slightly better using augmented lagrangian, but convergence is better without it.
contact
co15.feb
Two cylindrical tubes at 90 degrees are compressed against each other. The ends of one tube are constrained in x,y,z and the ends of the other tube have a prescribed displacement of 6.0 units in the xdirection.
This problem is similar to co13 except the ends of the moving tube are given a prescribed displacement instead of being attached to a moving rigid body.
contact
co16.feb
A smaller element with a prescribed zdisplacement of 0.1 units from time t=0 to t=0.5 and xdisplacement of 0.3 units from time t=0.5 to t=1 compresses a larger element with contact and friction.
contact
co17.feb
A block with friction has a zdisplacement of 0.4 units from time t=0 to t=0.2 and an xdisplacement of 3 units from time t=0.2 to t=1 along a rigid body.
contact
co18.feb
An element with friction has a zdisplacement of 0.1 units from time t=0 to t=0.2 and an xdisplacement of 3 units from time t=0.2 to t=1 along a rigid body.
Similar to problem c017, but with a greater coefficient of friction.
contact
co19.feb
A block with a prescribed zdisplacement of 1.5 units from time t=0 to t=7.5 slides with friction along a rigid inclined plane.
contact
co20.feb
Two cylindrical tubes at 90 degrees are compressed against each other. The ends of one tube are constrained in x,y,z and the ends of the other tube have a prescribed displacement of 6.0 units in the xdirection.
This problem is the same as co15 but with friction. Two vertical lines intersecting the plane y = 0 and the verticle tube are constrained in the ydirection and two horizontal lines intersecting the plane z = 5 and the horizontal tube are constrained in the zdirection.
contact
co21.feb
A stiff section of a cylinder has a prescribed displacement in the zdirection 0f 0.35 units from time t=0 to t=1 and then a displacement of 3 units in the xdirection from time t=1 to t=4. It presses down then slides with friction along a softer block.
This problem is similar to co07 but with friction.
contact
co22.feb
A 2dimensional billet block is compressed between two rigid bodies in plain strain which each have a prescribed zdisplacement (compressing the block) of 0.6 units from time t=0 to t=1.
Gets better convergence and less contact penetration using augmented Lagrangian.
contact
co26.feb
Two blocks have tied contact between them. The left end nodes of the left block are constrained in x,y,z and the right end nodes of the right block have a prescribed xdisplacement of 1 unit in time t=1.
contact
co27.feb
Cleveland Clinic 2d foot model
contact
co28.feb
A 2dimensional billet block is compressed between two rigid bodies which each have a prescribed zdisplacement (compressing the block) of 0.6 units from time t=0 to t=1.
This problem is similiar to co22 but with Ogden material and 100 time steps.
contact
co29.feb
A 2D quarter billet with 1/4th symmetrical representation is compressed by a rigid bar with prescribed ydisplacement of 0.3 units in time t=1.
This problem is similar to co28 but with quater symmetry. It has good penetration results and converges faster with augmented Lagrangian turned off
contact
co30.feb
A rigid block subject to gravity bounces off a deformable 25 x 25 element membrane, also subject to gravity.
contact
co31.feb
A 10 x 10 element planar block rises (with compression) in plane strain over a ramped rigid bump.
This model has significant pentration without using Augmented Lagrangian. The sliding_with_gaps contact algorithm performed significantly better than the other three contact algorithms.
contact
co32.feb
Two blocks undergo confined compression with nonconforming mesh contact between them. The bottom block is defined with a variable thickness mesh.
contact
co34.feb
Two incongruent cylindrical sections are compressed against each other in plain strain.
The penetration is actually a bit better using Augmented Lagrangian, but this demonstrates that a reasonable penetration can be obtained with trialanderror selection of the contact penalty parameter. The other contact methods either error terminated or had too significant penetration.
contact
co35.feb
Two 5 x 5 element blocks are compressed against each other with contact between them.
contact
co36.feb
A small hollow sphere rests inside a larger rigid hollow sphere and is allowed to contract due to the force of gravity.
contact
co37.feb
A semispherical rigid body in contact with a block of Ogden material is displaced 0.2 units in the zdirection (into the block) at time t=1.
This problem was changed to quarter symmetry to eliminate inconsistent convergence.
contact
co38.feb
A 2dimensional billet block is compressed between two rigid bodies which each have a prescribed zdisplacement (compressing the block) of 0.6 units from time t=0 to t=1.
This problem is similiar to co22 but with Ogden unconstrained material and 10 time steps.
contact
co39.feb
A 2D quarter billet with 1/4th symmetrical representation is compressed by a rigid bar with prescribed ydisplacement of 0.3 units in time t=1.
This problem is similar to co28 but with quater symmetry.
contact
co40.feb
A semispherical body in contact with a block, both consisting of Ogden unconstrained material. The semispherical body is displaced 0.2 units in the zdirection (into the block) at time t=1.
contact
co41.feb
A serpentine surface is compressed along the xdirection experiencing selfcontact and demonstrating the need of reducing the search_radius contact parameter from it's default value of 1.
contact
co42.feb
A strip is pulled off of a rigid surface using sticky contact.
contact
co43.feb
Patch contact of multiphasic rectangular strips in unconfined compression stress relaxation
contact
co44.feb
A smaller element with a prescribed zdisplacement of 0.1 units from time t=0 to t=0.5 and xdisplacement of 0.3 units from time t=0.5 to t=1 compresses a larger element with contact and friction.
contact
co45.feb
A block with friction has a zdisplacement of 0.4 units from time t=0 to t=0.2 and an xdisplacement of 3 units from time t=0.2 to t=1 along a rigid body.
contact
co46.feb
An element with friction has a zdisplacement of 0.1 units from time t=0 to t=0.2 and an xdisplacement of 3 units from time t=0.2 to t=1 along a rigid body.
Partial contact slip initiates at time t=0.23 and full contact slip at t=0.25.
contact
co47.feb
Two cylindrical tubes at 90 degrees are compressed against each other. The ends of one tube are constrained in x,y,z and the ends of the other tube have a prescribed displacement of 6.0 units in the xdirection.
This problem is the same as co15 but with friction. Two vertical lines intersecting the plane y = 0 and the verticle tube are constrained in the ydirection and two horizontal lines intersecting the plane z = 5 and the horizontal tube are constrained in the zdirection.
contact
co48.feb
A stiff section of a cylinder has a prescribed displacement in the z-direction 0f -0.35 units from time t=0 to t=1 and then a displacement of 3 units in the x-direction from time t=1 to t=4. It presses down then slides with friction along a softer block.
This problem is similar to co07 but with friction.
contact
co49.feb
Contact of two nonconforming multiphasic layers under displacement control
contact
co50.feb
A rigid noncircular cam rotates about a revolute joint and contacts a flat rigid follower, whose motion is restricted by a prismatic joint in parallel with a spring.
Lagrange augmentation cannot be used for contact between rigid bodies.
contact
co51.feb
Contact with top and bottom faces of shells.
There are two contact interfaces, one for each cylindrical shell contacting the rectangular shell. The contact area on the recangular shell was split to allow display of contact pressure and contact gap on each half correctly.
contact
co52.feb
Multistep analysis with tiedelastic contact defined in second step. Tied contact engages with no gap between surfaces.
contact
co53.feb
Multistep analysis with tiedelastic contact defined in second step. Tied contact engages with an initial gap between surfaces.
contact
co54.feb
2D frictional biphasic contact between biphasic square slab and rigid plane under cyclical loading and sliding
contact
co55.feb
3D frictional biphasic contact between biphasic cube and rigid plane under cyclical loading and sliding
contact
co56.feb
2D frictional biphasic contact between rigid cylinder and biphasic slab under prescribed load and reciprocal sliding
contact
co57.feb
2D frictional biphasic contact between elastic cylinder and biphasic slab under prescribed load and reciprocal sliding
contact
co58.feb
2D frictional biphasic contact between biphasic cylinder and biphasic slab under prescribed load and reciprocal sliding
FEBio Studio Model Repository