I am glad to share the developed computational codes during our research. Please cite the corresponding paper(s) if you use any of these codes. Feel free to contact me if you have any questions.

1. Cohesive zone model

Cohesive elements are commonly used in modelling adhesives, bonded interfaces, composite laminates. The constitutive response of these elements usually assume a traction-separation law of the interface. Just like assuming the cohesive element as a linear-elastic spring. It will be able to pick up loading in mode I (open) and mode II (shear). Damage initiates when the assumed criteria is reached (e.g. Maximum stress, Quadratic nominal stress criterion). The damage evolution will be governed by a defined fracture toughness value (e.g. critical energy release rate Gc).

Code link: Please contact me for more questions (


[1] W. Tan, B.G. Falzon, L.N.S. Chiu, M. Price, Predicting low velocity impact damage and Compression-After-Impact (CAI) behaviour of composite laminates, Compos. Part A 71 (2015) 212-226.

[2] W. Tan, and E. Martinez-Paneda, Phase field predictions of microscopic fracture and R-curve behaviour of fibre-reinforced composites. Compos. Sci. Technol. In press.

2. Crystal plasticity model for composite laminates

By analogy, we developed a micromechanical model to capture the matrix shearing and fibre rotation of CFRP under finite strain and different strain rates, inspired by the crystal plasticity theory. Strain rate dependency of the shear modulus and yield strength of matrix is modelled through scaling functions. We then validate the predictive capabilities of this micromechanical model against the measured stress-strain responses of unidirectional (UD) and cross-ply composite laminates.  

Codes will be shared for collaboration. Please contact me if you are interested (


[1] Tan, W. and Liu, B., 2020. A physically-based constitutive model for the shear-dominated response and strain rate effect of Carbon Fibre Reinforced composites. Composites Part B: Engineering, p.108032.