1
#多工况线性插值
2
#载荷文件格式.csv,每个值以","分隔，以"/n"换行
3
#第一列X坐标,第二列Y坐标,第三列Z坐标,第四列及以后为各个工况的载荷值,并且每一列对应一个工况
4
tk_messageBox -message "开始多工况线性插值-每个工况生成一个CSV文件"
5
#选择结果文件保存路径,并设置为全局变量
6
set ::ResultDirectoryFilePath [tk_chooseDirectory -title "选择结果文件保存路径"]
7
set file [tk_getOpenFile -title "选择工况文件excel文件" ]
8
set load_data [open "$file" r]
9
while {![eof $load_data]} {
10
    #gets 一行数据，遇到换行符停止
11
    gets $load_data load_data_line_all
12
    set line_load [split $load_data_line_all ","]
13
    set line_x [lindex $line_load 0]
14
    set line_y [lindex $line_load 1]
15
    set line_z [lindex $line_load 2]
16
    set m [llength $line_load]
17
    if {!($line_x eq "")} {
18
        for {set i 1} {$i <= $m - 3 } {incr i} {
19
        #a+一直打开这个文件，追加写入
20
        set aload_case [open "$::ResultDirectoryFilePath/case_$i.txt" a+]
21
        set line3 [lindex $line_load [expr $i+2]]
22
        #puts 写文件会自动添加换行符
23
        puts $aload_case "$line_x,$line_y,$line_z,$line3"
24
        flush $aload_case
25
        close $aload_case
26
        }
27
    }       else {
28
                continue
29
        }
30
    }
31
close $load_data
32
tk_messageBox -message "完成load文件生成多个txt文件"

1
proc start_interpolation {loadstepName loadstepID} {
2
    eval set txt_num  [llength [glob -nocomplain -directory $::ResultDirectoryFilePath *.txt]]
3
    for {set i 1 } {$i <= $txt_num } {incr i} {
4
    #加载单元，加载方向为单元法向，可以通过createmark的方式改变加载单元集合
5
    *createmark elements 1 "displayed"
6
    #线性插值依赖已经存在的load和载荷点位，需要空的load
7
    *createmark loads 2
8
    *createentity loadcols id=[expr $i + $loadstepID] name=[concat $loadstepName 1.0_ $i]
9
    #hypermesh中fill gap关掉，才能输入负值，当为负值时，所有单元都要进行插值（包括载荷点外的所有单元）
10
    *BCM elements 1 2 4 1 "$::ResultDirectoryFilePath/case_$i.txt" -50 1 30
11
    *createmark loads 1 all
12
    *hideentitybymark loads 1
13
    }
14
    tk_messageBox -message "完成多工况线性插值"
15
}