Приложение вариационно-разностного метода к расчету продольно-поперечного изгиба пластин силами инерции

Полный текст

Введение. Элементы тонких пластин, используемые в военно-морских и аэронавигационных структурах, часто подвергаются нормальным и сдвигающим силам, действующим в плоскости пластин [1]. Решения проблем устойчивости имеют практические значения при конструировании ракет, высокоскоростных самолетов, подвергнутых значительным ускорениям [2]. Следует отметить возможность приложения этих решений для проектирования баннеров, тонких ограждений в виде вертикальных консольных пластин, находящихся под действием собственного веса. Поэтому можно заявлять об актуальности рассматриваемой темы. Дифференциальная формулировка задачи устойчивости. Дифференциальная классическая формулировка задачи о продольно-поперечном изгибе изотропной пластинки в декартовой системе координат известна в литературе, она представлена уравнением равновесия [3] і д4 w 1XF + dx2 dy2 д4 w д4 w dy4 хг д2 w хг д2 w д2 w - Nx -- + Ny -- + 2S + q x dx2 y дy2 xy дxдy z (1) в котором D - Eh /12(1 - д2) - цилиндрическая жесткость; Е - модуль Юнга; д - коэффициент Пуассона; h - толщина пластинки; w - w(x, y) - функция прогиба. Мембранные усилия N - Nx ( x, y), Ny - Ny (x, y), SXy - Sy (x, y) правой части (1) выражаются через перемещения u - u(x,y), v - v(x,y) базисного слоя на основе закона Гука [4]: N- Eh 1 - д2 ды dv - + д- dx dy N- Eh 1 - д2 dv du - + д- dy dx Eh (du dv Sxy --1 -+- 2(1 + д) ^éy dx (2) Функции перемещений для (2) могут быть вычислены с помощью уравнений равновесия Сен-Венана [5]: д_ dx du dv 1-dx dy 1 - д ( д2u ' д2u 1 + д dx2 + dy2 д ( du dv ^ 1 - д ( д2 v д2 v dy ^ dx dy J 1 + д ^ dx2 dy2 2(1 - д) Eh 2(1 - д) ( Eh ' q,x, (3) (4) В (3), (4) для задач о действии ускорений qx ( x, y) - hX , qy ( x, y) - hY , qz ( x, y) - 0 , в которых объемные силы равны X - pax и Y - pay , зависят от плотности материала р и ускорений ax, ay по направлениям осей координат Ox и Oy соответственно. Таким образом, сформирована система трех уравнений (1)-(4) относительно трех неизвестных функций u, v, w. К этим уравнениям добавляются статические и геометрические граничные условия. Отметим, что в уравнения (3) и (4) не входят прогибы; задача распадается на две: в первой задаче вычисляются функции u и v по формулам (3), (4), а во второй -определяются функции N - Nx(x,y), Ny - Ny(x,y), Sxy - S y ( x, y ) по формулам (2) и реализуется проблема собственных чисел, т. е. проблема расчета устойчивости пластины (1). Обоснование предпочтения интегральной формулировки задачи устойчивости. Критерий устойчивости. При конечно-разностной аппроксимации дифференциальных уравнений (1), (3), (4), а именно, при рассмотрении свободных от закрепления углов пластинки, возникают алгоритмические трудности, y 103 Вестник СибГАУ. 2014. N 4(56) связанные с исключением законтурных узлов. Главной же особенностью является то, что аппроксимация вторых производных функций прогибов в (1) центральными конечными разностями дает несимметричную матрицу внутренних усилий, что затрудняет решение проблемы собственных чисел разрешающей дискретной краевой задачи. Этот факт явно прослеживается при формировании соответствующих матриц продольно-поперечного изгиба стержней [6]. Поэтому целесообразно к расчету пластин на устойчивость от действия сил инерции применить энергетические принципы. Здесь следует обратиться к основам теоретической механики и упомянуть теорему Лагранжа-Дирихле и Ляпунова, устанавливающую равновесие механической системы и определяющую критерии равновесных положений [7]. Для деформируемых конструкций, в частности для пластин, приложение энергетических принципов к проблеме устойчивости приводится в большом количестве научных разработок. Назовем здесь работы [1-5; 8-14], использующие как численные, так и аналитические методы расчета. Пусть ЭЛ = ЭЛ (^ w) - (5) функционал Лагранжа конечного числа переменных w1,wl, ...,wn краевой задачи. Приращение полной потенциальной энергии деформирования для вариаций w1 + 5wj,w2 + §w2,...,wn + §wn с удержанием членов не выше второго порядка имеет вид ДЭЛ = ЭЛ (Wj + 5wj, w2 + §w2,..., wn + §wn ) - -ЭЛ (WJ, W2,..., Wn ) = д ЭЛ + §2ЭЛ /2. (6) Здесь первая и вторая вариации (5) соответственно равны 5Эл = (дЭЛ / dwk )bwk, (7) 52Эл =(д2Эл / dwkdw, ))w, (8) с суммированием по повторяющимся индексам: k = 1, 2, ..., n , I = 1, 2, ..., n . Для любого равновесного состояния 5ЭЛ = 0 [8]. В задачах устойчивости о характере равновесия судят по знаку второй вариации функционала: устойчивое равновесие, если 52ЭЛ > 0 для всех допустимых виртуальных перемещений; неустойчивое равновесие, если хотя бы для одного допустимого виртуального перемещения 52ЭЛ < 0 . Но первая верхняя критическая нагрузка соответствует переходу от устойчивых равновесных положений к неустойчивым, следовательно, для этого значения нагрузки должно быть [10; 14; 15] (9) 52ЭЛ = 0. Введем в (8) операторы варьирования и 52. Теперь запишем (9) так: 52ЭЛ = 52(§1Эл (wJ, W2 , ..., Wn ; 51w1,w2, ..., wn; S2w1,§2w2, ..., §2wn)) = 0. (10) Условие (9) с учетом (10) эквивалентно следующему уравнению {§2 W1 §2W1 .. §2 w1} x д 2 Эл д 2 Эл д 2 Эл dw1dw1 дwJдw2 дwJдwn д 2Эл д 2Эл дwJдw2 дw2 дw2 д2 Эл д2 Эл д2 Эл дwn д^ дwn дw2 " ^n^n 01 w1 §1 W1 W1 = 0, (11) в котором, матрица n x n содержит коэффициенты при вариациях перемещений уравнения (10). Данная система будет иметь решение, отличное от нуля, лишь в том случае, если определитель, составленный из ее коэффициентов, равен нулю [10]. Группировка коэффициентов при независимых вариациях перемещений дает нулевой определитель д 2 Эл д 2 Эл д 2 Эл дwJдw2 дwJдwn д 2Эл д 2Эл дwJдw2 дw2 дw2 д2 Эл д2 Эл д2 Эл дwn дwJ дWnдW2 ' ' дwnдwn = 0. (12) Коэффициенты в (12) являются функциями параметров жесткости конструкции, соотношения ее размеров. Также они являются функциями внутренних усилий, зависящих в рассматриваемой работе от параметра ускорения. Параметр ускорения есть параметр нагрузки, который является единственной неизвестной характеристикой определителя (12); причем его значения, обращающие определитель в нуль, будут соответствовать критическим значениям нагрузки. Корни системы уравнений (12) определяют n собственных чисел, т. е. все критические нагрузки, в том числе наименьшую из них. Преобразования дифференциальных уравнений плоской задачи теории упругости в интегральные уравнения. Используя основы вариационных принципов механики [13], дифференциальные постановки рассматриваемой задачи (1), (3), (4) заменим интегральными формулировками. В начале преобразуем уравнения (3) и (4). Для этого составим уравнение возможной работы внешних и внутренних сил гг Гг ' дИ д£ 5 ff 1 -- + ^ + Xh §u + JJ S 1 дх ду J (дS дИу 5 ДУ + y + Yh §v dxdy = v дх ду J (13) из которого с учетом геометрических уравнений и закона Гука (2) для плоской задачи теории упругости, приемом интегрирования внутренних произведений по частям в прямоугольной области S : 0 < x < a и 0 < y < b получим: - функционал Лагранжа: x + 104 Математика, механика, информатика <Я E Ф л (u, v) =--- x 2(1 - ц2) duY (dv^ 2 du dv 1 -ц(du dv dx J ^ôy J dx dy 2 ^éy dx dxdy + - И (X u + Y v)dx dy + j ( + Txyv)dy + j (u + x, y =0 - его вариацию (14) y=0 5фл = 5,фл = - E (1 - ц2) xjj du d51u dv 35,v ( du 95,v dv d51u 1 +--- + ц I--- + 1 dx dx dy dy 1 - ц I du dv V db1u 55jv 2 1 ôy dx dx dx dy dy dx dxdy + - jj (X Su + Y Sv)dx dy j (Su+"CSv y y=0 j (y5v + tVxSu ) (15) y=0 - и вторую вариацию функционала E 52Фл = 82(81Эл ) = - (1 - ц2) JsJ jj 552u d51u dx dx dy dy d52v d5v ( 55,u d51v d5u 55,v +-2--- + ц I -2--^ + ^--2- | + dx dy +1 - ц I 952u + dS2v 11 95ju + dSjv dx dy dxdy. dx (16) du, i = 1,2, , p ; 51 u, = 1, при l = i; 0, при l Ф i. (17) Для формирования коэффициентов жесткости at j этой же системы уравнений составим вторую вариацию как вариацию первой вариации, реализуя (16) для дискретной задачи: ( p ъох ta s; a s; а\ Л -V t = V д V дФЛ (u, 51u, 52м ), aÿ ^ du ^ k =1 UUk V l=1 du, 52uk , (18) \ 1, при k = i; \ 1, при l = j; 52uk = ] n j ■ 51ul = ] n . [0, при k Ф i, [0, при l Ф j, i = 1,2, ...,p ; j = 1,2, ...,p . Решение сформированной системы уравнений ajjuj = bi дает вектор перемещений uj, по компонентам которого вычисляются внутренние усилия Nx, Ny, Sy согласно (2). Преобразования дифференциальных уравнений продольно-поперечного изгиба в интегральные уравнения. Для преобразования дифференциальной формулировки задачи (1) в интегральную формулировку воспользуемся уравнениями равновесия, полученными по деформируемой схеме в [3]: dQx , dQy + м d2w(x, y) d2w(x, y) + 1 " + Nx 2 dx dy dx У2 w( x, y) dxdy -(y + Syx У +Ny dy2 +q- ( x, y) = 0, Qx = dMx dH Qt =dM„ dH„ (19) dx dy y dy dx На основании (19) составим вариационное уравнение dy dx J\ dy Следует к (14)-(16) добавить главные граничные условия: u = u , v = v , 5iu = 5iu , 5,v = 5iv , i = 1,2 (заданные значения отмечены символом «*»). В (14)... * * * * (16) ax, Txy, Gy, tyx - заданные напряжения. Алгоритм формирования матриц жесткости и правых частей плоской задачи теории упругости в перемещениях. Таким образом, для решения краевой задачи (3), (4) вариационно-разностным методом применяем 5ФЛ = 51ФЛ = 0 , что сообщает функционалу (14) минимальное значение. Рассмотрим приложение формул (15) и (16) для формирования системы разрешающих уравнений с целью вычисления перемещений u = u( x, y), v = v( x, y) пластины. Пусть u = (u1, u2,...,up) - вектор перемещений для p переменных, а 51u = (51u1,51u2,...,51up) и 52u = (52u1, 52u2, ...,52up)-вариации вектора перемещений. Тогда для вычисления коэффициентов bi вектора правой части системы разрешающих уравнений запишем (15) в виде - c- - dФл (u, 51u) b, = 51Фл (u, 51u) = ^---51ul • 11 d2 M,, + dX + d2 H„, d2 M „ d2 w . . . + N--+ dx2 dxdy dxdy dy2 x dx2 d2 w t \ d2 w +N +( +S-y+q- 5wdS = 0, (20) где Mx = D ( æx + цю y ), My = D (æ y + ), Hxy = Hx = = (1 - ц) D xxy, - изгибающие и крутящие моменты; d2 w æx = т~г x dx2 d2 w X xy =■ d2 w 3x = dw / dx, y dy2 ' dxdy 3y =dw / dy - кривизны кривой линии и углы поворота; 5w( x, y) - возможное перемещение базисной поверхности пластины; dS = dxdy. Считая внутренние силы Nx, Ny, Sxy известными функциями, выполнив процедуры интегрирования внутренних произведений в интеграле (20) по частям, получим вариационное уравнение Л-U J’-P' 5ЭЛ = j j D[æx5æx + ^æx5æy + æy5æx) + æy5æy + x = 0 y =0 x=a y=b +2(1 -ц)Xxy5Xxy]dxdy- j j q-5wdxdyx = 0 y = o y=b x=a -j(i?x 5w - Mx 53x )dy y=0 S S X S S l =1 x=0 105 Вестник СибГАУ. 2014. № 4(56) - j (y§w -MyЪЗу )dx y-0 + (Hw + H )5w 0 + v x yx ' lx-0 'y -0 x-0 y-0 Л-и y-u + j j (a(x + NySy53y + Sy(3y53x + 3x53y))dxdy- -b x-a - j (A + S3 )§wdy y=0 - j (Syx3x + Ny3y )5wdx x=0 y-b - 0, (21) y-0 где 5 - оператор варьирования; Rx - Qx + dHxy / dy , Ry - Qy +dHx / dx - реакции. Из (21) следуют функционал Лагранжа 1 x-a y-b Эл - 2 j jD[æX+2джxæy + æy +2У1 - д)xXy]dxdyx-0 y-0 x-a y- y -0 y-b + (Hxy + Hyx )w|x -0 |y-0 + - j j qzwdxdy -j (,w - M,3, )dy x-0 y-0 x-o - j (yw - My3y )dx x=0 +2 j 1 (X + Ny32y + 2Sxy3y3x)dxdy 2 x-0 y-0 y-b - j (x3x + Sxy3y )wdy y-0 - j (yx3x + Ny3y )wdx x-0 y-b (22) y-0 и вторая вариация этого функционала x-0 y-b 52ЭЛ - j j D [52æx5tæx + д(52æx5tæy + 52æy5tæx) + x-0 y-0 + 52 æ y 51æy + 2 У1 - д)5 2 Xxy 51Xxy ] dxdy + x-o y-b + j jy N523, 5Д + Ny 52 3y 5,3y + x-0 y-0 +Sxy (523y513x + 523x513y ) ) dxdy ■ (23) где 5i, i = 1, 2, - операторы варьирования; 5i æ x - d2 5,. w d25i.w 5æ y - -~dr j d2 5,w 5 х - i- ; i xy dxdy dx2 1 y dy2 53 - d5i w / dx, 53y - d5i w / dy - вариации кривизны линий и вариации углов поворота. Формирование матриц жесткости и матриц внутренних усилий проблемы устойчивости. На основании (9) выражение (23) приравняем к нулю, что приводит к постановке задачи устойчивости дискретной задачи. Пусть w - (w1,w2, ...,wp) - вектор перемещений для p переменных, а 5tw - (5^, 5tw2, ...,5jwp) и 52w - (52w1,52w2, ...,52wp) - вариации вектора перемещений. Сформируем для дискретной задачи систему (12) подобно (18): p д Ґ р /vt, £ ,7, s; ,7\ Л " f-1 dw k 3 дЭЛ (w, 5jw, 52w) _ f --51w i dw, 52 #k - 0, (24) 52 wk - 5! wl - 1, при l - j; 1, при k - i; [0, при k Ф i, 1 ‘ [0, при l Ф j, i - 1, 2, ..., p ; j - 1, 2, ..., p . Здесь каждый коэффициент Cÿ - Bÿ - sAij, согласно (23), содержит характеристики жесткости ( Atj ) и внутренних усилий ( Btj ); s - корни уравнения (12). Умножив найденные из (24) Cij на вектор wj, получим Cywj - 0, (25) что представляет формулировку обобщенной проблемы собственных значений [ B ]{w} - s [ A]{w}, (26) где [B] назовем матрицей внутренних усилий; [A] -матрица жесткости; {w} - (wt,w2, ...,wp) - собственный вектор для p переменных; s - собственные числа. Дискретизация дифференциальных операторов конечно-разностными аналогами и аппроксимация континуальной области сеточной. Выполним дискретизацию уравнений (15), (16), (23). Выберем на области пластинки прямоугольную равномерную сетку œÿ - {( X - il x, yj - jl y ), i - 0,1,..., m, j - 0,1, ..., n} на отрезках [0, lx ] и [0, ly ]. Здесь x - xi и y - yj -узлы сетки; lx - lx / m и l - ly / n - шаг сетки, а lx и ly - размеры пластинки по направлениям осей координат x и y. Эту сетку (рис. 1) с узлами i, j назовем основной сеткой. Введем дополнительную сетку с узлами 4, П: - {(,4 - lx /2+ ilx, yy - ly /2 + jly), i - 0, 1,...,m -1, j - 0, 1, ...,n -1}. Интегралы по области в (15) и (16) и второй интеграл от внутренних усилий в (23) заменим численным интегрированием m-1 n-1 jj( )dxdy - ffy )S4n (здесь S4n - площадка S 4-1 n-1 интегрирования, равная lx l во внутренних узлах области пластины), а дифференциальные операторы заменим конечно-разностными аналогами [16]: d5ku dx d5kv dy 5kui+1, j+1 + 5kui+1, j 5kui, j+1 5kui, j 5kvi+1, j+1 + 5 У j+1 - 5kvi+1, j - 5kvi, j d5ku d5kv dy dx 2 X„ 5 kui+1, j+1 + 5 kui, j+1 5 kui+1, j 5 kui, j 2X„ 5kvi+1,j+1 + 5kv>+1,j - 5kvi, j+1 - 5kvi,j 2! x x-0 x-0 x-0 106 Математика, механика, информатика Далее, в первом интеграле (23) интегрирование для членов, содержащих вариации кривизны §k æ x и §k æ , заменим суммированием по узлам основной ГГ m n сетки fj( )dxdy=XS( )sj и применим централь- S i=1 j=1 ные разности. Здесь площадки интегрирования St j равны: XxX - во внутренних узлах области пластины; XxХу /2 - в узлах, расположенных на контуре; XxXy /4 - в узлах, расположенных в углах пластинки: §kwm,j - 2§kwi.j + §kWi-J.j i д 2Skw 5 v Sx2 ,.. 4 X'. J X i д 2Skw 5 ду 2 §kwu j+1 - 2§kw;, j + §kw,, j-1 x 2, si, j "y Кривизну кручения § k X xy вычислим в узлах типа £, J : І д2§^ 5 дхду §,, w. k i+1, j+1 - §kUL j+1 + §Л j-1 - §kU i+1, j и в узлах типа i, n: Ід2§kW'5 = §kWi+J.j+1 -§kui-1,j+1 + §kui-1,j -§kмi+J.j дхду , Здесь k = 1,2. 2 LX. Рис. 1. Конечно-разностная сетка: сплошными линиями обозначены линии основной сетки; штриховыми линиями обозначены линии дополнительной Примеры расчета. Для расчета продольнопоперечного изгиба тонких пластин используем математический пакет Maple [17], позволяющий вычислять собственные числа и собственные значения матриц. Приведем два примера расчета устойчивости от сил инерции. Пример 1. Пластина жестко защемлена по всему контуру; объемные силы равны X = рах = 1 и Y = рау = 0 . Размеры в плане: 1х = 0,8 м; I = 0,6 м; толщина h = 0,01 м. Материал пластинки имеет характеристики жесткости: E = 1010 Па и ц = 0,45 . Плотность материала р = 3 г/см3. Пусть допускаемое ускорение [а ] равно 50 g. Первый этап расчета включает в себя вычисление внутренних мембранных усилий (2), эпюры которых приведены на рис. 2. По найденным перемещениям базисной поверхности пластины (рис. 3) построим ее деформированный вид (рис. 3, а). Решение обобщенной проблемы собственных значений (26) (это второй этап расчета) представляет первая форма поперечного изгиба (рис. 3, б). ,Ш 12 14 16 1$ 2016 а Рис. 2. Эпюры внутренних усилий: а - продольные Nx; б - продольные Ny; в - сдвиговые Sxy (размерность внутренних усилий Н/м) / -1- 1 \ \ \ \ \ \ \ \ \ 10 20 30 40 50 а б Рис. 3. Деформированный вид поверхности пластины: а - перемещения линии сетки срединного слоя; б - выпучивание - первая равновесная форма сетки 107 Вестник СибГАУ. 2014. № 4(56) Инерционная нагрузка здесь такова, что формируемая матрица внутренних усилий имеет на главной диагонали нули и эта матрица не имеет обратной матрицы. Найдем ускорение, отвечающее первому собственному значению, равному 304 Н/м3: ax =304 = Р = 304, = 1000 м/с2. Ускорение слишком велико 3 -10-3 и, видимо, близко к ударным нагрузкам. Заключаем, что заданная пластина не потеряет устойчивость, так как полученное критическое ускорение, равное 100 g, меньше допускаемого [a] = 50 g . Пример 2. В качестве второго примера рассмотрим пластину, закрепленную только с одного торца, когда три других края свободны и не закреплены. В ее плоскости действует постоянное ускорение. Аналогичная задача расчета продольно-поперечного изгиба была приведена в [2], где изучалось действие на вертикальную консольную пластинку постоянного ускорения; для решения задачи применен метод Ритца. Пластина имеет размеры в плане 0,6*0,8 м и толщину 0,1-10-3 м; модуль Юнга материала равен 2 1011 Па; коэффициент Пуассона -0,25; плотность 780 кг/м3. Ускорение ay (это параметр) направлено по оси абсцисс. Требуется определить критические ускорения. На рис. 4 отображен деформированный вид, полученный на первом этапе решении задачи до возникновения критических ускорений. На рис. 5 приведен вид функций внутренних мембранных усилий, приводящих к потере устойчивости. На рис. 6 продемонстрированы полученные формы потери устойчивости пластины, там же отмечены четыре первых критических ускорения, соответствующих этим формам. Таким образом, чтобы достичь при разгоне ускорения, равного, к примеру, 41 м/с2, следует не допустить предыдущих критических ускорений 30 и 40 м/с2. Следует обратить внимание, что система Maple первое и четвертое собственные числа подобрала отрицательными и выдала соответствующие им формы (рис. 6, а, г) в противоположном направлении. На первый взгляд, отрицательное значение собственного числа противоречит критерию потери устойчивости, что если минимальное собственное значение положительно, тогда конструкция устойчива. Здесь рассматривалась обобщенная проблема [ A]{w } = s* [B]{w} (s* = 1/s ). Если всегда [A] > 0, тогда больше нуля должно быть произведение s * [ B ] > 0 , а не собственное значение s*. Дополнительно отметим в этой задаче особенность приложения вариационно-разностного подхода. На свободных краях пластины при формировании матриц жесткости записывались уравнения равновесия как в узлах на контуре, так и в законтурных узлах, поэтому размерности сформированных матриц [A] и [B] в (26) не совпадали. Дополнительные уравнения о равенстве нулю изгибающих моментов и обобщенных реакций на контуре понижали размерность матрицы [A] до размерности матрицы [B]. Эти дополнительные процедуры сложения соответствующих строк между собой и столбцов матриц были автоматизированы для сеток различной густоты, что несколько усложняло общий алгоритм. Реализация же дифференциальной формулировки задачи (1) методом конечных разностей для пластины со свободными краями затруднительна. Рис. 4. Деформированный вид пластинки ,-4 ,1-4 , 2 ,. -Г'г « « 4 Г1Тл1210 8 10і2147;,~Уг7 б ■ І-т- -п Т'"1 H fi s 1о9"Ггт-~,- --'ГЛі-Ий“ Рис. 5. Пространственный вид эпюр внутренних усилий: а - Nx; б - Ny; в - Sxy в 108 Математика, механика, информатика Рис. 6. Формы потери устойчивости при ускорениях: a -29 м/с2; б - 30 м/с2; в - 40 м/с2; г--43,3 м/с2 Заключение. Таким образом, рассмотренный подход к расчету пластин может быть использован для подбора и анализа геометрических и параметров жесткости панелей при заданном ускорении или торможении, применяемых в аэрокосмической технике. Матрица жесткости и матрицы внутренних усилий формировались с помощью прикладной программы интерпретатора Maple, что ограничило возможности сгущения конечно-разностной сетки и численного решения задачи большой размерности, исследования сходимости результатов. Поэтому для дальнейших исследований желательно использование пакетов программ, основанных на алгоритмических языках высокого уровня. Необходимо сказать, что на поиск всего спектра собственных чисел и форм программа Maple затрачивала на два порядка больше времени, по сравнению с временем формирования матриц с применением предложенных здесь методов. Для обеспечения достоверности был применен метод скалярных произведений, позволяющий вычислять максимальное собственное число. Однако это требовало вычисления обратной матрицы жесткости, тогда как матрицы внутренних усилий обратных матриц не имели.

Об авторах

Рашид Альтавович Сабиров

Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнева

Email: aaa@mail.sibsau.ru; rashidsab@mail.ru
кандидат технических наук, доцент, доцент кафедры технической механики и кафедры летательных аппаратов

Список литературы

Дополнительные файлы