Auto-adaptivity fails in 3D after first mesh refinement
Affects | Status | Importance | Assigned to | Milestone | |
---|---|---|---|---|---|
DOLFIN |
Fix Released
|
High
|
Marie Rognes |
Bug Description
from dolfin import *
# Load mesh and define function space
ou = 1
mesh = Box(0,-
mesh.order()
mesh.init(2)
V = VectorFunctionS
# Sub domain for clamp at left end
def left(x, on_boundary):
tol = 1.0e-14
return abs(x[0]) < tol and on_boundary
# Sub domain for shear at right end
class RightBoundary(
def inside(self, x, on_boundary):
tol = 1.0e-14
return abs(x[0]-48) < tol and on_boundary
# Mark boundary parts
domain_parts = MeshFunction(
boundary_parts = MeshFunction(
dOmega_t = RightBoundary()
dOmega_
# Define variational problem
v = TestFunction(V)
u = TrialFunction(V)
f = Constant((0.0, 0.0, 0.0))
tg = Expression(
E = 30000.0
nu = 0.25
rho = 1
mu = E / (2.0*(1.0 + nu))
lmbda = E*nu / ((1+nu)*(1.0 - 2*nu))
def sigma(v):
return 2.0*mu*sym(grad(v)) + lmbda*tr(
a = inner(grad(v), sigma(u))*dx
b = inner(rho*v,u)*dx
L = inner(v,
L.exterior_
# Define goal
u = Function(V)
g = u[1]*ds(0)
g.exterior_
u.name = "Displacements"
# Set up boundary condition at left end
c = Constant((0.0, 0.0, 0.0))
bcl = DirichletBC(V, c, left)
# Set up PDE and solve
eps = 1.e-2
problem = VariationalProb
problem.
Changed in dolfin: | |
status: | Confirmed → Fix Committed |
Changed in dolfin: | |
status: | Fix Committed → Fix Released |