Файл: Синицын А.П. Расчет балок и плит на упругом основании за пределом упругости пособие для проектировщиков.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 27.06.2024
Просмотров: 97
Скачиваний: 0
№
блоксхемы
1
4
Программа на языке АЛГОЛ-60 для выполнения вычислений на Э ВМ М-220а
JÖ |
|
|
|
|
Операции и команды |
|
|
|||
і.п. |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
1 |
’BEGIN’ REAL’L l, Н, Т, Е, |
DL, |
DH, Р1, |
N10, |
C l, H l, L2, N1, |
|||||
2 |
N2, N5, N6; |
|
|
|
|
|
|
|
||
3 |
’INTEGER’ I, J, К, N, F, Р, М, |
L, N3, |
N4; |
|
||||||
4 |
Р 0 0 4 2 |
(L2, HI, Т, Е, DL, DH, PI, NIO, N1, N2, N5, N6); |
||||||||
5 |
Р1041 |
(L2, |
HI, Т, Е, DL, DH, PI, N10, |
N1, N2, |
N5, N6); |
|||||
6 |
N3 : = N 5 + 0 ; N4 : = N 6 + 0 ; |
|
|
|
|
|||||
7 |
’BEGIN’ ’INTEGER' ’ARRAY’ А [I : 6, 1 : N31; |
|
||||||||
8 |
’BEGIN’ ’ARRAY’ R [1 : 6, 1 : N3]; P 0 0 4 2 |
(R ); |
|
|||||||
9 |
’FOR’ J: = |
1’STEP’ l ’UNTIL’ 6’ DO’ BEGIN’ |
|
|||||||
1 0 |
'FOR’K: = |
1’STEP’ l ’UNTIL’ N3’ DO’ |
|
|
|
|||||
11 |
A[J, K] : = R [J , K ] + 0 ’END’ ’END’; |
|
|
|
||||||
1 2 |
P1041 |
(A); |
|
|
|
|
|
|
|
|
13 |
’BEGIN’ ’ARRAY’ B6, |
B4 [1 :3 , |
1 :6 ], |
B5 [1 :6 , 1 :3 ], |
||||||
14 |
S[1 ; 3, |
1 : I], |
D [ 1 : 3, |
1 :3 ], |
U |
[1 : N4, |
1 |
: 1], |
|
|
15 |
С [1 : 6, 1 : 6], |
В [1 : N 4+1, |
1 : N 4+ 1], |
U l |
[1 : 6, 1 : 1]; |
|||||
16 |
’F 0R ’H: = |
H1’STEP’DH’UNTIL’ |
10XH1 |
’DO’ |
|
|||||
17 |
’BEGIN’ P1041 |
(H); |
|
|
|
|
|
|
18 ’FOR’ L 1:= L 2’STEP’ DL ’UNTIL’ 10XL2 ’DO’
19’BEGIN’ P1041 (LI)
20’FOR’ I : = 1 ’STEP’ 1’ UNTIL’ 3’ DO’ ’BEGIN’
21 |
’FOR’ J: = |
1 ’STEP’ 1’ UNTIL’ 6’ DO’ |
|||
22 |
B4 |
[I, J] |
= |
0 ’END’ |
|
23 |
B4 |
[1, 1] |
= |
—H /Nl; B4 [1, 5]: = H/N1; |
|
24 |
“ '.LI. " |
= |
L1/N2; B4[2, 4]: = |
—L1/N2; |
|
B4[2, 2] |
|||||
25 |
B4[3, 1] |
=L 1/N 2; B4[3, 2]: = |
—H /Nl; |
||
26 |
B4[3, 3] |
= |
—L1/N2; B 4[3, 6]: = H/N1; |
|
27 |
’FOR’ I: = |
1 |
’STEP’ |
Г UNTIL’3’D 0 ’ ’BEGIN’ |
||||
5 |
28 |
’FOR’ J: = |
1 ’STEP’ 1 ’UNTIL’ 6 ’DO’ |
|
|
||||
|
29 |
B5 [J, I]: = |
|
B4 [I, J] |
’END’; |
|
|
||
|
30 |
D [1- 1 |
= |
D[2, 2]: = |
1; |
|
|
|
|
6 |
31 |
D T, 2 |
= |
D[2, 1]: = |
NK); |
|
|
|
|
32 |
D '3, 3] |
= |
(1— NI0)/2; |
|
|
|
|||
|
33 |
D .1. 3] |
= D [2 , 3 ]:= D [3 , |
1 ]:= D [3 , |
2]: = |
0; |
|||
7 |
34 |
P I067 (3, 6, 3, D, B4, B6) ; |
|
|
|
||||
35 |
P1067 (3, 6, 6, B5, B6, С); |
P1024 (1, C); |
|
||||||
|
|
||||||||
|
36 |
’FORT: = |
1 |
’STEP’ |
1 ’UNTIL’ N 4+1 |
’DO’ |
’BEGIN’ |
||
|
37 |
’FOR’ J: = |
|
1 ’STEP’ 1 ’UNTIL’ N 4+1 |
’DO’ |
|
|||
8 |
38 |
В [I, J]: = 0 |
’END’; |
|
|
|
|
||
39 |
’FOR’ N: = |
l ’STEP’ 1 ’UNTIL’ N3 ’DO’ |
|
||||||
|
|
||||||||
|
40 |
’BEGIN’ |
’FOR’ K: = |
l ’STEP' 1 ’UNTIL’ 6 ’DO’ |
|||||
|
41 |
’BEGIN’ TF’ A[IC, N] = 0 |
’THEN’ ’GO TO’ FI; |
132
К* |
№ |
|
блок- |
п.п. |
|
схемы |
||
|
42
43
44
g45
46
47
48
49
950
51
10 52
Продолжение
|
|
Операции и комапды |
М := А [К , N1; F: = |
K; |
|
’FOR’ I: = |
1 ’ST EP’ 1 ’UNTIL’ 6 ’DO’ |
|
’BEG IN ' |
T F’ A [I, |
N 1 = 0 ’THEN’ ’GO TO’ F2; |
P : = I; L := A [I , N ]; |
В [M, L]: = B[M , L ]4 -C [F , P ]; F 2 : ’END’;
F I: ’END’ ’END’
B [l, N 4 4 - 1 ]:= —P1X H X L 1 (1— N K 0f2)/E/T/N l/N 2; P I024(2, B);
P I052 (N44-1, N4, B);
53 |
’FOR’ I: = |
1 ’ST EP’ 1 ’UNTIL’ N4 ’DO’ |
54 |
U [I, 1]: = |
B [I, N44-1]; P1041 (U ); |
55
56
57
58
59
60
1261
62
63
64
65
66
67
13 68
14 69
70
71 ю 72
73
C1: = E X N 1X N 2/H /L 1/(1—N ro f2 );
’FOR’ I: = |
l ’ST E P’ 1 ’UNTIL’ 3 ’DO’ ’BEGIN’ |
’FOR’ J: = |
l ’ST EP’ 1 ’UNTIL’ 6 ’DO’ |
B 6[I, J]: = |
B 6[I, J]X C 1, ’END’ |
’FOR’ N: = |
l ’ST EP’ 1 ’UNTIL’ N3 ’DO’ |
’BEG IN ' ’FO R’ K : = l ’ST EP’ 1 ’UNTIL’ 6 ’DO’
’B EG IN ’ TF’ A[K, N] = 0 |
’THEN’ ’BEGIN’ |
|
U1[K. 1] := 0; ’GO TO’ F3 ’END’ ’ELSE ’ |
||
I:= A [K , N ]; U 1[K, 1]: = |
U [I, 1]; |
|
F3 : ’END’; |
|
|
P I 067 |
(6, 1, 3, B6, U l, S ); |
|
PI 041 |
(S) |
|
’END’
’END ’
’END’
’END’ ’EN D ’ ’EN D ’
Значения вертикальных перемещений ѵ* |
и горизонтальных перемещений гр |
узла |
і: |
Оз |
= 4 ------- |
03 |
46409328 |
о4 |
= 4 ------- |
02 11348957; |
||
о, |
= |
4------- |
03 |
|
51850627 |
И9 |
= + 4 — оз 23552335; |
|
«14 = н —і— 03 |
|
13149515 |
U,5= + 4 — 03 22402670; |
|||||
Од |
= |
4------- |
02 |
|
18544202 |
«г |
= 4 ------- |
03 20516689; |
и6 = + Н — 03 |
|
13348722 |
о 6 |
= 4 ---------- |
03 45603356; |
|||
цп^ + Н — 03 |
13763200 |
i> u = 4 -'----- |
03 25825119; |
|||||
и5 = |
4 4 — 03 |
|
14231248 |
о5 = 4------- |
03 76149626; |
|||
о9 |
= |
4-------- |
03 |
|
39171004 |
о,о = 4 ------- |
03 31158627; |
|
о2 |
= |
4------- |
03 |
|
84074321 |
«з |
= 4 ---------- |
03 21322661; |
о7 |
= - | ------- |
03 |
|
66449605 |
и7 = 4-4— 03 16116591; |
|||
Ui2 = -} ------- |
03 |
|
2255893І |
о,2 = 4 ------- |
03 21679617. |
133
Выданные на печать цифры расшифровываются следующим образом:
"I- |
|
— |
|
|
— |
|
02 |
|
18544202 |
|
t |
|
|
t |
|
t |
числа |
I |
|
|
знак числа |
знак порядка |
порядок |
мантисса |
|||||
следовательно, получим: ѵ8 = —0,00185 |
см. |
|
|
|
|||||
|
Полученные в кгс/см2значения напряжений Oxj, o Vj, x XVj , |
||||||||
|
|
|
отнесенные к центру тяжести элемента / |
|
|||||
№ |
|
|
|
|
|
|
V |
|
Ѵѵу |
элемента |
|
|
а ѵ/ |
|
|
|
|
||
1 |
4------- |
1-01 64654843 |
|
-)------- |
[-02 59501609 |
+ + + 0 2 |
11489213 |
||
2 |
4------- |
[-02 20126679 |
|
-[------- |
[-02 12377760 |
4— 1— [-02 20498390 |
|||
3 |
4------- |
1-01 28658766 |
|
-1------- |
[-01 71995192 |
4— + 0 1 |
12928208 |
||
4 |
-1------- |
[-00 92111070 |
|
-[------- |
[-00 92111070 |
+ + + 0 0 |
92111070 |
||
5 |
+ + + 0 1 |
17622867 |
|
-[------- |
1-02 37103287 |
+ + + 0 1 |
44919932 |
||
6 |
+ + + 0 1 |
61024601 |
|
-1------- |
[-02 17608461 |
+ + + 0 2 |
10909108 |
||
7 |
+ 4 — |
00 12839087 |
■ |
+ — |
[-02 19400682 |
+ + + 0 1 |
33213167 |
||
8 |
-1------- |
[-01 24722915 |
|
4------- |
[-01 58875687 |
+ + + 0 1 |
62592789 |
||
9 |
+ + + 0 1 |
27920060 |
|
-1------- |
[-02 27395180 |
+ + + 0 1 |
16410793 |
||
10 |
+ + + 0 1 |
73045615 |
|
-1------- |
[-02 18629038 |
+ + + 0 1 |
52161135 |
||
11 |
+ + + 0 0 |
86864631 |
|
4------- |
1-02 20559812 |
+ + + 0 1 |
19996592 |
||
12 |
+ + + 0 1 |
19238049 |
|
Н------- |
[-02 13415968 |
+ + + 0 1 |
42070830 |
||
13 |
+ + + 0 1 |
33423533 |
|
4------- |
1-02 23924195 |
+ - — 00 00000000 |
|||
14 |
+ + + 0 1 |
52884964 |
|
4------- |
[-02 19073546 |
+ + + 0 1 |
18299056 |
||
15 |
+ + + 0 1 |
13236214 |
|
4------- |
[-02 20263008 |
+ 4 — 00 18882625 |
6.3. Определение пластических областей —в основании
Метод конечных элементов с успехом применяется для ре шения задач об исследовании упругопластического напряжен ного состояния различных систем [62, 45]. Используем этот ме тод для изучения образования пластических областей в упру гом основании, представляющем слой конечной толщины.
Для треугольной схемы сетки напряжения внутри каждого треугольника, примыкающего к данному узлу, различны, по этому за предел упругости перейдут сначала отдельные тре угольники, другие же останутся в упругой стадии. В п. 1.4 бы ли приведены общие формулы [46], позволяющие определить условия перехода за предел упругости участков упругого осно вания, которое в общем виде рассматривается как упругопла стический сжимаемый слой конечной толщины. Однако для уп рощения этой новой задачи сначала будет использован обыч ный критерий текучести:
^2 — 26,
где схі и 0 2 — главные напряжения; /г — предельное касательное напряжение.
134
Вычислим главные напряжения для единичной силы Д0= = 20 тс, приложенной на поверхности основания. Схема на грузки и полученные эпюры напряжений показаны на рис. 6 .5 — 6.7. Расчет ведется на 1 пог. см в направлении, перпендикуляр ном чертежу. Для подсчета напряжений используем формулы плоской задачи. Для треугольника 145 (см. рис. 6.5) получим:
стп = 0 Г.Ѵ+ О,, |
7 7 h - » , |
) ’ + 4 ; |
°2 і = |
+ і / К - ' М |
’+ Ч , - |
Определяем разницу <ji— 0 2 и сравниваем с величиной 2 k, принятой для данного грунтового основания. Заметим, что тре угольник 145 при данной нагрузке первым переходит в пласти ческую стадию. Затем переходит в пластическое состояние тре угольник 125 и т. д. По данным этой схемы определяются гра ницы пластической области, которая возникает в упругом ос новании, однако это будет только первым приближением, так как при определении поля напряжений для всех треугольников принимался одинаковый упругий модуль деформации.
Для пластической области модуль деформации снижается, поэтому при определении реакций в узлах, возникающих от еди ничных смещений, следует это учесть и в формулы (6.4) — (6 .6 )
вместо упругого |
модуля Е подставить модуль |
деформации |
£ P= ß £ , который |
соответствует пластическому |
участку били |
135
нейной диаграммы. Тогда величины реакций k uluj, возникаю
щих в узле і от смещения узла / на единицу, будут в ß раз меньше, поэтому для треугольников, находящихся в пластичес кой стадии, значения реакций будут также в ß раз меньше по сравнению с упругими треугольниками; их можно вычислить по формуле
|
иІ )п л |
= |
ß { К Uj |
)у п р ' |
|
|
|
|
|
|
|
|
|
|
Эп^зрз. 6Х |
||
|
|
|
|
|
|
15 2} |
-7,9 |
0.92 |
|
|
|
|
»* |
|
|
|
|
|
|
|
|
|
|
- ^ 0 |
^ |
-2.2 |
|
|
|
|
|
t |
|
|
- г |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
•4.0 |
•S,0 |
-Vi j |
|
|
|
|
|
J |
|
<■ |
|
|
|
|
|
|
-J.3 |
‘V |
||
|
|
|
|
|
|
|
||
|
|
|
|
<* |
|
*3.3 |
|
< |
|
|
|
|
|
|
*3,3 |
*tß |
|
|
|
|
|
|
Рис. |
6.7 |
|
|
Для нашего примера |
в |
пластическое состояние |
переходят |
|||||
треугольники |
145 и 125, |
пластическая |
область |
заштрихована |
||||
(см. рис. 6.5). |
|
|
|
|
|
|
|
|
6.4. Высокая |
фундаментная |
балка |
|
|
|
|
||
на упругопластическом основании |
|
|
|
|
||||
Образование пластических областей в основании |
приводит |
|||||||
к перераспределению реакций, |
возникающих между |
балкой |
и основанием. Наметим ход решения для этой задачи, исполь зуя результаты, полученные в предыдущих разделах. Но для упрощения решения прикрепим балку к упругому основанию несколькими стержнями, как это показано на рис. 6 .8 , а. Для определения реакций в стержнях составим систему уравнений, приведенную на рис. 6 .8 , б. Сначала получим решение в упру гой стадии. Однако в отличие от задач, решенных в главе 2, коэффициенты матрицы жесткости вычисляем с учетом резуль
татов, полученных в п. 6.2. Так, например, |
от Х2= 1 в формуле |
|||||
(2 .2 ) коэффициент 6 2 2 будет равен: |
|
|
|
|||
§22 = Уг-> + |
^22 = |
— я— — |
F ik |
+ І'ІЬ |
w lk = |
|
J i . |
|
Е ^ с |
|
I К |
|
|
|
, п |
. , |
|
1 |
—Ho |
|
|
= (Fm+hkawik) |
|
IZQCJi • |
|
a= с*Е0п/6(1 — (XQ) EJ.
136
и2 2 представляет собой прогиб высокой консольной балки от силы Хг=1, приложенной на ее конце. Этот прогиб был вычис лен в п. 6.2 и равнялся: оц=0,0117 м при Р о = 20 тс, I— 1 м, /г = 0,25 м; Ь — 0,1 м и £ '= 2 -10 6 тс/м2. Если эту балку рассмат ривать как тонкий стержень и использовать формулы, указан ные в п. 2.2, то получим 0^=0,016 м. Для высокой балки про
гиб будет несколько меньше.
Вычислим отношение этих прогибов и получим переходный
коэффициент Kih, с помощью которого |
можно учесть влияние |
высоты балки на прогиб. Используя |
формулы, полученные |
в главе 2 для тонких балок, |
|
. __0,0117
0,731.
~0,016
Формула для подсчета хшіь. будет теперь записана так:
-- ( ? ) '( * S |
- О** |
|
ш2 2 = |
16 • 0,731 = |
9,596. |
Следует отметить, что |
зависит от положения точки і по |
отношению к точке k. Для всех перемещений, входящих в мат рицу жесткости, эти коэффициенты были вычислены. Опреде ление реакций основания в упругой стадии является первым приближением при решении упругопластической задачи, поэто му можно учесть тот факт, что Хіь. изменяется мало. При перехо де от одной точки к другой этот параметр для всех точек можно считать постоянным, тогда учет влияния высоты
137
■балки сводится к замене |
во всех вычислениях а на |
а'=аК, |
||
и формула для |
вычисления б,-л будет иметь вид: |
|
||
|
|
. |
2 |
|
|
|
1 |
—Йо |
|
■Увеличение а |
приводит к |
уменьшению |
концентрации |
реакции |
к краю балки, |
а уменьшение а — к ее возрастанию. |
|
6.5. Предельная нагрузка балки, определяемая из условия образования пластических областей в основании
Образование пластических областей в основании |
приводит |
к перераспределению реакций, возникающих между |
балкой и |
основанием. Изменение эпюры реакций, как было |
указано |
в главе 2 , в свою очередь, влечет за собой изменение величины предельной нагрузки. Для определения границ пластической области в основании применим способ последовательных при ближений. Первым приближением будет расчет в упругой ста дии. Выполняя этот расчет с использованием конечных элемен тов, получим следующие значения сил Xf.
Силы Х і представлены в безразмерной форме для удобства дальнейших сравнительных подсчетов. Отметим, что концен трация реакций к краю балки получилась меньше, чем для бес конечно жесткой балки. Метод конечных элементов позволяет более точно определить концентрацию реакций к краю балки для высоких массивных фундаментных балок, которые обычно рассматриваются как бесконечно жесткие.
С помощью метода, конечных элементов в результате реше ния системы линейных уравнений определяются перемещения узлов сетки. Гіо этим перемещениям вычисляются напряжения, относящиеся к треугольным элементам. Однако к каждой уз ловой точке сетки примыкает несколько треугольников, поэто му для получения напряжений в узлах сетки приходится ус реднять полученные для треугольников напряжения.
Применяются разные способы усреднения в зависимости от желательной степени точности расчета. Для инженерных рас четов, связанных с грунтовыми основаниями, высокая степень точности вряд ли будет оправдана, так как реальные грунты имеют неоднородные свойства, что снижает теоретическую точ ность. Очень часто оказывается вполне достаточно использовать среднее арифметическое значение напряжений треугольников, сходящихся в узле.
Таким образом, полученное в упругой стадии решение явля ется первым приближением для выявления последовательности образования пластических областей в упругом основании. По
138