FEAP输出matlab刚度矩阵的代码

snowwave02

最近接触FEAP，发现有时需要将FEAP输出matlab格式的刚度矩阵。
其实在FEAP的官网上有类似的代码，但可能只能在FEAP中使用，可惜论述上的FEAP没法编译通过，所以只能用FEAPpv
自己在FEAPpv加了一段代码。
如下：
c       output tang to file
        if (1) then
            if(isw.eq.3) then
                open(unit = 12,file = 'tang',status = 'replace')
                rewind(12)
                do i = 1,nst ! {
                  do j = 1,nst ! {
                    if(s(i,j).ne.0.0d0) then
                        write(12, 2001) i,j,s(i,j)
                    endif
                  end do ! j  }
                end do ! i  }
                close(unit = 12, status = 'keep')
            endif
        endif
2001    format(2i10,1p,1d25.15)
把这段代码加到你需要输出的单元文件的fortran中就行了
譬如3d壳单元文件shell3d.f的91行。
运行后就会得到tang文件，里面包括matlab认识的稀疏矩阵形式。
然后在matlab中使用下面两句命令就能得到matlab标准的矩阵形式的刚度
load tang
FEAPStiff = full(spconvert(tang));
得到下面类似的24X24的矩阵。