2017
Математика и механика
№ 46
МАТЕМАТИКА
УДК 519.711.3, 514.181.22 Б01 10.17223/19988621/46/1
М.С. Бухтяк
СОСТАВНАЯ ПОВЕРХНОСТЬ, БЛИЗКАЯ К ПСЕВДОМИНИМАЛЬНОЙ
Данная публикация продолжает серию работ автора о моделировании деформированного лепестка осесимметричного параболического рефлектора. Схема моделирования формы деформированного лепестка сетеполотна осесимметричного рефлектора, предложенная В.М. Гряником и В.И. Ломаном, разработана для изотропного упругого материала, прикрепляемого к параболическим жестким ребрам, и не видно способов адаптировать её для орто-тропного сетеполотна с иными способами закрепления. Автором в публикациях 2016 года предложена методика моделирования формы ортотропного упругого материала, основанная на использовании поверхности, для которой отношение главных кривизн есть величина постоянная (выражающаяся через отношение коэффициентов растяжения материала в двух ортогональных направлениях. Такая поверхность названа псевдоминимальной. В указанных публикациях методика адаптирована под конкретную ситуацию, описанную Гряником и Ломаном (лепесток сетеполотна осесимметричного рефлектора). Автором решены вопросы и более общего характера. Доказана теорема существования (широта класса псевдоминимальных поверхностей оказалась такая же, как и для минимальных поверхностей - две функции одного аргумента, то есть - в принципе - псевдоминимальная поверхность указанного веса определяется граничной линией). Данные о широте класса допускают иное истолкование: возможность (теоретическая!) построить составную псевдоминимальную поверхность, присоединяя к одной плоской линии семейство других плоских линий. Пример построения такой поверхности (она задается вектор-функцией) строится в данной статье. Для этой же поверхности строится явное задание. Это открывает возможности для моделирования поверхности ортотропного материала при выборе иных способов закрепления.
Ключевые слова: гладкая поверхность, главные кривизны, ортотропия, псевдосредняя кривизна, задание поверхности посредством семейства образующих кривых, явное задание поверхности.
1. Постановка задачи
Автор продолжает исследование о моделировании деформированного лепестка осесимметричного параболического рефлектора [1-6]. Схема моделирования формы деформированного лепестка сетеполотна осесимметричного рефлектора, предложенная в [7], разработана для изотропного упругого материала, и не видно способов адаптировать её для ортотропного сетеполотна. Автором в [1, 6] предложена методика моделирования формы ортотропного упругого материала, осно-
ванная на использовании поверхности, для которой отношение главных кривизн есть величина постоянная (выражающаяся через отношение коэффициентов растяжения материала в двух ортогональных направлениях. Такая поверхность названа псевдоминимальной. В [1] методика адаптирована под конкретную ситуацию (лепесток сетеполотна осесимметричного рефлектора). Работа [6] посвящена вопросам более общего характера. Доказана теорема существования (широта класса псевдоминимальных поверхностей оказалась такая же, как и для минимальных поверхностей - две функции одного аргумента, то есть - в принципе -псевдоминимальная поверхность указанного веса [6] определяется граничной линией). Данные о широте класса допускают иное истолкование: возможность (теоретическая!) построить составную псевдоминимальную поверхность, присоединяя к одной плоской линии семейство других плоских линий. Пример построения такой поверхности (она задается вектор-функцией) строится в данной статье. Для этой же поверхности строится явное задание.
В работе [1] применен способ моделирования ортотропных свойств лепестка сетеполотна осесимметричного рефлектора, основанный на зависимости отношения главных кривизн поверхности сетеполотна от коэффициентов растяжения в указанных направлениях.
На этом пути автор пришел к понятию псевдосредней кривизны и псевдоминимальной поверхности [6]. А именно, если главные кривизны поверхности
Z : r = r (u, v) e C2
упорядочены номерами (k1 и k2), то псевдосредней кривизной веса L называется величина
HL = kj + Lk2, (L = const, L Ф 0), а поверхность, для которой
Hl = 0, (1)
называется псевдоминимальной поверхностью веса L. Ясно, что при L = 1 получаем минимальную поверхность.
Если отсутствует правило упорядочения главных кривизн, то, имея в виду условность нумерации, говорим об ослабленном условии псевдоминимальности в виде
(k1 + Lk2)(Lk1 + k2) = 0 ,
которое, очевидно, от перемены номеров не меняется.
Пусть gn, g12, g22 - коэффициенты первой квадратичной формы поверхности Z , а bu ,b12 ,b22 - коэффициенты второй квадратичной формы. Тогда ослабленное условие псевдоминимальности веса а запишется [4] в виде
L(2b12 g12 - g11b22 - b11 g 22 )2 + (1 - L)2 (b11b22 - ^2X^22 - gfo = 0 . (2)
Соответственно безразмерный показатель псевдоминимальности с весом L
5(u v )= (2b12 g12 - g11b22 - b11g22)2 (1 - L )2 (b122 - b11b22 )(g11 g22 - g122) L
2. Построение поверхности
В [1, 2] составлено дифференциальное уравнение, решение которого позволяет построить псевдоминимальную поверхность, проходящую через сечение параболоида вращения двумя плоскостями, проходящими через ось его симметрии. Уравнение в общем случае решается приближенно. Указанная поверхность может быть отнесена к составным поверхностям: она образована однопараметрическим семейством парабол, «подклеенных» своими вершинами к определенной линии О («гребневая линия») - она и определяется упомянутым дифференциальным уравнением.
Возникает вопрос о моделировании псевдоминимальной поверхности (хотя бы приближенно), при этом можно обойтись более скромными средствами: обычно это упрощает исследование. Именно, пусть «гребневая линия» О - заданная парабола и к её текущей точке «подклеивается» парабола, зависящая от точки прикрепления так, что условие (1) выполняется вдоль линии О (рис. 1). Вопрос в том, насколько точно при этом выполняется хотя бы ослабленное условие (1.6) уже вдоль поверхности.
Гребневая линия
Рис. 1. Схема построения составной поверхности Fig. 1. Scheme of constructing the combined surface
Схема построения пояснена на рис. 1. Заштрихована часть нормальной плоскости гребневой линии. Векторы ei составляют ортогональный репер. Пунктирная линия - «подклеенная парабола».
Гребневая линия пусть задана вектор-функцией
r = {t,0, at2}, a = const, a > 0.
Кривизна равна
k = 2a
g = V1 + 4a212 '
Базисные векторы
e = r'(t), e = {0,1,0}, ¿2 =[e3,e] Радиус-вектор текущей точки присоединенной параболы P (t, u) = r (t)- h(t)u2E2 (t) + uE3 (t) .
(3)
(4)
Таким образом,
р(, и) = - 2аи2к(/), и, ^2 + и2к(/)}. Условие (2), примененное вдоль гребневой линии, дает нам уравнение Ь2 ка + ка + 4 ка3Г2 + Ьк2 + Ьа2 + 96 ¿А2 а4/4 + 256Ьк2а6Г6 + +16 Ьк2а2Г2 + 256Ьк2 а8/8 + 4 Ь2 ка3/2 + 8Ька3/2 + 32 Ька5Г4 = 0. Уравнение второй степени относительно к(/). Дискриминант равен а2 (1 + 4а2/2 )2 (Ь -1)2 (Ь2 +16Ьа2/2 + 2Ь +1) .
Совпадение корней равносильно тому, что вдоль гребневой линии сумма главных кривизн равна нулю.
В качестве первого корня уравнения (4) примем
h(t)
(-8La2t2 -1 - L2 +416La2t2 +16lla2t2 +1 - 2L2 + L4 - 32a2L2t2 ) a 2 L (1 + 8a2t2 + 16a 4t4 )( + 4 a 2t2 )
Тогда при а = 1 для различных значений Ь получаем поверхности (3), изображенные на рис. 2.
z
0.4 0.2
-0.2
-0.4
Рис. 2. Поверхности, приближенно псевдоминимальные при L = 0.4, 0.6, 0.8, 1.0 (очередность сверху вниз) Fig. 2. Surfaces approximately pseudominimal at L = 0.4, 0.6, 0.8, and 1.0 (from top to bottom)
0
Вектор-функция (3) слишком сложна для дальнейшего анализа, и для неё строится дробно-рациональное приближение
р , У, ^ } . (5)
Здесь
t
Xs =-
12a2 L2t2 + 12La 2t2 + 48L2 a4t4 + 48 La4t4 + 64 L2a6t6 + + 64 La6t6 + L2 + L + 16a4u 2 Lt2 + 2a2u2 L + 2a2u 2
Zs =-
(L + 1)L (1 + 4a 2t2 ) Y = u,
Г 64L2t8 a6 + 64Lt8 a6 + 48L2t6 a4 + 48Lt6 a4 +12 L214 a 2 + ' +12Lt4a2 + L2t2 + Lt2 - 8u2La2t2 - u2L - u2
(6)
(L + 1)L (1 + 4a 2t2 )
Относительная погрешность построенного приближения оценивается выражением
s =Р - Psl
Значения погрешности при a = 1 и L e {0.6, 0.7, 0.8,0.9,1.0} приведены на рис. 3.
S -
2.0-10-3: 1.5-10-3: 1.0-10-3: 0.5-10-30:
Рис. 3. Относительная погрешность s для L = 0.6,0.7,0.8,0.9,1.0 (очередность
снизу вверх). Для L e {1.0,1.1,1.2,1.3,1.4} результат аналогичный
Fig. 3. Relative error е for L = 0.6, 0.7, 0.8, 0.9, and 1.0 (from bottom to top). For L = {1.0, 1.1, 1.2, 1.3, 1.4}, the result is similar.
Безразмерный индекс псевдоминимальности для а = 1 и Ь е {1.00,1.05,1.10, 1.15,1.20} представлен на рис. 4. Этот же индекс для а = 1 и Ь е {0.6,0.7,0.8, 0.9,1.0} приведен на рис. 5. Использование второго корня уравнения (4) в качестве функции к(/) принципиально картину не меняет.
8(u,t) 0
0.510-3
1.010-3
1.510-3
2.010-3
2.5-10-3 -0
2
Рис. 4. Безразмерный индекс псевдоминимальности для L = 1.00,1.05,1.10,1.15,1.20 (очередность сверху вниз) Fig. 4. Dimensionless index of pseudominimality for L = 1.00, 1.05, 1.10, 1.15, and 1.20 (from top to bottom)
8(u,t) 0
-0.01 -0.02 -0.03
2
1.5
Рис. 5. Безразмерный индекс псевдоминимальности для L = 0.80, 0.85,0.90, 0.95,1.00 (очередность снизу вверх)
Fig. 5. Dimensionless index of pseudominimality for L = 0.80, 0.85, 0.90, 0.95, and 1.00 (from bottom to top)
t
u
u
t
3. Задание поверхности явным уравннием
В некоторых случаях задание поверхности как годографа вектор-функции не столь удобно как использование уравнения вида
^ = / (X, у). (7)
Наличие уравнений (6) позволяет решить этот вопрос локально и приближенно. В окрестности начала координат приближенное разрешение первых двух уравнений приводит к соотношению
Л2ху2 а4 (Ь + 3)х3у2 ха4у4 х3а6 (Ь + 3)у4 г и Т = х - 2-— + 8—у' + 4-4---64- v
ь Ь + ь Ь Ь (ь+1)
Уравнение вида (7) получаем, внеся в третье из уравнений (6) Т вместо г и у вместо и . Результат подстановки не приводим ввиду крайней громоздкости. Для
а = 1, Ь = 6 приводим изображение поверхности (5) и поверхности (7).
Рис. 6. Изображение поверхности (5) (крупная клетка) и поверхности (7) (мелкая клетка) Fig. 6. Image of surface (5) (large squares) and surface (7) (small squares)
4. Заключение
Автор полагает, что предложенная работа (наряду с предшествующими) служит прояснению понятия псевдоминимальной поверхности, способов её аналитического построения (приближенно) и использования для моделирования формы поверхности упругого ортотропного материала.
ЛИТЕРАТУРА
1. БухтякМ.С. Геометрическое моделирование деформации сетеполотна параболического рефлектора // Математическое моделирование. 2016. Т. 39. № 1. С. 97-106.
2. Бухтяк М.С. Дефект отображения для деформированного лепестка сетеполотна // Вестник Томского государственного университета. Математика и механика. 2016. № 2(40). С. 5-17.
3. БухтякМ.С., Соломина А.В. Об одном инварианте пары поверхностей применительно к раскрою сетеполотна // Вестник Томского государственного университета. Математика и механика. 2016. № 1(39). С. 13-24.
4. Бухтяк М.С., Соломина А.В. Геометрическое моделирование раскроя сетеполотна для осесимметричного рефлектора. Часть 1 // Вестник Томского государственного университета. Математика и механика. 2015. № 2(34). С. 5-17.
5. Бухтяк М.С., Соломина А.В. Геометрическое моделирование раскроя сетеполотна для осесимметричного рефлектора. Часть 2 // Вестник Томского государственного университета. Математика и механика. 2015. № 4(29). С. 5-14.
6. Бухтяк М.С. Обобщение минимальных поверхностей и моделирование формы конструкции из ортотропного материала // Вестник Томского государственного университета. Математика и механика. 2016. № 6(44). С. 5-17.
7. Гряник М.В., Ломан В.И. Развертываемые зеркальные антенны зонтичного типа. М.: Радио и связь, 1987. 72 с.
Статья поступила 21.11.2016 г.
Bukhtyak M.S. (2017) A COMPOSITE SURFACE CLOSE TO PSEUDO-MINIMAL. Tomsk State University Journal of Mathematics and Mechanics. 46. pp. 5-13
DOI 10.17223/19988621/46/1
This publication continues the series of the author's papers on simulation of a deformed petal of an axisymmetric parabolic reflector. V.M. Gryanik and V.I. Loman proposed a simulation scheme for the deformed petal shape of the metallic mesh axisymmetric parabolic reflector. The scheme was developed for an isotropic elastic material attached to parabolic hard edges and there seem to be no method to adapt it for an orthotopic metallic mesh with other methods of fixing. The author's papers of 2016 propose a shape simulation method for an orthotopic elastic material. The method is based on using a surface for which the ratio of principal curvatures is constant (expressed by the ratio of the stretching factors of the material in two orthogonal directions). Such a surface is called pseudo-minimal. The method is adapted in these papers to a specific situation described by Gryanik and Loman (metallic mesh petal of an axisymmetric reflector). The author has also solved more general problems. The existence theorem has been proven (the latitude of the class of pseudo-minimal surfaces is the same as for minimal surfaces, namely, two functions of one argument; basically, a pseudo-minimal surface with a specified weight is determined by the boundary line). The data about the latitude class data allow other interpretation: a composite pseudo-minimal surface can be (theoretically!) constructed by attaching a family of flat lines to one flat line. An example of such a surface (given by a vector function) is constructed in this paper. The same surface is constructed in an explicit form. This opens possibilities for simulating an orthotropic material surface when choosing different methods of fixing.
Keywords: smooth surface, principal curvatures, orthotropy, pseudomean curvature, plotting a surface by a family of generating curves, explicit surface plotting.
BUKHTYAK Mikhail Stepanovych (Candidate of Physics and Mathematics, Tomsk State University, Tomsk, Russian Federation)
REFERENCES
1. Bukhtyak M.S. (2016) Geometrical modeling of a metallic mesh deformation of the parabolic reflektor. Mathematical Models and Computer Simulations. 8(4). pp. 453-461. DOI: 10.1134/S2070048216040050.
2. Bukhtyak M.S. (2016) Defekt otobrazheniya dlya deformirovannogo lepestka setepolotna [Defect of mapping for deformed segment of metallic mesh]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 2(40). pp. 5-17. DOI: 10.17223/19988621/40/1.
3. Bukhtyak M.S., Solomina A.V. (2016) Ob odnom invariante otobrazheniya poverkhnostey primenitel'no k raskroyu setepolotna [On an invariant of surface mapping as applied to metallic mesh tailoring]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i
mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 1(39). pp. 13-24. DOI: 10.17223/19988621/39/2.
4. Bukhtyak M.S., Solomina A.V. (2015) Geometricheskoe modelirovanie raskroya setepolotna dlya osesimmetrichnogo reflektora. Chast' 1 [Geometric modeling of metallic mesh sheet tailoring for an axissymmetric reflector. Part 1]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 2(34). pp. 5-17. DOI: 10.17223/19988621/34/1.
5. Bukhtyak M.S., Solomina A.V. (2015) Geometricheskoe modelirovanie raskroya setepolotna dlya osesimmetrichnogo reflektora. Chast' 2 [Geometric modelling of metallic mesh tailoring for axisymetric refltor. Part 2]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 4(36). pp. 5-14. DOI: 10.17223/19988621/36/1.
6. Bukhtyak M.S. (2017) Obobshchenie minimal'nykh poverkhnostey i modelirovanie formy konstruktsii iz ortotropnogo materiala [Generalization of minimal surfaces and simulation of the shape of an orthotopic material construction]. Tomsk State University Journal of Mathematics and Mechanics. 45. pp. 5-24. DOI: 10.17223/19988621/45/1.
7. Gryanik M.V., Loman V.I. (1987) Razvertyvaemye zerkal'nye antenny zontichnogo tipa [Deployable Mirror Umbrella-Type Antennas]. Moscow: Radio i svyaz'.