ОПРЕДЕЛЕНИЕ МАКСИМАЛЬНОГО ПРОГИБА ТРЕУГОЛЬНЫХ ПЛАСТИНОК НА УПРУГОМ ОСНОВАНИИ
С КОМБИНИРОВАННЫМИ ГРАНИЧНЫМИ УСЛОВИЯМИ
В.И. КОРОБКО, доктор технических наук, профессор, А.А. АКТУГАНОВ, студент
ФБГОУ ВПО «Государственный университет - учебно- научно- производственный комплекс», г. Орел
Задача поперечного изгиба треугольных пластинок с комбинированными граничными условиями (все возможные комбинации шарнирного опирания и жесткого защемления по сторонам) на упругом основании, нагруженных равномерно распределенной по всей площади нагрузкой решается методом интерполяции по коэффициенту формы (МИКФ). Для каждой комбинации граничных условий построены аппроксимирующие функции для подсчета двух вспомогательных коэффициентов, необходимых для выбора опорных решений при использовании методики МИКФ.
КЛЮЧЕВЫЕ СЛОВА: треугольные пластинки, комбинированные граничные условия, максимальный прогиб, упругое основание, коэффициент постели, коэффициент формы, метод интерполяции по коэффициенту формы.
Пластинки на упругом основании широко используются в строительстве. К ним относятся фундаменты зданий и сооружений, настилы, покрытия автомобильных дорог и взлетно-посадочных полос аэродромов, днища резервуаров и емкостей. Упругое основание оказывает существенное влияние на напряженно-деформированное состояние всего сооружения.
Расчет пластинок сложного вида с различными комбинациями граничных условий является весьма трудной задачей, которая решается в настоящее время в основном с помощью метода конечных элементов (МКЭ). Однако по-прежнему в строительной механике актуальной остается проблема развития аналитических методов расчета таких пластинок, позволяющих с помощью достаточно простых аналитических зависимостей проводить анализ напряженно-деформированного состояния пластинок на упругом основании. Аналитические решения задачи по определению максимального прогиба треугольных пластинок на упругом основании в научной литературе не приводятся.
Метод интерполяции по коэффициенту формы [1] до настоящего времени к расчету пластинок на упругом основании не применялся. Однако в статье [2] приводятся исследования, позволяющие сделать вывод о возможности и перспективности его применения к таким задачам. Используя этот метод, можно получать удовлетворительные оценки интегральных физических характеристик пластинок (таких, как максимальный прогиб) без решения сложных дифференциальных уравнений и использования мощных ЭВМ.
Для расчета пластинок сложных форм с помощью МИКФ необходимо предварительно получить граничные кривые «максимальный прогиб - коэффициент формы» пластинки, на которых будут лежать так называемые «опорные» решения. Подробные сведения о коэффициенте формы и методике МИКФ приводятся в монографии [1]. Для четырехугольных пластинок (ив частности, тра-пециидальных) одну из граничных кривых образуют решения, полученные для треугольных пластинок. Построению этой границы для данного типа пластинок с различными комбинациями граничных условий (рис. 1) и посвящена предлагаемая статья.
--/7777777777777777777 77777777777777777777
Рис. 1. Возможные комбинации граничных условий для треугольной пластинки (жесткое защемление - шарнирное опирание)
В статье [2] при исследовании вопроса о возможности использования МИКФ для определения максимального прогиба пластинок на упругом основании была получена формула:
1 _ 1 Ч
у0
ВВ+кС А2
(1)
где ч - интенсивность равномерно распределенной нагрузки, к - коэффициент постели; А - площадь пластинки; В - её цилиндрическая жесткость;
В _ Ф& • К} 4\%(р)сСр и с _ 2\§2(р)сСр |%(р)сСр - коэффициенты пропорциональ-10 0 10 ности, зависящие от формы пластинки и вида граничных условий; К/ - коэффициент формы пластинки; ф^ _ | ^"2р + 2g"g' + g'2р^р, g(P) - приближенная
0
функция прогибов пластинки, представленная в однопараметрическом виде; р -г/г(ф) (г - полярная координата произвольной точки пластинки, г(ф) - уравнение контура пластинки в полярных координатах с центром системы координат в точке с максимальным прогибом).
Коэффициенты В и С могут быть найдены, если известен максимальный прогиб пластинки при двух значениях коэффициента постели, причём один из
них может быть равен нулю, то есть соответствовать решению для пластинки без упругого основания, а второй может иметь произвольное значение, не равное нулю:
1
1
V0 (k=0)
D
qA2
B ^ B
1
qA
V0 (k=0)
D
w
0(k *0)
D
qA2
• B+k • С q
w
V0{k * 0) w0 k = 0)
- С ^ С = 1
J_
0 (k=0)
л
+ k • c, q
w0(k * 0) w0 (k = 0)
(2)
(3)
Значения параметров В и С не зависят от площади, коэффициента постели, цилиндрической жесткости и интенсивности равномерно распределенной нагрузки, их величина определяется лишь значением коэффициента формы пластинки и видом граничных условий. Величины w0(k ф 0) и w0(k = 0) могут быть найдены с помощью приближенных методов строительной механики.
Для определения параметров В и С были рассчитаны пластинки на упругом основании с использованием МКЭ в программном комплексе «SCAD». Расчет производился при следующих параметрах: материал пластинок - сталь (Е = 2,1-10 Па, коэффициент Пуассона и = 0,3), толщина h = 15 мм, А = 1 м , интенсивность равномерно распределенной нагрузки q = 40 кН/м2, число конечных элементов 225. Значение максимального прогиба w0(k ф 0) определялось при коэффициенте постели k = 20 МПа/м. По данным расчёта были построены кривые С - Kf (рис. 2) и В - Kf (рис. 3).
Анализ рис. 2 показывает, что несмотря на относительную простоту и универсальность предложенного подхода, при подборе аппроксимирующих функций С - Kf возникают трудности применения коэффициента С при использовании методики МИКФ из-за сложного вида этих функций и неоднозначности их поведения при некоторых видах граничных условий. Указанные трудности можно преодолеть, если вместо С ввести коэффициент С :
c=C=
B
D
kA2
V0 (k=0)
-1
v0(k *0)
Тогда формула (1) примет следующий вид:
где
± = 1В
v0
q
D+kc
A2
С = 8jg2(p)dp Ф& • Kf.
(4)
(5)
Полученные решения при определении коэффициентов В и С, сведены в таблицу, в которой приведено их сопоставление с результатами, найденными по аппроксимирующим функциям, построенным по этим численным решениям.
Используя программу ТаЫеСигуе 2D, по данным численного расчёта, построены аппроксимирующие функции В - К (рис. 3) и С - К (рис. 4). Как видно из рис. 4, замена параметра С на С упорядочили кривые С - К, которые стали соответствовать логике МИКФ.
Пластинки в виде равнобедренных треугольников с углом при вершине а > 60о При 1-й комбинации граничных условий:
B
(a+b• Kff
где a = -0,2138; b = 0,03458;
(6)
1
1
1
1
k
Таблица
а 20 40 60 80 90 100 120 140 160
К/ 16,07 11,20 10,39 10,96 11,66 12,67 16,09 23,42 46,10
Значения параметра В-10-3
«1» МКЭ 0,640 0,368 0,326 0,355 0,390 0,439 0,610 1,08 3,30
(6) 0,640 0,368 0,329 0,352 0,381 0,425 0,593 1,05 3,27
д, % 0 0 0,92 -0,85 -2,31 -3,19 -2,79 -2,78 -0,91
«2» МКЭ 1,22 0,611 0,520 0,536 0,580 0,642 0,893 1,59 4,95
(8) 1,22 0,607 0,512 0,541 0,580 0,642 0,892 1,59 4,95
д, % 0 -0,66 -1,54 0,93 0 0 -0,11 0 0
«3» МКЭ 0,846 0,517 0,520 0,602 0,678 0,790 1,20 2,32 7,45
(10) 0,846 0,517 0,520 0,597 0,676 0,792 1,21 2,32 7,45
д, % 0 0 0 -0,83 -0,30 0,25 0,83 0 0
«4» МКЭ 2,17 0,975 0,771 0,784 0,835 0,929 1,32 2,42 7,94
(12) 2,17 0,963 0,760 0,795 0,845 0,931 1,31 2,43 7,94
д, % 0 -1,23 -1,43 1,15 1,20 0,22 -0,76 0,41 0
«5» МКЭ 1,54 0,854 0,756 0,872 0,979 1,14 1,71 3,04 9,28
(14) 1,54 0,846 0,770 0,866 0,981 1,15 1,71 3,04 9,28
д, % 0 -0,94 1,85 -0,69 0,20 0,88 0 0 0
«6» МКЭ 2,63 1,32 1,11 1,27 1,41 1,62 2,35 4,47 15,38
(16) 2,63 1,33 1,13 1,26 1,41 1,62 2,36 4,48 15,38
д, % 0 0,76 1,80 -0,79 0 0 0,43 0,22 0
Значения параметра С
«1» МКЭ 1,031 1,781 1,981 1,844 1,694 1,512 1,090 0,627 0,209
(7) 1,028 1,782 1,980 1,843 1,697 1,515 1,087 0,628 0,207
д, % -0,29 0,06 -0,05 -0,05 0,18 0,20 -0,28 0,16 -0,96
«2» МКЭ 0,548 1,064 1,255 1,215 1,131 1,021 0,741 0,428 0,141
(9) 0,545 1,073 1,275 1,204 1,124 1,018 0,746 0,427 0,141
д, % -0,55 0,85 0,16 -0,91 -0,62 -0,29 0,68 -0,23 0
«3» МКЭ 0,788 1,247 1,255 1,079 0,958 0,827 0,558 0,303 0,094
(11) 0,788 1,246 1,255 1,072 0,955 0,827 0,558 0,302 0,094
д, % 0 -0,08 0 -0,65 -0,31 0 0 -0,33 0
«4» МКЭ 0,311 0,669 0,863 0,824 0,775 0,702 0,508 0,287 0,090
(13) 0,311 0,669 0,861 0,818 0,767 0,698 0,512 0,286 0,090
д, % 0 0 -0,23 -0,73 -1,03 -0,57 0,79 -0,35 0
«5» МКЭ 0,428 0,768 0,863 0,741 0,667 0,582 0,400 0,224 0,075
(15) 0,429 0,768 0,856 0,759 0,681 0,591 0,397 0,218 0,076
д, % 0,23 0 -0,81 2,43 2,10 1,55 -0,75 -0,27 1,33
«6» МКЭ 0,249 0,495 0,588 0,521 0,471 0,416 0,286 0,153 0,045
(17) 0,252 0,495 0,586 0,520 0,470 0,415 0,286 0,153 0,045
д, % 1,21 0 -0,34 -0,19 -0,21 -0,24 0 0 0
Примечание - В таблице для сокращения объема статьи приведены лишь часть исследуемых данных.
— Ь С
С = а+-+-, где а = -0,00005959; Ь = 0,09728; с = -0,09334; (7)
К/ К/
при 2-й комбинации граничных условий:
В = а+Ь-К/5 + с• Кг + с• К/+е- К/, (8)
где а = 1,6699; Ь = -1,0722; с = 0,2285; С = -0,0009256; е = 2,036-10"5;
а)
0,000675 0,000672 0,000669 0,000666 0,000663 0,00066 0,000657 0,000654 0,000651 0,000648 0,000645 0,000642
10 10,5 11 11,5 12 12,5 13 13,5 14 14,5 15 15,5 16 16,5
С
0,00072 0,00071 0,0007 0,00069 0,00068 0,00067 0,00066 0,00065 0,00064
Рис. 2. Кривые С - К/ для пластинок в виде равнобедренных треугольников на упругом основании с различными комбинациями граничных условий (1.. .6 - номера схем граничных условий по рис. 1): а) пластинки с углом при вершине а < 60°, б) то же при а > 60°
С
1
а+Ь ■ К22+с ■ К2/5
, где а = 352,0919; Ь = 4,748; с = -0,2319; (9)
при 3-й комбинации граничных условий:
В = а+Ь Кг+& К/ + ё ■ К/ + е^ К/ +/ ■ К/, (10)
где а = -1,6227; Ь = 0,4039; с = -0,03214; ё = 0,001604; е = -0,0000349; / = 0,2845-10"6;
С
В
2.7
2.4 2,1
1.8
1.5 1,2 0,9 0,6 0,3
В
10 10,5 11 11,5 12 12,5 13 13,5 14 14,5 15 15,5 16 16,5 Кг
б)
. Кг = 10.392
Кг
Рис. 3. Кривые В - К для пластинок в виде равнобедренных треугольников на упругом основании с различными комбинациями граничных условий (1.6 - номера схем граничных условий по рис. 1): а) пластинки с углом при вершине а < 60°, б) то же при а > 60°
С
а+Ь К/+ с К
1 + & Кг+е^ К/ + /■ КЗ
(11)
0
Рис. 4. Кривые С - К для пластинок в виде равнобедренных треугольников на упругом основании с комбинациями граничных условий по рисунку 1: а) пластинки с углом при вершине а < 60°, б) тоже при а > 60°
где а = -0,003647; Ь = 0,0003714; с = -0,1241-10"5; ё = -0,3929; е = 0,01942; / = 0,0009338;
при 4-й комбинации граничных условий:
В = а + Ь ■ К/ + с ■ К1/ + ё ■ К2/ + е/Кг , (12)
где а = -2,1112; Ь = 0,2022; с = -0,01405; С = 0,0003382; е = 11,674;
_ а+ Ь• Кг
С=-(13)
1 + с• К+С- К2
где а = 0,001207; Ь = 0,115710"4; с = -0,04797; С = 0,009631; при 5-й комбинации граничных условий:
В = а+Ь-К/+с- К2 + С • К3 + е- К4 + / • К5, (14)
где а = -1,587; Ь = 0,334; с = -0,01602; С = 0,0006613; е = -0,115710"4; / = 0,8162-10-7;
С =-1—тт, где а = -225,79; Ь = 42,529; (15)
а+Ь-К1/5
при 6-й комбинации граничных условий:
В = а+Ь-К/5 + с- К+С • К/5 + е• К/ + / • К2/, (16)
где а = -76,2907; Ь = 82,8033; с = -35,3718; С = 7,4666; е = -0,7699; / = 0,03182; - а+Ь-К+с- К2 + С - К/
С=---/2--, (17)
1 + е - К/ + / - K2/+g - К/
где а = 0,001042; Ь = -0,758-10"4; с = -0,3242-10"5; С = 0,2811-10"7; е = -0,1649;/= 0,02632; g = 0,001994;
Пластинки в виде равнобедренных треугольников с углом при вершине а < 60о При 1-й комбинации граничных условий:
В = а+Ь- К1/5, где а = -0,01151; Ь = 0,01012; (18)
V
К, '
С=а+— где а = -0,0007089; Ь = 0,0279; (19)
при 2-й комбинации граничных условий:
В =-1—~, где а = 0,04123; Ь = 201,695; (20)
а+Ь- К-2
С = а+—, где а = 0,4621-10"4; Ь = 0,1288; (21)
К2
при 3-й комбинации граничных условий:
а+Ь-К+ с- К2
В =-/-(22)
1 + С- К/+ е- К}
где а = -0,3028; Ь = 0,09734; с = -0,00681; С = -0,1055; е = 0,0004201;
— а+Ь- К+с- К2
С=-/-2, (23)
1 + С- К/+е- К} У '
где а = 0,004796; Ь = -0,0006176; с = 0,1422-10"4; С = 0,03339; е = -0,01311; при 4-й комбинации граничных условий:
В = а+Ь-К ^, где а = -0,1741; Ь = 0,009063; (24)
_ а+Ь-К2 + с- К4
С=-2-т, (25)
1 + С- К2+е- К4
где а = -0,01690; Ь = 0,0001607; с = -0,2171-10"7; ё = -0,2166; е = 0,00194; при 5-й комбинации граничных условий:
В =-1--, где а = -2,0573; Ь = 10,8419; (26)
а+Ь К/'5
С=а+— где а = -0,0003511; Ь = 0,01254; (27)
К/
при 6-й комбинации граничных условий:
V
В = а+Ь К1;5, где а = -0,4925; Ь = 0,04853; (28)
С =-1—гт, где а = -680,474; Ь = 72,072. (29)
а+Ь■ Ку
Как видно из таблицы, погрешность результатов расчёта, полученных с помощью функций (6).. .(29), не превышает 3 %.
Несколько сложный вид аппроксимирующих функций (6).. .(29) и значения коэффициентов а и Ь до четырех значащих цифр после запятой объясняются тем, что нами ведётся работа по составлению программного комплекса для решения рассматриваемых задач с помощью МИКФ.
Таким образом, зная коэффициент формы заданной треугольной пластинки, можно, используя аппроксимирующие функции В - К/ и С - К/, соответствующие заданным граничным условиям, найти значения параметров В и С , и, воспользовавшись формулой (5), определить её максимальной прогиб.
Построенные кривые (6).. .(29) образуют верхнюю границу максимального прогиба для всего множества четырехугольных пластинок с выпуклым контуром и комбинированными граничными условиями. Нижнюю границу для таких пластинок образуют решения, соответствующие прямоугольным пластинкам. Задача по построению этой границы нами решена в статье [3]. Используя эти обе границы, можно определять прогиб любой четырехугольной пластинки, применяя методику МИКФ [1].
Выводы:
1. Для пластинок в виде равнобедренных треугольников, лежащих на упругом основании, с различными комбинациями граничных условий жесткого защемления и шарнирного опирания по сторонам и нагруженных равномерно распределённой нагрузкой, построены аппроксимирующие функции В - К/ и С - К/, входящие в выражение (1), для определения максимального прогиба треугольных пластинок.
2. Все множество значений параметров В и С для всех комбинаций граничных условий, представленное в координатах В - К/ и С - К/ при любых значениях коэффициента постели ограничено с двух сторон соответствующими значениями этих параметров для пластинок с однородными граничными условиями (либо шарнирное опирание по контуру, либо жесткое защемление).
3. Полученные аппроксимирующие функции могут использоваться для непосредственного расчета треугольных пластинок с комбинированными граничными условиями, лежащих на упругом основании, а также в качестве граничных кривых при расчете четырёхугольных пластинок с помощью МИКФ.
Л и т е р а т у р а
1. Коробко, А.В. Геометрическое моделирование формы области в двумерных задачах теории упругости [Текст] / А.В. Коробко. - М.: Изд-во АСВ, 1999. - 320 с.
2. Коробко, В.И. Расчет пластинок на упругом основании методом интерполяции по коэффициенту формы [Текст] / В.И. Коробко, А.А. Актуганов. // Фундаментальные исследования РААСН по научному обеспечению развития архитектуры, градостроительства и строительной отрасли Российской Федерации в 2010 году. - Москва-Орел, 2011. - С. 261-267.
3. Коробко, В.И. Определение максимального прогиба прямоугольных пластинок на упругом основании с комбинированными граничными условиями [Текст]/ В.И. Коробко, А.А. Актуганов// Строительство и реконструкция, №4, 2011. - Орел, 2011. - С. 15-20.
DETERMINATION OF MAXIMAL BENDING OF THREE-CORNERED PLATES ON RESILIENT FOUNDATION WITH THE COMBINED SCOPE TERMS
V.I. Korobko, A.A. Aktuganov Task of cross-bending of three-cornered plates with the combined scope terms (all possible combinations of the joint leaning and gesture-who jammings on sides) on resilient foundation, loaded with the equipartition on all area loading decides the method of interpolation on the coefficient of form (MIKF). For every combination of scope terms approximating functions are built for the count of two auxiliary coefficients, necessary for the choice of supporting decisions at the use of method of MIKF.
KEYWORDS: three-cornered plates, combined scope terms, maximal bending, resilient foundation, coefficient of bed, coefficient of form, method of interpolation on the coefficient of form.
-0- -0- Hh