Methodology of determination of balancing weights mounting places inside spacecraft compartments

Cover Page

Cite item

Abstract

The paper presents a methodology for determining the mounting places of balancing weights inside spacecraft compartments, based on the use of a combination of methods of analytical and computational geometry, mathematical programming and computer graphics. The use of balancing weights is necessary to ensure the required position of the center of gravity of the compartment and the product as a whole. When using the methodology, the problems of ensuring a minimum mass of balancing weights and reducing the labor intensity of developing options for their installation are solved to speed up the preparation and approval of design documentation. The balancing weight placement zone is considered as a set of spatial regions free from compartment structural elements and other component parts. To minimize the overall mass of the balancing weights by determining their placement locations on a coordinate grid, the balancing weight placement problem is proposed to be represented as a linear programming problem. For testing, a conical compartment of a product with a spherical bottom was used as an example. It was determined that the balancing weight placement zone should be located near the junction of the bottom and the hull shell. The configuration of the placement zone was identified, taking into account the surrounding structural elements. The coordinates for placing the balancing weights were determined, and their masses were selected. Testing has showed the performance of the proposed methodology and the algorithm based on it. Effective use of the methodology is possible with the availability of a specialized calculation software package.

Full Text

Introduction

When developing design documentation (DD) for spacecraft, at the conceptual or technical design stage, mass-centering and inertial requirements for the product and permissible deviations from them are determined. Ideally, these requirements should be met through a rational layout of the spacecraft's compartments. In practice, numerous errors (design, computational, manufacturing, metrological, etc.) lead to unacceptable deviations in the mass-centering and inertial characteristics (MCIC) of the product [1; 2]. To eliminate these, partial or complete reconfiguration of the compartments is carried out, and when this is impossible for technical and economic reasons, balancing weights (BW) are used.

Various balancing rigs are used to weigh and balance products and their components [2 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@ 11]. During spacecraft testing, they provide data on the actual values of mass, center of mass or center of gravity (COG), and moments of inertia. If unacceptable deviations from MCIC are detected, the layout or arrangement of BW is adjusted according to change notices.

There may be cases where the design documentation does not provide for the installation of a BW. During the working documentation stage, during its development or during prototype testing, the need for a BW may be identified. In such cases, a new design group is added to the working documentation, and technically, the number of BWs must be determined and locations for their installation on the spacecraft must be found. The task of locating payloads under such conditions increases the complexity of the design documentation development. A heuristic solution to this problem does not guarantee the use of the minimum number (mass) of BWs.

The objective of this work is to develop and refine a methodology for determining BW installation locations within spacecraft compartments. The objectives of this work include developing computational mathematical models, identifying sources and methods for processing initial data, developing an algorithm for a software package to be developed in the future, and testing (refining) the methodology using a spacecraft compartment as an example.

The relevance of this research topic, in addition to reducing the labor intensity of design documentation development and minimizing the weight of the BW, is determined by the fact that spacecraft balancing using additional BWs will be carried out on both existing and newly developed products. Therefore, a reliable mathematical apparatus and calculation software package for solving the problem of determining BW installation locations will be in demand for quite a long time.

This paper examines the static balancing of spacecraft. For questions on dynamic balancing of spacecraft, it is recommended to refer to works [2; 3; 5; 7 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@ 11]. Although determining the actual values of the maximum permissible coefficient of performance of products using dynamic balancing methods can be more accurate, as noted in [11], for large-sized and outsized products, where BW placement zones can cover large areas and, therefore, hundreds of small payloads of various configurations can be used, static balancing methods appear more rational in terms of the test base. This provision does not exclude the fact that dynamic balancing of the product must also be carried out on par with static balancing in order to ensure not only the required position of the COG, but also the required values of the moments of inertia.

This article presents mathematical models for calculating the COG and mass of a BW installation, determining the configuration of a BW placement zone, calculating the BW mass separately, presents a block diagram of the algorithm of actions according to the proposed method, provides a list of the necessary initial data indicating their possible sources, and presents the results of developing this method.

Center of Mass for Balancing Weight Installation

First, it is necessary to determine the coordinates of the COG for the BW installation, assuming that only one BW is required. This point is formed at the intersection of the centering line and the BW placement surface. The centering line is a line that passes through the specified COG of the spacecraft and the current COG of the spacecraft without the BW. The balancing weight placement surface is a theoretical surface on which the median number of balancing weight centers of mass is located. It is situated at a distance from the shell surface of the spacecraft compartment equal to half the thickness h BWavg MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamiAamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbGaamyyaiaadAhacaWGNbaa paqabaaaaa@3BD1@  of the standard balancing weights from the standard nomenclature, which are in the form of plates of equal length and width l BWavg MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbGaamyyaiaadAhacaWGNbaa paqabaaaaa@3BD5@ . In general, the surface is not smooth. The calculation scheme is shown in Fig. 1.

 

Рис. 1. Расчётная схема для определения ЦМ установки БГ

Fig. 1. Computational model for balancing weights (BW) center of gravity (COG) determination

 

The standard form of the equation of the centring line is

x BW x C x 0 x C = y BW y C y 0 y C = z BW z C z 0 z C x BW x C Δ r x = y BW y C Δ r y = z BW z C Δ r z , MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WG4bWaaSbaaSqaamaaBaaameaaqaaaaaaaaaWdbiaadkeacaWGxbaa paqabaaaleqaaOGaeyOeI0IaamiEamaaBaaaleaacaWGdbaabeaaaO qaaiaadIhadaWgaaWcbaGaaGimaaqabaGccqGHsislcaWG4bWaaSba aSqaaiaadoeaaeqaaaaakiabg2da9maalaaabaGaamyEamaaBaaale aadaWgaaadbaWdbiaadkeacaWGxbaapaqabaaaleqaaOGaeyOeI0Ia amyEamaaBaaaleaacaWGdbaabeaaaOqaaiaadMhadaWgaaWcbaGaaG imaaqabaGccqGHsislcaWG5bWaaSbaaSqaaiaadoeaaeqaaaaakiab g2da9maalaaabaGaamOEamaaBaaaleaadaWgaaadbaWdbiaadkeaca WGxbaapaqabaaaleqaaOGaeyOeI0IaamOEamaaBaaaleaacaWGdbaa beaaaOqaaiaadQhadaWgaaWcbaGaaGimaaqabaGccqGHsislcaWG6b WaaSbaaSqaaiaadoeaaeqaaaaakiabgkDiEpaalaaabaGaamiEamaa BaaaleaadaWgaaadbaWdbiaadkeacaWGxbaapaqabaaaleqaaOGaey OeI0IaamiEamaaBaaaleaacaWGdbaabeaaaOqaaiabgs5aejaadkha daWgaaWcbaGaamiEaaqabaaaaOGaeyypa0ZaaSaaaeaacaWG5bWaaS baaSqaamaaBaaameaapeGaamOqaiaadEfaa8aabeaaaSqabaGccqGH sislcaWG5bWaaSbaaSqaaiaadoeaaeqaaaGcbaGaeyiLdqKaamOCam aaBaaaleaacaWG5baabeaaaaGccqGH9aqpdaWcaaqaaiaadQhadaWg aaWcbaWaaSbaaWqaa8qacaWGcbGaam4vaaWdaeqaaaWcbeaakiabgk HiTiaadQhadaWgaaWcbaGaam4qaaqabaaakeaacqGHuoarcaWGYbWa aSbaaSqaaiaadQhaaeqaaaaakiaacYcaaaa@7B73@  (1)

