Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimal start”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Microwave multispectral investigations of snow

A long term observational program on microwave emission and scatter behavior under controlled conditions, was started at a high altitude alpine test site. All stages of development of snow cover during the whole season are under investigation. The study was done to achieve required knowledge on microwave radiative properties of snow for optimization of the microwave payloads of air and spaceborne snow sensors and for the interpretation of large scale snow maps obtained by these sensors. Preliminary results of the first month obtained with radiometers at 4.9, 10.5, 21 and 36 GHz are presented.

Schanda, E.↗
NASA NTRS1989

Fuel-optimal trajectories of aeroassisted orbital transfer with plane change

The problem of minimization of fuel consumption during the atmospheric portion of an aeroassisted, orbital transfer with plane change is addressed. The complete mission has required three characteristic velocities, a deorbit impulse at high earth orbit (HEO), a boost impulse at the atmospheric exit, and a reorbit impulse at low earth orbit (LEO). A performance index has been formulated as the sum of these three impulses. Application of optimal control principles has led to a nonlinear, two-point, boundary value problem which was solved by using a multiple shooting algorithm. The strategy for the atmospheric portion of the minimum-fuel transfer is to start initially with the maximum positive lift in order to recover from the downward plunge, and then to fly with a gradually decreasing lift such that the vehicle skips out of the atmosphere with a flight path angle near zero degrees.

Naidu, Desineni Subbaramaiah↗
NASA NTRS1989

Fuel-optimal trajectories of aeroassisted orbital transfer with plane change

The problem of minimization of fuel consumption during the atmospheric portion of an aeroassisted, orbital transfer with plane change is addressed. The complete mission has required three characteristic velocities, a deorbit impulse at high earth orbit (HEO), a boost impulse at the atmospheric exit, and a reorbit impulse at low earth orbit (LEO). A performance index has been formulated as the sum of these three impulses. Application of optimal control principles has led to a nonlinear, two-point, boundary value problem which was solved by using a multiple shooting algorithm. The strategy for the atmospheric portion of the minimum-fuel transfer is to start initially with the maximum positive lift in order to recover from the downward plunge, and then to fly with a gradually decreasing lift such that the vehicle skips out of the atmosphere with a flight path angle near zero degrees.

Naidu, D. S.↗
NASA NTRS1999

Investigation of the Influence of Microgravity on Transport Mechanism in a Virtual Spaceflight Chamber: A Flight Definition Program

A need exists for understanding precisely how particles move and interact in a fluid in the absence of gravity. Such understanding is required, for example, for modeling and predicting crystal growth in space where crystals grow from solution around nucleation sites as well as for any study of particles or bubbles in liquids or in experiments where particles are used as tracers for mapping microconvection. We have produced an exact solution to the general equation of motion of particles at extremely low Reynolds number in microgravity that covers a wide range of interesting conditions. We have also developed diagnostic tools and experimental techniques to test the validity of the general equation . This program, which started in May, 1998, will produce the flight definition for an experiment in a microgravity environment of space to validate the theoretical model. We will design an experiment with the help of the theoretical model that is optimized for testing the model, measuring g, g-jitter, and other microgravity phenomena. This paper describes the goals, rational, and approach for the flight definition program. The first objective of this research is to understand the physics of particle interactions with fluids and other particles in low Reynolds number flows in microgravity. Secondary objectives are to (1) observe and quantify g-jitter effects and microconvection on particles in fluids, (2) validate an exact solution to the general equation of motion of a particle in a fluid, and (3) to characterize the ability of isolation tables to isolate experiments containing particle in liquids. The objectives will be achieved by recording a large number of holograms of particle fields in microgravity under controlled conditions, extracting the precise three-dimensional position of all of the particles as a function of time and examining the effects of all parameters on the motion of the particles. The feasibility for achieving these results has already been established in the ongoing ground-based NRA, which led to the "virtual spaceflight chamber" concept.

Trolinger, James D.↗
NASA NTRS2018

Information Theory Applied to Decision Making Structures

