3.2.3 Elasticity problem

Assume that we want to solve an elasticity problem with an isotropic elasticity tensor determined by the Lame parameters \lambda=\mu=1 on the unit square \Omega=(0,1)^2 with right-hand side \vecf(x,y)\equiv (1,0)^t and Dirichlet boundary conditions \vecu(x,y)=0 for (x,y) \in \partial \Omega. The following command solves this equation approximately on a uniformly refined mesh using as termination criterion that the time for approximating the solution has increased beyond 20 seconds.

     (let* ((problem
             (standard-elasticity-problem
              (n-cube-domain 2) :lambda 1.0 :mu 1.0
              :force (constant-coefficient (vector #m(1.0) #m(0.0)))))
            (blackboard
             (blackboard :problem problem :output t
                         :success-if '(> :time 20.0))))
       (solve blackboard)
       (plot (getbb blackboard :solution) :component 0)
       (plot (getbb blackboard :solution) :component 1))