Namespace Example_commented

const arch = require('arch')


const model = new arch.Model()
model.setMaterial(0.25, 1, 1000)

// --------------------------------------

const vertices  = [...]
const triangles = [...]
const surface = new arch.Surface(vertices, triangles)
model.addSurface(surface)
surface.setBC("dip",    "free", 0)
surface.setBC("strike", "free", 0)
surface.setBC("normal", "free", (x,y,z) => -2000*9.81*z)

// --------------------------------------

// Add a "pre-defined" Traction-Inequality-constraint (Coulomb friction)
// NOTE: Arch for node.js does not support yet the inequalities.
// If needed, just let me (Frantz) know...
//
const c1 = arch.Coulomb()
c1.setFriction(0.2)
c1.setCohesion(0.01)
surface.addConstraint(c1)

// --------------------------------------

// Add a "user-defined" Displacement Inequality Constraint
// Make dip displacement to be >= 0
//
surface.addConstraint(
     new arch.UserDic( u => [u[0], u[1]<0?0:u[1], u[2]] )
)

// --------------------------------------

// Add another "user-defined" Displacement Inequality Constraint
// Make:  0 < dx < 5
surface.addConstraint(
     new arch.UserDic( u => {
         if (u[0] > 5) u[0] = 5
         if (u[0] < 0) u[0] = 0
         return u
     })
)

// --------------------------------------

// Add a remote stress
const rho = 2200 // Density between 0 and 3000
const Rh  = 0.1  // Normalized Sh according to Sv, between 0 and 2 and <= SH
const RH  = 0.6  // Normalized SH according to Sv, between 0 and 2 and >= Sh
const Sv  = rho*9.81*z
const remote = new arch.UserRemote( (x,y,z) => [Rh*Sv, 0, 0, RH*Sv, 0, Sv] )
model.addRemote( remote )

// --------------------------------------

// Solve the problem
const solver   = new arch.Forward(model)
solver.select("seidel")
solver.setMaxIter(2000)
solver.setEps(1e-9)
solver.run()

// --------------------------------------

// Post-process at some random points (with z=0)
const solution = new arch.Solution(model)
const stresses = new Array(100)
     .fill([0,0,0])
     .map( v => [Math.random(), Math.random(), 0])
     .map( p => solution.stressAt(p[0], p[1], p[2]))  // compute the stress at p
     .map (s => console.log(s))