This paper presents a method for the efficient simulation of a large number of plant leaves, realistically showing their weathering process in autumn time, with various deformed shapes. In modeling various shapes of leaves for weathering, the authors introduce a botanically and physically based method is introduced for simulating variations of leaf shape by using a double layered model (DLM). This model takes advantage of the key observation that the withered shapes of leaves are caused by different shrinkages between the mesophyll cells and vein cells during the aging process of leaves. By setting up the mass-spring system in two interconnecting layers, different kinetic manners of mesophyll and veins can be simulated respectively. Meanwhile the double layers affect each other and determine together the deformed shape of leaves.