Calcule la deformación de una tensión de plano 

Defina un operador de tensión de plano.

In[1]:=

X planeStrainOperator[ Y_, \[Nu]_] := {Inactive[ Div][({{0, -((Y \[Nu])/((1 - 2 \[Nu]) (1 + \[Nu])))}, {-(Y/( 2 (1 + \[Nu]))), 0}}.Inactive[Grad][v[x, y], {x, y}]), {x, y}] + Inactive[ Div][({{-((Y (1 - \[Nu]))/((1 - 2 \[Nu]) (1 + \[Nu]))), 0}, {0, -(Y/(2 (1 + \[Nu])))}}.Inactive[Grad][ u[x, y], {x, y}]), {x, y}], Inactive[ Div][({{0, -(Y/(2 (1 + \[Nu])))}, {-(( Y \[Nu])/((1 - 2 \[Nu]) (1 + \[Nu]))), 0}}.Inactive[Grad][ u[x, y], {x, y}]), {x, y}] + Inactive[ Div][({{-(Y/(2 (1 + \[Nu]))), 0}, {0, -(( Y (1 - \[Nu]))/((1 - 2 \[Nu]) (1 + \[Nu])))}}.Inactive[ Grad][v[x, y], {x, y}]), {x, y}]};

Resuelva la ecuación en una región rectangular. En el lado izquierdo ambas variables dependientes se mantienen fijas con desplazamiento 0 y en el lado derecho con una carga límite de una unidad en la dirección de ejes y negativos.

In[2]:=

X {uif, vif} = NDSolveValue[{planeStrainOperator[10^3, 33/100] == {0, NeumannValue[-1., x == 5]}, DirichletCondition[{u[x, y] == 0., v[x, y] == 0.}, x == 0]}, {u, v}, {x, 0, 5}, {y, 0, 1}];

Visualice la deformación en las direcciones x y y usando el paquete de elemento finito.

muestre la entrada completa de Wolfram Languageoculte la entrada

In[3]:=

X Needs["NDSolve`FEM`"] mesh = uif["ElementMesh"]; Show[{ mesh["Wireframe"[ "MeshElement" -> "BoundaryElements"]], ElementMeshDeformation[mesh, {uif, vif}][ "Wireframe"[ "ElementMeshDirective" -> Directive[EdgeForm[Red], FaceForm[]]]]}, ImageSize -> 300]

Out[3]=