Tag Info

New answers tagged

2

NDSolve is complaining, because v[0, t] == u[1, t] appears to link boundary conditions at x == 0 and x == 1. This can be circumvented by redefining x as 1 - x for v, so that the code becomes, {su, sv} = NDSolveValue[{ D[u[x, t], t] == D[u[x, t], x, x] - u[x, t], D[v[x, t], t] == D[v[x, t], x, x] - v[x, t], u[x, 0] == 1, v[x, 0] == ...


3

Here is something to get you started: The idea is to use a 2D cross section of the orifice plate in the xz-direction (not the xy-direction, as then you could not apply the surface force). a = 10*10^-3; b = 5*10^-3; \[Nu] = 1/3; p0 = 0.1*10^6; Ey = 200*10^9; h = 1*10^-3; De = (Ey h^3)/(12 (1 - \[Nu]^2)); planeStrain = {Inactive[ Div][{{0, ...


1

This works. a = 10*10^-3; b = 5*10^-3; ν = 1/3; p0 = 0.1*10^6; Ey = 200 *10^9; h = 1*10^-3; De = (Ey h^3)/(12 (1 - ν^2)); eqn = w''''[r] + (2/r) w'''[r] - (1/(r^2)) w''[r] + (1/(r^3)) w'[ r] == -p0/De w1 = NDSolveValue[{eqn, DirichletCondition[w[r] == 0, r == a], DirichletCondition[w'[r] == 0, r == a], DirichletCondition[ ...



Top 50 recent answers are included