Initial conditions for a mixed element

1

I am trying to set initial conditions (initial displacement) for a mixed element (from the von-karman demo). The following gives me a warning " *** FFC warning: evaluate_dof(s) for enriched element not implemented." but the rest of my script still runs.
Although there are no errors, the initial conditions are not actually applied when I print or visualize my results. Any help would be appreciated.
Thanks

P1 = FiniteElement("Lagrange", triangle, degree = 1)
P2 = FiniteElement("Lagrange", triangle, degree = 2)
bubble = FiniteElement("B", triangle, degree = 3)
enriched = P1 + bubble

element = MixedElement([VectorElement(P2, dim=2), P2, VectorElement(enriched, dim=2)])

e_v0 = Expression(("0.36*x[0]", "0.36*x[1]") ,degree=1)
e_w0 = Expression(("0*x[0]") ,degree=1)
e_theta0 = Expression(("0*x[0]", "0*x[1]") ,degree=1)
v0 = interpolate(e_v0,Q.sub(0).collapse())
w0 = interpolate(e_w0,Q.sub(1).collapse())
theta0 = interpolate(e_theta0,Q.sub(2).collapse())
q_, q, q_t = Function(Q), TrialFunction(Q), TestFunction(Q)
assign(q_,[v0,w0,theta0])
v_, w_, theta_ = split(q_)
q_old = Function(Q)
assign(q_old,[v0,w0,theta0])
v_old, w_old, theta_old = split(q_old)
2

Try projecting your initial condition instead of interpolating it.

3

Thanks, this solved the warning but the initial conditions are still not actually applied when I see the results. Could this be an issue somewhere else in my script?

4

As you have not supplied a minimal working example, i cant really tell you if you are doing anything wrong.