Code covered by the BSD License  

Highlights from
An LCP solution of the linear elastic contact mechanics problem

4.0

4.0 | 1 rating Rate this file 16 Downloads (last 30 days) File Size: 2.73 KB File ID: #43216
image thumbnail

An LCP solution of the linear elastic contact mechanics problem

by

 

An LCP solution of the one-dimensional linear elastic contact mechanics problem

| Watch this File

File Information
Description

A script that solves the linear complementarity problem associated with
the Boussinesq contact mechanics problem

  h = K*p+g, h'*p=0, h,p >=0,

where the unknown h is the gap between two elastic (half-space) bodies and the unknown p is the corresponding contact pressure.

Acknowledgements

A Pivoting Algorithm Solving Linear Complementarity Problems inspired this file.

MATLAB release MATLAB 8.0 (R2012b)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (1)
19 Jun 2014 Daniel

Very useful, well written and fast code. However, it lacks dimensional definitions/explanations for the input parameters (Radii, E, iv, g0). It seems there is a mix of meters and millimeter but it is hard to judge the quantities.

Typically:
2/E' = (1-nu1^2)/E1+(1-nu2^2)/E2 is also defined as:
1/E' = (1-nu1^2)/E1+(1-nu2^2)/E2 by Johnson but the factor 2 is compensated for in your calculations of ph and b, but I cannot see whether it is compensated in K as well.

Contact us