This series of examples shows how to scale the penalty stiffness in contact problems. The influence of this scaling on the solution is shown by the example of a soft ring in contact with a rigid wall. In this first example the penalty factor obviously is to low.
*CONTACT_ONE_WAY_SURFACE_TO_SURFACE *CONTROL_CPU *CONTROL_ENERGY *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_BINARY_D3THDT *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_NODOUT *DATABASE_RCFORC *ELEMENT_SHELL *END *INITIAL_VELOCITY_GENERATION *KEYWORD *MAT_ELASTIC *NODE *PART *SECTION_SHELL *TITLE
*KEYWORD
*TITLE
Soft Ring Impacting Rigid Wall (SF=0.01)
*CONTROL_TERMINATION
0.0300000 0 0.0000000 0 0.0000000
*CONTROL_TIMESTEP
0.0000000 0.0000000 0 0.0000000 0.0000000 0 0 0
*CONTROL_ENERGY
1 2 1 1
*DATABASE_BINARY_D3PLOT
5.00000-4 0
*DATABASE_BINARY_D3THDT
1.00000-4
*DATABASE_NODOUT
1.00000-4
*DATABASE_GLSTAT
1.00000-4
*DATABASE_MATSUM
1.00000-4
*DATABASE_RCFORC
1.00000-4
*CONTROL_CPU
0.0000000
*MAT_ELASTIC
1 7.85000-8 100.00000 0.0000000
*MAT_ELASTIC
2 7.85000-8 100.00000 0.0000000
*SECTION_SHELL
1 0 0.0000000 3.0000000 0.0000000 0.0000000
1.0000000 1.0000000 1.0000000 1.0000000 0.0000000
*SECTION_SHELL
2 0 0.0000000 3.0000000 0.0000000 0.0000000
1.0000000 1.0000000 1.0000000 1.0000000 0.0000000
*PART
material type # 1 (elastic)
1 1 1 0 0 0 0 0
*PART
material type # 1 (elastic)
2 2 2 0 0 0 0 0
*NODE
$
*ELEMENT_SHELL
$
*INITIAL_VELOCITY_GENERATION
$ ID STYP OMEGA VX VY VZ
1 2 0.0 2000.0 -1000.0 0.0
$ XC YC ZC NX NY NZ PHASE
0.0 0.0 0.0 0.0 0.0 0.0 0
$ INTERFACE NAME: 1 $$$
*CONTACT_ONE_WAY_SURFACE_TO_SURFACE
1 2 3 3 0 0 0 0
0.0 0.0 0.0 0.0 0.0 0 0.0 0.0
$ SFS SFM SST MST SFST SFMT FSF VSF
0.01 0.01 0.0 0.0 0.0 0.0 0.0 0.0
*END
