Hay que tener en cuenta las propiedades elásticas del contacto para determinar la distribución de fuerzas. En realidad, lo que te interesa es la distribución de presiones, ya que sumas las presiones para obtener la fuerza, y sumas el momento de presión para obtener el par. Imagina un área de contacto con el centro en el punto de contacto, el z -eje normal al contacto, y el xy plano a lo largo del contacto.
$$ \begin{aligned} F_z & = \int P(x,y) \, {\rm d} A \\ M_x & = \int (-y) P(x,y) \,{\rm d} A \\ M_y & = \int (+x) P(x,y)\, {\rm d} A \end{aligned} $$
Así que dada una distribución de presión $P(x,y)$ se pueden encontrar las fuerzas, pero ¿cómo obtener la distribución de la presión?
El ejemplo más básico es un objeto que empuja sobre un plano. Consideremos el caso en el que hay una cantidad de solapamiento conocida $a(x,y)$ en función de la ubicación en el contacto
Como resultado, la superficie del objeto (curva roja) debe deformarse en $\delta(x,y)$ para ajustarse al solapamiento. Para ello, debe aplicarse al objeto una presión superficial (curva amarilla)
La relación de desviación de la presión fue desarrollada por Boussinesq y se ve así
$$ \delta(x,y) = \frac{1-\nu^2}{\pi E} \iint \frac{P(u,v)}{r} {\rm d} A $$
donde <span class="math-container">$r=\sqrt{(x-u)^2+(y-v)^2}$</span> y <span class="math-container">$\nu$</span> es el <a href="https://en.wikipedia.org/wiki/Poisson%27s_ratio" rel="nofollow noreferrer">Relación de Poisson </a>del material y <span class="math-container">$E$</span> es el <a href="https://en.wikipedia.org/wiki/Elastic_modulus" rel="nofollow noreferrer">Módulo elástico </a>.
Un caso especial para lo anterior es el Contacto hertziano que trata de una esfera sobre un plano, una esfera sobre una esfera, un cilindro sobre un cilindro (contacto 2D) y un contacto general con forma de fútbol (dos radios).
El caso general sólo puede ser tratado numéricamente con el contacto se dividen en una cuadrícula y se desarrolla un sistema gigante de ecuaciones de la forma $\delta_i = C_{i j} P_j$ (Solución Hartnett).
Los pasos son los siguientes:
- Asumir una cierta cantidad de solapamiento $a(x,y)$ derivado de la cantidad máxima de penetración.
- Establece las deflexiones iguales al solapamiento, y resuelve las presiones $P(x,y)$ .
- Calcule la carga total y los momentos, y ajuste el solapamiento en consecuencia. Si la carga total es inferior a la carga aplicada, aumente el solapamiento, y si el momento total es inferior al par aplicado, incline el objeto
- Si las fuerzas y los momentos siguen estando desequilibrados, vaya al paso 2 y repita la operación hasta lograr la convergencia.
Ver este documento detallado para un ejemplo de análisis numérico de contactos. Y aquí hay un documento del NIST con algo de teoría sobre el problema general de contacto.