Coupling a simple nonlinear membrane to point displacement Michael Alletto

• contributor: Michael Alletto
• affiliation: WRD MS
• contact: click here for email address
• OpenFOAM versions: v2112
• Published under: CC BY-NC-ND license (creative commons licenses)

Go back to Collection by authors.

You can download the case https://gitlab.com/mAlletto/openfoamtutorials/-/tree/master/membranBCSend] here.

Introduction

In the tutorial https://wiki.openfoam.com/Coupling_a_membrane_with_pretension_to_point_displacement_Michael_Alletto we saw how to couple different physical models. In the before mentioned tutorial we coupled the equations of motion of a prestressed membrane with the displacement of the mesh. The membrane is deformed by the action of the pressure on the solid surface. The equations of motion for a prestressed membrane assume that the stress of the membrane is constant and does not change when the membrane is deformed. The question which arises is what is if we want to simulate a membrane which is in an initial undeformed and unstressed configuration and deforms by the action of the fluid pressure? How should we proceed to figure out what are the correct equations to be solved, what kind of assumptions has to be made to simplify the equation and how do we verify if our assumptions are correct? In this tutorial we will answer this question.

Equation of motion

The equations of motion for a curved thick membrane under the action of external forces can be found in Reddy [3] p. 421:

The above equation is described in a curvilinear coordinate system. x₁ and x₁ represent curved coordinates which follow the surface. N₁₁ and N₁₁ are the normal (membrane) forces per unit length. N₁₂ represents the shear force parallel to the curved surface and Q₁ and Q₂ represent the shear forces normal to curved surface. R₁ and R₂ are the curvature radii and f₃ is the external force acting normal to the plane. u₁ and u₂ are the displacements parallel to the curved surface and $w$ is the displacement normal to the curved surface. I₀ = &rho h is the mass per unit surface. a₁ and a₂ are the scale factors (see e.g. https://mathworld.wolfram.com/ScaleFactor.html).

In our case we want only to consider the effect of the pressure on the membrane motion. Since the pressure always acts in surface normal direction, it reasonable to assume that the displacements parallel to the surface u₁ and u₂ can be neglected. For thin membranes the transverse shear forces Q₁ and Q₂ are usually set to zero (see \cite{reddy2006theory}). Since the pressure acts in surface normal direction is it also reasonable to assume that the membrane is not rotated parallel to the surface and hence also the shear force N₁₂ can be neglected. With this assumption we obtain the following equation:

References

[1] Fadl Moukalled, L Mangani, Marwan Darwish, et al. The finite volume method in computational fluid dynamics, volume 113. Springer, 2016.

[2] Draga Pihler-Puzović, Anne Juel, Gunnar G Peng, John R Lister, and Matthias Heil. Displacement flows under elastic membranes. part 1. ex- periments and direct numerical simulations. Journal of Fluid Mechanics, 784:487–511, 2015.

[3] Junuthula Narasimha Reddy. Theory and analysis of elastic plates and shells. CRC press, 2006.

[4] Stephen Timoshenko, Sergius Woinowsky-Krieger, et al. Theory of plates and shells, volume 2. McGraw-hill New York, 1959.