//-------------------------------------------------
// rectangle
//
// Geo: Plate [-0.5,+0.5]*[0.0,5.0]
//      Interface y=2.5
// Properties: E1= 1, v1=0. (down)
//             E2=10, v2=0. (up)
// Loading: displacement field v(y=0)=0, v(y=5)=1
// Results: strain1=0.364 (cst)
//          strain2=0.0364 (cst)
//          stress1=0.364 (cst)
//          stress2=0.364 (cst)
//          elastic energy=0.181818
//          total volume = 5.0 [L^2]
// Rq: elastic energy=0.5*Integ(Plate)(strain:stress)
//-------------------------------------------------

MeshFile rectangle.msh

// Material Domains
Domain 1 2 IsotropicElastic Data E Scalar  1. nu Scalar 0. EndData
Domain 2 2 IsotropicElastic Data E Scalar 10. nu Scalar 0. EndData

// Loading
BoundaryCondition Edge 1 Displacement 0 Constant Scalar 0.
BoundaryCondition Edge 1 Displacement 1 Constant Scalar 0.

BoundaryCondition Edge 2 Displacement 0 Constant Scalar 0.
BoundaryCondition Edge 2 Displacement 1 Constant Scalar 1.

Solver 3