where r 0 = x 0 y 0 z 0 T MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaacaaIWaaabeaakiabg2da9maabmaabaqbaeqabeWaaaqaaiaa dIhadaWgaaWcbaGaaGimaaqabaaakeaacaWG5bWaaSbaaSqaaiaaic daaeqaaaGcbaGaamOEamaaBaaaleaacaaIWaaabeaaaaaakiaawIca caGLPaaadaahaaWcbeqaaiaadsfaaaaaaa@416E@  is the coordinate vector of the specified position of the spacecraft’s COG; r C = x C y C z C T MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaacaWGdbaabeaakiabg2da9maabmaabaqbaeqabeWaaaqaaiaa dIhadaWgaaWcbaGaam4qaaqabaaakeaacaWG5bWaaSbaaSqaaiaado eaaeqaaaGcbaGaamOEamaaBaaaleaacaWGdbaabeaaaaaakiaawIca caGLPaaadaahaaWcbeqaaiaadsfaaaaaaa@41A6@  is the coordinate vector of the current position of the spacecraft’s COG without BW; r BW = x BW y BW z BW T MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGccqGH9aqpdaqa daqaauaabeqabmaaaeaacaWG4bWaaSbaaSqaa8qacaWGcbGaam4vaa WdaeqaaaGcbaGaamyEamaaBaaaleaapeGaamOqaiaadEfaa8aabeaa aOqaaiaadQhadaWgaaWcbaWdbiaadkeacaWGxbaapaqabaaaaaGcca GLOaGaayzkaaWaaWbaaSqabeaacaWGubaaaaaa@459E@  is the coordinate vector of the COG of the BW installation; Δr= Δ r x Δ r y Δ r z T = r 0 r C MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam OCaiabg2da9maabmaabaqbaeqabeWaaaqaaiabgs5aejaadkhadaWg aaWcbaGaamiEaaqabaaakeaacqGHuoarcaWGYbWaaSbaaSqaaiaadM haaeqaaaGcbaGaeyiLdqKaamOCamaaBaaaleaacaWG6baabeaaaaaa kiaawIcacaGLPaaadaahaaWcbeqaaiaadsfaaaGccqGH9aqpcaWGYb WaaSbaaSqaaiaaicdaaeqaaOGaeyOeI0IaamOCamaaBaaaleaacaWG dbaabeaaaaa@4CA0@  is the vector of deviations of the coordinates of the spacecraft’s COG without BW from the specified position.

Geometrically, the spacecraft’s compartment hull may take the form of a hypersurface MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@  a combination of several segments of various elementary surfaces. For design purposes, the shapes of such segments of the spacecraft’s hull can be adequately described by surfaces of no higher than second order [12; 13]. In general form, the surface equation takes the form

U T AU+2bU+ a 44 =0, MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyvamaaCa aaleqabaGaamivaaaakiaadgeacaWGvbGaey4kaSIaaGOmaiaadkga caWGvbGaey4kaSIaamyyamaaBaaaleaacaaI0aGaaGinaaqabaGccq GH9aqpcaaIWaGaaiilaaaa@42EB@  (2)

where U= u x u y u z T MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyvaiabg2 da9maabmaabaqbaeqabeWaaaqaaiaadwhadaWgaaWcbaGaamiEaaqa baaakeaacaWG1bWaaSbaaSqaaiaadMhaaeqaaaGcbaGaamyDamaaBa aaleaacaWG6baabeaaaaaakiaawIcacaGLPaaadaahaaWcbeqaaiaa dsfaaaaaaa@4120@  is the vector of coordinates of points on a segment of the spacecraft’s hull surface;

A= a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyqaiabg2 da9maabmaabaqbaeqabmWaaaqaaiaadggadaWgaaWcbaGaaGymaiaa igdaaeqaaaGcbaGaamyyamaaBaaaleaacaaIXaGaaGOmaaqabaaake aacaWGHbWaaSbaaSqaaiaaigdacaaIZaaabeaaaOqaaiaadggadaWg aaWcbaGaaGOmaiaaigdaaeqaaaGcbaGaamyyamaaBaaaleaacaaIYa GaaGOmaaqabaaakeaacaWGHbWaaSbaaSqaaiaaikdacaaIZaaabeaa aOqaaiaadggadaWgaaWcbaGaaG4maiaaigdaaeqaaaGcbaGaamyyam aaBaaaleaacaaIZaGaaGOmaaqabaaakeaacaWGHbWaaSbaaSqaaiaa iodacaaIZaaabeaaaaaakiaawIcacaGLPaaaaaa@50B8@  is the affine matrix of coefficients for the quadratic term of the equation;

b= a 14 a 24 a 34 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOyaiabg2 da9maabmaabaqbaeqabeWaaaqaaiaadggadaWgaaWcbaGaaGymaiaa isdaaeqaaaGcbaGaamyyamaaBaaaleaacaaIYaGaaGinaaqabaaake aacaWGHbWaaSbaaSqaaiaaiodacaaI0aaabeaaaaaakiaawIcacaGL Paaaaaa@415F@  is the vector of coefficients for the linear term of the equation;

a jk :j,k= 1,4 ¯ MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGQbGaam4AaaqabaGccaGG6aGaamOAaiaacYcacaWGRbGa eyypa0Zaa0aaaeaacaaIXaGaaiilaiaaisdaaaaaaa@3FA0@  is a set of coefficients for the equation of the spacecraft’s hull surface segment.

The coefficients of equation (2) are determined according to Table 1 [14], depending on the type of surface segment.

 

Table 1

Coefficients of the equation for a spacecraft hull surface segment

Surface

a11a22a33a14a24a34a44

Plane

0

0

0

nxnynznxxCOMPT+nyyCOMPT+nzzCOMPT

Sphere

1

1

1

xCOMPTyCOMPTzCOMPTxCOMPT2+yCOMPT2+zCOMPT2RCOMPT2

Cylinder

0

1

1

0

yCOMPTzCOMPTyCOMPT2+zCOMPT2RCOMPT2

Cone

tg2φ

1

1

xCOMPTtg2φyCOMPTzCOMPTyCOMPT2+zCOMPT2xCOMPT2tg2φ

 

In Table 1, n x n y n z MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaeWaaeaafa qabeqadaaabaGaamOBamaaBaaaleaacaWG4baabeaaaOqaaiaad6ga daWgaaWcbaGaamyEaaqabaaakeaacaWGUbWaaSbaaSqaaiaadQhaae qaaaaaaOGaayjkaiaawMcaaaaa@3E26@  represents the cosines of the normal to the plane; r COMPT = x COMPT y COMPT z COMPT MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaaqaaaaaaaaaWdbiaadoeacaWGpbGaamytaiaadcfacaWGubaa paqabaGccqGH9aqpdaqadaqaauaabeqabmaaaeaacaWG4bWaaSbaaS qaa8qacaWGdbGaam4taiaad2eacaWGqbGaamivaaWdaeqaaaGcbaGa amyEamaaBaaaleaapeGaam4qaiaad+eacaWGnbGaamiuaiaadsfaa8 aabeaaaOqaaiaadQhadaWgaaWcbaWdbiaadoeacaWGpbGaamytaiaa dcfacaWGubaapaqabaaaaaGccaGLOaGaayzkaaaaaa@4E7C@  represents the coordinates of a point on the plane, the centre of the sphere, the cylinder or the cone; R COMPT MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOuamaaBa aaleaaqaaaaaaaaaWdbiaadoeacaWGpbGaamytaiaadcfacaWGubaa paqabaaaaa@3B67@  represents the radius of the spacecraft’s compartment; φ MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOXdaaa@3734@  represents the angle of the cone’s semi-apex.

The coordinates of the center of gravity (COG) of the BW assembly can be found using a vector transformation of the canonical equation of the centering line (1):

r BW = r 0 +rΔ r BW , MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGccqGH9aqpcaWG YbWaaSbaaSqaaiaaicdaaeqaaOGaey4kaSIaeyOaIyRaamOCaiabgs 5aejaadkhadaWgaaWcbaWdbiaadkeacaWGxbaapaqabaGccaGGSaaa aa@444A@  (3)

where r= λ χ γ MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeyOaIyRaam OCaiabg2da9maabmaabaqbaeqabeWaaaqaaiaabU7aaeaacaqGhpaa baGaae4SdaaaaiaawIcacaGLPaaaaaa@3EDA@  represents certain parameters of the centering line equation, depending on the shape of the spacecraft body surface segment; Δ r BW = ξ ζ ς T MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam OCamaaBaaaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGccqGH 9aqpdaqadaqaauaabeqabmaaaeaacaqG+oaabaGaaeOTdaqaaiaabk 8aaaaacaGLOaGaayzkaaWaaWbaaSqabeaacaWGubaaaaaa@41EA@  represents the deviations of the spacecraft’s COG coordinates (without BW) from the specified position, depending on the shape of the spacecraft body surface segment.

Substituting equation (3) into equation (2) produces quadratic equations in the unknown parameters of the centering line:

ξ 2 λ 2 +2 r 0 r COMPT Т ξλ+ R 2 =0  sphere, ζ 2 χ 2 +2 r COMPT Т ζχ+ Р 2 =0  cylinder, ς 2 A γ 2 2ς r COMPT Т AE γ+ Q 2 =0  cone. MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaab67adaahaaWcbeqaaiaaikdaaaGccaqG7oWaaWbaaSqabeaa caaIYaaaaOGaey4kaSIaaGOmamaabmaabaGaamOCamaaBaaaleaaca aIWaaabeaakiabgkHiTiaadkhadaWgaaWcbaaeaaaaaaaaa8qacaWG dbGaam4taiaad2eacaWGqbGaamivaaWdaeqaaaGccaGLOaGaayzkaa WaaWbaaSqabeaacaWGIqaaaOGaaeOVdiaabU7acqGHRaWkcaWGsbWa aWbaaSqabeaacaaIYaaaaOGaeyypa0JaaGimaiaabccacqGHsislca qGGaWdbiaabohacaqGWbGaaeiAaiaabwgacaqGYbGaaeyza8aacaqG SaaabaGaaeOTdmaaCaaaleqabaGaaGOmaaaakiaabE8adaahaaWcbe qaaiaaikdaaaGccqGHRaWkcaaIYaGaamOCamaaDaaaleaadaWgaaad baWdbiaadoeacaWGpbGaamytaiaadcfacaWGubaapaqabaaaleaaca WGIqaaaOGaaeOTdiaabE8acqGHRaWkcaWGGqWaaWbaaSqabeaacaaI YaaaaOGaeyypa0JaaGimaiaabccacqGHsislcaqGGaWdbiaabogaca qG5bGaaeiBaiaabMgacaqGUbGaaeizaiaabwgacaqGYbWdaiaabYca aeaacaqGcpWaaWbaaSqabeaacaaIYaaaaOGaamyqaiaabo7adaahaa WcbeqaaiaaikdaaaGccqGHsislcaaIYaGaaeOWdmaabmaabaGaamOC amaaDaaaleaadaWgaaadbaWdbiaadoeacaWGpbGaamytaiaadcfaca WGubaapaqabaaaleaacaWGIqaaaOGaamyqaiabgkHiTiaadweaaiaa wIcacaGLPaaacaqGZoGaey4kaSIaamyuamaaCaaaleqabaGaaGOmaa aakiabg2da9iaaicdacaqGGaGaeyOeI0Iaaeiia8qacaqGJbGaae4B aiaab6gacaqGLbWdaiaab6caaaGaay5Eaaaaaa@93E0@  (4)

In system (4), the following convolutions are used for the equation on the sphere:

ξ x 2 + ξ y 2 + ξ z 2 = ξ Т ξ= ξ 2 , x 0 x COMPT ξ x + y 0 y COMPT ξ y + z 0 z COMPT ξ z = r 0 r COMPT T ξ, x 0 x COMPT 2 + y 0 y COMPT 2 + z 0 z COMPT 2 R COMPT h BWavg 2 = r 0 r COMPT 2 R COMPT h BWavg 2 = R 2 . MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaab67adaqhaaWcbaGaamiEaaqaaiaaikdaaaGccqGHRaWkcaqG +oWaa0baaSqaaiaadMhaaeaacaaIYaaaaOGaey4kaSIaaeOVdmaaDa aaleaacaWG6baabaGaaGOmaaaakiabg2da9iaab67adaahaaWcbeqa aiaadkcbaaGccaqG+oGaeyypa0JaaeOVdmaaCaaaleqabaGaaGOmaa aakiaacYcaaeaadaqadaqaaiaadIhadaWgaaWcbaGaaGimaaqabaGc cqGHsislcaWG4bWaaSbaaSqaaabaaaaaaaaapeGaam4qaiaad+eaca WGnbGaamiuaiaadsfaa8aabeaaaOGaayjkaiaawMcaaiaab67adaWg aaWcbaGaamiEaaqabaGccqGHRaWkdaqadaqaaiaadMhadaWgaaWcba GaaGimaaqabaGccqGHsislcaWG5bWaaSbaaSqaa8qacaWGdbGaam4t aiaad2eacaWGqbGaamivaaWdaeqaaaGccaGLOaGaayzkaaGaaeOVdm aaBaaaleaacaWG5baabeaakiabgUcaRmaabmaabaGaamOEamaaBaaa leaacaaIWaaabeaakiabgkHiTiaadQhadaWgaaWcbaWdbiaadoeaca WGpbGaamytaiaadcfacaWGubaapaqabaaakiaawIcacaGLPaaacaqG +oWaaSbaaSqaaiaadQhaaeqaaOGaeyypa0ZaaeWaaeaacaWGYbWaaS baaSqaaiaaicdaaeqaaOGaeyOeI0IaamOCamaaBaaaleaapeGaam4q aiaad+eacaWGnbGaamiuaiaadsfaa8aabeaaaOGaayjkaiaawMcaam aaCaaaleqabaGaamivaaaakiaab67acaqGSaaabaWaaeWaaeaacaWG 4bWaaSbaaSqaaiaaicdaaeqaaOGaeyOeI0IaamiEamaaBaaaleaape Gaam4qaiaad+eacaWGnbGaamiuaiaadsfaa8aabeaaaOGaayjkaiaa wMcaamaaCaaaleqabaGaaGOmaaaakiabgUcaRmaabmaabaGaamyEam aaBaaaleaacaaIWaaabeaakiabgkHiTiaadMhadaWgaaWcbaWdbiaa doeacaWGpbGaamytaiaadcfacaWGubaapaqabaaakiaawIcacaGLPa aadaahaaWcbeqaaiaaikdaaaGccqGHRaWkdaqadaqaaiaadQhadaWg aaWcbaGaaGimaaqabaGccqGHsislcaWG6bWaaSbaaSqaa8qacaWGdb Gaam4taiaad2eacaWGqbGaamivaaWdaeqaaaGccaGLOaGaayzkaaWa aWbaaSqabeaacaaIYaaaaOGaeyOeI0YaaeWaaeaacaWGsbWaaSbaaS qaa8qacaWGdbGaam4taiaad2eacaWGqbGaamivaaWdaeqaaOGaeyOe I0IaamiAamaaBaaaleaapeGaamOqaiaadEfacaWGHbGaamODaiaadE gaa8aabeaaaOGaayjkaiaawMcaamaaCaaaleqabaGaaGOmaaaakiab g2da9maabmaabaGaamOCamaaBaaaleaacaaIWaaabeaakiabgkHiTi aadkhadaWgaaWcbaWdbiaadoeacaWGpbGaamytaiaadcfacaWGubaa paqabaaakiaawIcacaGLPaaadaahaaWcbeqaaiaaikdaaaGccqGHsi sldaqadaqaaiaadkfadaWgaaWcbaWdbiaadoeacaWGpbGaamytaiaa dcfacaWGubaapaqabaGccqGHsislcaWGObWaaSbaaSqaa8qacaWGcb Gaam4vaiaadggacaWG2bGaam4zaaWdaeqaaaGccaGLOaGaayzkaaWa aWbaaSqabeaacaaIYaaaaOGaeyypa0JaamOuamaaCaaaleqabaGaaG Omaaaakiaac6caaaGaay5Eaaaaaa@CF73@

In system (4), the following convolutions are used for the equation on the cylinder:

0+ ζ y 2 + ζ z 2 = ζ T ζ= ζ 2 , 0+ y COMPT ζ y + z COMPT ζ z = r COMPT T ζ, 0+ y COMPT 2 + z COMPT 2 R COMPT 0,5 l BWavg 2 = r COMPT R COMPT 0,5 l BWavg 2 = P 2 . MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaaicdacqGHRaWkcaqG2oWaa0baaSqaaiaadMhaaeaacaaIYaaa aOGaey4kaSIaaeOTdmaaDaaaleaacaWG6baabaGaaGOmaaaakiabg2 da9iaabA7adaahaaWcbeqaaiaadsfaaaGccaqG2oGaeyypa0JaaeOT dmaaCaaaleqabaGaaGOmaaaakiaacYcaaeaacaaIWaGaey4kaSIaam yEamaaBaaaleaaqaaaaaaaaaWdbiaadoeacaWGpbGaamytaiaadcfa caWGubaapaqabaGccaqG2oWaaSbaaSqaaiaadMhaaeqaaOGaey4kaS IaamOEamaaBaaaleaapeGaam4qaiaad+eacaWGnbGaamiuaiaadsfa a8aabeaakiaabA7adaWgaaWcbaGaamOEaaqabaGccqGH9aqpcaWGYb Waa0baaSqaamaaBaaameaapeGaam4qaiaad+eacaWGnbGaamiuaiaa dsfaa8aabeaaaSqaaiaadsfaaaGccaqG2oGaaeilaaqaaiaaicdacq GHRaWkcaWG5bWaa0baaSqaamaaBaaameaapeGaam4qaiaad+eacaWG nbGaamiuaiaadsfaa8aabeaaaSqaaiaaikdaaaGccqGHRaWkcaWG6b Waa0baaSqaamaaBaaameaapeGaam4qaiaad+eacaWGnbGaamiuaiaa dsfaa8aabeaaaSqaaiaaikdaaaGccqGHsisldaqadaqaaiaadkfada WgaaWcbaWdbiaadoeacaWGpbGaamytaiaadcfacaWGubaapaqabaGc cqGHsislcaaIWaGaaiilaiaaiwdacaWGSbWaaSbaaSqaa8qacaWGcb Gaam4vaiaadggacaWG2bGaam4zaaWdaeqaaaGccaGLOaGaayzkaaWa aWbaaSqabeaacaaIYaaaaOGaeyypa0JaamOCamaaBaaaleaapeGaam 4qaiaad+eacaWGnbGaamiuaiaadsfaa8aabeaakiabgkHiTmaabmaa baGaamOuamaaBaaaleaapeGaam4qaiaad+eacaWGnbGaamiuaiaads faa8aabeaakiabgkHiTiaaicdacaGGSaGaaGynaiaadYgadaWgaaWc baWdbiaadkeacaWGxbGaamyyaiaadAhacaWGNbaapaqabaaakiaawI cacaGLPaaadaahaaWcbeqaaiaaikdaaaGccqGH9aqpcaWGqbWaaWba aSqabeaacaaIYaaaaOGaaiOlaaaacaGL7baaaaa@A06F@

In system (4), the following convolutions are used for the equation on the cone:

ς x 2 + ς y 2 + ς z 2 = ς T ς= ς 2 , tg 2 φ 1 1 T =A, ς x x COMPT tg 2 φ1 + ς y y COMPT 1 + ς z z COMPT 1 =ς r COMPT T AE , x 0 x COMPT 2 tg 2 φ+ y 0 y COMPT 2 + z 0 z COMPT 2 2 x COMPT 2 tg 2 φ= r 0 r COMPT 2 A2 R COMPT 2 = Q 2 . MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaabk8adaqhaaWcbaGaamiEaaqaaiaaikdaaaGccqGHRaWkcaqG cpWaa0baaSqaaiaadMhaaeaacaaIYaaaaOGaey4kaSIaaeOWdmaaDa aaleaacaWG6baabaGaaGOmaaaakiabg2da9iaabk8adaahaaWcbeqa aiaadsfaaaGccaqGcpGaeyypa0JaaeOWdmaaCaaaleqabaGaaGOmaa aakiaacYcaaeaadaqadaqaauaabeqabmaaaeaacaqG0bGaae4zamaa CaaaleqabaGaaGOmaaaakiaabA8aaeaacaaIXaaabaGaaGymaaaaai aawIcacaGLPaaadaahaaWcbeqaaiaadsfaaaGccqGH9aqpcaWGbbGa aiilaaqaaiaabk8adaWgaaWcbaGaamiEaaqabaGcdaqadaqaaiaadI hadaWgaaWcbaaeaaaaaaaaa8qacaWGdbGaam4taiaad2eacaWGqbGa amivaaWdaeqaaOGaaeiDaiaabEgadaahaaWcbeqaaiaaikdaaaGcca qGgpGaeyOeI0IaaGymaaGaayjkaiaawMcaaiabgUcaRiaabk8adaWg aaWcbaGaamyEaaqabaGcdaqadaqaaiaadMhadaWgaaWcbaWdbiaado eacaWGpbGaamytaiaadcfacaWGubaapaqabaGccqGHsislcaaIXaaa caGLOaGaayzkaaGaey4kaSIaaeOWdmaaBaaaleaacaWG6baabeaakm aabmaabaGaamOEamaaBaaaleaapeGaam4qaiaad+eacaWGnbGaamiu aiaadsfaa8aabeaakiabgkHiTiaaigdaaiaawIcacaGLPaaacqGH9a qpcaqGcpWaaeWaaeaacaWGYbWaa0baaSqaamaaBaaameaapeGaam4q aiaad+eacaWGnbGaamiuaiaadsfaa8aabeaaaSqaaiaadsfaaaGcca WGbbGaeyOeI0IaamyraaGaayjkaiaawMcaaiaacYcaaeaadaqadaqa aiaadIhadaWgaaWcbaGaaGimaaqabaGccqGHsislcaWG4bWaaSbaaS qaa8qacaWGdbGaam4taiaad2eacaWGqbGaamivaaWdaeqaaaGccaGL OaGaayzkaaWaaWbaaSqabeaacaaIYaaaaOGaaeiDaiaabEgadaahaa WcbeqaaiaaikdaaaGccaqGgpGaey4kaSYaaeWaaeaacaWG5bWaaSba aSqaaiaaicdaaeqaaOGaeyOeI0IaamyEamaaBaaaleaapeGaam4qai aad+eacaWGnbGaamiuaiaadsfaa8aabeaaaOGaayjkaiaawMcaamaa CaaaleqabaGaaGOmaaaakiabgUcaRmaabmaabaGaamOEamaaBaaale aacaaIWaaabeaakiabgkHiTiaadQhadaWgaaWcbaWdbiaadoeacaWG pbGaamytaiaadcfacaWGubaapaqabaaakiaawIcacaGLPaaadaahaa WcbeqaaiaaikdaaaGccqGHsislcaaIYaGaamiEamaaDaaaleaadaWg aaadbaWdbiaadoeacaWGpbGaamytaiaadcfacaWGubaapaqabaaale aacaaIYaaaaOGaaeiDaiaabEgadaahaaWcbeqaaiaaikdaaaGccaqG gpGaeyypa0ZaaeWaaeaacaWGYbWaaSbaaSqaaiaaicdaaeqaaOGaey OeI0IaamOCamaaBaaaleaadaWgaaadbaWdbiaadoeacaWGpbGaamyt aiaadcfacaWGubaapaqabaaaleqaaaGccaGLOaGaayzkaaWaaWbaaS qabeaacaaIYaaaaOGaamyqaiabgkHiTiaaikdacaWGsbWaa0baaSqa amaaBaaameaapeGaam4qaiaad+eacaWGnbGaamiuaiaadsfaa8aabe aaaSqaaiaaikdaaaGccqGH9aqpcaWGrbWaaWbaaSqabeaacaaIYaaa aOGaaiOlaaaacaGL7baaaaa@D547@

The roots of the equations in system (4) are

λ= r 0 r COMPT T ξ± D sph ξ 2   sphere, χ= r COMPT Т ζ± D cyl ζ 2   cylinder, γ= ς r COMPT Т AE ± D cone ς 2 A   cone. MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaabU7acqGH9aqpdaWcaaqaaiabgkHiTmaabmaabaGaamOCamaa BaaaleaacaaIWaaabeaakiabgkHiTiaadkhadaWgaaWcbaaeaaaaaa aaa8qacaWGdbGaam4taiaad2eacaWGqbGaamivaaWdaeqaaaGccaGL OaGaayzkaaWaaWbaaSqabeaacaWGubaaaOGaaeOVdiabgglaXoaaka aabaGaamiramaaBaaaleaapeGaam4CaiaadchacaWGObaapaqabaaa beaaaOqaaiaab67adaahaaWcbeqaaiaaikdaaaaaaOGaaeiiaiabgk HiTiaabccapeGaae4CaiaabchacaqGObGaaeyzaiaabkhacaqGLbWd aiaabYcaaeaacaqGhpGaeyypa0ZaaSaaaeaacqGHsislcaWGYbWaa0 baaSqaamaaBaaameaapeGaam4qaiaad+eacaWGnbGaamiuaiaadsfa a8aabeaaaSqaaiaadkcbaaGccaqG2oGaeyySae7aaOaaaeaacaWGeb WaaSbaaSqaa8qacaWGJbGaamyEaiaadYgaa8aabeaaaeqaaaGcbaGa aeOTdmaaCaaaleqabaGaaGOmaaaaaaGccaqGGaGaeyOeI0Iaaeiia8 qacaqGJbGaaeyEaiaabYgacaqGPbGaaeOBaiaabsgacaqGLbGaaeOC a8aacaqGSaaabaGaae4Sdiabg2da9maalaaabaGaaeOWdmaabmaaba GaamOCamaaDaaaleaadaWgaaadbaWdbiaadoeacaWGpbGaamytaiaa dcfacaWGubaapaqabaaaleaacaWGIqaaaOGaamyqaiabgkHiTiaadw eaaiaawIcacaGLPaaacqGHXcqSdaGcaaqaaiaadseadaWgaaWcbaWd biaadogacaWGVbGaamOBaiaadwgaa8aabeaaaeqaaaGcbaGaaeOWdm aaCaaaleqabaGaaGOmaaaakiaadgeaaaGaaeiiaiabgkHiTiaabcca peGaae4yaiaab+gacaqGUbGaaeyza8aacaqGUaaaaiaawUhaaaaa@934E@  (5)

The following convolution operations are used in system (5):

D sph = r 0 r COMPT Т ξ 2 ξ 2 R 2   sphere, D cyl = r COMPT Т ζ 2 ζ 2 Р 2   cylinder, D cone = ς r COMPT Т AE 2 ς 2 A Q 2   cone. MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaadseadaWgaaWcbaaeaaaaaaaaa8qacaWGZbGaamiCaiaadIga a8aabeaakiabg2da9maadmaabaWaaeWaaeaacaWGYbWaaSbaaSqaai aaicdaaeqaaOGaeyOeI0IaamOCamaaBaaaleaapeGaam4qaiaad+ea caWGnbGaamiuaiaadsfaa8aabeaaaOGaayjkaiaawMcaamaaCaaale qabaGaamOieaaakiaab67aaiaawUfacaGLDbaadaahaaWcbeqaaiaa ikdaaaGccqGHsislcaqG+oWaaWbaaSqabeaacaaIYaaaaOGaamOuam aaCaaaleqabaGaaGOmaaaakiaabccacqGHsislcaqGGaWdbiaaboha caqGWbGaaeiAaiaabwgacaqGYbGaaeyza8aacaqGSaaabaGaamiram aaBaaaleaapeGaam4yaiaadMhacaWGSbaapaqabaGccqGH9aqpdaWa daqaaiaadkhadaqhaaWcbaWaaSbaaWqaa8qacaWGdbGaam4taiaad2 eacaWGqbGaamivaaWdaeqaaaWcbaGaamOieaaakiaabA7aaiaawUfa caGLDbaadaahaaWcbeqaaiaaikdaaaGccqGHsislcaqG2oWaaWbaaS qabeaacaaIYaaaaOGaamiiemaaCaaaleqabaGaaGOmaaaakiaabcca cqGHsislcaqGGaWdbiaabogacaqG5bGaaeiBaiaabMgacaqGUbGaae izaiaabwgacaqGYbWdaiaabYcaaeaacaWGebWaaSbaaSqaa8qacaWG JbGaam4Baiaad6gacaWGLbaapaqabaGccqGH9aqpdaWadaqaaiabgk HiTiaabk8adaqadaqaaiaadkhadaqhaaWcbaWaaSbaaWqaa8qacaWG dbGaam4taiaad2eacaWGqbGaamivaaWdaeqaaaWcbaGaamOieaaaki aadgeacqGHsislcaWGfbaacaGLOaGaayzkaaaacaGLBbGaayzxaaWa aWbaaSqabeaacaaIYaaaaOGaeyOeI0IaaeOWdmaaCaaaleqabaGaaG OmaaaakiaadgeacaWGrbWaaWbaaSqabeaacaaIYaaaaOGaaeiiaiab gkHiTiaabccapeGaae4yaiaab+gacaqGUbGaaeyza8aacaqGUaaaai aawUhaaaaa@98E7@  

In each case, one of the two possible solutions in system (5) is selected based on the following criterion:

r BW min. MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGccqGHsgIRciGG TbGaaiyAaiaac6gacaGGUaaaaa@3E89@  (6)

From the known equation of the COG of a spacecraft, we can derive a formula for estimating the total mass of the BW:

M BW = r C r 0 r 0 r BW M SC , MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGccqGH9aqpdaWc aaqaaiaadkhadaWgaaWcbaGaam4qaaqabaGccqGHsislcaWGYbWaaS baaSqaaiaaicdaaeqaaaGcbaGaamOCamaaBaaaleaacaaIWaaabeaa kiabgkHiTiaadkhadaWgaaWcbaWdbiaadkeacaWGxbaapaqabaaaaO GaamytamaaBaaaleaapeGaam4uaiaadoeaa8aabeaakiaacYcaaaa@480C@  (7)

where M SC MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaaBa aaleaaqaaaaaaaaaWdbiaadofacaWGdbaapaqabaaaaa@38E6@  is the mass of the spacecraft without BW.

The proposed approach to determining the coordinates of the center of gravity (COG) of the BW installation using Equation (3), with the coefficient values from (5) based on criterion (6) and the total mass of the BW calculated using equation (7), gives only an approximate result. As mentioned earlier, in reality, the surface on which the BW is placed is not smooth when BWs of different configurations are used. Therefore, the COG of the BW assembly will shift in the neighborhood along the centering line from the wall of the spacecraft compartment closer to the specified COG of the spacecraft without the BW. Consequently, the BW placement must be modeled in order to derive the final coordinates of the BW center points from the resulting electronic geometric model of the BW installation and then recalculate the spacecraft’s MCIC. To do this, first of all, one must identify available spaces in the neighborhood of the previously determined center of gravity of the BW installation.

Identification of the Configuration of the Balancing Weight Accommodation Zone

The balancing weight accommodation zone is a collection of spatial regions free from compartment structural elements and other components of the spacecraft. Generally, it can be separable and contain discontinuities.

During automated balancing weight placement, the configuration identification of the accommodation zone must also be automated to ensure the processing and use of model information for further actions. During this procedure, the boundaries of free regions along the surface of the spacecraft compartment body are determined, and the coordinates of points within these boundaries are stored at a specified interval.

The points of the structural elements on the surface of the spacecraft compartment body form a layer. Within this layer, points belonging to the boundaries of the BW accommodation zone are determined. They are connected by straight lines, thereby forming closed exclusion zones. The number of such zones is equal to the number of structural elements in a given area surrounding the BW installation's center of gravity, as the points in question possess attributes of belonging to specific electronic geometric models. By logically subtracting the exclusion zones from the surface of the BW accommodation zone, the available space for the payload installation is determined. Built-in procedures of CAD systems are used to check the space of the BW installation's intersections with its surroundings.

The sequence of steps described above is known as the Jarvis algorithm or method [15], according to which the points of interest are determined based on the following criterion:

cos v i+1 = x i+1 x i x i1 x i + y i+1 y i y i1 y i + z i+1 z i z i1 z i min, MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaci4yaiaac+ gacaGGZbGaaeODamaaBaaaleaacaWGPbGaey4kaSIaaGymaaqabaGc cqGH9aqpdaqadaqaaiaadIhadaWgaaWcbaGaamyAaiabgUcaRiaaig daaeqaaOGaeyOeI0IaamiEamaaBaaaleaacaWGPbaabeaaaOGaayjk aiaawMcaamaabmaabaGaamiEamaaBaaaleaacaWGPbGaeyOeI0IaaG ymaaqabaGccqGHsislcaWG4bWaaSbaaSqaaiaadMgaaeqaaaGccaGL OaGaayzkaaGaey4kaSYaaeWaaeaacaWG5bWaaSbaaSqaaiaadMgacq GHRaWkcaaIXaaabeaakiabgkHiTiaadMhadaWgaaWcbaGaamyAaaqa baaakiaawIcacaGLPaaadaqadaqaaiaadMhadaWgaaWcbaGaamyAai abgkHiTiaaigdaaeqaaOGaeyOeI0IaamyEamaaBaaaleaacaWGPbaa beaaaOGaayjkaiaawMcaaiabgUcaRmaabmaabaGaamOEamaaBaaale aacaWGPbGaey4kaSIaaGymaaqabaGccqGHsislcaWG6bWaaSbaaSqa aiaadMgaaeqaaaGccaGLOaGaayzkaaWaaeWaaeaacaWG6bWaaSbaaS qaaiaadMgacqGHsislcaaIXaaabeaakiabgkHiTiaadQhadaWgaaWc baGaamyAaaqabaaakiaawIcacaGLPaaacqGHsgIRciGGTbGaaiyAai aac6gacaGGSaaaaa@7709@  (8)

where x i1 y i1 z i1 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaeWaaeaafa qabeqadaaabaGaamiEamaaBaaaleaacaWGPbGaeyOeI0IaaGymaaqa baaakeaacaWG5bWaaSbaaSqaaiaadMgacqGHsislcaaIXaaabeaaaO qaaiaadQhadaWgaaWcbaGaamyAaiabgkHiTiaaigdaaeqaaaaaaOGa ayjkaiaawMcaaaaa@430F@  is the coordinates of the second-to-last identified point; x i y i z i MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaeWaaeaafa qabeqadaaabaGaamiEamaaBaaaleaacaWGPbaabeaaaOqaaiaadMha daWgaaWcbaGaamyAaaqabaaakeaacaWG6bWaaSbaaSqaaiaadMgaae qaaaaaaOGaayjkaiaawMcaaaaa@3E17@  is the coordinates of the last identified point; x i+1 y i+1 z i+1 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaeWaaeaafa qabeqadaaabaGaamiEamaaBaaaleaacaWGPbGaey4kaSIaaGymaaqa baaakeaacaWG5bWaaSbaaSqaaiaadMgacqGHRaWkcaaIXaaabeaaaO qaaiaadQhadaWgaaWcbaGaamyAaiabgUcaRiaaigdaaeqaaaaaaOGa ayjkaiaawMcaaaaa@42EE@  is the coordinates of the next point to be identified; v i+1 MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaaeODamaaBa aaleaacaWGPbGaey4kaSIaaGymaaqabaaaaa@39C8@  is the angle between the points in consideration.

The advantage of the Jarvis method, for example, over the Graham fast shell scan, is that it requires fewer calculations and is more convenient to use in three-dimensional space [16].

Coordinates of the placement of balancing weights

The identified BW placement area is delineated with a coordinate grid, the spacing of which is equal to the dimension l BWavg MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbGaamyyaiaadAhacaWGNbaa paqabaaaaa@3BD5@  of the main BWs from the nomenclature. The markings should be made on the BW placement surface.

To minimize the overall BW mass by determining BW installation locations using the coordinate grid, the BW placement problem can be represented as a linear programming problem, where the objective function is

M BW = j=1 N BW m BW j min, MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGccqGH9aqpdaae Wbqaaiaad2gadaWgaaWcbaWdbiaadkeacaWGxbaapaqabaGcdaWgaa WcbaGaamOAaaqabaaabaGaamOAaiabg2da9iaaigdaaeaacaWGobWa aSbaaWqaa8qacaWGcbGaam4vaaWdaeqaaaqdcqGHris5aOGaeyOKH4 QaciyBaiaacMgacaGGUbGaaiilaaaa@4B1F@  (9)

where j= 1, N BW ¯ MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOAaiabg2 da9maanaaabaGaaGymaiaacYcacaWGobWaaSbaaSqaaabaaaaaaaaa peGaamOqaiaadEfaa8aabeaaaaaaaa@3C5B@  is the ordinal number of the BW; N BW MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOtamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaaaaa@38EA@  is the number of all BWs; m BW j MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaGa amOAaaqabaaaaa@3A2E@  is the mass of the BW.

This objective function is subject to mass constraints depending on the available range of BWs used:

m BW j 0, m BW j max m BW j , MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaad2gadaWgaaWcbaaeaaaaaaaaa8qacaWGcbGaam4vaaWdaeqa aOWaaSbaaSqaaiaadQgaaeqaaOGaeyyzImRaaGimaiaacYcaaeaaca WGTbWaaSbaaSqaa8qacaWGcbGaam4vaaWdaeqaaOWaaSbaaSqaaiaa dQgaaeqaaOGaeyizImQaciyBaiaacggacaGG4bGaamyBamaaBaaale aapeGaamOqaiaadEfaa8aabeaakmaaBaaaleaacaWGQbaabeaakiaa cYcaaaGaay5Eaaaaaa@4BDF@  (10)

where max m BW j MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaciyBaiaacg gacaGG4bGaamyBamaaBaaaleaaqaaaaaaaaaWdbiaadkeacaWGxbaa paqabaGcdaWgaaWcbaGaamOAaaqabaaaaa@3D01@  is the maximum possible mass of BW according to the nomenclature.

In addition, the following centering constraints are imposed on the objective function (9):

M SC r 0 min r 0 + j=1 N BW m BW j r BW j min r 0 0, M SC r 0 max r 0 + j=1 N BW m BW j r BW j max r 0 0, MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaiqaaqaabe qaaiaad2eadaWgaaWcbaaeaaaaaaaaa8qacaWGtbGaam4qaaWdaeqa aOWaaeWaaeaacaWGYbWaaSbaaSqaaiaaicdaaeqaaOGaeyOeI0Iaci yBaiaacMgacaGGUbGaamOCamaaBaaaleaacaaIWaaabeaaaOGaayjk aiaawMcaaiabgUcaRmaaqahabaGaamyBamaaBaaaleaapeGaamOqai aadEfaa8aabeaakmaaBaaaleaacaWGQbaabeaaaeaacaWGQbGaeyyp a0JaaGymaaqaaiaad6eadaWgaaadbaWdbiaadkeacaWGxbaapaqaba aaniabggHiLdGcdaqadaqaaiaadkhadaWgaaWcbaWdbiaadkeacaWG xbaapaqabaGcdaWgaaWcbaGaamOAaaqabaGccqGHsislciGGTbGaai yAaiaac6gacaWGYbWaaSbaaSqaaiaaicdaaeqaaaGccaGLOaGaayzk aaGaeyyzImRaaGimaiaacYcaaeaacaWGnbWaaSbaaSqaa8qacaWGtb Gaam4qaaWdaeqaaOWaaeWaaeaacaWGYbWaaSbaaSqaaiaaicdaaeqa aOGaeyOeI0IaciyBaiaacggacaGG4bGaamOCamaaBaaaleaacaaIWa aabeaaaOGaayjkaiaawMcaaiabgUcaRmaaqahabaGaamyBamaaBaaa leaapeGaamOqaiaadEfaa8aabeaakmaaBaaaleaacaWGQbaabeaaae aacaWGQbGaeyypa0JaaGymaaqaaiaad6eadaWgaaadbaWdbiaadkea caWGxbaapaqabaaaniabggHiLdGcdaqadaqaaiaadkhadaWgaaWcba WdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaGaamOAaaqabaGccqGH sislciGGTbGaaiyyaiaacIhacaWGYbWaaSbaaSqaaiaaicdaaeqaaa GccaGLOaGaayzkaaGaeyizImQaaGimaiaacYcaaaGaay5Eaaaaaa@8510@  (11)

where min r 0 , max r 0 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaciyBaiaacM gacaGGUbGaamOCamaaBaaaleaacaaIWaaabeaakiaacYcacaqGGaGa ciyBaiaacggacaGG4bGaamOCamaaBaaaleaacaaIWaaabeaaaaa@40D6@  are the vectors of the boundary admissible values of the spacecraft's center of mass coordinates.

A linear programming problem with an objective function of the form (9) and constraints of the form (10) and (11) can be solved, for example, using the simplex method or the multiparametric Newton method [17]:

m BW j k+1 T = m BW j k T H k 1 M BW k T M BW , MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaamWaaeaaca WGTbWaaSbaaSqaaabaaaaaaaaapeGaamOqaiaadEfaa8aabeaakmaa BaaaleaacaWGQbaabeaaaOGaay5waiaaw2faamaaDaaaleaacaWGRb Gaey4kaSIaaGymaaqaaiaadsfaaaGccqGH9aqpdaWadaqaaiaad2ga daWgaaWcbaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaGaamOAaa qabaaakiaawUfacaGLDbaadaqhaaWcbaGaam4AaaqaaiaadsfaaaGc cqGHsislcaWGibWaa0baaSqaaiaadUgaaeaacqGHsislcaaIXaaaaO WaaeWaaeaacaWGnbWaaSbaaSqaa8qacaWGcbGaam4vaaWdaeqaaaGc caGLOaGaayzkaaGaey4bIe9aa0baaSqaaiaadUgaaeaacaWGubaaaO WaaeWaaeaacaWGnbWaaSbaaSqaa8qacaWGcbGaam4vaaWdaeqaaaGc caGLOaGaayzkaaGaaiilaaaa@5A32@  (12)

where k MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4Aaaaa@3709@  is the ordinal number of the computational iteration; H k 1 M BW MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamisamaaDa aaleaacaWGRbaabaGaeyOeI0IaaGymaaaakmaabmaabaGaamytamaa BaaaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaaakiaawIcaca GLPaaaaaa@3E17@  is the inverse Hessian matrix of the objective function (9); k T M BW MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaey4bIe9aa0 baaSqaaiaadUgaaeaacaWGubaaaOWaaeWaaeaacaWGnbWaaSbaaSqa aabaaaaaaaaapeGaamOqaiaadEfaa8aabeaaaOGaayjkaiaawMcaaa aa@3E01@  is the Nabla operator of the objective function (9).

For technical reasons, the obtained BW masses may be rounded to integer values. The grid point numbers at which the BW mass is nonzero are used to construct an electronic geometric model of the BW installation. It should be noted that BW units of different masses but identical overall dimensions l BW j MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaGa amOAaaqabaaaaa@3A2D@  will have different thicknesses h BW j MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamiAamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaGa amOAaaqabaaaaa@3A29@ . Therefore, if geometric intersections with structural elements occur, some BW units will need to be repositioned.

After that, the coordinates of the BW unit’s COG are recalculated to account for the fasteners using the formula

r BW = j m BW j r BW j + m BW fast r BW fast j m BW j + m BW fast , MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGccqGH9aqpdaWc aaqaamaaqafabaGaamyBamaaBaaaleaapeGaamOqaiaadEfaa8aabe aakmaaBaaaleaacaWGQbaabeaakiaadkhadaWgaaWcbaWdbiaadkea caWGxbaapaqabaGcdaWgaaWcbaGaamOAaaqabaaabaGaamOAaaqab0 GaeyyeIuoakiabgUcaRiaad2gadaWgaaWcbaWdbiaadkeacaWGxbaa paqabaGcdaWgaaWcbaWdbiaadAgacaWGHbGaam4Caiaadshaa8aabe aakiaadkhadaWgaaWcbaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWc baWdbiaadAgacaWGHbGaam4Caiaadshaa8aabeaaaOqaamaaqafaba GaamyBamaaBaaaleaapeGaamOqaiaadEfaa8aabeaakmaaBaaaleaa caWGQbaabeaaaeaacaWGQbaabeqdcqGHris5aOGaey4kaSIaamyBam aaBaaaleaapeGaamOqaiaadEfaa8aabeaakmaaBaaaleaapeGaamOz aiaadggacaWGZbGaamiDaaWdaeqaaaaakiaacYcaaaa@63E9@  (13)

where m BW fast MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyBamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaWd biaadAgacaWGHbGaam4Caiaadshaa8aabeaaaaa@3D20@  is the mass of the fastener; r BW fast = x BW fast y BW fast z BW fast T MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFy0Jg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaaqaaaaaaaaaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaWd biaadAgacaWGHbGaam4Caiaadshaa8aabeaakiabg2da9maabmaaba qbaeqabeWaaaqaaiaadIhadaWgaaWcbaWdbiaadkeacaWGxbaapaqa baGcdaWgaaWcbaWdbiaadAgacaWGHbGaam4Caiaadshaa8aabeaaaO qaaiaadMhadaWgaaWcbaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWc baWdbiaadAgacaWGHbGaam4Caiaadshaa8aabeaaaOqaaiaadQhada WgaaWcbaWdbiaadkeacaWGxbaapaqabaGcdaWgaaWcbaWdbiaadAga caWGHbGaam4Caiaadshaa8aabeaaaaaakiaawIcacaGLPaaadaahaa Wcbeqaaiaadsfaaaaaaa@55FA@  is the coordinate vector of the fastener's COG.

When using adhesive and other permanent types of joints, the calculation using formula (13) is performed under the assumption that the layer of fastening material is distributed uniformly across the surface of the spacecraft compartment wall. The coordinates of the fastener’s center of mass are determined during simulation.

A general flowchart of the algorithm for placing BWs using formulas (1) MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@ (13) is shown in Fig. 2. The list of required input data and their possible sources is presented in Table 2.

 

Рис. 2. Общая блок-схема алгоритма определения мест установки БГ

Fig. 2. Control flow chart for BW mounting places determination

 

Table 2

List of input data and their sources for the algorithm

Data Module

Parameter

Sources

Collection options

Processing methods

Coordinates of the spacecraft’s COG at a specified position

r0

Explanatory note on spacecraft

General types of spacecrafts

Retrieving or entering information from DD

Entering

Retrieving

Coordinates of the spacecraft’s COG at the current position, excluding the BW

rC

Calculation of MCIC of spacecraft

Model of the spacecraft compartment

Balancing rig

Retrieving or entering information from DD

Retrieving or entering from the model

Retrieving or entering information from the rig

Retrieving Calculation

Coordinates of the geometric center of the compartment hull

rCOMPT

General types of spacecrafts

Model of the spacecraft compartment

Retrieving or entering information from DD

Retrieving or entering from the model

Entering

Retrieving

Inner radius of the spacecraft compartment hull

RCOMPT

General types of spacecrafts

Model of the spacecraft compartment

Retrieving or entering information from DD

Retrieving or entering from the model

Entering

Retrieving

Mass of the spacecraft, excluding the BW

MSC

Explanatory note on spacecraft

General types of spacecrafts

Calculation of MCIC of spacecraft

Model of the spacecraft compartment

Balancing rig

Retrieving or entering information from DD

Retrieving or entering from the model

Retrieving or entering information from the rig

Retrieving Calculation

Average thickness of the BW from the nomenclature

hBWavg

DD on BW

Model of the spacecraft compartment

Retrieving or entering information from DD

Retrieving or entering from the model

Entering

Retrieving

Average dimensions of the BW from the nomenclature

lBWavg

DD on BW

Model of the spacecraft compartment

Retrieving or entering information from DD

Retrieving or entering from the model

Entering

Retrieving

Mass Constraints

maxmBWj

Explanatory note on spacecraft

Retrieving or entering information from DD

Entering

Retrieving

Centering Constraints

minr0maxr0

Explanatory note on spacecraft

Retrieving or entering information from DD

Entering

Retrieving

 

The method is implemented according to the algorithm presented in Fig. 2, in the following order of actions:

  1. Determine (specify) the surface type of the spacecraft compartment body according to equation (2).
  2. Calculate the coordinates of the COG of the BW installation using formulas (3) MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@ (6).
  3. Calculate the total mass of the BW using formula (7).
  4. Construct the boundaries of the BW placement zone and the exclusion zones according to criterion (8) using the Jarvis method.
  5. Construct a coordinate grid on the surface of the BW placement zone with a given step equal to the average BW size from the nomenclature.
  6. Record the coordinates of the grid nodes as possible coordinates for the placement of the BW.
  7. Set the maximum masses of the BW at all nodes of the coordinate grid.
  8. Calculate the optimal masses of the payload at the coordinate grid nodes by solving the linear programming problem (9) MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@ (11) using Newton's method (12) or the simplex method.
  9. Taking into account the results of optimization and the mass of the fasteners, calculate the coordinates of the COG of the BW installation using formula (13).
  10. Calculate the MCIC of the spacecraft layout.
  11. Check that the conditions for the permissible deviations of the COG and the moments of inertia are met.
  12. If the conditions in point 11 are not met, repeat points 6 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@ 11 of this procedure, taking into account the current optimization results.
  13. If the conditions in point 11 are met, take the balancing result into account in subsequent calculations of the MCIC of the spacecraft compartment.

Test Results

A conical spacecraft compartment with a spherical bottom was used as an example for testing. Calculations using formulas (2) MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbaqcLbwaqa aaaaaaaaWdbiaa=nbiaaa@37A4@ (6) revealed that the BW placement area should be located near the junction of the bottom and the hull shell. The configuration of the placement area was identified, taking into account the surrounding structural elements. The coordinates for the placement of the BWs were determined and their masses were selected. 9 of the 40 applied BWs needed to be moved to eliminate intersections with the structure. After this, the determination of the BW installation locations was completed, the coordinates of the COG of the BW installation center were recalculated, the MCIC of the spacecraft with the BW were also recalculated, and the fulfillment of the typical mass-centering and inertial requirements was ensured. To improve the assembly efficiency, 24 adjacent BWs were combined into monolithic massive BWs with complex geometry (topology). The simplified modeling results are shown in Fig. 3.

 

Рис. 3. Результаты определения мест установки БГ:

1 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  внутренние стенки корпуса; 2 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  сечение корпуса; 3 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  конструкция; 4 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  БГ тип 1 (массивные); 5 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  БГ тип 2; 6 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  БГ тип 3; 7 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  БГ тип 4

Fig. 3. Results of BW mounting places determination:

1 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  hull surface inside; 2 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  hull sectional view; 3 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  construction; 4 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  BW type 1 (massive); 5 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  BW type 2; 6 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  BW type 3; 7 MathType@MTEF@5@5@+= feaahGart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaacbeqcLbtaqa aaaaaaaaWdbiaa=nbiaaa@3775@  BW type 4

 

