Визначення критичного значення стискаючих зусиль, рівномірно розподілених по двох протилежних поперечних торцях прямокутної ортотропної вільно опертої пластини, що відповідає моменту втрати її стійкості.
А. С. Вольмир, Устойчивость деформируемых систем, Москва, Наука, 1967, стр. 374.
Прямокутна ортотропна вільно оперта пластина піддається впливу стискаючих зусиль σ, рівномірно розподілених по двох протилежних поперечних торцях. Визначити критичне значення стискаючих зусиль σcrort, що відповідає моменту втрати стійкості прямокутної ортотропної пластини.
Розрахункова схема - система загального вигляду. Розглядаються дві розрахункові схеми із співвідношеннями сторін пластини a/b = 1,0; 4,0. Сітка скінченних елементів розбита (вздовж осей X та Y загальної системи координат) з кроком 0,075 м. На одному з двох протилежних поперечних торців пластини, схильних до впливу стискаючих зусиль, задається рівномірно розподілене по лінії навантаження з початковим значенням p = σ*h = 1000 кН/м.
Розмір сторони прямокутної пластини, вільної від навантаження (вздовж осі X загальної системи координат) a = 0,6; 2,4 м;
Розмір сторони прямокутної пластини, схильної до впливу стискаючих зусиль (вздовж осі Y загальної системи координат) b = 0,6 м;
Товщина прямокутної пластини h = 0,01 м.
Модуль пружності матеріалу пластини, що відповідає поздовжнім деформаціям осі X загальної системи координат E1 = 5,6*108 кН/м2;
Коефіцієнт Пуассона, що відповідає поперечним деформаціям по осі Y загальної системи координат v21 = 0,3;
Модуль пружності матеріалу пластини, що відповідає поздовжнім деформаціям осі Y загальної системи координат E2 = 2,123*108 кН/м2;
Коефіцієнт Пуассона, що відповідає поперечним деформаціям по осі X загальної системи координат v12 = 0,114;
Модуль зсуву матеріалу пластини G12 = 0,769*108 кН/м2.
Забезпечення граничних умов досягається за рахунок накладання в'язів на вузли опорного контуру пластини у напрямку ступеня вільності Z. На вузли протилежного торця, де задається навантаження, накладаються в'язі за відповідним напрямком (вздовж осі X загальної системи координат). З метою забезпечення геометричної незмінності розрахункової схеми на вузли одного з двох протилежних поздовжніх торців пластини, вільних від навантаження, накладаються в'язі по нормальному до нього напрямку (вздовж осі Y загальної системи координат). У цих же цілях на вузол одного з кутів пластини накладається в'язь у напрямку UZ загальної системи координат.
Початкове значення стискаючих зусиль σ = 1*105 кН/м2.
Задача розв’язується у просторовій постановці (ознака схеми 5).
Для опису роботи пластини використаний СЕ 44 – універсальний чотирикутний СЕ оболонки.
Кількість вузлів: 81; 297. Кількість елементів: 64; 256.
|
Розрахункова схема 1 із співвідношеннями сторін a/b = 1,0
|
|
1-ша форма втрати стійкості для схеми 1
|
2-га форма втрати стійкості для схеми 1
|
3-тя форма втрати стійкості для схеми 1
|
|
Розрахункова схема 2 із співвідношеннями сторін a/b = 4,0
|
|
1-ша форма втрати стійкості для схеми 2
|
2-га форма втрати стійкості для схеми 2
|
3-тя форма втрати стійкості для схеми 2
|
При аналітичному рішенні критичне значення стискаючих зусиль σкрорт, що відповідає моменту втрати стійкості прямокутної ортотропної пластини, визначається за такою формулою:
n, m = 1, 2, 3 … – число напівхвиль форми втрати стійкості у напрямках поперечному та поздовжньому відносно стиснення пластини.
Жорсткісні характеристики ортотропної пластини приймаються з умов еквівалентності жорстким характеристикам підкріпленої пластини:
і визначаються за такими формулами:
Критичні значення стискаючих зусиль σкр для підкріпленої пластини повинні бути зменшені з коефіцієнтом k відносно критичних значень стискаючих зусиль σкр для ортотропної пластини, оскільки при визначенні останніх не враховується складова, що діє на ребра першої:
де F – площа ребра жорсткості, s – число ребер.
Без додаткових вузлів на сторонах:
| Критичне значення стискаючих зусиль σкрорт, кН/м2 | |||||
| Розрахункова схема | Форма втрати стійкості | Кількість напівхвиль у поперечному n та поздовжньому m напрямках | Теорія | LIRA-FEM | Похибка, % |
| a/b=1,0 | 1 | 1; 1 | 283093 | 2,77828*105=277828 | 1,86 |
| 2 | 1; 2 | 642810 | 6,27352*105=627352 | 2,40 | |
| 3 | 2; 2 | 1132373 | 10,54700*105=1054700 | 6,86 | |
| a/b=4,0 | 1 | 1; 3 | 264196 | 2,60524*105=260524 | 1,39 |
| 2 | 1; 4 | 283093 | 2,77841*105=277841 | 1,86 | |
| 3 | 1; 2 | 334385 | 3,31868*105=331868 | 0,75 | |
| Критичне значення стискаючих зусиль σкр, кН/м2 | |||||
| Розрахункова схема | Форма втрати стійкості | Кількість напівхвиль у поперечному n та поздовжньому m напрямках | Теорія | LIRA-FEM | Похибка, % |
| a/b=1,0 | 1 | 1; 1 | 235900 | 277828*0,869565=241590 | 2,41 |
| 2 | 1; 2 | 533934 | 627352*0,869565=545523 | 2,17 | |
| 3 | 2; 2 | 942681 | 1054700*0,869565=917130 | 2,71 | |
| a/b=4,0 | 1 | 1; 3 | 220165 | 260524*0,869565=226543 | 2,89 |
| 2 | 1; 4 | 235900 | 277841*0,869565=241601 | 2,42 | |
| 3 | 1; 2 | 278652 | 331868*0,869565=288581 | 3,56 | |
З використанням додаткових вузлів на сторонах:
| Критичне значення стискаючих зусиль σкрорт, кН/м2 | |||||
| Розрахункова схема | Форма втрати стійкості | Кількість напівхвиль у поперечному n та поздовжньому m напрямках | Теорія | LIRA-FEM | Похибка, % |
| a/b=1,0 | 1 | 1; 1 | 283093 | 2,83157*105=283157 | 0,02 |
| 2 | 1; 2 | 642810 | 6,42733*105=642733 | 0,01 | |
| 3 | 2; 2 | 1132373 | 11,34010*105=1134010 | 0,15 | |
| a/b=4,0 | 1 | 1; 3 | 264196 | 2,64266*105=264266 | 0,03 |
| 2 | 1; 4 | 283093 | 2,83181*105=283181 | 0,03 | |
| 3 | 1; 2 | 334385 | 3,34456*105=334456 | 0,02 | |
| Критичне значення стискаючих зусиль σкр, кН/м2 | |||||
| Розрахункова схема | Форма втрати стійкості | Кількість напівхвиль у поперечному n та поздовжньому m напрямках | Теорія | LIRA-FEM | Похибка, % |
| a/b=1,0 | 1 | 1; 1 | 235900 | 283157*0,869565=246223 | 4,38 |
| 2 | 1; 2 | 533934 | 642733*0,869565=558898 | 4,68 | |
| 3 | 2; 2 | 942681 | 1134010*0,869565=986095 | 4,60 | |
| a/b=4,0 | 1 | 1; 3 | 220165 | 264266*0,869565=229796 | 4,38 |
| 2 | 1; 4 | 235900 | 283181*0,869565=246244 | 4,39 | |
| 3 | 1; 2 | 278652 | 334456*0,869565=290831 | 4,37 | |
Помилка в тексті? Виділіть її та натисніть Ctrl + Enter, щоб повідомити нам.
Коментарі