Abstract
The problem of axisymmetric deformation of a peripherally fixed and uniformly loaded circular membrane under deflection restrictions (by a frictionless horizontal rigid plate) was analytically solved, where the assumption of constant membrane stress adopted in the existing work was given up, and a closed-form solution of this problem was presented for the first time. The numerical analysis shows that the closed-form solution presented here has higher calculation accuracy than the existing approximate solution.
1. Introduction
Elastic membrane structures and components are widely used in many fields [1–7]. The large deflection phenomena of membrane problem usually give rise to nonlinear differential equations [8–11]. These nonlinear equations generally present serious analytical difficulties when applied to boundary-value problems. Due to these somewhat intractable nonlinear equations, the analytical solutions of membrane problems are available in a few cases, but in practice they are often found to be necessary.
Hencky [12] originally dealt with the problem of axisymmetric deformation of the circular membrane fixed at the outer edge under the action of a uniformly distributed transverse loads and presented the power series solution of the problem, as shown in Figure 1, where is the uniformly distributed transverse loads. A calculation error in [12] was corrected by Chien [13] and Alekseev [14], respectively. This problem and its solution are usually called well-known Hencky problem and well-known Hencky solution for short, which are often referred to or cited in a number of related studies [15–22]. However, if we use a frictionless horizontal rigid plate to restrict the deflection of the membrane in the well-known Hencky problem, as shown in Figure 2, then such a problem will probably become somewhat complicated, where is the radial coordinate, is the transversal displacement, is the radius of the membrane contacting with the frictionless rigid plate, and is the gap between the frictionless rigid plate and the initially flat membrane. So, Xu and Liechti had to use the following four assumptions to deal with this problem [23]: the membrane has negligible flexural rigidity and only membrane stresses are considered; the slope angle of membrane is so small that the condition could approximately hold, that is, the so-called small-rotation-angle assumption; a constant radial stress is assumed, that is, the radial stress in the membrane under loads has nothing to do with the radial coordinate; the contact between the membrane and the rigid plate is frictionless. Obviously, assumption above seems to be too harsh, and, in fact, it can be given up, as will be seen later.
In the next section, the problem shown in Figure 2 was analytically solved, where the assumption of constant membrane stress adopted in the existing work was given up, and the closed-form solution of this problem was presented. In Section 3, based on numerical calculations, the reliability of the presented closed-form solution is verified and the beneficial effect of giving up the assumption of constant membrane stress was discussed. Section 4 is the concluding remarks.
2. Membrane Equation and Its Solution
Suppose that, an initially flat, linearly elastic, rotationally symmetric, taut circular membrane with Young’s modulus of elasticity , Poisson’s ratio , thickness , and radius is fixed at the outer edge and uniformly distributed transverse loads are quasi-statically applied onto the membrane surface (as shown in Figure 1, the well-known Hencky problem), and the applied loads especially have reached a value large enough so that the deflected membrane has a contact with a frictionless rigid plate being parallel to the initially flat membrane (as shown in Figure 2). Such a problem can be viewed as consisting of the two local problems in the central portion of and in the annular portion of , which are connected by the continuity conditions at , where the problem in may be simplified as a plane tension or compression problem of membrane, while the problem in is still a problem of membrane deflection. The following assumptions are here made in order to reach a closed-form solution: during deformation the variation in thickness of membrane could be ignored, as is seen usually in membrane problems; the so-called small-rotation-angle assumption reported in the existing literatures [18–20] is still adopted here, that is, the slope angle of membrane is so small that the condition could approximately hold; the contact between the deflected membrane and the rigid plate is frictionless. Here, we give up the assumption of the constant radial stress reported in [23]. The continuous condition at and the boundary conditions at will be applied during the process of solution. As for the case prior to the contact between the deflected membrane and the frictionless rigid plate, that is, the well-known Hencky problem, it has been dealt with and its solution may be found in [12, 17].
Let us consider firstly the problem in the annular portion and take a piece of the central portion of the annular membrane whose radius is ; with a view of studying the static problem of equilibrium of this membrane under the joint action of the uniformly distributed loads , the membrane force acted on the boundary and the reaction force from the rigid plate, just as it is shown in Figure 3, where is the radial stress and is the slope angle of membrane. Right here there are three vertical forces, that is, the total force (in which ) of the uniformly distributed loads , the total reaction force from the rigid plate, and the total vertical force which is produced by the membrane force .
The out-plane equilibrium equation iswhereSubstituting (2) into (1), one hasThe in-plane equilibrium equation iswhere is the circumferential stress. The relations of strain and displacement of the large deflection problem arewhere is the radial strain, is the circumferential strain, and is the radial displacement. The relations of stress and strain areSubstituting (5a) and (5b) into (6a) and (6b),By means of (4), (7a), and (7b), one hasIf we substitute the of (8) into (7a), the compatibility equation may be written asThe detailed derivation from (4) to (9) may be obtained from any general theory of plates and shells. It is not necessary to discuss this problem here.
Obviously, in the central portion of , the membrane is always in a plane state of radial tensile or compression; in other words, it is in fact only a plane problem with a salient characteristic of . After substituting into (5a) and (5b), it may be found thatSubstituting (10a) and (10b) into (6a) and (6b),Substituting (11a) and (11b) into (4),So, considering the conditionsthe solution of (12) may be written asSubstituting (14) into (10a), (10b), (11a), and (11b), it may finally be found that
Let us introduce the following nondimensional variables:and transform (3), (9), (4), and (8) intoThe boundary conditions, under which (17) and (18) may be solved, are
where the subscripts and denote the values of various variables on two sides of the interconnecting circle . The side of region is under the plane state of radial tensile or compression, while the side of region is under the deflection state of membrane. Eliminating from (17) and (18), we obtain a differential equation which contains only :Letting and expanding into the power series of , one hasAfter substituting (24) into (23), the coefficients can be expressed into the polynomial of the undetermined constants and , which are shown in Appendix A. Further, expanding into the power series of ,and substituting (24) and (25) into (17), the coefficients can also be expressed into the polynomial of and , which are shown in Appendix B, while is left as another undetermined constant, except and . With the given and all the undetermined constants can be determined by applying the continuous condition at , as follows.
From (20) and (24), (21a) and (22a) giveFrom (25), (21b) and (22b) giveFrom (24), (22c) givesEliminating the from (27) and (30), one hasHence, for the given problem where , , , , , and are known in advance, the undetermined constants and can be determined by (26) and (31), then the coefficients and can also be determined, and consequently the undetermined constant can be determined by (28). After this, the nondimensional variable can be determined by (29), that is, the applied loads q, corresponding to this , can be determined. Thus, the problem dealt with can be solved.
3. Results and Discussions
The undetermined constants , , and depend, in fact, on only and , as seen above. For the convenience of application, we here present the variation of the undetermined constants , , and with , as shown in Figures 4–6, where takes 0.1, 0.2, 0.3, 0.4, and 0.5, respectively.
From the process of solution above we may see that, after giving up the assumption of the constant radial stress reported in [23], the solution to the problem dealt with becomes really complicated. Generally, the structural response is derived from an action of structure, but here the action, the applied loads , has to be derived from the response, the contact radius of the membrane with the frictionless rigid plate.
For verifying the reliability of the closed-form solution presented here, we made a comparison between the solution presented here and the well-known Hencky solution [12, 17], based on the numerical example of a polyethylene terephthalate thin-film (considered in [23]) with 10 mm, 0.003 mm, 0.8 mm, 4.65 GPa, and . When the contact radius is equal to 0.0001 mm, the action, the applied loads , is about 2.69 KPa. Figure 7 shows the variation of with , where the solid line represents the result obtained by the solution presented here, and the dashed line represents the result obtained by the well-known Hencky solution. From Figure 7 it may be seen that the two profiles are very close to each other; this shows that the closed-form solution presented here is, to some extent, reliable.
Further loading, when the contact radius is equal to 1 mm, the applied loads are about 2.998 KPa. Figure 8 shows the variation of the radial stress with the radial coordinate , where the solid line represents the result obtained by the solution presented here, and the dashed dotted line represents the result obtained by the solution presented in [23] (the case without residual stress ). This survey found a wide variation in the radial stresses , which indicates that the assumption of the constant radial stress adopted in [23] should be given up.
Figure 9 shows the variation of the contact radius with the applied loads , where the solid line represents the result obtained by the solution presented here, and the dashed dotted line represents the result obtained by the solution presented in [23] (the case without residual stress ). As we can see from Figure 9, the dashed dotted line shows that the deflected membrane will start touching the frictionless rigid plate at KPa, while it starts touching the solid line at KPa. This difference is caused mainly by the assumption of the constant radial stress adopted in [23].
4. Concluding Remarks
In this study, the problem of axisymmetric deformation of a peripherally fixed circular membrane under uniformly distributed transverse loads was analytically dealt with under the condition of deflection restrictions, where the assumption of constant membrane stress adopted in [23] was given up, and the closed-form solution of this problem was presented for the first time. As far as well-known Hencky problem is concerned, the problem dealt with here can be understood as a new circular membrane problem under the combined loads, that is, the uniformly distributed transverse loads and the reaction force produced by the frictionless horizontal rigid plate.
The work presented here makes a significant and new contribution to thin-film mechanics. It ought to have significant implication in the mechanical characterization of film/substrate surface and interface and the interpretation of film/substrate delamination experiment; it could be incorporated especially into the study work reported in [23]. Further studies are expected to focus on giving up the so-called small-rotation-angle assumption adopted usually in membrane problems, in order to reach a closed-form solution with higher calculation accuracy and wider application scope.
Appendix
A. The Coefficients
Letting , then the coefficients may be written as
B. The Coefficients
Letting , then the coefficients may be written as
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was funded by the National Natural Science Foundation of China (Grants nos. 11772072 and 11572061).