# Flexure of a viscous lithosphere - derivation

I have been trying to figure out solution for the fourth order differential equation of flexure for a viscous plate.

The fourth order differential equation for a viscous plate is given as

$$Dv*\dfrac{\partial^4}{\partial x^4}\bigg(\dfrac{\partial w}{\partial t}\bigg) + kw = 0,$$

here, $Dv$ is a function of viscosity and elastic thickness, and $k$ is the specific weight. My solution so far is $$w = w_0*\exp(t/b)*\exp(-x/\alpha)*\cos(x/\alpha),$$ where, $w_0$ is the maximum deflection, $b$ is function of $Dv$, $k$ and $\alpha$; and $\alpha$ is function of $Dv$ and $k$.

The problem with my derivation is I am not sure if it correct and I do not have a way to check it. Moreover, the $b$ term is an issue as it resolves to constant integer, which is not physically likely. Anyone out there, who could help me get on the right path? Any papers, any links?

• Please review my changes to your markup to make sure I haven't botched up the equations as you intended. – casey May 12 '15 at 13:28
• Can you post your derivation? – Antonio Mar 13 '16 at 1:17

For example, as you have written your solution, the $b$ and the $\alpha$ are constant over time ($t$) and location ($x)$. First substitute your solution into the differential equation, and carry out the differentiation to confirm that it is in fact a solution. If it can be verified as a solution, turn your attention to the dimensions of b and alpha. The dimension of b must be time and the dimension of alpha must be length. You are modeling the system as a viscoelastic solid, so the the initial condition ($t=0$) for your displacements $w(0,x)$ are just the elastic solution, and the time-dependent term is unity. Use the elastic solution to solve for alpha (the characteristic wavelength) because this does not change over time. Name the purely elastic displacement '$e$' and it has dimension of length.
You know you want to include viscosity somewhere, so let's try the kinematic viscosity $v$ (dimensions $M*M/s$). So try $b = -(L * L)/v$, or $b = -(L * 1)/v$, because this is a 2D problem and can be one unit deep. You'll need the negative sign to make the displacement grow faster for a smaller viscosity.
What could $L$ be? I would try $L=e$, or $L=e/2$. Because you really don't know the kinematic viscosity of the lithosphere, the factor is not critical.
Substitute $b$ and $\alpha$ into your solution. Does the result seem to make physical and geological sense?