Уравнения первого приближения теории устойчивости

Подчеркнем также значение теорем Н. Г. Четаева, дающих признаки устойчивости и неустойчивости для того случая, когда система дифференциальных уравнений первого приближения неправильная [c.130]

Ссылка автора на теорему Ляпунова ошибочна, а его точка зрения на значение метода малых колебаний при рассмотрении частных практических вопросов может ввести читателя в заблуждение. Метод малых колебаний приводит к исчерпывающему ответу, если все корпи характеристического уравнения имеют действительные отрицательные части или в том случае, когда хотя бы один из них имеет положительную вещественную часть. Если же имеются корни, действительные части которых равны нулю, то нельзя судить об устойчивости и неустойчивости по первому приближению, так как все будет зависеть от членов более высокого порядка в уравнениях возмущенного движения. Если псе корпи чисто мнимые, то требуется дополнительное исследование. Обычно это встречается при исследовании устойчивости консервативных систем, по в этих случаях можно вывести необходимое заключение из анализа интеграла энергии. Если в рассмотрение входят диссипативные силы, что обычно и бывает при решении технических проблем, то можно потребовать, чтобы все корни характеристического уравнения имели отрицательные действительные части. В тех случаях, когда все же нельзя удовлетворить этому условию и когда входит, например, один нулевой корень, следует обратиться к исследованиям особых случаев" Ляпунова или изменить постановку задачи, что иногда бывает возможно. [c.425]

Сделаем по поводу полученных результатов два замечания. Во-первых, устойчивость по первому приближению еще не означает устойчивости при рассмотрении точных уравнений (гл. XIX). Кроме того, в этом случае мы лишены возможности вывести суждение об устойчивости из интеграла энергии, как это мы делали в теории малых колебаний (гл. IX). Во-вторых, если система устойчива при рассмотрении точных уравнений, а также в первом приближении, то это связано с влиянием линейных членов Ti в выражении для L. Благодаря им в уравнениях движения появляются гироскопические члены. При отсутствии слагаемых мы имели бы задачу о движении в поле консервативных сил, а для такого поля потенциальная функция в точках Ni и имеет максимум, и эти точки являются положениями неустойчивого равновесия. [c.570]

Признаки устойчивости или неустойчивости движения, в зависимости от знаков корней характеристического уравнения нормальной системы первого приближения, составляют содержание теорем Ляпунова об устойчивости автономных систем по первому приближению. [c.74]

Однородные линеаризованные уравнения теории упругой устойчивости — основной рабочий инструмент этой теории — относятся к разделу математики, называемому задачи на собственные значения (см. приложение I). Кроме однородных линеаризованных уравнений, служащих для определения точек бифуркации, в теории упругой устойчивости широко применяют неоднородные линеаризованные уравнения для приближенного описания поведения систем с начальными неправильностями при малых, но конечных значениях отклонений. Такие уравнения достаточно полно характеризуют поведение систем вблизи точек бифуркаций первого типа (см., например, 18). [c.26]

Значительное расхождение результатов линейной теории устойчивости и эксперимента вызвало необходимость обратиться к решению нелинейных уравнений. Большинство решений нелинейных уравнений было получено в первом приближении при аппроксимации функции прогиба тремя-четырьмя членами тригонометрического ряда [c.119]

