Dear all,

Is there any way of modeling the Payne effect (also known as Fletcher-Gent effect) by using Abaqus connector elements in the frequency domain SSD analysis?

The Payne effect can be defined as the dependency of the dynamic stiffness (k_dyn) on the relative dynamic amplitude (u_dyn) which is a harmonic amplitude.

I’m trying to model the Coulumb friction related behavior damper, not a rubber related behavior damper. Fortunately, in my theoretical model, dynamic stiffness can be assumed as only just the function of dynamic amplitude, since I don’t deal with a rubber damper.

In summary, I have a lookup table defining the k_dyn = f(u_dyn) relationship. Thus, I would like to embed this lookup table in the Abaqus “direct-SSD” frequency domain solution procedure, by iterating within each frequency value, until it converges at each frequency value.

Note-1: In Abaqus documentation, it asserts that “in doing perturbation analysis, (such as frequency-domain steady-state dynamic analysis) with Abaqus, temperature and field variable variations are not permitted within an analysis step”. Thus, I understand that “predefined field variable” option is not applicable to the problem.

Note-2: In the following webinar link (https://www.youtube.com/watch?v=ZX9MsX1Lt48&ab_channel=MSCSoftware), the capability of implementing a “nonlinear harmonic analysis of rubber components” by using the MSC.Marc software is mentioned. It is mentioned that, this capability is developed to analyze the 3D continuum elements.

Inspired by this webinar, I succeed to develop a Matlab script (for a single degree of freedom system with harmonic base excitation) that can model the Payne effect using this iterative kind of “nonlinear harmonic analysis” approach. As a result, the Payne effect was clearly observed on the frequency domain transmissibility function.

Thanks in advance.