deltaC_x = np.zeros((n_c,)) deltaC_r = np.zeros((n_c,)) deltaC_r[-1] = 10 #second to last element, set to 10 deltaC = np.array(zip(deltaC_x,deltaC_r)) geom.deform(centerbody=deltaC) print "Run Time: ", time.time()-start_time start_time = time.time() geom.writeSTL('new.stl', ascii=False) print "STL Write Time: ", time.time()-start_time start_time = time.time() geom.writeFEPOINT('deform.dat') print "FEPoint Write Time: ", time.time()-start_time #geometry.writeSTL('new.stl',body.coords.cartesian,ascii=False) #centerbody.writeFEPOINT('deform.dat',body.coords.cartesian)
deltaC_c = np.array(zip(deltaC_cx,deltaC_cr)) deltaC_tx = np.zeros((n_c,)) deltaC_tr = np.zeros((n_c,)) # deltaC_tr[:-1] = 1 #second to last element, set to 10 deltaC_t = np.array(zip(deltaC_tx,deltaC_tr)) geom.deform(cowl=(deltaC_c,deltaC_t)) print "Run Time: ", time.time()-start_time start_time = time.time() geom.writeSTL('new.stl', ascii=False) print "STL Write Time: ", time.time()-start_time start_time = time.time() geom.writeFEPOINT('new.fepoint') import pylab as p profile = geom.project_profile() for point_set in profile: X = point_set[:,0] Y = point_set[:,1] p.plot(X,Y) p.show()