increase stiffness in contact between shell and solid - LSDYNA

WasifWasif Member Posts: 35

Hello everyone,

I am doing a three point bending test. In the results below, you will see that when I use mid-plane surface in the shell, the contact stiffness is lower. However, when I use NLOC, the stiffness is higher but there seem to vibrations in results. I tried to set VDC to 100% percent but the vibrations still exist. Is there a workaround to this problem where we can increase the stiffness without going to material properties?

PS. I tried changing penalty factors in the contact but that does not work.


  • AndreasKoutrasAndreasKoutras Posts: 11Member

    Hello Wasif,

    You may tune the contact penalty stiffness through the parameters SLSFAC (in CONTROL_CONTACT), SFS and SFM (in *CONTACT), SOFSCL (in *CONTACT optional card A) depending on what contact you are using. See Sections 29.7.1, 29.7.2, and 29.7.3 of the Theory Manual for how the penalty stiffness is determined in SOFT=0, SOFT=1, and SOFT=2 contacts. For a general purpose contact in explicit analysis we recommend *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE with SOFT=1. The penalty contact stiffness of SOFT=1 is determined from:

    k = max(SLSFAC*SFi*k0, SOFSCL*k1)


     k is the penalty stiffness

     SLSFAC is user input in *CONTROL_CONTACT

     SFi is either SFS or SFM in *CONTACT card 3

     SOFSCL is user input in *CONTACT optional card A

     k0 is the stiffness calculated from the bulk modulus of the material and the element dimensions, as described in Section 29.7.1.

     k1 is the stiffness calculated from nodal masses and the solution time step, as described in Section 29.7.2.

    If you increase the contact penalty stiffness to a value higher than the default one, you may need to specify a lower time step through TSSFAC in *CONTROL_TIMESTEP to ensure stability of the explicit analysis.

    For implicit analysis, the recommended contact is *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR. The contact pressure -versus- penetration constitutive relation can be tuned through the parameters SLSFAC, SFS, IGAP, and PENMAX as described in Section 29.21.2 of the Theory Manual.

    If you want to tie the shells to the solids, you will need to look at the CONTACT_..._TIED and CONTACT_..._TIEBREAK options.

    I hope this helps.



Sign In or Register to comment.