Decision making structures, such as control boards, process information on many topics as they select options for the system and project. The decisions are based on the information known to the board members and presenters (subject matter experts) and shared at the board during discussion. This information flow through this process can be modelled using information theory. Information theory provides a mathematical basis to understand the flow of information through the decision making process and the information needed for a particular decision. Information theory also provides the mathematical relationships on which to base the optimal decision making structure for a specific system development and organizational structure. Since decision making bodies provide control for the system or project, control theory can be used to construct a decision making model. This provides a starting point for adding cognitive science models. Information processing by each individual board participant can be represented through cognitive processes which are integrated across the board participants through information theory relationship. The set theory view of information theory provides a structure in which to look at the relationships between the participants in a decision making structure.

Watson, Michael D.↗
NASA NTRS2018

Information Theory Applied to Decision Making Structures

Decision making structures, such as control boards, process information on many topics as they select options for the system and project. The decisions are based on the information known to the board members and presenters (subject matter experts) and shared at the board during discussion. This information flow through this process can be modelled using information theory. Information theory provides a mathematical basis to understand the flow of information through the decision making process and the information needed for a particular decision. Information theory also provides the mathematical relationships on which to base the optimal decision making structure for a specific system development and organizational structure. Since decision making bodies provide control for the system or project, control theory can be used to construct a decision making model. This provides a starting point for adding cognitive science models. Information processing by each individual board participant can be represented through cognitive processes which are integrated across the board participants through information theory relationship. The set theory view of information theory provides a structure in which to look at the relationships between the participants in a decision making structure.

Watson, Michael D.↗
NASA NTRS1991

Flow analysis and design optimization methods for nozzle afterbody of a hypersonic vehicle

This report summarizes the methods developed for the aerodynamic analysis and the shape optimization of the nozzle-afterbody section of a hypersonic vehicle. Initially, exhaust gases were assumed to be air. Internal-external flows around a single scramjet module were analyzed by solving the three dimensional Navier-Stokes equations. Then, exhaust gases were simulated by a cold mixture of Freon and Argon. Two different models were used to compute these multispecies flows as they mixed with the hypersonic airflow. Surface and off-surface properties were successfully compared with the experimental data. In the second phase of this project, the Aerodynamic Design Optimization with Sensitivity analysis (ADOS) was developed. Pre and post optimization sensitivity coefficients were derived and used in this quasi-analytical method. These coefficients were also used to predict inexpensively the flow field around a changed shape when the flow field of an unchanged shape was given. Starting with totally arbitrary initial afterbody shapes, independent computations were converged to the same optimum shape, which rendered the maximum axial thrust.

Baysal, Oktay↗
NASA NTRSDate not supplied

Baseline Medical System Translation for the Impact Medical Database

NASA has developed a new evidence-based data-driven probabilistic risk assessment and tradespace analysis tool as a successor to the Integrated Medical Model (IMM). This updated decision support tool is known as IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces). Whereas IMM focuses on the resources and risks associated with International Space Station (ISS) and Low Earth Orbit (LEO) missions, IMPACT estimates the frequency and consequences of medical conditions that might arise during exploration missions. One of the services offered by IMM is a series of generic or baseline medical systems associated with typical mission types or DRMs (Design Reference Missions), such that requestors may prioritize questions pertaining to the mission itself over the medical supplies indicated by the model outputs for that DRM. In a mission focused request, the appropriate baseline medical system is used in place of a prepared or optimized medical kit. In preparation, an effort was undertaken to create baseline systems of medical resources within IMPACT suitable for typical DRMs. Using an existing baseline medical system within IMM’s Integrated Medical Evidence Database (iMED) as a starting point, the resources found within the medical kit were compared to and substituted for equivalent resources available within the IMPACT MD (Medical Database). In consultation with clinicians with knowledge of IMPACT’s MD, each resource was matched as closely as possible to a similarly purposed resource in IMPACT, seeking to preserve the treatment capabilities and procedures offered by the IMM medical system while reconciling the differences in modeled resources and medical conditions between the models. To demonstrate the efficacy of this work, a prototype ISS medical system in IMPACT was translated from the baseline used in the IMM. The appropriate medical system was then run through its associated medical model to compare and validate the resultant risks and risk mitigation provided by each baseline ISS medical kit.

