pyNastranでop2を読む

以下、CQUAD4のFORCEを読むサンプルです。

from pyNastran.op2.op2 import OP2

op2 = OP2(debug = False)
op2.read_op2('test.op2')

vectors = getattr(op2, 'cquad4_force')

icase = 0
for subcase, vector in vectors.items():
    eids = getattr(vector, 'element', None)
    print(len(eids))
    forces = getattr(vector, 'data', None)
    print('SUBCASE', subcase)
    ielem = 0
    print('ID,FX,FY,FXY,MX,MY,MXY,QX,QY')
    for e in eids:
        print(e, end=' ')
        for val in forces[icase][ielem]:
            print(',', val, end=' ')
        print()
        ielem += 1
    icase += 1

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA