The present invention refers to a system for locating the barycenter of at least one object orbiting in space and to a related process of physical and mechanical characterization of the identified object.
It is known that operations like recovery, capture or deviation of the orbit of unknown and non-collaborative space objects are affected by inertial and gyroscopic effects which depend on quantities like: mass, inertial tensor, barycenter position, orientation of main inertial axes.
These quantities, in general, cannot be directly measured. A possible way for measuring them provides for observing the object to be identified by means of remote sensors, typically video cameras, which can be placed on board the observing satellite.
Known image analysis techniques allow observing particular points of the object to be identified, and measuring the followed trajectories by defining them in a reference system integral with an observer or with respect to any other reference system.
EP2340998 discloses a control system of an aerospace vehicle based on an on-line estimation of inertia. The system comprises a processor on board the aerospace vehicle adapted to on-line estimate the inertia of the aerospace vehicle and to determine the angular position and the angular speed of the aerospace vehicle. Through the measure of the attitude of the aerospace vehicle and the on-line estimation of inertia, it is possible to check movement and orientation of the aerospace vehicle. In particular, the processor allows generating a matrix of inertia with respect to a center of mass of the aerospace vehicle by using an estimation of the angular speed and an estimated measure of the moment of an actuator placed on board the aerospace vehicle.
EP2340998 however does not solve the problem of determining the barycenter of a satellite or body lacking the necessary means for generating the matrix of inertia.
The problem is particularly felt in the aerospace sector for capturing and removing orbiting debris by means of systems for recognizing, capturing and removing orbiting objects which populate the space surrounding the Earth, and which can be configured as potential sources of risks for present and future aerospace activities.
From 4 Oct. 1957, day in which the first Sputnik artificial satellite was launched by the Soviet Union, more than 4800 launches have been performed, which have orbited about 15000 satellites, of which, so far, only 800 are still operating. This means that a great number of objects which have concluded their missions are still orbiting around the Earth. This set of objects built by men is identified as “Space Debris”.
Apart from the objects orbiting at the highest heights of terrestrial atmosphere, all other debris are aimed to re-enter the Earth due to the natural decay process effect.
As an average, every day, an object catalogued as bigger than 10 cm of diameter falls onto the Earth. The risk of re-enter is given not only by the mechanical impact, but also by the chemical or radiological environmental contamination. Moreover, space debris can compromise the
operation of active satellites, by damaging these latter ones due to collision or reducing their performance by depositing onto the surfaces of optical systems, degrading solar panels and antennas, thereby reducing their transmission capacity and generating interferences with the signals.
Therefore, there is the need of developing the necessary technology for capturing and removing big-sized space debris.
The problem of identifying the barycenter of an orbiting object occurs also in relation to the prevention of risks linked to the so-called NEOs (Near Earth Objects).
CN-A-101 435 732 discloses a system for the detection of the center of gravity according to the preamble of Claim 1.
Object of the present invention is providing a system capable of identifying the barycenter of an orbiting body or anyway a moving body in space.
A further object is providing a process adapted to allow identifying the barycenter.
The above and other objects and advantages of the invention, as will appear from the following description, are obtained with a system for locating the barycenter of at least one object orbiting in space, as claimed in claim 1.
Moreover, the above and other objects and advantages of the invention are obtained with a process of physical and mechanical characterization of the identified object, as claimed in claim 5.
Preferred embodiments and non-trivial variations of the present invention are the subject matter of the dependent claims.
It is intended that all enclosed claims are an integral part of the present description.
It will be immediately obvious that numerous variations and modifications (for example related to shape, sizes, arrangements and parts with equivalent functionality) could be made to what is described, without departing from the scope of the invention as appears from the enclosed claims.
The present invention will be better described by some preferred embodiments thereof, provided as a non-limiting example, with reference to the enclosed drawings, in which:
With reference to
second means adapted to determine the instantaneous rotation axes of the identified object associated with the trajectory, determine a segment perpendicular to each pair of instantaneous rotation axes in a sequence and locate the mean point of the segment.
The system for locating the barycenter of at least one object orbiting in space, moreover, comprises third means adapted to compute a discrete function d(tk) of the length of the segments, compute an envelope curve c of the local maxima of the function d(tk) and determine the minimum of the envelope curve c for locating the barycenter G(tk) of the identified object.
The remote sensor 1 can be an optical system or the like.
The station 2 can be a space craft adapted to follow the identified object. Alternatively, the station 2 can be an earth base.
A process of physical and mechanical characterization of the orbiting object identified through the system is based on the consideration that typical motions of an orbiting body can be decomposed into barycenter motion along the orbit and into body rotation motions around the barycenter. Assuming that the body does not exchange with the outside other forces in addition to the mass ones, it is not possible that rotary motions occur around axes which do not pass by the barycenter.
With reference to
For each one of the observed points P1, P2 and P3, the coordinates must be available in a suitable reference system at following instants tk and the sampling frequency must be high enough.
Data are considered which are related to two following sampling intervals, [tk−1,tk] and [tk,tk+1]. In each interval, the motion can be described as rotary-translating, and its axis of instantaneous rotation is certain by using, for example, the process discovered by Ebherarter-Ravani which is known in the art.
By analyzing the two following intervals [tk−1,tk] and [tk,tk+1], two instantaneous rotation axes are obtained.
If the two instantaneous rotation axes are coplanar and not parallel, it is possible to locate the intersection point. It is moreover possible to state that the body motion is a pure rotation and that the body barycenter coincides with the located intersection point; it is not necessary, in this case, to perform other operations.
If, instead, the two axes are skewed, the body has a generic rotary-translating motion. Moreover, there is no intersection between these axes, but it is anyway possible to locate the segment which is perpendicular to both, whose length is the distance between the axes.
Taking into account a high enough sampling frequency, the distance is small and the position of the mean point of the segment perpendicular to both axes is an estimation of the position of the body barycenter G(tk) at instant tk.
With reference to
The discrete function d(tk) has several local maxima and the corresponding minima are next to zero.
In order to define the barycenter position, the envelope curve c is constructed for the local maxima, and its minimum point is located. The estimation computed next to the local minimum is assumed as best estimation of the barycenter position.
With reference to
F1) acquiring, at a suitable sampling frequency, data related to space coordinates of the positions assumed in time by certain points Pi of the identified object through the sensor 1 of the system according to the present invention;
F2) verifying the availability of trajectory data of at least three of the points Pi at three following instants tk−1, tk, tk+1 and, if available, reconstructing the trajectory followed by at least three Pi. Should this condition not be verified, wait for new trajectory data of three points Pi at following instants tk, tk+1, tk+2;
F3) locating the instantaneous rotation axes, typically two, associated with the trajectory of the at least three points Pi through the second means of the system according to the present invention;
IF-A) verifying the co-planarity condition of axes located in step F3 through the second means of the system according to the present invention. If a co-planarity is found, proceed to step F8), otherwise to step F4);
F4) locating the segment perpendicular to each pair of instantaneous rotation axes in a sequence and computing its length d(tk) through the second means of the system according to the present invention;
F5) determining the mean point G(tk) of the segment located through the second means of the system according to the present invention;
IF B) verifying whether conditions occur for proceeding with the process;
F6) constructing the envelope curve c of local maxima of the discrete function d(tk) through the third means of the system according to the present invention;
F7) searching for the minimum of the envelope curve c and locating the instant tk* in which the minimum occurs, and locating the best estimation of the object barycenter as value of the discrete function G(tk) obtained in step F5) at instant tk*, for example through the third means of the system according to the present invention;
F8) setting the position of the barycenter G(tk*) available, for example through the third means of the system according to the present invention.
The barycenter position, therefore, is expressed in three reference systems: inertial, integral with the observer, integral with the body.
The sequence of positions assumed by the barycenter, when can be measured with a satisfactory accuracy, is part of the orbit on which the observed body moves. These data then allow, with known processes, computing the actual orbit on which the observed body moves.
It is advisable to state that, in the trajectory of each monitored point, there are instants in which the coordinates of the i-th point cannot be measured, because the point itself is hidden to the observer.
The algorithm can be applied at a given instant tk when the coordinates of at least three different points Pi are available, at instant tk, at previous instant tk−1 and at following instant tk+1. If this condition is not verified, the chance can be verified of applying it at the immediately following instant tk+1, thereby verifying the availability of data at instants tk, tk+1, tk+2. This iteration will have to be repeated till three instants in which required data are available can be identified.
A first sequence of steps is described, belonging to the process dealing with the location of instantaneous rotation axes associated with the trajectory of the at least three points Pi. Starting data are composed of the coordinates of at least three points P1, P2 and P3 at two generic following time instants tk and tk+1.
Assuming that the instants are near enough to be able to approximate the rotary-translating motion with a rotary motion characterized by an instantaneous rotation axis, the algorithm determines the axis.
Firstly, vectors are built which describe the displacement of each point Pi in the time interval tk, tk+1. The knowledge of vectors designated as gi allows locating a plane ε perpendicular to the body rotation axis.
The projections ri are built for two of the three vectors gi on plane ε, and their mean points are located.
Then, perpendicular lines are traced to vectors ri passing by their mean points, and their intersection point Q is located. The point is the trace of the rotation axis on plane ε.
The rotation axis is then determined depending on the knowledge of a plane perpendicular thereto and of the intersection point Q with this plane. Assuming that data related to the coordinates of points Pi are affected by noise, it is possible to determine the searched plane ε by means of a minimum square estimation of results which can be obtained, taking into account a number of points which exceeds three.
The search for point Q, trace of the rotation axis on plane ε, is performed by projecting on plane ε all available vectors, locating for each pair of projections rm, rn the related intersection point Qm,n of the perpendicular lines passing by the mean points, and finally computing the geometric barycenter of the set of located points.
With reference to
F3.1) computing the vectors gi adapted to connect homologous points Pi(tk) and Pi(tk+1);
F3.2) computing the plane ε perpendicular to the body rotation axis starting from vectors gi;
F3.3) projecting each pair of points Pi(tk) and Pi(tk+1) on piano ε;
F3.4) constructing the axes of the segments connecting homologous points Pi(tk) and Pi(tk+1) projected on plane ε;
F3.5) intersecting all possible pairs of axes of the segments;
F3.6) estimating a point belonging the rotation axis by computing the geometric barycenter of all intersection points;
F3.7) constructing the rotation axis starting from collected information.
It is wholly clear for a skilled person in the art that step F3) can be performed through any other algorithm suitable for the object of the present invention. A second sequence of steps belonging to the process dealing with the definition of the position of the barycenter is described. The envelope curve c is constructed for the local maxima of the discrete function d(tk) and its minimum point is located. The estimation computed next to the local minimum point tk* is assumed as best estimation of the barycenter G(tk*) position.
With reference to
F7.1) locating the local maxima of said function d(tk);
F7.2) evaluating the curve c, as envelope of local maxima;
F7.3) evaluating the time instant tk* in which the curve c has an absolute minimum;
F7.4) evaluating the function G(tk) in the found instant tk*.
Number | Date | Country | Kind |
---|---|---|---|
TO2014A000550 | Jul 2014 | IT | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IT2015/000177 | 7/7/2015 | WO | 00 |