S Schwartz↗
NASA NTRS1989

An algorithm for the solution of dynamic linear programs

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation scheme.

Psiaki, Mark L.↗
NASA NTRS2011

Post-Fragmentation Whole Genome Amplification-Based Method

This innovation is derived from a proprietary amplification scheme that is based upon random fragmentation of the genome into a series of short, overlapping templates. The resulting shorter DNA strands (<400 bp) constitute a library of DNA fragments with defined 3 and 5 termini. Specific primers to these termini are then used to isothermally amplify this library into potentially unlimited quantities that can be used immediately for multiple downstream applications including gel eletrophoresis, quantitative polymerase chain reaction (QPCR), comparative genomic hybridization microarray, SNP analysis, and sequencing. The standard reaction can be performed with minimal hands-on time, and can produce amplified DNA in as little as three hours. Post-fragmentation whole genome amplification-based technology provides a robust and accurate method of amplifying femtogram levels of starting material into microgram yields with no detectable allele bias. The amplified DNA also facilitates the preservation of samples (spacecraft samples) by amplifying scarce amounts of template DNA into microgram concentrations in just a few hours. Based on further optimization of this technology, this could be a feasible technology to use in sample preservation for potential future sample return missions. The research and technology development described here can be pivotal in dealing with backward/forward biological contamination from planetary missions. Such efforts rely heavily on an increasing understanding of the burden and diversity of microorganisms present on spacecraft surfaces throughout assembly and testing. The development and implementation of these technologies could significantly improve the comprehensiveness and resolving power of spacecraft-associated microbial population censuses, and are important to the continued evolution and advancement of planetary protection capabilities. Current molecular procedures for assaying spacecraft-associated microbial burden and diversity have inherent sample loss issues at practically every step, particularly nucleic acid extraction. In engineering a molecular means of amplifying nucleic acids directly from single cells in their native state within the sample matrix, this innovation has circumvented entirely the need for DNA extraction regimes in the sample processing scheme.

Benardini, James↗
NASA NTRS2014

Quantification of In-flight Physical Changes: Anthropometry and Neutral Body Posture

Currently, NASA does not have sufficient in-flight anthropometric data gathered to assess the impact of physical body shape and size changes on suit sizing. For developing future planetary and reduced gravity suits, NASA needs to quantify the impacts of microgravity on anthropometry, body posture, and neutral body postures (NBP) to ensure optimal crew performance, fit, and comfort. To obtain these impacts, anthropometric data, circumference, length, height, breadth, and depth for body segments (i.e. chest, waist, bicep, thigh, calf) from astronauts for pre, in-, and postflight conditions needs to be collected. Once this data has been collected, a comparison between pre, in-, and postflight anthropometric values will be analyzed, yielding microgravity factors. The NBP will be used to determined body posture (joint angle) changes between subjects throughout the duration of a mission. Data collection, starting with Increments 37/38, is still in progress with the completion of 3 out of 12 subjects. NASA suit engineers and NASA's Extravehicular Activity (EVA) Project Office have identified that suit fit in microgravity could become an issue. It has been noted that crewmembers often need to adjust their suit sizing once they are in orbit. This adjustment could be due to microgravity effects on anthropometry and postural changes, and is necessary to ensure optimal crew performance, fit, and comfort in space. To date, the only data collected to determine the effects of microgravity on physical human changes have been during Skylab, STS-57, and a recent HRP study on seated height changes due to spinal elongation (Spinal Elongation, Master Task List [MTL] #221). The Skylab and the STS-57 studies found that there is a distinct neutral body posture (NBP) based on photographs. The still photographs showed that there is a distinguishable posture with the arms raised and the shoulder abducted; and, in addition, the knees were flexed with noticeable hip flexion and the foot plantar flexed [1,2]. This is the one standard set of body joint angles for a NBP in microgravity. A recent simulated microgravity NBP study [3] has shown an individual variability and inconsistencies in defining NBP. This variation may be influenced by spinal growth, the type of suit fit, and other potential anthropometry factors such as spinal curvature, age, and gender. The variation aspect of this essential data is required for all kinds of space device designs (e.g. suits, habitat, mobility aids, etc.). The method proposed considers the dynamic nature of body movement and will use a measurement technique to continually monitor posture and develop a probability likelihood of the natural posture and how the NBP postures are affected by anthropometry. Additionally, Skylab studies found that crewmembers experienced a stature growth of up to 3%. The data included 3 crewmembers that showed that there is a bi-phasic stature growth once the crew enters into weightlessness. However, the Spinal Elongation study identified that the crewmembers could experience about a 6% growth in seated height and a 3% stature growth, when exposed to microgravity. The results prove that not all anthropometric measurements have the same microgravity percent growth factor. For EVA and suit engineers to properly update the sizing protocol for microgravity, they need additional anthropometric data from space missions. Hence, this study is aimed to gather additional in-flight anthropometric measurements, such as length, depth, breadth, and circumference, to determine the changes to body shape and size due to microgravity effects. It is anticipated that by recording the potential changes to body shape and size, a better suit sizing protocol will be developed for ISS and other space missions. In essence, this study will help NASA quantify the impacts of microgravity on anthropometry to ensure optimal crew performance, fit, and comfort. This study will use simplistic data collection techniques, 3D laser scanning, digital still, and video data, and perform photogrammetric analyses to determine the changes that occur to the body shape, size, and NBP when exposed to a microgravity environment.

Young, K. S.↗
NASA NTRS2005

A Parallel Particle Swarm Optimization Algorithm Accelerated by Asynchronous Evaluations

A parallel Particle Swarm Optimization (PSO) algorithm is presented. Particle swarm optimization is a fairly recent addition to the family of non-gradient based, probabilistic search algorithms that is based on a simplified social model and is closely tied to swarming theory. Although PSO algorithms present several attractive properties to the designer, they are plagued by high computational cost as measured by elapsed time. One approach to reduce the elapsed time is to make use of coarse-grained parallelization to evaluate the design points. Previous parallel PSO algorithms were mostly implemented in a synchronous manner, where all design points within a design iteration are evaluated before the next iteration is started. This approach leads to poor parallel speedup in cases where a heterogeneous parallel environment is used and/or where the analysis time depends on the design point being analyzed. This paper introduces an asynchronous parallel PSO algorithm that greatly improves the parallel e ciency. The asynchronous algorithm is benchmarked on a cluster assembled of Apple Macintosh G5 desktop computers, using the multi-disciplinary optimization of a typical transport aircraft wing as an example.

Venter, Gerhard↗
NASA NTRS1978

Stirling engine design manual

This manual is intended to serve both as an introduction to Stirling engine analysis methods and as a key to the open literature on Stirling engines. Over 800 references are listed and these are cross referenced by date of publication, author and subject. Engine analysis is treated starting from elementary principles and working through cycles analysis. Analysis methodologies are classified as first, second or third order depending upon degree of complexity and probable application; first order for preliminary engine studies, second order for performance prediction and engine optimization, and third order for detailed hardware evaluation and engine research. A few comparisons between theory and experiment are made. A second order design procedure is documented step by step with calculation sheets and a worked out example to follow. Current high power engines are briefly described and a directory of companies and individuals who are active in Stirling engine development is included. Much remains to be done. Some of the more complicated and potentially very useful design procedures are now only referred to. Future support will enable a more thorough job of comparing all available design procedures against experimental data which should soon be available.

Martini, W. R.↗
NASA NTRS1997

Thermal History and Mantle Dynamics of Venus

