def box_spar():
" test wing models "
Wing.sparModel = BoxSpar
W = Wing()
W.substitutions[W.W] = 50
fs = FlightState()
perf = W.flight_model(W, fs)
loading = [W.spar.loading(W)]
loading[0].substitutions["W"] = 100
loading.append(W.spar.gustloading(W))
loading[1].substitutions["W"] = 100
from gpkit import settings
if settings["default_solver"] == "cvxopt":
for l in loading:
for v in ["\\bar{M}_{tip}", "\\bar{S}_{tip}",
"\\bar{\\delta}_{root}", "\\theta_{root}"]:
l.substitutions[v] = 1e-3
m = Model(perf.Cd, [
loading[1].v == fs.V,
loading[1].cl == perf.CL,
loading[1].Ww == W.W,
loading[1].Ww <= 0.5*fs.rho*fs.V**2*perf.CL*W.planform.S,
W, fs, perf, loading])
m.solve(verbosity=0)
def box_spar():
" test wing models "
Wing.sparModel = BoxSpar
W = Wing()
W.substitutions[W.W] = 50
W.substitutions[W.planform.tau] = 0.115
fs = FlightState()
perf = W.flight_model(W, fs)
loading = [W.spar.loading(W, fs)]
loading[0].substitutions["W"] = 100
loading.append(W.spar.gustloading(W, fs))
loading[1].substitutions["W"] = 100
from gpkit import settings
if settings["default_solver"] == "cvxopt":
for l in loading:
for v in ["Mtip", "Stip", "wroot", "throot"]:
l.substitutions[v] = 1e-2
m = Model(perf.Cd, [
loading[1].v == fs.V, loading[1].cl == perf.CL, loading[1].Ww == W.W,
loading[1].Ww <= fs.qne * perf.CL * W.planform.S, W, fs, perf, loading
])
m.solve(verbosity=0)