Обычно вихри одного ряда располагаются не посередине между вихрями другого ряда. Все вихревые дорожки, которые удовлетворяют этому уравнению, являются неустойчивыми во втором приближении, в то время как все другие вихревые системы неустойчивы уже в первом приближении. По фотографиям, полученным различными исследователями, числовые значения кЦ не одинаковы, поскольку кЦ зависит от времени [26—28]. При больших дозвуковых скоростях образовавшиеся вихри быстро затухают и дорожка становится визуально ненаблюдаемой. Тем не менее происходит периодический отрыв потока. Измерения поля скоростей с помощью термоанемометров и приближенные вычисления показали, что данные, полученные с помощью термоанемометров, недостаточны для характеристики вихревой дорожки 129, 30]. Было установлено, что метод расчета, предложенный в работе 129], может дать более подробную информацию о вихрях [301. Так как результаты не согласуются друг с другом, можно сказать, что в настоящем виде теория устойчивости вихревой дорожки не удовлетворительна. Теория устойчивости первого приближения достаточно точно описывает физические явления, но математический анализ предсказывает неустойчивость, указывая, что упорядоченное расположение вихрей не может сохраняться. [c.90]

Случай Л =0. В этом случае диск представляет собою консервативную систему. Корни характеристического уравнения (3.26), не считая двух нулевых корней, оказываются либо действительными (и тогда движение диска неустойчиво), либо чисто мнимыми сопряженными. В последнем случае обычно считают, что система обладает консервативной устойчивостью. Однако с точки зрения теории устойчивости по первому приближению (см. 1) последний случай является так называемым сомнительным случаем Ляпунова и поэтому требует дальнейшего исследования. Это исследование было проведено в 2 гл. II. [c.310]

Пользуясь своим вторым методом, А. М, Ляпунов решил задачу об устойчивости по первому приближению, независимо от членов выше первого порядка в функциях Хд", в решении этой задачи он видел свое главное достижение. Случаи, когда первое приближение не решает вопроса об устойчивости, названы Ляпуновым критическими. В некоторых из критических случаев установившихся движений, а именно, в случаях одного нулевого корня, пары чисто мнимых корней и двух нулевых корней характеристического уравнения, а также в некоторых случаях периодических движений Ляпунов дал решение задачи об устойчивости. В замечательной работе Ляпунова общая теория дифференциальных уравнений получила существенное развитие. [c.10]

Следует отметить, что положительное решение вопроса о существовании функций.Ляпунова не только обосновало универсальность второго метода Ляпунова, но и позволило развить теорию устойчивости движений по первому приближению, при постоянно действующих возмущениях, при вариациях параметров, при наличии запаздываний и т. п. Это объясняется тем, то наличие функций Ляпунова обычно позволяет доказать сохранение соответствующих свойств при малых изменениях правых частей уравнений (1.1). [c.20]

Мы не будем останавливаться на построении решений уравнений в частных производных, которые выполнил Ляпунов и которые также играют. большую роль в теории устойчивости. Что же на основе этого первого метода Ляпунов получил при других предположениях относительно характеристических чисел первого приближения (7) Те случаи, когда первое приближение (7) имеет только нулевые и отрицательные характеристические числа, Ляпунов назвал сомнительными. В этих случаях характер поведения общего решения или всей совокупности интегральных кривых вблизи начала координат (точки покоя) определяется коэффициентами при нелинейных членах правых частей дифференциальных уравнений. Ляпунов рассмотрел тот случай системы (4), когда коэффициенты рм и. ....постоянные и когда одно характеристическое число матрицы Р = р г нулевое, а все остальные [c.71]

Во второй — пятой главах рассмотрены задачи теории гамильтоновых систем и ее приложений. Вторая глава посвящена линейным гамильтоновым системам. Приводятся результаты Ляпунова об устойчивости линейных гамильтоновых систем с постоянными или периодическими коэффициентами. Для устойчивых систем в случае простых корней характеристического уравнения строятся конструктивные алгоритмы приведения системы к нормальным координатам. Тут же приводится теорема Ляпунова — Пуанкаре о характеристическом уравнении гамильтоновых систем и рассматривается задача о параметрическом резонансе в гамильтоновых системах, содержащих малые периодические возмущения. В последнем параграфе второй главы получены области параметрического резонанса в первом приближении по малому параметру и приведены необходимые расчетные формулы. [c.11]

Решения задачи (114), (1.16), (1.17) мы будем по аналогии с теорией обыкновенных линейных дифференциальных уравнений называть решениями Флоке, а и — показателем Флоке. Как окажется в дальнейшем, используя уравнения лучевого метода в малом, эту задачу в некотором смысле можно решить, точнее свести к решению обыкновенных дифференциальных уравнений. Для возможности такого решения геодезическая I должна удовлетворять условию устойчивости в первом приближении (ср. гл. 4, 7). При этом оказывается, что существует счетное множество показателей Флоке и, которым соответствуют решения задачи (1.14), (1.16), (1.17). Каждому и отвечает конечное множество решений Флоке. [c.235]

В связи с этим в практических инженерных расчетах, в част-рости, в теории автоматического регулирования, большое распространение получили приближенные методы, одним из основоположников которых стал профессор Петербургского Технологического института И. А. Вышнеградский (1831—1895). В 1876 г. Ц. А. Вышнеградский впервые применил свой приближенный метод к задаче об устойчивости регуляторов прямого действия. Основной предпосылкой метода Вышнеградского было допущение, что свойства системы в отношении устойчивости установившегося ее движения обнаруживаются уже в тех малых возмущенных движениях, которые возникают около невозмущенного движения в течение небольшого промежутка времени вслед за моментом сообщения системе достаточно малого начального возмущения. На этом основании при решении вопросов об устойчивости движения в уравнениях возмущенного движения отбрасывались все члены выше первого порядка (относительно координат и скоростей) и по форме интегралов линеаризованных уравнений делались заключения об устойчивости невозмущенного движения. Совокупность методов исследования устойчивости на основании линеаризованных уравнений составляет содержание теории первого приближения. [c.425]

Оригинальный приближенный метод интегрирования дифференциальных уравнений теории упругости был разработан профессором Петербургского политехнического института и Морской академии И. Г. Бубновым (1872— 1919). Впервые этот метод Бубнов описал в 1911 г. в отзыве на только что упомянутое сочинение Тимошенко, представленное на премию имени Журавского. Затем Бубнов использовал свой метод для решения задач на устойчивость пластин, важных в расчетах обшивки корабельного корпуса. Такие задачи разобраны в известном курсе Бубнова Строительная механика корабля (СПб., 1912). Бубнову, как и А. Н. Крылову, принадлежат очень большие заслуги в теории и практике кораблестроения. В частности, он явился в России пионером строительства подводных лодок, первая из которых была спу-ш ена на воду в 1903 г. [c.263]

Для анализа местной потери устойчивости могут быть использованы технические теории пологих ортотропных или многослойных цилиндрических оболочек [18,62]. В первом случае критические параметры нагрузок можно найти из характеристического уравнения, например (5.2) гл. 2, а для приближенного анализа воспользоваться аналитическими выражениями (5.3), (5.11), [c.226]

Дастся изложение основ теории усхойчпвоети движения, базирующееся на общем курсе высшей математики для втузов. Основное внимание уделено наиболее эффективным методам иссл< дова-ния — прямому методу Ляпунова, исследованию устойчивости по уравнениям первого приближения и частотным методам. Отдельные главы посвящены исследованию устойчивости движения но стру -туре действующих сил, устойчивости неавтономных систем, в тол числе систем с периодическими коэффициентами, и систем автоматическою регулирования. [c.2]

Для эффективного построения приближенного решения необходимо предварительно решить уравнения первого или второго приближения (усредненные уравнения). Однако эти уравнения (так же, как и точные) являются дифференциальными, что накладывает определенные ограничения на возможность применения изложенного метода. В большинстве случаев усредненные уравнения, в особенности уравнения первого приближения, более простые и поддаются исследованию. Во многих случаях, в которых общее решение не удается получить, можно найти важные частные решения, например, соответствующие установившимся колебательным процессам. При п = 1 уравнения первою приближения (125) интегрируются в квадратурах при п = 2 для их исследования может быть использована известная теория Пуанкаре. При любом п, если Хо ( ) обращается в нуль в некоторой точке = о, можем рассматривать квазистатическое решение j = уравнений первого приближения. Для исследования устойчивости этого решения можно поступать обычным образом, составив уравнения для малых отклонений (уравнения в вариациях) [c.86]

Строгие методы теории устойчивости движения могут быть распространены на распределенные системы. При этом, например, вместо функций Ляпунова вводят функционалы Ляпунова, производные от которых по времени в силу уравнений движения обладают определен-Егыми свойствами. По этим свойствам судят об устойчивости (неустойчивости) невозмущенного движения. Если модель распределенной системы линейна или если для выводов об устойчивости используют уравнения первого приближения (уравнения в вариациях), то анализ устойчивости приводит к некоторым обобщенным задачам о собственных значениях. [c.461]

Соотношения (П.381Ь)—искомые уравнения в вариациях. Они соответствуют уравнениям (П.331Ь) первого приближения В теории устойчивости движения А. М. Ляпунова. Вместе с уравнениями (11.379) уравнения (П.381Ь) составляют систему 2д уравнений с неизвестными функциями х, Х2, х , 6x1, [c.382]

Ляпунов дал строгое решение вопроса о том, когда при исследовании задачи об устойчивости движения можно ограничиваться рассмотрением первого приближения. Он установил особые слу 1аи, при которых использование первого приближения не решает задачу об устойчивости. Большой заслугой его явилось подробное исследование уравнений, в которых коэффициентами являются периодические функции с одним и тем же периодом. Оп указал признаки устойчивости и неустойчивости для периодических двингений. Отметим еще, что он впервые доказал теорему, согласно которой положение равновесия при некоторых дополнительных условиях неустойчиво, если в положении равновесия потенциальная энергия не мини- [c.248]

В инженерной практике широко распространены конструкции, элементы которых имеют полости или отсеки, содержащие жидкость, иапример, объекты авиационной и ракетно-космической техники, танкеры и плавучие топливозаправочные станции, суда для перевозки сжиженных газов и стационарные резервуары, предназначенные для хранения нефтепродуктов и сжиженных газов, ректификационные колонны и т. д. В большинстве случаев жидкость-заполняет соответствующие полостн или отсеки лишь частично, так что имеется свободная поверхность, являющаяся границей раздела между жидкостью и находящимся над ней газом (в частности, воздухом). Обычно можно считать (за исключением особых случаев движения тела с жидкостью в условиях, близких к невесомости, которые здесь не рассматриваются), что колебания жидкости происходят в поле массовых сил, гравитационных и инерционных, связанных с некоторым невозмущенным движением. Как правило, это поле можно в первом приближении считать потенциальным, а само возмущенное движение отсека и жидкости — носящим характер малых колебаний, что Оправдывает линеаризацию уравнений возмущенного движения. Ряд актуальных для практики случаев возмущенного движения жидкости характеризуется большими числами Рейнольдса, что позволяет использовать при описании этого движения концепцию пограничного слоя, считая, кроме того, жидкость несжимаемой. Эти гипотезы лежат в основе теории, излагаемой ниже [23, 28, 32, 34, 45, 54J. Учету нелинейности немалых колебаний жидкости посвящены, например, работы [15, 26, 29, 30]. Взаимное влияние колебаний отсека и жидкости при ее волновых движениях может сильно изменять устойчивость системы, а иногда порождать неустойчивость, невозможную при отсутствии подвижности жидкости. В качестве примера можно привести резкое ухудшение остойчивости корабля при наличии жидких грузов и Динамическую неустойчивость автоматически управляемых ракет-носителей и космических аппаратов с жидкостными ракетными двигателями при неправильном выборе структуры или параметров автомата стабилизации. Поэтому одной из основных Задач при проектировании всех этих объектов является обеспечение их динамической устойчивости [9, 10, 39, 43]. Для гражданских и промышленных сооружений с отсеками, содержащими жидкость, центр тяжести при исследовании их динамики смещается в область определения дополнительных гидродинамических нагрузок, например при сейсмических колебаниях сооружения [31]. [c.61]

Первое систематическое рассмотрение устойчивости равновесия упругих тел принадлежит Дж. Брайану Он выяснил пределы применимости теоремы Кирхгофа и показал, что при условии малых деформаций она отпадает, если только один или два размера тела можно считать малыми. При этом явление неустойчивости может иметь место в пределах упругости, если произведение модуля упругости Е на квадрат отношения малого размера к конечному будет того же порядка, что и предел упругости материала. Дальнейшая разработка общей теории устойчивости равновесия упругих тел принадлежит Р. Саусвеллу Он устраняет ограничение относительно малости деформаций и оперирует с идеальным телом бесконечно большой прочности. При этих условиях и тела, у которых все размеры одного порядка, могут оказаться в состоянии неустойчивого равновесия. Исходя из однородного напряженного состояния тела, Р. Саусвелл дает точкам тела весьма малые перемещения и, v, w ) и для этой отклоненной формы пишет дифференциальные уравнения нейтрального равновесия, причем считает начальные деформации конечными. То соотношение между внешними силами и размерами тела, при котором полученные уравнения дают для и, у и w решения, удовлетворяющие условиям на поверхности, определяет критическое значение нагрузки в рассматриваемом случае. Применяя свой общий метод к тонким стержням и пластинкам, Р. Саусвелл нашел, что имеющееся решения задач устойчивости являются лишь первыми приближениями, хотя и вполне достаточными для практических приложений. Мы в дальнейшем ограничимся этими приближенными решениями, отсылая интересующихся теорией вопроса к работе Р. Саусвелла. [c.258]

Для н-гексана константы в (9.52) найдены Тимпаном по значениям критических параметров и по удельному объему жидкости при 20 °С и атмосферном давлении. Граница термодинамической устойчивости жидкости по (9.52)-вполне удовлетворительно согласуется с другими определениями (см. рис. 80, 81). Формулу Фюрта и уравнение Гимпана можно рекомендовать в качестве первого приближения для оценки положения спинодали жидкости. Уравнение состояния по ячеечной теории, протабу.чиро-ванное в [221], и уравнение по дырочной модели в приближении Оно [39] приводят к заметно более низким [c.271]

Советская научная литература по устойчивости чрезвычайно обширна и весьма богата результатами как в области развития теории, так л в области ее практических приложений (см. А. М. Ляпунов. Библиография . Составила А, М. Лукомская, под редакцией В. И. Смирнова, М.—Л., 1953). Разработка идей Ляпунова ведется по многим направлениям. Здесь надо отметить развитие и применение первого и, особенно, второго методов Ляпунова, установление новых теорем, расширяющих ж углубляющих эти методы анализ существования функций Ляпунова и их эффективного построения исследования устойчивости по первому приближению и в критических случаях, а также при постоянно действу-лопщх возмущениях исследования устойчивости не установившихся и периодических движений, а также уртойчивости на конечном интервале времени развитие теории приводимых и правильных систем, а также качественной теории дифференциальных уравнений распространение методов Ляпунова на механические системы, описываемые аппаратом, отличным от обыкновенных дифференциальных уравнений (в особенности на сплошные среды), и многие другие. В последние годы выяснилось, что метод функций Ляпунова можно с успехом применять и в получении оценок приближенных интегрирований, и в теории оптимального управления (см. обзор Н. Н, Красовского в настоящем сборнике, стр. 179— 243), и в теории нелинейных колебаний и во многих других разделах науки. По теории устойчивости движения опубликован ряд прекрасных монографий. [c.11]

Задачи устойчивости типичны для тонких и тонкостенных тел. Решения этих задач для стержней, пластин и оболочек строятся обычно на основе приближенных уравнений, в которых используются некоторые кинематические и динамические гипотезы. Имеется несколько путей для получения этих уравнений. Первый, наиболее ранний способ состоит в непосредственном рассмотрении форм движения (равновесия), смежных с невозмущенным. При этом ищется некоторая приведенная нагрузка, которая вводится в уравнение невозмущенного движения. Все рассуждения носят наглядный характер однако в достаточно сложных задачах эта наглядность оказывается обманчивой. Другой путь состоит в использовании нелинейных уравнений соответствующих прикладных теорий. Линеаризуя последние в окрестности невозмущенного движения, получим искомые уравнения. В теории оболочек этот путь использовался X. М. Муштари (1939), Н. А. Алумяэ (1949), X. М. Муштари и К. 3. Галимовым (1957), Н. А. Кильчевским (1963), В. М. Даревским (1963) и другими авторами. Однако в нелинейной теории имеется еще меньше единства взглядов на то, как должны записываться основные уравнения. Следо вательно, идя по этому пути, мы лишь смещаем все трудности в другую, еще менее согласованную область. Третий путь состоит в использовании общих уравнений теории упругой устойчивости (В. В. Новожилов, 1940, 1948). Метод, основанный на соответствующем вариационном принципе, был применен [c.332]

Разгружающие сводики. На основании теории сыпучих тел в первом приближении можно принять, что при формировании осадка над порами фильтрующего материала образуются разгружающие сводики. Причем параболическая форма сводика по М. М. Протодьяконову является наиболее устойчивой (рис. 35). В данном случае осадок, заключенный в разгружающем своди-ке, попадает в фильтрат. Форма наиболее устойчивого свода сыпучего материала описывается уравнением параболы [c.116]

В практическом использовании второй метод Ляпунова значитель но сложнее, чем способы исследования устойчивости по первому при ближению, ибо общих рецептов построения функций Ляпунова не су ществует- Теория устойчивости по первому приближению сводит вопр об устойчивости к чисто алгебраической задаче - к анализу расположени корней характеристического уравнения в комплексной плоскости корне а для этой цели разработаны различные стандартные приемы [14, 33]. [c.44]

В первых пяти главах учебника рассматриваются общие вопросы теории упругости (теория напряжений и деформаций, основные соотношения и теоремы, постановка и лгетоды решения задач теории упругости, плоская задача в декартовых координатах, плоская задача в полярных координатах). В шестой и седьмой главах излагаются основные уравнения теории тонких пластин (гибких и жестких) и некоторые задачи изгиба и устойчивости пластин. Восьмая глава учебника посвящена рассмотрению приближенных методов решения задач прикладной теории упругости (вариационных, конечных разностей, конечных элементов). В девятой главе рассматриваются основы расчета тонких упругих оболочек, причем основное внимание уделено вопросам расчета безмоментных и пологих оболочек. В десятой главе изучаются основы теории пластичности. Здесь рассмотрена и теория расчета конструкций по предельнол1у состоянию. [c.6]

На протяжении XVIII—XIX столетий были составлены различные варианты дифференциальных уравнений задачи п тел, установлены методы их редукции посредством найденных первых интегралов, изучены частные решения, указаны приближенные методы решения задачи с помощью рядов, рассмотрены вопросы устойчивости соответствующих систем дифференциальных уравнений, предложены качественные методы изучения задачи. Созданная на рубеже XIX и XX столетий качественная теория дифференциальных уравнений открыла новые возможности в исследовании проблемы п тел, которые в значительной степени реализованы в трудах ученых XX в. [c.87]

При исследовании ползучести тонких оболочек и решении вопросов устойчивости может иметь значение учет нелинейных слагаемых (квадратов углов поворота) в выражениях для деформаций. Одна из первых работ в этом направлении была выполнена А. С. Вольмиром и П., Г. Зыкиным [31, 32]. Здесь рассматривалась квадратная цилиндрическая панель с начальным прогибом при продольном сжатии. Для решения задачи о прощелкивании панели в условиях ползучести используется. приближенное решение нелинейной упругой задачи панели с начальным прогибом. В процессе ползучести этот начальный прогиб растет и рассчитывается с помощью некоторого приближенного приема, не учитывающего перераспределения напряжений в процессе ползучести. За счет переменного начального прогиба меняется значение верхней критической нагрузки, определяемой уравнениям-и упругой задачи, соответствующее ее прощелкиванию. Когда ве-,личина прогиба достигает значения, при котором соответствующая верхняя критическая нагрузка для упругой панели станет равной действующей нагрузке, произойдет прощелки-вание панели. Существенным результатом этой работы явилось определение критического времени, по истечении которого оболочка скачком перейдет в новое состояние. Учет перераспределения напряжений в процессе ползучести в этой схеме при использовании, как и в [32], теории старения проводился в работе [79]. Аналогичные задачи для сжатой цилин- дрической панели при нелинейной ползучести рассматривались в [60, 95]. [c.272]

Во-первых, общие уравнения нелинейной теории упругости используются для обоснованного вывода уравнений устойчивости для тонких и тонкостенных тел. Работы этого направления (В. В. Новожилов, 1940, 1948 В. В. Болотин, 1956, 1965 А. И. Лурье, 1966, и др.) уже обсуждались в 3. Во-вторых, решения задач, полученные на основе теории упругости, могут быть использованы для оценки точности и установления границ применения известных приближенных решений. К этому направлению относятся работы Л. С. Лейбензона (1917) и А. Ю. Ишлинского (1954). Заметим, что в этих работах в качестве уравнений для описания форм равновесия, смежных с невозмущенной формой, предлагалось использовать классические уравнения теории упругости внешние силы входили при этом только в возмущенные граничные условия. Этот подход обсуждался недавно А. Н. Гузем (1967). В-третьих, необходимость в привлечении уравнений теории упругости возникает в задачах об устойчивости пластин и оболочек, находящихся в контакте с упругим материалом пониженной жесткости. Применительно к слоистым пластинам с мягким наполнителем этот подход развивался А. П. Вороновичем (1948), В. Н. Москаленко (1964) и другими. Устойчивость цилиндрических оболочек с мягким упругим ядром рассматривалась А. П. Варваком (1966). Типичным для этих задач является применение теории пластин и оболочек к несущим слоям и трехмерной теории упругости — к заполнителю. [c.346]

Энциклопедия по машиностроению XXL

Оборудование, материаловедение, механика и ...