📄 138-143.html
字号:
<HTML>
<HEAD>
<META name=vsisbn content="0849398010">
<META name=vstitle content="Industrial Applications of Genetic Algorithms">
<META name=vsauthor content="Charles Karr; L. Michael Freeman">
<META name=vsimprint content="CRC Press">
<META name=vspublisher content="CRC Press LLC">
<META name=vspubdate content="12/01/98">
<META name=vscategory content="Web and Software Development: Artificial Intelligence: Other">
<TITLE>Industrial Applications of Genetic Algorithms:Optimized Non-Coplanar Orbital Transfers Using Genetic Algorithms</TITLE>
<!-- HEADER -->
<STYLE type="text/css">
<!--
A:hover {
color : Red;
}
-->
</STYLE>
<META NAME="ROBOTS" CONTENT="NOINDEX, NOFOLLOW">
<!--ISBN=0849398010//-->
<!--TITLE=Industrial Applications of Genetic Algorithms//-->
<!--AUTHOR=Charles Karr//-->
<!--AUTHOR=L. Michael Freeman//-->
<!--PUBLISHER=CRC Press LLC//-->
<!--IMPRINT=CRC Press//-->
<!--CHAPTER=8//-->
<!--PAGES=138-143//-->
<!--UNASSIGNED1//-->
<!--UNASSIGNED2//-->
<CENTER>
<TABLE BORDER>
<TR>
<TD><A HREF="135-138.html">Previous</A></TD>
<TD><A HREF="../ewtoc.html">Table of Contents</A></TD>
<TD><A HREF="143-148.html">Next</A></TD>
</TR>
</TABLE>
</CENTER>
<P><BR></P>
<P>Reichert [5] used a GA to minimize the total ΔV in three coplanar orbit transfer cases. The first was the general transfer between circular orbits studied by Hohmann in 1925, which can typically be solved by hand. The second and third were the two cases studied by Lawden. The classical solution to these cases was used to verify the accuracy of the genetic algorithm.
</P>
<P>In 1995, Pinon and Fowler [6] used a GA to generate two-dimensional lunar trajectories. The goal was to analyze orbital transfers from the earth to the moon. Their results were comparable to the trajectories used by the Apollo missions. They concluded that the two-dimensional trajectories produced by the GA could be used as a starting point for three-dimensional transfers.</P>
<P>The above-mentioned studies were simplified to have either coplanar or circular initial and final orbits or both. The purpose of the current effort is to examine the use of a GA in non-coplanar orbit transfers. A GA is effective in finding near optimum solutions, but does not guarantee to find the optimum solution. Because of the GA’s capability, test cases will include both circular and elliptical problems. This situation is encountered more often in orbital maneuvers. However, this problem, like the coplanar transfer, is not easily solved. In this case, not only must an optimum transfer orbit be found, but also the amount of plane change at each ΔV.</P>
<P><FONT SIZE="+1"><B>PROBLEM STATEMENT</B></FONT></P>
<P>In this chapter, orbital perturbations will be ignored, impulsive maneuvers will be assumed, and test cases will be limited to bodies orbiting the earth. In a real world analysis, these issues would need to be addressed. However, neglecting them will be sufficient for this study. This investigation will limit the number of velocity changes to two, the first when leaving the initial orbit and the second at the intersection of the transfer orbit and the final orbit. It will also be assumed that the semi-major axes of the initial, final and transfer orbits are aligned with one another.
</P>
<P>The analysis will begin by user specification of the initial and final orbits, which will be limited to elliptical or circular. Parameters used to define an orbit will be the eccentricity (e), the argument of perigee (ω), the semi-major axis (a), the inclination with respect to the equatorial plane (i), and the longitude of the ascending node (Ω). The GA will be used to search for defining parameters for the transfer orbit and an initial plane change (α<SUB><SMALL>1</SMALL></SUB>) that minimizes the velocity change needed.</P>
<P>This study assumes the first velocity change will occur at the perigee point of the initial orbit and of the transfer orbit. Therefore, the perigee radius of the initial and transfer orbits will correspond. Since the perigee radius of the transfer orbit is defined by the initial orbit, the semi-major axis can be calculated (Equation 8.1) and need not be one of the search parameters for the GA.</P>
<P ALIGN="CENTER"><IMG SRC="images/08-01d.jpg"></P>
<P>where
</P>
<P ALIGN="CENTER"><IMG SRC="images/08-02d.jpg"></P>
<P>The perigee velocity on both the initial and transfer orbits can then be calculated using Equation 8.3.
</P>
<P ALIGN="CENTER"><IMG SRC="images/08-03d.jpg"></P>
<P>where
</P>
<P ALIGN="CENTER"><IMG SRC="images/08-04d.jpg"></P>
<P>The velocity change required to go from the initial orbit to the transfer orbit (ΔV<SUB><SMALL>1</SMALL></SUB>) can be found by applying the law of cosines (Figure 8.5, Equation 8.5).</P>
<P><A NAME="Fig5"></A><A HREF="javascript:displayWindow('images/08-05.jpg',350,176)"><IMG SRC="images/08-05t.jpg"></A>
<BR><A HREF="javascript:displayWindow('images/08-05.jpg',350,176)"><FONT COLOR="#000077"><B>Figure 8.5</B></FONT></A> Physical representation of the law of cosines at the first velocity change.</P>
<P ALIGN="CENTER"><IMG SRC="images/08-05d.jpg"></P>
<P>where
</P>
<DL>
<DD><IMG SRC="images/08-13i.jpg"> = perigee velocity on the transfer orbit<DD><IMG SRC="images/08-14i.jpg"> = perigee velocity on the initial orbit</DL>
<P>In order for the transfer orbit to be sufficient, it must intersect the final orbit (Figure 8.6). The next step will be to determine the intersection location. The initial, random population will undoubtedly contain transfer orbits that do not intersect the final orbit and are not possible solutions. In these situations, the objective function will be penalized according to how badly the transfer orbit misses the final orbit (details will be discussed later).
</P>
<P><A NAME="Fig6"></A><A HREF="javascript:displayWindow('images/08-06.jpg',350,459)"><IMG SRC="images/08-06t.jpg"></A>
<BR><A HREF="javascript:displayWindow('images/08-06.jpg',350,459)"><FONT COLOR="#000077"><B>Figure 8.6</B></FONT></A> (a) transfer orbit does not intersect final orbit, (b) transfer orbit intersects final orbit at two locations, (c) transfer orbit intersects final orbit at one location.</P>
<P>The process of finding the intersection point or points of two non-coplanar elliptical orbits is not a trivial task. The point will be defined by the following two values: a true anomaly angle on the transfer orbit, v<SUB><SMALL>t</SMALL></SUB>, and a true anomaly angle on the final orbit, v<SUB><SMALL>f</SMALL></SUB> (Figure 8.7).</P>
<P><A NAME="Fig7"></A><A HREF="javascript:displayWindow('images/08-07.jpg',150,209)"><IMG SRC="images/08-07t.jpg"></A>
<BR><A HREF="javascript:displayWindow('images/08-07.jpg',150,209)"><FONT COLOR="#000077"><B>Figure 8.7</B></FONT></A> Representation of intersection points.</P>
<P>The first step is to represent the radius vector to the intersection point on both the transfer orbit and the final orbit in a “perifocal coordinate system” (Figure 8.8). The perifocal coordinate system is defined as having a fundamental plane in the plane of the orbit. Positive <I>x</I> points toward perigee of the orbit and positive <I>y</I> is rotated 90° in the direction of the orbital motion. This is one of the most convenient coordinate systems for describing an orbit because the <I>z</I> component is always zero (Equation 8.6).</P>
<P><A NAME="Fig8"></A><A HREF="javascript:displayWindow('images/08-08.jpg',350,480)"><IMG SRC="images/08-08t.jpg"></A>
<BR><A HREF="javascript:displayWindow('images/08-08.jpg',350,480)"><FONT COLOR="#000077"><B>Figure 8.8</B></FONT></A> Perifocal coordinate system.</P>
<P ALIGN="CENTER"><IMG SRC="images/08-06d.jpg"></P>
<P>where
</P>
<P ALIGN="CENTER"><IMG SRC="images/08-07d.jpg"></P>
<P><BR></P>
<CENTER>
<TABLE BORDER>
<TR>
<TD><A HREF="135-138.html">Previous</A></TD>
<TD><A HREF="../ewtoc.html">Table of Contents</A></TD>
<TD><A HREF="143-148.html">Next</A></TD>
</TR>
</TABLE>
</CENTER>
<hr width="90%" size="1" noshade>
<div align="center">
<font face="Verdana,sans-serif" size="1">Copyright © <a href="/reference/crc00001.html">CRC Press LLC</a></font>
</div>
<!-- all of the reference materials (books) have the footer and subfoot reveresed -->
<!-- reference_subfoot = footer -->
<!-- reference_footer = subfoot -->
</BODY>
</HTML>
<!-- END FOOTER -->
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -