The purpose of this work is to use level set topology optimization to improve the design of a representative wing box structure for the NASA common research model. The objective is to minimize the total compliance of the structure under aerodynamic and body force loading, where the aerodynamic loading is coupled to the structural deformation. A taxi bump case was also considered, where only body force loads were applied. The trim condition that aerodynamic lift must balance the total weight of the aircraft is enforced by allowing the root angle of attack to change. The level set optimization method is implemented on an unstructured three-dimensional grid, so that the method can optimize a wing box with arbitrary geometry. Fast matching and upwind schemes are developed for an unstructured grid, which make the level set method robust and efficient. The adjoint method is used to obtain the coupled shape sensitivities required to perform aerostructural optimization of the wing box structure.
NomenclatureA = sensitivity factor for angle of attack a = vector of Doublet Lattice Method (DLM) box areas b = vector defining influence of wing deformed shape on lift C = compliance of the structure c p = pressure coefficient vector D = aerodynamic influence coefficient matrix E = material property tensor e = number of elements attached to a node f a = aerodynamic load vector f g = body force load vector f t = total load vector g = acceleration due to gravity h = element edge length i, j = indices K = global structural stiffness matrix K c = stiffness matrix of an element cut by the boundary K E = stiffness matrix of a finite element k = iteration number 1 Research Scholar, peter.d.dunning@nasa.gov, AIAA Member. 2 Research Aerospace Engineer, Aeroelasticity Branch, bret.k.stanford@nasa.gov, AIAA Member. 3 Senior Lecturer, Department of Mechanical Engineering, h.a.kim@bath.ac.uk, AIAA Senior Member. Downloaded by PURDUE UNIVERSITY on July 30, 2015 | http://arc.aiaa.org | 2 L = total lift force L c = lift force from built-in twist and camber Lα = lift force from unit angle of attack N = load factor n = unit normal vector p = adjoint state vector q = dynamic pressure Q = aerodynamic stiffness matrix S = force transfer matrix T = displacement transfer matrix t = fictitious time variable u = displacement field or vector V n = velocity function v = virtual displacement W b = wing box weight W c = fixed aircraft weight w = downwash dependent on deformed wing shape w c = constant downwash from built-in camber x = point in the design domain z = column vector of 1's α = angle of attack β c = volume of a cut element that lies inside the structure β E = volume of an element γ = small number Γ = structural boundary Γ D = part of boundary subject to displacement boundary conditions Γ N = part of boundary subject to aerodynamic loads Γ 0 = part of boundary free from boundary conditions and aerodynamic loads Δt = time step ε = strain tensor θ = arbitrary vector (shape derivative auxiliary variable) ρ = material density φ = implicit function Ω ...