One objective of this research proposal is to develop a 3-D thermal history model for Venus. The basis of our study is a finite-element computer model to simulate thermal convection of fluids with highly temperature- and pressure-dependent viscosities in a three-dimensional spherical shell. A three-dimensional model for thermal history studies is necessary for the following reasons. To study planetary thermal evolution, one needs to consider global heat budgets of a planet throughout its evolution history. Hence, three-dimensional models are necessary. This is in contrasts to studies of some local phenomena or local structures where models of lower dimensions may be sufficient. There are different approaches to treat three-dimensional thermal convection problems. Each approach has its own advantages and disadvantages. Therefore, the choice of the various approaches is subjective and dependent on the problem addressed. In our case, we are interested in the effects of viscosities that are highly temperature dependent and that their magnitudes within the computing domain can vary over many orders of magnitude. In order to resolve the rapid change of viscosities, small grid spacings are often necessary. To optimize the amount of computing, variable grids become desirable. Thus, the finite-element numerical approach is chosen for its ability to place grid elements of different sizes over the complete computational domain. For this research proposal, we did not start from scratch and develop the finite element codes from the beginning. Instead, we adopted a finite-element model developed by Baumgardner, a collaborator of this research proposal, for three-dimensional thermal convection with constant viscosity. Over the duration supported by this research proposal, a significant amount of advancements have been accomplished.

Hsui, Albert T.↗
NASA NTRS2018

Static Controls Performance Tool for Lunar Landers

This document presents a static analysis tool used to evaluate the controllability of Lunar landers. This was created as part of the NASA Lunar Cargo Transportation and Landing by Soft Touchdown (Lunar CATALYST) program. This tool is capable of accepting typical design information such as location and direction of thrusters, maximum thruster forces, gravity vectors, and center of mass locations. The tool evaluates how far the center of gravity can move from its starting position while still maintaining control. This type of analysis is intended to support results produced by time domain simulations. The code created for this project was implemented in Python, and it was designed to be integrated into systems level optimization tools to yield first-cut results on optimal thruster placement.

Aretskin-Hariton, Eliot D.↗
NASA NTRSDate not supplied

OSIRIS-REx Earth Return & Entry: Targeting Strategy and Maneuver Performance

The OSIRIS-REx sample return capsule landed safely in the Utah desert on September 24, 2023. Starting with the departure from asteroid Bennu in May 2021 through the release of the sample capsule and its arrival at the edge of Earth’s atmosphere hours later, this phase of the mission was planned and analyzed thoroughly to optimize the probability of a successful sample return in a large number of nominal and off-nominal scenarios, all while emphasizing spacecraft safety. This paper discusses the spacecraft and sample capsule Earth return and atmospheric entry trajectory design, maneuver analysis, and flight performance, which ultimately helped lead to a successful sample return and initiation of the OSIRIS-APEX extended mission.

Daniel R. Wibben↗
NASA NTRS2004

Starting Vortex Identified as Key to Unsteady Ejector Performance

