-
Notifications
You must be signed in to change notification settings - Fork 11
/
sparloading.py
100 lines (86 loc) · 3.69 KB
/
sparloading.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
" spar loading "
from gpkit import Model, parse_variables, SignomialsEnabled
from gpkitmodels.GP.beam.beam import Beam
from gpkitmodels.GP.aircraft.tail.tail_boom import TailBoomState
from numpy import pi
#pylint: disable=no-member, unused-argument, exec-used, invalid-name
#pylint: disable=undefined-variable, attribute-defined-outside-init
class SparLoading(Model):
""" Spar Loading Model
Variables
---------
Nmax 5 [-] max loading
Nsafety 1.0 [-] safety load factor
kappa 0.2 [-] max tip deflection ratio
W [lbf] loading weight
N [-] loading factor
twmax 30.*pi/180 [-] max tip twist
Stip 1e-10 [N] tip loading
Mtip 1e-10 [N*m] tip moment
throot 1e-10 [-] root deflection angle
wroot 1e-10 [m] root deflection
Variables of length wing.N
--------------------------
q [N/m] distributed wing loading
S [N] shear along wing
Mbar [-] wing section root moment
M [N*m] wing section root moment
th [-] deflection angle
w [m] wing deflection
Variables of length wing.N-1
----------------------------
Mtw [N*m] local moment due to twisting
theta [-] twist deflection
EIbar [-] EIbar
Sout [N] outboard weights
LaTex Strings
-------------
Nmax N_{\\mathrm{max}}
kappa \\kappa
Mr M_r
"""
def new_qbarFun(self, c):
" define qbar model for chord loading "
barc = self.wing.planform.cbar
return [f(c) for f in self.wing.substitutions[barc]]
new_SbarFun = None
def setup(self, wing, state, sp=False):
self.wing = wing
exec parse_variables(SparLoading.__doc__)
b = self.b = self.wing.planform.b
I = self.I = self.wing.spar.I
Sy = self.Sy = self.wing.spar.Sy
cave = self.cave = self.wing.planform.cave
cbar = self.cbar = self.wing.planform.cbar
E = self.wing.spar.material.E
sigma = self.wing.spar.material.sigma
deta = self.wing.planform.deta
if sp:
print "here"
with SignomialsEnabled():
shear = S[:-1] >= (S[1:] + 0.5*deta*(b/2.)*(q[:-1] + q[1:])
- Sout)
else:
shear = S[:-1] >= S[1:] + 0.5*deta*(b/2.)*(q[:-1] + q[1:])
constraints = [
N == Nsafety*Nmax, q >= N*W/b*cbar,
shear, S[-1] >= Stip,
M[:-1] >= M[1:] + 0.5*deta*(b/2)*(S[:-1] + S[1:]), M[-1] >= Mtip,
th[0] >= throot,
th[1:] >= th[:-1] + 0.5*deta*(b/2)*(M[1:] + M[:-1])/E/I,
w[0] >= wroot, w[1:] >= w[:-1] + 0.5*deta*(b/2)*(th[1:] + th[:-1]),
sigma >= M[:-1]/Sy, w[-1]/(b/2) <= kappa,
]
self.wingSparJ = hasattr(self.wing.spar, "J")
if self.wingSparJ:
qne = self.qne = state.qne
J = self.J = self.wing.spar.J
G = self.wing.spar.shearMaterial.G
cm = self.wing.planform.CM
constraints.extend([
Mtw >= cm*cave**2*qne*deta*b/2*Nsafety,
theta[0] >= Mtw[0]/G/J[0]*deta[0]*b/2,
theta[1:] >= theta[:-1] + Mtw[1:]/G/J[1:]*deta[1:]*b/2,
twmax >= theta[-1]
])
return constraints