OS-V: 0760 MacNeal-Harder Solid Patch Test

MacNeal-Harder TestThe patch test is a classical benchmark problem for the element, if it produces correct results for the test, the result for any problem solved with the element will converge toward the correct solution. The intended purpose of proposed problem set is to ascertain the accuracy of finite element in various applications.



Figure 1. Cube subjected to uniform displacement

Benchmark Model

The outer dimension have unit cube of 1mm size. A mesh of cube having node location as mentioned in the table with first order CHEXA element. The eight corners of the cube are constrained in all three translational direction and free in all three rotational directions. Displacement is enforced using SPCD on the eight nodes of cylinder in X, Y and Z translation directions of the cube.

The material properties are:
Material Properties
Value
Young's Modulus
1 x 106 Pa
Poisson's Ratio
0.25


Figure 2. Patch test for solids
Table 1. Location of Inner Nodes
x y z
1 0.249 0.342 0.192
2 0.826 0.288 0.288
3 0.850 0.649 0.263
4 0.273 0.750 0.230
5 0.320 0.186 0.643
6 0.677 0.305 0.683
7 0.788 0.693 0.644
8 0.165 0.745 0.702

The arbitrarily distorted element shapes are an essential part of the test. The principal virtue of a patch test is that if an element produces correct results for the test, the results for any problem solved with the element will converge toward the correct solution as the elements are subdivided. On the other hand, passing the patch test does not guarantee satisfaction, since the rate of convergence may be too slow for practical use. The above patch test is an extension of Robinson’s patch test to three dimensions.

Displacement boundary conditions for the test are:
u
10−3(2x+y+z)/2 MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8qacaaIXaGaaGima8aadaahaaWcbeqaa8qacqGHsislcaaIZaaaaOWdamaabmaabaWdbiaaikdacaWG4bGaey4kaSIaamyEaiabgUcaRiaadQhaa8aacaGLOaGaayzkaaWdbiaac+cacaaIYaaaaa@421B@
v
10−3(x+2y+z)/2 MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8qacaaIXaGaaGima8aadaahaaWcbeqaa8qacqGHsislcaaIZaaaaOWdamaabmaabaWdbiaadIhacqGHRaWkcaaIYaGaamyEaiabgUcaRiaadQhaa8aacaGLOaGaayzkaaWdbiaac+cacaaIYaaaaa@421B@
w
10−3(x+y+2z)/2 MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8qacaaIXaGaaGima8aadaahaaWcbeqaa8qacqGHsislcaaIZaaaaOWdamaabmaabaWdbiaadIhacqGHRaWkcaWG5bGaey4kaSIaaGOmaiaadQhaa8aacaGLOaGaayzkaaWdbiaac+cacaaIYaaaaa@421B@

Results

εx= εy= εz= γxy= γyz= γzx= 10−3 MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8qacqaH1oqzpaWaaSbaaSqaa8qacaWG4baapaqabaGcpeGaeyypa0Jaaeiiaiabew7aL9aadaWgaaWcbaWdbiaadMhaa8aabeaak8qacqGH9aqpcaqGGaGaeqyTdu2damaaBaaaleaapeGaamOEaaWdaeqaaOWdbiabg2da9iaabccacqaHZoWzpaWaaSbaaSqaa8qacaWG4bGaamyEaaWdaeqaaOWdbiabg2da9iaabccacqaHZoWzpaWaaSbaaSqaa8qacaWG5bGaamOEaaWdaeqaaOWdbiabg2da9iaabccacqaHZoWzpaWaaSbaaSqaa8qacaWG6bGaamiEaaWdaeqaaOWdbiabg2da9iaabccacaaIXaGaaGima8aadaahaaWcbeqaa8qacqGHsislcaaIZaaaaaaa@58E7@

σx= σy= σz= 2000 MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8qacqaHdpWCpaWaaSbaaSqaa8qacaWG4baapaqabaGcpeGaeyypa0Jaaeiiaiabeo8aZ9aadaWgaaWcbaWdbiaadMhaa8aabeaak8qacqGH9aqpcaqGGaGaeq4Wdm3damaaBaaaleaapeGaamOEaaWdaeqaaOWdbiabg2da9iaabccacaaIYaGaaGimaiaaicdacaaIWaGaai4oaaaa@4839@

τxy= τyz= τzx= 400 MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8qacqaHepaDpaWaaSbaaSqaa8qacaWG4bGaamyEaaWdaeqaaOWdbiabg2da9iaabccacqaHepaDpaWaaSbaaSqaa8qacaWG5bGaamOEaaWdaeqaaOWdbiabg2da9iaabccacqaHepaDpaWaaSbaaSqaa8qacaWG6bGaamiEaaWdaeqaaOWdbiabg2da9iaabccacaaI0aGaaGimaiaaicdaaaa@49C2@


Figure 3. Elemental strains in all 6 direction plot


Figure 4. Elemental stresses in all 6 direction plot

The results CHEXA elements agree with the reference results.

Model Files

Refer to Access the Model Files to download the required model file(s).

The model files used in this problem include:

chexa_mh_patch_test.fem

Reference

MacNeal, R.H., and Harder, R.L., A Proposed Standard Set of Problems to Test Finite Element Accuracy, Finite Elements in Analysis and Design, 1 (1985) 3-20