We study the morphology of lipid membranes subjected to intra-membrane viscous flows and interactions with elliptical cylinder substrates. From the non-linear theory of elastic surfaces, a linearized shape equation and admissible boundary conditions are formulated in elliptical coordinates via the Monge representation of a surface. In particular, the intra-membrane viscosity terms are linearized and mapped into elliptic coordinates in order to accommodate more general forms of viscous flow. The assimilated viscous flow is characterized by potential functions which satisfies the continuity condition. A complete solution in terms of Mathieu function is then obtained within the prescription of incremental deformations superposed on large. The results describe smooth morphological transitions over the domain of interest and, more importantly, predicts wrinkle formations in the presence of intramembrane viscous flow in the surface. Lastly, the obtained solution accommodates the results from the circular cases in the limit of vanishing eccentricity and intra-membrane viscous flow.The mechanics of lipid membranes has consistently been the subject of intense research for its importance in the understanding of a wide variety of essential cellular processes 1-5 . Traditionally, it was believed that cells are surrounded by a thin oil-based barrier, yet the structure of this membrane was not well understood. In 1920s, E. Gorter and F. Grendel. 6 found that the cell membrane is composed of lipid molecules (phospholipids) which are generally divided into two important groups: the hydrophilic head parts and hydrophobic tail groups. When dispersed into aqueous solutions, lipid molecules are driven by hydrophobic effects to form a unique bilayer structure (a lipid bilayer) with opposing orientations, that maintain symmetry about a mid-surface. In fact, with the advances in electron microscopy, the bilayer structure was identified as a characteristic of all biological membranes (biomembranes) 7 . Since they are negligibly thin (typically 5-10 nm), and fragile, the study of various aspects of lipid bilayers is often achieved by employing mathematical models in order to overcome the formidable difficulties of experimental studies. Also, from a mechanical perspective, the responses of a lipid membrane can be idealized as a thin elastic film. Within this context, the development of theoretical models describing the behavior of lipid bilayers has greatly benefited from the differential geometry of a surface and the theory of an elastic surface such that the deformation of energy of a thin membrane can be expressed by the mean and Gaussian curvatures of a surface 8,9 . In particular, Helfrich proposed a well-known Helfrich energy potential 1 which addresses the symmetry of lipid bilayers and further ensures the resulting equilibrium state to be energy minimizing 10 . This, together with the use of variational principle, and the virtual-work statement, furnish the Euler-Lagrange equations also known as "membrane shape equations", which h...