Unsteady ejectors are currently under investigation for use in some pulse-detonation-engine-based propulsion systems. Experimental measurements made in the past, and recently at the NASA Glenn Research Center, have demonstrated that thrust augmentation can be enhanced considerably when the driver is unsteady. In ejector systems, thrust augmentation is defined as = T(sup Total)/T(sup j), where T(sup Total) is the total thrust of the combined ejector and driving jet and T(sup j) is the thrust due to the driving jet alone. There are three images in this figure, one for each of the named thrust sources. The images are color contours of measured instantaneous vorticity. Each image is an ensemble average of at least 150 phase-locked measurements. The flow is from right to left, and the shape and location of each driver is shown on the far right of each image. The emitted vortex is a clearly defined "doughnut" of highly vortical (spinning) flow. In these planar images, the vortex appears as two distorted circles, one above, and one below the axis of symmetry. Because they are spinning in the opposite direction, the two circles have vorticity of opposite sign and thus are different colors. There is also a rectangle shown in each image. Its width represents the ejector diameter that was found experimentally to yield the highest thrust augmentation. It is apparent that the optimal ejector diameter is that which just "captures" the vortex: that is, the diameter bounding the outermost edge of the vortex structure. The exact mechanism behind the enhanced performance is unclear; however, it is believed to be related to the powerful vortex emitted with each pulse of the unsteady driver. As such, particle imaging velocimetry (PIV) measurements were obtained for three unsteady drivers: a pulsejet, a resonance tube, and a speaker-driven jet. All the drivers were tested with ejectors, and all exhibited performance enhancement over similarly sized steady drivers. The characteristic starting vortices of each driver are shown in these images. The images are color contours of measured instantaneous vorticity. Each image is an ensemble average of at least 150 phase-locked measurements. The flow is from right to left. The shape and location of each driver is shown on the far right of each image. The rectangle shown in each image represents the ejector diameter that was found experimentally to yield the highest thrust augmentation. It is apparent that the optimal ejector diameter is that which just "captures" the vortex: that is, the diameter bounding the outermost edge of the vortex structure. Although not shown, it was observed that the emitted vortex spread as it traveled downstream. The spreading rate for the pulsejet is shown as the dashed lines in the top image. A tapered ejector was fabricated that matched this shape. When tested, the ejector demonstrated superior performance to all those previously tested at Glenn (which were essentially of straight, cylindrical form), achieving a remarkable thrust augmentation of 2. The measured thrust augmentation is shown as a function of ejector length. Also shown are the thrust augmentation values achieved with the straight, cylindrical ejectors of varying diameters. Here, thrust augmentation is plotted as a function of ejector length for several families of ejector diameters. It can be seen that large thrust augmentation values are indeed obtained and that they are sensitive to both ejector length and diameter, particularly the latter. Five curves are shown. Four correspond to straight ejector diameters of 2.2, 3.0, 4.0, and 6.0 in. The fifth curve corresponds to the tapered ejector contoured to bound the emitted vortex. For each curve, there are several data points corresponding to different lengths. The largest value of thrust augmentation is 2.0 for the tapered ejector and 1.81 for the straight ejectors. Regardless of their diameters, all the ejectors trend toward peak performance at a particular leng. That the cross-sectional dimensions of optimal ejectors scaled precisely with the vortex dimensions on three separate pulsed thrust sources demonstrates that the action of the vortex is responsible for the enhanced ejector performance. The result also suggests that, in the absence of a complete understanding of the entrainment and augmentation mechanisms, methods of characterizing starting vortices may be useful for correlating and predicting unsteady ejector performance.

Paxson, Daniel E.↗
NASA NTRS1990

Mathematical theory of a relaxed design problem in structural optimization

Various attempts have been made to construct a rigorous mathematical theory of optimization for size, shape, and topology (i.e. layout) of an elastic structure. If these are represented by a finite number of parametric functions, as Armand described, it is possible to construct an existence theory of the optimum design using compactness argument in a finite dimensional design space or a closed admissible set of a finite dimensional design space. However, if the admissible design set is a subset of non-reflexive Banach space such as L(sup infinity)(Omega), construction of the existence theory of the optimum design becomes suddenly difficult and requires to extend (i.e. generalize) the design problem to much more wider class of design that is compatible to mechanics of structures in the sense of variational principle. Starting from the study by Cheng and Olhoff, Lurie, Cherkaev, and Fedorov introduced a new concept of convergence of design variables in a generalized sense and construct the 'G-Closure' theory of an extended (relaxed) optimum design problem. A similar attempt, but independent in large extent, can also be found in Kohn and Strang in which the shape and topology optimization problem is relaxed to allow to use of perforated composites rather than restricting it to usual solid structures. An identical idea is also stated in Murat and Tartar using the notion of the homogenization theory. That is, introducing possibility of micro-scale perforation together with the theory of homogenization, the optimum design problem is relaxed to construct its mathematical theory. It is also noted that this type of relaxed design problem is perfectly matched to the variational principle in structural mechanics.

Kikuchi, Noboru↗