In origami structures, creases have been modeled as discrete hinges [1]}, [2]}, [3]}, [4]}, [5]}, [6]}, smooth folds with \(C^1\)  continuity [7]}, [8]}, [9]} and as thinner structural elements [10]}, which all need a careful specification of matching conditions between the creases and the joining facets. Recently, Jules et al. used the Heaviside feature of a hyperbolic tangent function to describe the local geometry of creases as \(C^{\infty }\)  continuity and studied the mechanics of creased elastica [11]}.
