[研究]MESAP的单元更新_MCFT单元理论的实现

最近很少上网，专心做一个有限元程序员。

研究剪力墙和连梁构件的非线性抗剪行为是我的博士论文的另一个主要内容，由于MVLEM在分析剪压耦合效应的构件存在不足，于是我在博士论文中的研究手段引入MCFT单元理论，即Modified Compression Field Theory。这个单元的好处在于可以通过简单的应变转换求出双轴应变应力关系。程序也比较容易编制了，由于水平有限，所以MCFT是比较好的选择，主要参考多伦多大学Vecchio教授多年的研究成果。引入MCFT到MESAP，使MESAP又增加了一些分析的功能，这部分功能对于构件分析还是有很大的意义的。
====================================================================================
MESAP增加功能如下：
（1）单元：增加MCFT单元（矩形一阶）
（2）材料：Concrete03（不考虑双压强化），Concrete04（双压强化，Kupfer二维本构相关曲线）
（3）求解器：基于割线刚度的力控制的迭代算法
（4）后处理：显示MCFT单元的应变角，应变分布图
====================================================================================
相关理论与单元适用性的论文可能往些时间整理及发表，先给几个截图给大家看看。作为山塞有限元程序，请大家给点支持，拍砖。

（1）MESAP中MCFT的模型图，以后开发接口直接读取SAP2000模型，现在这个是半自动读SAP2000模型的，未全自动化。


（2）计算后显示变形与裂缝。

 
（3）显示应变图（由于一阶单元，应变单元内等值），（4）为方便研究，建立MESAP的模型Auto Builder，这是一种新型的开发程序的功能，说白了，就是傻瓜式自动生成常用模型，如剪力墙构件。这个功能在SAP2000也体现了。




在未调整参数下能够较为准确地描述剪力墙剪坏时的荷载与变形，结果算比较满意的。
师夷之长以自强，希望一步步完善MESAP，能够实现一些实用的弹塑性分析与理论研究。