Hydrogels are nowadays widely used in various biomedical applications, and show great potential for the making of devices such as biosensors, drug- delivery vectors, carriers, or matrices for cell cultures in tissue engineering, etc. In these applications, due to the irregular complex surface of the human body or its organs/structures, the devices are often designed with a small thickness, and are required to be flexible when attached to biological surfaces. The devices will deform as driven by human motion and under external loading. In terms of mechanical modeling, most of these devices can be abstracted as shells. In this paper, we propose a mixed graph-finite element method (FEM) phase field approach to model the fracture of curved shells composed of hydrogels, for biomedical applications. We present herein examples for the fracture of a wearable biosensor, a membrane-coated drug, and a matrix for a cell culture, each made of a hydrogel. Used in combination with experimental material testing, our method opens a new pathway to the efficient modeling of fracture in biomedical devices with surfaces of arbitrary curvature, helping in the design of devices with tunable fracture properties.

The unique, hierarchical patterns of leaf veins have attracted extensive attention in recent years. However, it remains unclear how biological and mechanical factors influence the topology of leaf veins. In this paper, we investigate the optimization mechanisms of leaf veins through a combination of experimental measurements and numerical simulations. The topological details of three types of representative plant leaves are measured. The experimental results show that the vein patterns are insensitive to leaf shapes and curvature. The numbers of secondary veins are independent of the length of the main vein, and the total length of veins increases linearly with the leaf perimeter. By integrating biomechanical mechanisms into the topology optimization process, a transdisciplinary computational method is developed to optimize leaf structures. The numerical results show that improving the efficiency of nutrient transport plays a critical role in the morphogenesis of leaf veins. Contrary to the popular belief in the literature, this study shows that the structural performance is not a key factor in determining the venation patterns. The findings provide a deep understanding of the optimization mechanism of leaf veins, which is useful for the design of high-performance shell structures.