//-------------------------------------------------
// Bimaterial bar
// cf. Sukumar et al_2001_Modeling holes and inclusions by level sets in the extended finite-element.pdf
//
// Geo: Plate 2*2,
//      Interface position 0.1
// Properties: E1=1, v1=0. (left)
//             E2=10, v2=0. (right)
// Loading: displacement field u(x=-1)=0, u(x=1)=1
//-------------------------------------------------

Dim 2
MeshFile bimaterialBar_cut.msh

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

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

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

BoundaryCondition Edge 16 Displacement 1 Constant Scalar 0.
BoundaryCondition Edge 19 Displacement 1 Constant Scalar 0.
BoundaryCondition Edge 17 Displacement 1 Constant Scalar 0.
BoundaryCondition Edge 18 Displacement 1 Constant Scalar 0.

//BoundaryCondition Face 44 Displacement 2 Constant Scalar 0.
//BoundaryCondition Face 43 Displacement 2 Constant Scalar 0.

// Interface Domains
LevelSet 24 1 Function Scalar 1.*x+0.*y+0.*z-0.1 [x,y,z]
Domain 24 1 Bimaterial Data DomainInt Integer 44 DomainExt Integer 43 EndData

Solver 3