One of the forces acting on an artificial satellite is the Earth’s magnetic field (Roy, 1982) . If the satellite has metal in its construction the Earth’s magnetic field induces eddy currents in the satellite; in addition, a slight retardation acts on the satellite. The changes in the orbit due to this force are small; in any real situation the effect of the small force is to cause small departures from Keplerian motion, but these deviations can hardly give rise to large-scale changes of the effect of another force, like as in our work the force arising from gravitational waves. To gain concept of the concentration of charged particles in space that might lead to electromagnetic forces on a satellite, it is useful to estimate the flux and energy of those particles in the field of the Earth associated the aurorally zone, solar flares, solar winds, cosmic rays, and the inner and outer Van Allen radiation belts, and recently the gravitational waves. The question of the effect of such particles and the associated magnetic fields on satellite dynamics is yet to be resolved. As noted by (Bourdeau et al., 1961) the motion of a satellite through a magnetic field B, produces an induced potential that is a function of position on the satellite surface; that is, a function of a vector R from the satellite center to any point on the satellite surface . The magnetic field strongly affects the motion of charged particles and the gravitational waves interact with magnetic field producing other effects on the motion of an orbiter in the Earth’s gravitational field. The influence of the geo-magnetic field in satellite charging is quite important and should be taken into account in any detailed treatment  - . In this work we investigate the perturbation on the elements of the satellite’s orbit moving with velocity V through a magnetic field and existing of gravitational waves, considering the propagation of the gravitational waves in the same direction of the magnetic field.
2. The Acceleration Components of Magnetic field
A force arises from the interaction of the Earth’s magnetic field and any electric charge of an orbiting body given by (Gelying and Westerman, 1971) 
where the charge is acquired by the satellite, its velocity and is the magnetic induction of the Earth’s field. In MKS units
where the vector potential is
is the magnetic moment of the Earth, and is the permeability of free space, carrying out the operations indicated in (1) through (3), we obtain
where m is the mass of satellite. As a first step to analyze effect of the magnetic field on the Kaplerian elements (a, e, i, , , ), it is necessary to represent the acceleration exerted by the magnetic field on a satellite as a function of the radial, orthogonal, and normal, perturbing components (Baker, 1967) . This can be done in terms of either the true anomaly f, or eccentric anomaly E. The resolving components of the disturbing force in the direction of the satellite radius-vector S, tangential to the orbit T and perpendicular to it W, are given as
where is the parameter, a is the semi-major axis, e is the eccentricity of the orbit, is the Earth’s gravitational constant and f is the true anomaly.
3. The Change of the Orbital Elements
According to Gauss form of Lagrange planetary equations under the action of the geomagnetic field and normal incident of gravitational waves the changes of the osculating elements for an elliptical orbit have the form
where ( , , ) and ( , , ) are referred to the acceleration components for the magnetic field and gravitational waves in the radial direction S, the transverse direction T at right angle to S in the orbital plane and perpendicular to the orbital plane W as in Figure 1.
Where , , and i are the angles of longitude of ascending node, argument of perigee, true anomaly and inclination respectively. ( , , ) are the resolving acceleration components of geomagnetic force which derived in Equation (5). The resolving acceleration components of gravitational waves force are
Figure 1. The disturbing force of gravitational waves in (S, T, W) directions.
, and are the components of the acceleration of normal ancient gravitational waves in (x, y, z) coordinates and
where is the frequency of the wave, and are the phase difference, and are the amplitude of the wave in the two orthogonal directions in the transverse plane . Therefore
We have and are the unit vectors in the direction of r, normal to r in the orbital plane and normal to the orbital plane respectively,
It follows that
Therefore from Equations ((8) to (13)) Equation (7) yield to
For the approximate integration of Equations (6), we expand all the functions of the orbital coordinates on the right sides in series in powers of e. The coefficients in these series will be trigonometric functions of the mean anomaly
where n the mean angular velocity of satellite, t the initial time of the motion, and is the time of perigee passage. Without loss of generality, we can take the satellite lies on the line of nodes at the initial time, so and Equations (6) become in the form
Because of the smallness of the eccentricity e, it is of interest to consider such effects only for small powers of e. Therefore Equation (17) is represented the first order approximations of Equation (6) and the effect of geomagnetic and GW on the orbital elements. A general analysis of Equation (17) shows that the semi-major axis a and the position of the orbital plane in space determined by the angles and i, change for a wave frequency, wave amplitude and mass ratio of satellite due to the induction of Earth’s magnetic field. For , the plane of orbit of which is parallel to the direction of the gravitational waves and the direction of magnetic field given by the quantities and , in this case only the longitude of ascending node changes and the shape of the orbit is constant. If the initial data for the orbit and the waves are defined by , , and then the position of the orbital plane is constant in space and only the semi-major axis a changes without a change in the other parameters. For the initial data , and there is only a deviation of the orbital plane from the wave direction, determined by the angle i. For the case and only eccentricity changes. In case and , this time only the angular distance changes.
4. Numerical Simulation
A numerical simulation is developed to test several of the above results. The simulation is a Runge-Kutta 4th order integration of system defined by Equations (17) with respect to the mean anomaly M performed by MATHEMATICA V10. The simulation is valid for any orbit for charged satellite in a non-tilted dipole field with GW. The results of two different situations are presented here for polar and equatorial orbits. Considering the period of satellite is four hours and with semi-major 12,600 km. The effects of the magnetic field and the gravitational waves are proportional to the charge to mass ratio and the frequency of the gravitational waves. Considering the wave’s frequency of the same order of magnitude of the magnetic field 10−9, wave’s amplitude of order 10−21 as coming from the bursts sources, and the mass ratio for satellite has mass 300 kg and . This means that with two parameters we determined the perturbations on the orbital motion and we can control on it. Also due to the interaction of GW with magnetic field the variation on the orbital element will periodic with the time.
Figure 2 displays the variation in the orbital plane due the perturbation of in radian for polar orbits with inclination , or 180° or 360° and .
Figure 3 displays the variation of the semi-major a in radians for the polar circular orbit with , and .
Figure 4 displays the variation of the inclination i in radian for the polar circular orbit with , and .
Figure 2. The variations of the orbital plane by angle Ω for polar orbit.
Figure 3. The variations of the semi-major α for polar orbit.
Figure 5 displays the variation of the eccentricity e for equatorial orbit with the case and .
Figure 6 displays the variation of the angular distance in radian for equatorial orbit with the case and .
Figure 4. The variations of the inclination i for polar orbit.
Figure 5. The variations of the eccentricity e for equatorial orbit.
Figure 6. The variations of the angular distance ω for equatorial orbit.
5. Discussion and Conclusion
Summarizing the results we conclude that there is secular effect due to geomagnetic and gravitational waves on the orbital plane of polar orbits (the longitude of node Ω and inclination i) and on the shape (the semi-major axis a) according to phase and phase difference of gravitational waves. For equatorial orbits the secular effect will on the size on the orbit (the eccentricity e and the argument of perigee ω). The amount of the variation of the orbital elements depends on the direction of propagation of gravitational waves and on the frequency and its phase difference which changes the amount of charges on the satellite and consequently the perturbations on the orbital elements. This effect is small but it is important for studying the effect of GW with the magnetic field of the Earth during interval of time and their effects on artificial satellites.
 Vokrouhlicky, D. (1989) The Geomagnetic Effects on the Motion of an Electrically Charged Artificial Satellite. Celestial Mechanics and Dynamical Astronomy, 46, 85-104.
 Streetman, B. and Peck, M.A. (2006) New Synchronous Orbits Using the Geomagnetic Lorentz Force. AIAA Guidance, Navigation, and Control Conference, Keystone, 21-24 August 2006, 1677.