Conclusion

As a result of testing, the methodology for determining BW installation locations within the spacecraft compartment was refined. Electronic dimensional models of the simulated compartment and its structural elements were used for modeling. The degree of detail of the models will presumably impact the speed of the algorithm and the quality of the solution. Intersections with structural elements may occur when their detailed models are replaced by dimensional models, as well as if the surrounding elements do not belong to the payload placement surface. It is also necessary to shift BWs whose thickness does not equal the average value from the nomenclature.

Based on the placement results, the user can, at his own discretion, combine adjacent BWs into more massive BWs of any shape in order to reduce the labor intensity of their manufacture and installation. These actions may lead to deviations in the mass of the BW, therefore, after completing all transformations, it is recommended to recalculate the MCIC of the spacecraft. If the enterprise does not have a nomenclature of standardized BWs, then the most suitable dimensions of square BWs should be specified as initial data, guided by analogies from statistics. When modeling, it is necessary to take into account the variation in the mass of the BWs, so their models must be parameterized by thickness.

To determine the COG of the BW installation, identify the configuration of the BW placement zone, calculate the BW placement coordinates, and verify the requirements for the spacecraft's MCIC, one has to perform many mathematical operations. This method is advisable to use only when a specialized calculation software package is available, including modules for static and dynamic balancing.

Furthermore, it must be taken into account that the model of the BW installation obtained by the method must be supplemented with models of fasteners and other components. For issuing the DD, it will still be necessary to go through approvals regarding the surrounding structure and strength. Including auxiliary sections in the method to automate these stages makes no sense, because these simulation results are performed using special software products. Nevertheless, they are very useful when it is necessary to work out and discuss several options for BW placement before starting to prepare the DD.

Thus, as a result of the work performed, a method for determining the installation locations of BWs inside spacecraft compartments was developed and tested. Using the proposed algorithm, it is possible to ensure the minimum total mass of the BWs used. Simulation showed that in order to meet the requirements for the spacecraft's MCIC regarding allowable values of moments of inertia, dynamic balancing should be performed. In addition, some comments were identified regarding the particular procedures, the introduced assumptions, and the order of using this method. Nevertheless, the proposed approach allows reducing the labor intensity of developing the DD for BW installation by accelerated elaboration of various BW placement options, assessment of possible intersections with structural elements of the compartment and other surrounding equipment at the stages of preliminary design and working DD. Also, the possibility of organizing the process of transferring initial data from the balancing stand software is not excluded.

×

About the authors

Andrey А. Belyakov

S. P. Korolev Rocket and Space Corporation “Energia”; Samara National Research University

Author for correspondence.
Email: post@rsce.ru
ORCID iD: 0000-0002-5789-8048

Postgraduate, Samara National Research University; Design Engineer of II Category, S. P. Korolev Rocket and Space Corporation “Energia”

Russian Federation, 4a, Lenin street, Moscow region, Korolev, 141070; 34, Moskovskoe shosse, Samara, 443086

Alexander I. Shulepov

Samara National Research University

Email: shulepov-al@mail.ru

Cand. Sc., Distinguished Designer of Russian Federation, Assistant Professor of Department of Space Engineering of Institute of Aircraft and Spacecraft

Russian Federation, 34, Moskovskoe shosse, Samara, 443086

Vladimir М. Papazov

S. P. Korolev Rocket and Space Corporation “Energia”

Email: post@rsce.ru

Cand. Sc., Leading Researcher

Russian Federation, 4a, Lenin street, Moscow region, Korolev, 141070

References

  1. Kamalov V. S. Proizvodstvo kosmicheskih apparatov [Spacecraft Production]. Moscow, Mashinostroenie Publ., 1982, 280 p.
  2. Belyakov I. Т., Zernov I. А., Antonov Е. F. et al. Tehnologija sborki i ispytanij kosmicheskih apparatov [Spacecraft Assembly and Testing Technology]. Moscow, Mashinostroenie Publ., 1990, 352 p.
  3. Goryachev А. S. Balansirovka agregatov letatel'nyh apparatov [Balancing Aircraft Components]. Kuibyshev, KuAI Publ., 1982, 70 p.
  4. Kamaldinov А. M., Aksyonov Е. G., Ermolov М. I. [Using a Gyroplatform Static Balancing Stand to Determine the Position of the Center of Mass of “Blits-М” Spacecraft]. Vestnik NPO “Tehnomash”. 2017, No. 1, P. 38 (In Russ.).
  5. Videnkin N. А. [Automated Stands for Measuring Control of Mass Geometry Parameters]. Vestnik NPO “Tehnomash”. 2017, No. 1, P. 41–42 (In Russ.).
  6. Kochkin Е. V., Medar А. V., Matveev Е. V., Medenkov V. I., Makarov I. I. [Stand for Static Balancing and Determination of Product Mass]. Vestnik NPO “Tehnomash”. 2017, No. 1, P. 44 (In Russ.).
  7. Videnkin N. А., Kochkin Е. V., Matveev Е. V. [Method for Determining the Components of the Spacecraft Inertia Tensor on Combined Automated Stand]. Vestnik NPO “Tehnomash”. 2017, No. 3, P. 44–46 (In Russ.).
  8. Kuznetsov V. S., Kamaldinov А. М. [Development of a Combined Technology for Adjusting Masses When Balancing the Engine Rotor-Flywheel]. Vestnik NPO “Tehnomash”. 2018, No. 6, P. 38–40 (In Russ.).
  9. Shepetkov N. А., Makarov I. I., Kochkin Е. V. [Modern Technical Means for Monitoring the Characteristics of the Mass Geometry of Aircraft]. Vestnik NPO “Tehnomash”. 2018, No. 6, P. 81–82 (In Russ.).
  10. Klyuchnikov А. V. [Developing and Refining the Algorithm o Single-Plane Dynamic Balancing of a High-Speed Flying Prototype]. Vestnik SibGAU. 2015, Vol. 16, No. 2, P. 411–416.
  11. Klyuchnikov А. V. [Numerical Algorithm for Optimization of a Process of Conical Flying Model Equilibration on a Vertical Balancing Stand]. Vestnik SibGAU. 2016, Vol. 17, No. 2, P. 309–317.
  12. Stoyan Yu. G., Kulish Е. N. Avtomatizacija proektirovanija komponovki oborudovanija letatel'nyh apparatov [Automatizatoin of Aircraft Equipment Layout Engineering]. Moscow, Mashinostroenie Publ., 1984, 192 p.
  13. Gavrilov V. N. Avtomatizirovannaja komponovka pribornyh otsekov letatel'nyh apparatov [Automated Layout of Aircraft Instrument Compartments]. Moscow, Mashinostroenie Publ., 1988, 137 p.
  14. Shulepov А. I., Lu Tszya [About a Problem of Cargo Arrangement in Landing Module]. Sb. Trudov XVII Vserossiyskogo seminara po upravleniyu dvizheniem i navigatsii letatelnyh apparatov. Chast’ 1. [Proc. XVII Russ. Seminar of flight control and navigation of aircraft. Part 1]. Samara, 2015, No. 4, P. 187–190 (In Russ.).
  15. Thomas H. Cormen, Charles E. Leiserson, Ronald L. Rivest, Clifford Stein. Algoritmy: postroenie i analiz [Introduction to Algorithms]. Moscow, Williams Publ., 2013, 1328 p.
  16. Belyakov А. А., Shulepov А. I. [Topological Methods for Determining the Regular Zone of Placement of Balancing Weights in a Spacecraft Compartment]. Vestnik AmurSU. 2022, No. 97, P. 59–65 (In Russ.).
  17. Belyakov А. А. [An Automated Solution of the Problem of Arranging Balancing Weights in Spacecraft Compartments]. Kosmicheskaja tehnika i tehnologii. 2024, No. 46, P. 14–27 (In Russ.).

Supplementary files

Supplementary Files
Action
1. JATS XML
2. Fig. 1. Computational model for balancing weights (BW) center of gravity (COG) determination

Download (153KB)
3. Fig. 2. Control flow chart for BW mounting places determination

Download (241KB)
4. Fig. 3. Results of BW mounting places determination: 1 – hull surface inside; 2 – hull sectional view; 3 – construction; 4 – BW type 1 (massive); 5 – BW type 2; 6 – BW type 3; 7 – BW type 4

Download (251KB)

Copyright (c) 2026 Belyakov A.А., Shulepov A.I., Papazov V.М.

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.