CSI API编程训练第二课，这次做基于虚功原理的伸臂桁架几何优化，主要训练如何通过编程控制SAP2000修改节点位置，自动提交计算分析，并提取构件的内力等。
对于桁架结构，根据虚功原理，结合单位荷载法可知，结构任意一点在指定方向的位移可按以下公式表示：
其中，
为杆件的虚拟轴力，
为杆件的真实轴力，
为杆件的弹性模量，
为杆件的截面面积，
为杆件的长度。
根据
Baker 的研究可知，对于静定桁架结构，各杆件处于等应力状态时，结构杆件是最优的。这个最优说的是，对于给定挠度，当所有杆件均处于等应力状态时，所需结构材料用量最小；或者说对于给定材料用量的结构，当所有杆件均处于等应力状态时，结构挠度最小。
假设各杆件的应力水平均达到同一个值，设为
，此时结构任意一点在指定方向的位移公式变为以下：
由上述公式可见，当
一样时，
越小，结构的位移
越小，结构的刚度越大，此时
可以看做是一个能量指标，杆件虚拟力与杆件长度的乘积。
我们不妨把上述公式中的虚拟轴力换位真实轴力，则
变为
，这就是 Baker 在另外一篇文章中提到的传力路径，也常称为力流。关于力流，这个文章也介绍过： 弘构大讲堂 | 基于力流(∑N ×L)的结构效率评价方法 。
对于不同的结构，我们可以利用
评估各个结构的几何构型，
越小，表明这个几何构型传力越直接，传力路径越短，在构件达到满应力或同等应力水平的条件下，刚度最大。
于是可以利用这点，找到静定桁架的最优几何构型。
下面我们用一个X型伸臂桁架的几何优化为例进行演示。在SAP2000中进行X型桁架建模，如下图。
其中X型伸臂桁架的宽、高均为8000，我们通过迭代寻找中间X点的最优位置。
打开
软件，设置控制节点、坐标变化范围、监测点、监测位移和迭代步等参数。
设置好参数后，点击运行分析，软件控制SAP2000进行迭代计算，可以在软件界面看到实时的迭代过程。
下图是最终分析结果。其中第一个图是迭代坐标和监测点竖向位移的关系，第二个图是迭代坐标和力流
的关系。
由第一图可以看到，位移最小的X坐标是1000，第二图力流最小的X坐标是2000。
其中第二个图的
X=2000才是最优几何构型，如下图所示，得到了我们经常看到的3/4位置。
[01] Baker WF. Energy-Based Design of Lateral Systems [J]. Structural Engineering International, 1992, 2:99-102.
[02] Baker WF, Beghini LL, Mazurek A, Carrion J, Beghini A . Maxwell’s reciprocal diagrams and discrete Michell frames [J]. Struct Multidisc Optim. (2013)
全部回复（1 ）
只看楼主 我来说两句 抢板凳