Plainly, this work's application very nice and usefully for several civil, mechanical etc.. engineering areas.
I think, this file documantation, description and program-running-algorithm's flow not-regular. I know epoxi-laminate composite plate system with analytical and finite element solution technique but all-user don't know this theory and method.
If you use one-main program instead of more-sub programs (with input and output veriables) than this program and work be more likeable.
I was tested your program running algorithms to several thin and thick plate model for true-solution and shear locking effect. But I don't see shear locking program module or procedure in main-program.
Analysis results, different the Reissner-Mindlin thick plate model.
I think this files need more improvement but I know this area's works are very hard and require the more works-time.
Good works, Best Regard