solution
    case 'buckling'
        buckling
          nmodes 10
          shift = -1e5
end

file
    geometry_file cylinder_solid.g
end

parameters
 # eig_tol = 1.0e-8
end

boundary
    sideset = bottom
        fixed
end

loads
    #sideset = top
    #    # force = {force = 1.0}
    #    # area = {area = 2.654*2}
    #    pressure = {force/area}
    nodeset = top_ns
       force = 0 0 {-1./2313}      # we have to divide by number of nodes in the nodeset
end

outputs
    database name = 'can_solution.e'
    disp
end

block 1
    material metal
end
block 2
    material metal_stiff
end

material metal
    density 2.5e-4
    E 10.0e6
    nu 0.27
end

material metal_stiff
    density 2.5e-4
    E 10.0e8
    nu 0.27
end