【想法】

各位老师，同学，大家下午好！

从学期中到现在，我主要做了一下两方面的工作。

首先是水哥的简化单自由度系统计算程序的开发。

实践表明，在熟悉原文计算部分的逻辑结构的前提下，这类程序的编写速度还是比较快的。

虽然一开始写错了程序，加之原文本身也有点问题，但整体所花费的时间不是很多，到12月初也就基本完成了。

第二是张老师的项目。

我主要负责极限强度和剩余强度校核部分。

即规范第五章第一节的前半段、第二节、第三节以及附录二。

其中附录二讲的是应用逐步破坏法计算极限和剩余承载力。

剩余部分主要包括横剖面定义、损伤定义、极限强度和剩余强度校核的载荷、系数和判据。

图中左下部分的曲率需要用到郭德松负责的第一节的后半部分的数据。

除核心程序以外的部分在12月初完成了初稿。

核心程序部分，主要有单元划分和逐步破坏法计算等部分。

根据规范要求，并考虑到与项目主程序的对接和自己的编程水平，最后决定采用以节点和板格为中心的编程思路。

程序的几大基本输入对象如左图所示，包括节点、加强筋、面板、板格和特别硬角单元。

而基本的结构单元则有硬角单元、加强筋单元和加筋板单元。

程序中，首先需要分别将节点与加强筋、面板和板格进行绑定。

其中，前两者的较为简单，而板格作为确定各单元范围的关键对象，处理方式较为复杂，由于时间关系，这里不做展开。

总而言之，经过对节点对象的遍历，可以得到右图中除加筋板单元和特别硬角单元以外的所有单元。

在此基础上，对板格进行遍历，得到加筋板单元。

特别硬角单元单独输入，单独处理。

得到以上三类六项单元，便完成了单元划分部分。

逐步破坏法计算部分，需要在每一曲率和中性轴高度下下进行，且需遍历所有结构单元。

根据前一步计算得到的单元重心位置、平均最小屈服应力，求得单元的相对应变。

再结合其他属性和条件，依照右图所示的单元失效模式和相关的一系列公式求得单元的应力、轴力和弯矩。

在满足截面力的平衡和中性轴高度的精确度要求后，得到当前曲率下的弯矩值。

连续输出可得到整条曲率－弯矩曲线，其峰值即为结构的极限承载力或剩余承载力。

计算完成。

这个项目还没有完成，还需要花费一段时间去处理各项问题。

因此，在后续规划中，首先需要完成核心程序的编写，修订除核心程序以外部分的需求文件初稿，并编制程序文档。

另一方面，加快论文进度，多看文献，继续编程，引入中性轴偏转算法，对比IMO和GOALDS等损伤定义下的剩余强度，并研究结构和材料参数对极限强度的影响。