Семинар NAS101 | 2006 | MSC.Software Corporation Постоянное представительство в СНГ Москва Раздел 7 Линейный анализ устойчивости
2 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Теория устойчивости Запишем уравнение равновесия конструкции, к которой приложена система постоянных сил, в следующем виде: [ K ] { u } = { P } Рассмотрим эффект дифференциальной (геометрической) жесткости. Дифференциальная жесткость [ K d ] появляется в результате учета членов высокого порядка в зависимостях деформация-перемещение. Эти зависимости подразумевают, что перемещения в конструкции не зависят от интенсивности нагрузки.
3 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Теория устойчивости (продолжение) Пусть произвольный скалярный множитель для другой "интенсивности" нагрузки. При нагружении конструкции данной силой с различной интенсивностью, могут быть найдены несколько положений неустойчивого равновесия. Эти положения равновесия являются решениями задачи на собственные значения.
4 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Решение задачи на собственные значения [ K – K d ] { } = 0(1) Решение нетривиально (отлично от нуля) только для определенных значений = i для i = 1, 2, 3,…, n, которые делают матрицу [ K – K d ] сингулярной.
5 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Каждому собственному значению i, соответствует единственный собственный вектор { i }. { i } может быть масштабирован с помощью любого скалярного множителя и по прежнему оставаться решением уравнения (1). Компоненты вектора { i } - вещественные числа. Решение задачи на собственные значения (продолжение)
6 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation SOL 105Линейная устойчивость SOL 106Нелинейная устойчивость Ограничения для SOL 105 Требования к конструкции до потери устойчивости: Перемещения должны быть малы. Напряжения должны быть в упругой области (и линейно зависеть от деформаций). Последовательности решений для задач устойчивости
7 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример:Три класса колонн (под центральной нагрузкой, материал без дефектов) Тонкая Средняя Короткая Теряет форму из-за потери устойчивости. Прогиб, при потере устойчивости мал и критическая нагрузка не достигает предела текучести материала. Это колонна Эйлера. Причина потери формы – комбинация текучести и потери устойчивости. Прогиб при потере устойчивости мал, но некоторые напряжения превышают предел упругости материала. Причина потери формы – текучесть (также, как при тестовом сжатии образца ) Последовательности решений для задач устойчивости (продолжение)
8 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Примечание: SOL 105 может применяться для конструкций с небольшими дефектами материала или с очень малым эксцентриситетом нагрузки (например, когда нагрузка направлена не строго в центр и приводит к небольшому изгибу). Здесь инженеру надо исходить из здравого смысла. Те же соображения применимы при анализе пластин. Последовательности решений для задач устойчивости (продолжение)
9 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Колонна под нагрузкой с большим эксцентриситетом Прощелкивание тонкой оболочки (подобно дну бака) Перед потерей устойчивости возникают большие изгибающие напряжения, такие же как и сжимающие осевые напряжения. Перед потерей устойчивости возникают большие прогибы, и возможно не упругое поведение конструкции. Примеры нелинейного анализа устойчивости
10 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation (Для справки необходимо смотреть раздел 13, MSC NASTRAN Linear Statics Users Guide ) Секция CASE CONTROL должна содержать не менее двух SUBCASE. Запросы на вывод, которые относятся только к решению статической задачи должны быть помещены в первом SUBCASE. METHOD должен появиться в отдельном SUBCASE для выбора записей EIGB или EIGRL из секции BULK DATA для решения задач устойчивости. Если имеется несколько статических решений, тогда используйте команду STATSUB для выбора варианта одного из статических решений для дальнейшего решения задачи устойчивости. Правила для анализа устойчивости SOL 105
11 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Правила для анализа устойчивости SOL 105 (продолжение) Если необходимо, могут использоваться различные условия в SPC узлах в SUBCASE статического решения и SUBCASE решения задачи устойчивости. Запросы на вывод могут быть помещены в любом выбранном SUBCASE. Запросы на вывод, действующие одновременно и в статическом расчете и в анализе устойчивости могут быть помещены выше уровня SUBCASE.
12 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Секция EXECUTIVE CONTROL SOL 105 Секция CASE CONTROL SUBCASE 1 LOAD = M Определяет условия статического нагружения (LOAD, TEMP, DEFORM) SUBCASE 2 METHOD = N STATSUB = i Выбирает метод нахождения собственных значений Выбирает SUBCASE статического решения чтобы использовать его для решения задачи устойчивости (по умолчанию берется первый SUBCASE) Записи для линейного анализа устойчивости
13 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Секция CASE CONTROL должна содержать не менее двух SUBCASE. Секция BULK DATA Определение условий статического нагружения EIGB Данные для нахождения собственного значения или EIGRL Данные для нахождения собственного значения по методу Ланцоша. Записи для линейного анализа устойчивости (продолжение)
14 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation EIGRL Рекомендуемая запись для расчета устойчивости Определяет данные, необходимые для проведения решения задачи на собственные значения и анализа устойчивости методом Ланцоша EIGRLSIDV1V2NDMSGLVLMAXSETSHFSCLNORM EIGRL Запись EIGRL
15 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Поле Содержание SIDИдентификатор набора (уникальное целое число > 0) V1, V2Анализ вибрации: диапазон интересующих частот. Анализ устойчивости: диапазон интересующих l. Если необходимы все моды ниже определенной частоты, то установите V2 на это значение и оставьте V1 пустым. Не рекомендуется приравнивать V1 к нулю. Эффективнее использовать небольшое отрицательное значение или оставить его пустым. NDЧисло необходимых корней (целое > 0 или оставлять чистым) MSGLVLУровень диагностики (целое, от 0 до 3 или оставлять чистым) MAXSETЧисло векторов в блоке (целое, от 1 до 15 или оставлять чистым) Запись EIGRL (продолжение)
16 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример - простая колонна Эйлера Задача: Найти критическую нагрузку и соответствующую ей первую форму потери устойчивости цилиндрического стержня. 21 дюйм Свободно Зафиксировано Круглое сечение Диаметр = 0.25 дюйма E = 30.E+6 фунтов/дюйм 2 I = 1.917E-4 дюйм 4 Площадь = 4.909E-2 дюйм 2
17 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример - простая колонна Эйлера (продолжение) Теоретическое решение где L eff – эффективная длина колонны L eff = 2L – для свободно опертой колонны
18 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример - простая колонна Эйлера (продолжение) Модель MSC Nastran Решение MSC Nastran Значение силы в записи FORCE P cr = ×1.0 = фунта Собственное значение
19 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример – простая колонна Эйлера – входной файл
20 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример – простая колонна Эйлера – выходной файл Первое собственное значение: Pcr = 1 x 10 фунтов = фунта Первый собственный вектор (дает форму потери устойчивости)
21 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation MSC Seminar Notes, MSC NASTRAN Material and Geometric Nonlinear Analysis: MSC Nastran Linear Static Analysis Users Guide, Section 13. MSC Nastran Verification Problem Manual ( Version 64, January 1986 Edition ): Problem A, Lateral Buckling of a Cantilever Beam Problem A, Simple Frame Analysis with Buckling Problem S, Euler Buckling of a Simply Supported Beam Литература по анализу устойчивости
22 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation MSC Nastran Demonstration Problem Manual ( Version 64, March 1985 Edition ): Under Elastic Stability Analysis, see Demonstration Problem D0504A, Flexural Buckling of a Beam MSC Nastran Application Notes October 1978Buckling and Real Eigenvalue Analysis of Laminated Plates September 1979Static Stability of Structures with Nonlinear Differential Stiffness February 1982Elastic-Plastic Buckling of a Thin Spherical Shell November 1985Nonlinear Buckling Analysis Литература по анализу устойчивости (продолжение)
Семинар NAS101 | 2006 | MSC.Software Corporation Постоянное представительство в СНГ Москва Пример 9 Анализ устойчивости пластины
24 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример 9 (продолжение) Описание модели Та же самая модель панели, что и в примере 5, без подкреплений. В этой модели применены следующие граничные условия: Опирание на левом конце Ролики на правом конце Нулевые вертикальные перемещения на верхней и нижней гранях Приложим 100 фунт/дюйм 2 сжимающие нагрузки к правому краю пластины Общая нагрузка на правой стороне = (100) (8) (.01) = 8 Приложим 1 фунт в каждый из узлов 11 и 55 сетки Приложим 2 фунта в каждый из узлов 22, 33 и 44 сетки
25 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример 9 (продолжение) – Граничные условия Простое опирание Опора на роликах Защемление в вертикальном направлении
26 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Пример 9 (продолжение) – Приложенные нагрузки
27 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Жесткие элементы Лагранжа Новые жесткие элементы Лагранжа поддерживают дифференциальную жесткость Новые жесткие элементы Лагранжа поддерживают задачи на устойчивость В настоящее время решение SOL 105 (линейная задача устойчивости) поддерживает только метод исключений Лагранжа (error )
28 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Усовершенствования жестких элементов (продолжение) Пример: Анализ устойчивости модели двумя методами Жесткие элементы Лагранжа Жесткая балка x y K 5 = 100., K 6 =
29 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Усовершенствования жестких элементов (продолжение) Входной файл (элементы Лагранжа) $ $ lagrange1. dat $ SOL 105 CEND TITLE = BUCKLING ANALYSIS - RBAR SUBTI = LAGRANGE ELIMINATION METHOD DISP = ALL SPC = 10 RIGID = LGELIM SUBCASE 1 LABEL=STATIC PRELOAD CASE LOAD = 100 SUBCASE 2 LABEL = BUCKLING CASE METHOD = 10 BEGIN BULK EIGRL,10,,,10 FORCE,100,4,0,-100.0,1.0,0.0,0.0 CELAS2,101,100.0,3,5 CELAS2,102,120.0,3,6 GRID, 3,,2.0,0.0,0.0 GRID, 4,,4.0,0.0,0.0 RBAR, 3,3,4,123456,,, SPC1,10,1234,3 ENDDATA Входной файл (жесткие балки) $ $ lagrange2. dat $ SOL 105 CEND TITLE = BUCKLING ANALYSIS - USE STIFF BEAM DISP = ALL SPC = 10 SUBCASE 1 LABEL=STATIC PRELOAD CASE LOAD = 100 SUBCASE 2 LABEL = BUCKLING CASE METHOD = 10 BEGIN BULK EIGRL,10,,,10 FORCE,100,4,0,-100.0,1.0,0.0,0.0 CELAS2,101,100.0,3,5 CELAS2,102,120.0,3,6 CBEAM,100,100,3,4,0.,0.,1. PBEAM,100,1,100.,100.,100.,,100. MAT1,1,1.E7,,.32 GRID, 3,,2.0,0.0,0.0 GRID, 4,,4.0,0.0,0.0 SPC1,10,1234,3 ENDDATA
30 MSC Confidential Семинар NAS101 | 2006 | Раздел 7 | MSC.Software Corporation Усовершенствования жестких элементов (продолжение) Результаты (жесткие элементы Лагранжа) Результаты (жесткие балки)