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 25 Downloads (last 30 days) File Size: 2.73 KB File ID: #43216 Version: 1.0
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 (2)
30 Dec 2015 Krishnan S

Which point load solution has been used finding the kernel?

Comment only
19 Jun 2014 Daniel

Daniel (view profile)

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