Esempi in Python per Geometry.deform

Esempio n. 1

File: test_body.py Progetto: JustinSGray/pyBspline

geom = Geometry()
#geom.add(body0,name="cb0")
geom.add(body,name="centerbody")

#params = geom.get_params() #params['centerbody'] = [(0,0),]

print "Bspline Compute Time: ", time.time()-start_time
start_time = time.time()

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)

Esempio n. 2

File: test_shell.py Progetto: JustinSGray/pyBspline

#params = geom.get_params() #params['centerbody'] = [(0,0),]

print "Bspline Compute Time: ", time.time()-start_time
start_time = time.time()

#calculate new P's

deltaC_x = [0,0,0]
deltaC_r = [0,0,0]
deltaC_c = np.array(zip(deltaC_x,deltaC_r))

deltaC_x = [0,0,0]
deltaC_r = [0,1,0]
deltaC_t = np.array(zip(deltaC_x,deltaC_r))

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('deform.dat')

print "FEPoint Write Time: ", time.time()-start_time

Esempio n. 3

File: test_combined.py Progetto: JustinSGray/pyBspline

shroud = Shell(outer_shroud, inner_shroud, n_C, n_C, r_ref=100, x_ref=1000)

geom = Geometry()
geom.add(centerbody, name="centerbody")
geom.add(cowl, name="cowl")
geom.add(shroud, name="shroud")

print "Geom Calc Time: ", time.time()-start_time
start_time = time.time()

deltaC_x = np.zeros((n_C,))
deltaC_x[3:] = 1
deltaC_r = np.zeros((n_C,))
deltaC_r[3] = 1 #second to last element, set to 10
deltaC = np.array(zip(deltaC_x,deltaC_r))
geom.deform(centerbody=deltaC)


deltaC_cx = np.zeros((n_C,))
deltaC_cr = np.zeros((n_C,))
deltaC_cx[2:] = 1 #second to last element, set to 10
deltaC_c = np.array(zip(deltaC_cx,deltaC_cr))
deltaC_tx = np.zeros((n_C,))
deltaC_tr = np.zeros((n_C,))
deltaC_t = np.array(zip(deltaC_tx,deltaC_tr))
geom.deform(cowl=(deltaC_c,deltaC_t))


import pylab as p
profile = geom.project_profile()
for point_set in profile: 

Esempio n. 4

File: simple_nozzle.py Progetto: JustinSGray/pyBspline

geom = Geometry()
geom.add(body,name="plug")
geom.add(shell,name="cowl")

orig_points = geom.points.copy()

print "Geometry Object Building: ", time.time()-start_time
start_time = time.time()

deltaC_x = np.zeros((n_c,))
deltaC_r = np.zeros((n_c,))
# deltaC_x[3:] = 2 #second to last element, set to 10
deltaC = np.array(zip(deltaC_x,deltaC_r))
geom.deform(plug=deltaC)


deltaC_cx = np.zeros((n_c,))
deltaC_cr = np.zeros((n_c,))
# deltaC_cr[:] = 3 #second to last element, set to 10
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()