Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Upload

Welcome to Scribd!

  • 86 views

    ASEN3200 LabO3 2005

    Uploaded by

    ChristianRomiSihombing
    The document describes an assignment to numerically integrate the equations of motion for a satellite orbit perturbed by Earth's oblateness (J2 term). Students are asked to: 1) Write code to integrate the equations of motion in MATLAB and plot the resulting orbit elements, energy, and angular momentum deviations over three revolutions. 2) Generate the same perturbed orbit in STK and plot the relative inclination-circular (RIC) differences between the perturbed and two-body orbits. 3) Discuss topics like verifying the differential equations, analyzing the orbital element variations, and comparing the MATLAB and STK results.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Download as PDF, TXT or read online from Scribd

    ASEN3200 LabO3 2005

    Uploaded by

    ChristianRomiSihombing
    86 views4 pages
    The document describes an assignment to numerically integrate the equations of motion for a satellite orbit perturbed by Earth's oblateness (J2 term). Students are asked to: 1) Write code to integrate the equations of motion in MATLAB and plot the resulting orbit elements, energy, and angular momentum deviations over three revolutions. 2) Generate the same perturbed orbit in STK and plot the relative inclination-circular (RIC) differences between the perturbed and two-body orbits. 3) Discuss topics like verifying the differential equations, analyzing the orbital element variations, and comparing the MATLAB and STK results.

    Original Title

    ASEN3200_LabO3_2005

    Copyright

    © Attribution Non-Commercial (BY-NC)

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document describes an assignment to numerically integrate the equations of motion for a satellite orbit perturbed by Earth's oblateness (J2 term). Students are asked to: 1) Write code to integrate the equations of motion in MATLAB and plot the resulting orbit elements, energy, and angular momentum deviations over three revolutions. 2) Generate the same perturbed orbit in STK and plot the relative inclination-circular (RIC) differences between the perturbed and two-body orbits. 3) Discuss topics like verifying the differential equations, analyzing the orbital element variations, and comparing the MATLAB and STK results.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    86 views4 pages

    ASEN3200 LabO3 2005

    Uploaded by

    ChristianRomiSihombing
    The document describes an assignment to numerically integrate the equations of motion for a satellite orbit perturbed by Earth's oblateness (J2 term). Students are asked to: 1) Write code to integrate the equations of motion in MATLAB and plot the resulting orbit elements, energy, and angular momentum deviations over three revolutions. 2) Generate the same perturbed orbit in STK and plot the relative inclination-circular (RIC) differences between the perturbed and two-body orbits. 3) Discuss topics like verifying the differential equations, analyzing the orbital element variations, and comparing the MATLAB and STK results.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 4

    ASEN 3200 Orbit Mechanics / Attitude Dynamics & Control Lab O-3 Assigned 2/18/05, Due 3/1/05

    The potential function and equations of motion for a satellite perturbed by Earth oblateness are given below. U X X U Y = [ U ]ECI = Y Z ECI U Z
    2 J2 R Z2 U = 1 + 1 3 2 , E = r 2 r r ECI

    v2 U 2

    = U = U I + U J + U K = XI + YJ + ZK r X Y Z where
    2 2 = X 1 3 J 2 R 5 Z 1 X 2 r3 2 r r Y Y= X X 2 2 = Z 1 3 J 2 R 5 Z 3 Z 2 3 r 2 r r J 2 = 0.00108263,

    = 398600.4418 km
    R = 6378.137km

    s2

    J2, , and R are parameters for the JGM-2 gravity field. E is the total energy and will be constant for a potential function dependent only on zonal harmonics. The initial conditions for an orbit are given as follows: P=6200.0 sec, e=0.001, i=45, (RAAN) = 10.0, = 45, =315. Use an epoch and initial time of 1 Jan 2004 00:00:00 and the J2000 coordinate system. Use a stop time that is three orbit periods later. Assume the orbit is perturbed only by Earth oblateness, J2. Perform the following study: 1. Write the appropriate Matlab code and numerically integrate the equations of motion given above in the ECI Cartesian frame for three revolutions. Compute and store the orbit elements and ECI Cartesian coordinates at 1 min intervals. Starting with the same 1

    initial conditions, generate the ECI Cartesian coordinates for a two-body orbit for the same time interval. a. Plot the orbit elements, the energy deviation from its epoch value, and the deviation of the Z component of angular momentum from its epoch value for the perturbed orbit. Both the energy and hz should deviate from their epoch values by around 10-7 if your integration is done properly. Use the equation for energy given above. b. Generate a file of Cartesian position coordinate differences between the perturbed and two body orbits ( rpeturbed r2 body ). c. Code the RIC transformation matrix and plot the RIC differences using the twobody orbit as the reference. 2. Generate the perturbed orbit using STK. You will need to choose HPOP (High Precision Orbit Propagator) from the satellite basic properties. After choosing HPOP, click on the force model and choose the JGM2 gravity model. Set the maximum degree to 2 and the maximum order to 0. This will select J2 as the perturbing force. Turn off all other perturbing forces (drag, solar pressure, lunar, and solar gravity). Generate the RIC plots between the perturbed orbit and the two-body orbit as the reference using STK. Your RIC plots should look identical to those from MATLAB. You may want to do the STK plots first so that you will know if your MATLAB results are correct. 3. Your report should contain the following plots a. Plots of the orbit elements (a, e, i, , ,Tp), and the energy and hz deviations for your perturbed orbit generated with MATLAB. Do not plot the true anomaly. b. RIC differences between your perturbed and two body orbits. c. RIC differences between STKs perturbed and two-body orbit. 4. In addition to observations you wish to make, discuss the following topics in your report: a. Verify the differential equations of motion by generating the partial derivative, U , by hand as requested in the pre-lab assignment. X b. Describe the variations of the orbit elements and note the frequencies that are present, i. e., secular, short period, and long period. Can you detect any long period perturbations? Why or why not? Is the energy and hz constant? c. Discuss the RIC result. What is the direction of the dominant variation from the two-body orbit? Why does the major perturbation occur in this direction? d. Do the STK results differ from yours? If so, why?

    The RIC transformation for part 1 is obtained as follows (see also the handout on coordinate systems). Given the position and velocity vector of the reference orbit in the ECI frame, we have j r = Xi + Y + Zk,

    j r = Xi + Y + Zk, h = r r, then, r R= , r h C= , h I =CR In matrix notation, R R R R i i y z x XYZ j j I = I x I y I z = [T ]RIC C C C y z k C x k Note that the elements of the RIC transformation matrix at each point are the components of the RIC unit vectors. Hence, R X I = T XYZ Y [ ]RIC C Z

    Pre-Lab assignment:

    U Verify the equation for X by forming . X

    PROGRAMING NOTES The ode45 integrator in matlab is a variable step integrator. This means that it chooses its own step size based on default tolerances built into the program. These tolerances are very loose and as a result ode45 gives very poor results when integrating the orbit equations of motion. To solve this problem you can reset the tolerances by replacing the [t,x] statement with the following. tol = 1e-10; options = odeset('RelTol',tol,'AbsTol',tol); 3

    [t,x] = ode45('YourDerivativeFunction', time, x0, options) the time array can be defined as time=(0:60:18600) Where 0 is the initial time, 60s is the time step at which the ephemeris will be returned and 18600 sec is the final time.

    You might also like