150 B. Bahr et al. dofi = float(dof[i]) nodeNum = int(dofi) dirval = int(round((dofi-nodeNum)*10)) force = fshape[i] if (nodeNum == 0 or abs(force) < 1.0e-15): continu e elif (dirval == 1): f.write(str(nodeNum)+’, 1, ’) f.write(str(force)+’\n’) elif (dirval == 2): f.write(str(nodeNum)+’, 2, ’) f.write(str(force)+’\n’) elif (dirval == 3): f.write(str(nodeNum)+’, 3, ’) f.write(str(force)+’\n’) elif (dirval == 4): f.write(str(nodeNum)+’, 4, ’) f.write(str(force)+’\n’) elif (dirval == 5): f.write(str(nodeNum)+’, 5, ’) f.write(str(force)+’\n’) elif (dirval == 6): f.write(str(nodeNum)+’, 6, ’) f.write(str(force)+’\n’) else : continu e f.write(’*END STEP\n’) with open(modelname + ’_staticforce.inp’,’w’) as f: f.write(’** Static Force\n’) f.write(’**\n**\n’) f.write(’*Restart, read, step=1\n’) f.write(’** ----------------------------------------------------------- -----\n’) f.write(’**\n’) f.write(’** ---------------------------------------------------\n’) f.write(’** ------ STEP data for SUBCASE 1\n’) f.write(’** ---------------------------------------------------\n’) f.write(’*INCLUDE, INPUT=’ + modelname + ’_force.inp’) mdic = {’alpha’:alpha, ’Mphi’: mPhi, ’fn’:freq, ’phi_max_ind’:phi_max_ind, ’phi_max’:phi_max} savemat("VariablesFromModal", mdic, True) numpy.save("Mphi",mPhi) toc = timeit.default_timer() print(toc - tic)
RkJQdWJsaXNoZXIy MTMzNzEzMQ==