Option to choose whether to include tip rod or not
@@ -1119,6 +1119,7 @@ class Pazy(FlyingWing):
u_inf, # flight cond
alpha,
rho=1.225,
+ tip_rod=True,
b_ref=2. * 0.55, # geometry
main_chord=0.1,
aspect_ratio=(2. * 0.55) / 0.1,
@@ -1131,6 +1132,7 @@ class Pazy(FlyingWing):
route='.',
case_name='pazy',
RollNodes=False):
+
super().__init__(M=M, N=N,
Mstar_fact=Mstar_fact,
u_inf=u_inf,
@@ -1158,6 +1160,8 @@ class Pazy(FlyingWing):
# other
self.c_ref = main_chord
+ self.tip_rod = tip_rod
+
def update_mass_stiff(self):
'''
This method can be substituted to produce different wing configs.
@@ -1237,6 +1241,7 @@ class Pazy(FlyingWing):
self.elem_stiffness = np.zeros((self.num_elem_tot,), dtype=int)
self.elem_mass = np.zeros((self.num_elem_tot,), dtype=int)
+ if self.tip_rod:
n_lumped_mass = 2
self.lumped_mass = np.zeros((n_lumped_mass))
self.lumped_mass_position = np.zeros((n_lumped_mass, 3))
@@ -1268,6 +1273,7 @@ class PazyControlSurface(Pazy):
alpha,
cs_deflection=[0],
rho=1.225,
+ tip_rod=True,
b_ref= 2. * 0.55, # geometry
main_chord= 0.1,
pct_flap= 0.2,
@@ -1287,6 +1293,7 @@ class PazyControlSurface(Pazy):
u_inf=u_inf,
alpha=alpha,
rho=rho,
+ tip_rod=tip_rod,
b_ref=b_ref,
main_chord=main_chord,
aspect_ratio=aspect_ratio,
@@ -1343,7 +1350,6 @@ class PazyControlSurface(Pazy):
# if inode >= num_node_surf // 2:
ws_elem = 0
for i_surf in range(2):
- print('Surface' + str(i_surf))
for i_elem in range(num_elem_surf):
for i_local_node in range(self.num_node_elem):
if i_elem >= int(num_elem_surf *(1- pct_flap)):
