1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: Adding Perturbing Accelerations to J2

  1. Apr 12, 2017 #1
    1. The problem statement, all variables and given/known data
    I'm numerically integrating Gaussian Variational Equations on MATLAB. I have a script which, given initial values for a, e, i, Ω (RAAN), ω (argument at perigee), and θ (true anomaly), will compute the r and v vectors through numerical integration. This is done with J2 and thrust (given a value Tr, Tv, Th, radial thrust, along-V thrust, and along-h thrust per unit mass, respectively.) I'm then plotting the orbital elements obtained via GVE's and the orbital elements obtained via conversion of state vector on the same plot to ensure that they are equivalent. This needs to be done with J2 and additional perturbing accelerations.
    2. Relevant equations
    Note that there is a typo in Eq. (12.88). The final line reads ps but should read pw. Also, the equation I used for da/dt is as follows:

    \frac{da}{dt} = \frac{2a^2esin(\theta)}{h}p_r + \frac{2a^3\sqrt{1-e^2}}{hr}p_s

    3. The attempt at a solution
    My scripts work when I run them without including the perturbing accelerations, but including J2. The orbital elements found from conversion of state vector are equivalent to those found from integrating GVE's. I was given a function that will apply the thrust to the state vector version. The only thrust that comes into play is the thrust acting along V. The author multiplied Tv by v, and that was all. I'm not exactly sure how to apply this thrust to the GVE's, since I am given orbital elements, and as such, the v vector is not explicitly defined. Also, I'm not sure at what point the v vector comes into play in this equation set? Regardless, I've tried taking the norm of the thrust elements and adding them to J2, and I've also tried adding each component to pr, ps, and pw. I think a major problem is that I don't exactly understand what the number (0.00108263) associated with J2 is. I think that it somehow represents the effect of the oblateness of earth, but I'm not sure exactly how. How does it relate to perturbing accelerations? How could one combine the effects of J2 and perturbing accelerations? If necessary, I can provide the MATLAB files associated with this problem. Thanks for reading!
    Last edited: Apr 12, 2017
  2. jcsd
  3. Apr 12, 2017 #2
    Not sure how to edit my post, but I think that the script that factors in the perturbing acceleration into the orbital elements found through conversion of state vector actually just uses one perturbing acceleration at a time. This is due to the problem statement, found below (#2).
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted