Notepad++中OpenSEES钢抗弯框架Tcl建模报错排查修复

建模背景 #

• 工具：使用Notepad++编辑器编写Tcl语言OpenSEES代码，该编辑器支持错误行快速定位、未闭合括号匹配，操作便捷
• 建模目标：参照钢抗弯框架（MRF）示意图建立梁端转动弹簧模型，建模过程不考虑RBS（梁翼缘削弱段）
• 建模逻辑：目标结构为10层大尺寸钢框架，采用程序化批量建模，代码中j代表楼层编号，i代表柱/跨编号
• 故障现象：代码运行触发域警告，多次调整后仍无法修复

原始问题代码 #

for {set j 1} {$j <= $nFlrs} {incr j} {
    
    for {set i 1} {$i <= $nBays} {incr i} { 
        set ID1 [expr 4*1000+($i*10+$j)*10+1]
        set nodeR1 [expr ($i*10+$j)*100+5]
        set nodeC1 [expr ($i*10+$j)*10+1]
    }   
    for {set i 1} {$i <= $nBays+1} {incr i} {
        set ID2 [expr 4*1000+($i*10+$j)*10+2]
        set nodeR2 [expr (($i+1)*10+$j)*100+10]
        set nodeC2 [expr (($i+1)*10+$j)*10+2]
    }
}
    uniaxialMaterial Bilin $ID1 $ke $alfah $alfah $my -$my $lambda 0\
        0 0 1 0 0 0 $tetap $tetap $tetapc $tetapc 0 0 $tetau $tetau 1 1 $nFactor
    uniaxialMaterial Bilin $ID2 $ke $alfah $alfah $my -$my $lambda 0\
        0 0 1 0 0 0 $tetap $tetap $tetapc $tetapc 0 0 $tetau $tetau 1 1 $nFactor
    
    # define beam springs   
    element zeroLength $ID1 $nodeR1 $nodeC1 -mat $ID1 -dir 6
    element zeroLength $ID2 $nodeR2 $nodeC2 -mat $ID2 -dir 6

    #Constrain the translational DOF with a multi-point constraint
    #                   retained constrained DOF_1 DOF_2 ... DOF_n
           equalDOF    $nodeR1     $nodeC1     1     2
           equalDOF    $nodeR2     $nodeC2     1     2
}
incr ID1
incr ID2

报错信息 #

warning: WARNING Domain::addElement - In element 4401 no Node 4005 exists in the domain
中文翻译：添加编号为4401的单元时，计算域内不存在编号为4005的节点

问题根因 #

• 循环作用域错误：材料定义、零长度单元定义、多点约束代码全部放在两层for循环外部，循环执行过程中只会反复覆盖变量值，不会执行建模操作，最终仅会在循环结束后用最后一组变量值创建1组弹簧，和逐层逐跨创建的预期逻辑完全不符。
• 节点编号规则不匹配：报错提到的4005节点由nodeR1 = ($i*10+$j)*100+5公式计算得出，属于千位级编号，若前期梁柱节点建模时未采用该规则生成对应编号的节点，自然会触发节点不存在的报错。
• 语法错误：代码末尾存在多余的闭合大括号}，会触发Tcl解析错误。
• ID冲突：材料编号与零长度单元编号使用完全相同的计算规则，会出现同ID对应不同类型对象的问题，极易引发OpenSEES对象索引错误。
• 冗余代码：循环外的incr ID1、incr ID2无实际作用，属于无效代码。

修正后代码框架 #

# 需提前定义全局参数：nFlrs、nBays、ke、alfah、my、lambda、tetap、tetapc、tetau、nFactor
for {set j 1} {$j <= $nFlrs} {incr j} {
    # 逐跨生成梁左端转动弹簧
    for {set i 1} {$i <= $nBays} {incr i} { 
        # 区分材料ID与单元ID，避免重名冲突：材料ID以4开头，单元ID以5开头
        set matID1 [expr 4*1000+($i*10+$j)*10+1]
        set eleID1 [expr 5*1000+($i*10+$j)*10+1]
        # 注意：此处nodeR1、nodeC1的计算公式需和前期梁柱节点定义的编号规则完全对齐
        set nodeR1 [expr ($i*10+$j)*100+5]
        set nodeC1 [expr ($i*10+$j)*10+1]

        # 建模操作移入循环内部，逐组生成对象
        uniaxialMaterial Bilin $matID1 $ke $alfah $alfah $my -$my $lambda 0\
            0 0 1 0 0 0 $tetap $tetap $tetapc $tetapc 0 0 $tetau $tetau 1 1 $nFactor
        element zeroLength $eleID1 $nodeR1 $nodeC1 -mat $matID1 -dir 6
        equalDOF $nodeR1 $nodeC1 1 2
    }   
    # 逐跨生成梁右端转动弹簧
    for {set i 1} {$i <= $nBays+1} {incr i} {
        set matID2 [expr 4*1000+($i*10+$j)*10+2]
        set eleID2 [expr 5*1000+($i*10+$j)*10+2]
        set nodeR2 [expr (($i+1)*10+$j)*100+10]
        set nodeC2 [expr (($i+1)*10+$j)*10+2]

        uniaxialMaterial Bilin $matID2 $ke $alfah $alfah $my -$my $lambda 0\
            0 0 1 0 0 0 $tetap $tetap $tetapc $tetapc 0 0 $tetau $tetau 1 1 $nFactor
        element zeroLength $eleID2 $nodeR2 $nodeC2 -mat $matID2 -dir 6
        equalDOF $nodeR2 $nodeC2 1 2
    }
}

内容的提问来源于stack exchange，提问作者sahar ebrahimi