3 семестр / Manchev_E_I_Posobie
.pdfЧасть II. Напряженно-деформированное состояние в точке. Тензор напряжений
Задача 5
1.Получить характеристическое уравнение для определения главных
значений (главных напряжений) и системы уравнений для нахождения направлений главных осей тензора напряжений .
2.Для проверки правильности найденных корней характеристического уравнения, вычислить инварианты тензора напряжений в исходных и главных осях.
3.Вычислить главные значения и направляющие косинусы главных осей
jk , МПа
4. По заданному критерию прочности вычислить эквивалентные напряжения.
Исходные данные
|
11 |
МПа |
|
22 |
МПа |
|
33 |
МПа |
|
|
12 |
= |
21 |
МПа |
|
12 |
= |
21 |
МПа |
|
23 |
= |
32 |
МПа |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
-70 |
|
-20 |
|
50 |
|
|
|
-10 |
|
|
|
-10 |
|
|
22 |
|
|||||||
Критерий прочности – Мизеса |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Запишем тензор в матричной форме: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
70 |
10 |
10 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
10 |
20 |
22 |
|
|
|
|
|
|
|
|
|
|
(5.1) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
10 |
22 |
50 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
NB! Главные напряжения – это нормальные напряжения, действующие по площадкам, где отсутствуют касательные напряжения. Координатные оси,
являющиеся нормалями к таким площадкам, называются главными осями тензора напряжений, а сами площадки – главными площадками.
Главные значения (главные напряжения) определяются из уравнения
|
3 |
I1 |
2 |
I |
2 I3 0 |
(5.2) |
|
|
|
|
|
где I1, I2 , I3 это инварианты напряженного состояния (самого тензора).
NB! Инвариант – это величина, которая не зависит от системы координат.
Выражения для нахождения инвариантов:
21
I |
1 |
|
x |
|
|
y |
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
||
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
2 |
2 |
2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
2 |
|
x |
|
y |
|
|
y |
|
|
z |
|
z |
|
x |
|
xy |
|
yz |
|
zx |
|
|
|
|
x |
|
xy |
|
xz |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
I |
3 |
|
yx |
|
|
y |
|
yz |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
zx |
|
zy |
|
z |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5.3)
Условием существования собственных значений матрицы является равенство нулю определителя, составленного из этой матрицы (в нашем случае из матрицы (5.1)), в которой к диагональным элементам прибавляем .
Получится:
70 |
10 |
10 |
|
10 |
20 |
22 |
0 |
10 |
22 |
50 |
|
(5.4)
Раскрывая этот определитель, получим характеристическое уравнение для нахождения собственных значений матрицы (5.1)
|
3 |
( 40) |
|
2 |
( 3784) 105280 0 |
|
|
|
|
Сравнивая (5.5) и (5.2), нетрудно догадаться, что
I1 40 I2 3784 I3 105280
Попробуем теперь найти инварианты по формулам (5.3)
I |
1 |
70 20 50 40 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
I |
|
70 ( 20) ( 20) (50) (50) ( 70) ( 10) |
2 |
(22) |
2 |
( 10) |
2 |
3784 |
|||
2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
70 |
10 |
10 |
|
|
|
|
|
|
|
I |
3 |
10 |
20 |
22 |
105280 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
22 |
50 |
|
|
|
|
|
|
|
(5.5)
Получается то же самое, следовательно, использование формул (5.3) позволяет проверить правильность получения характеристического уравнения (5.5), получаемого при раскрытии определителя (5.4).
Решим это характеристическое уравнение (5.5):
22
(Можно пользоваться программными комплексами, можно решить аналитически с помощью метода Кардано)
Выберем первый вариант. Кубическое уравнение (5.5) имеет 3 корня:
57.547; –72.212; –25.335
Индексация этих трех корней проводится по такому правилу:
|
|
|
1 |
|
max |
|
3 |
|
min |
|
|||
|
|
|
|
|
|
|
|
|
|||||
Значит: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
57.547 |
|
|
2 |
25.335 |
|
3 |
72.212 |
||||
|
|
|
|
|
|
|
|
|
|
|
Известно, что на главных площадках отсутствуют касательные напряжения, поэтому тензор в главных осях будет иметь такой вид:
|
|
57.547 |
0 |
0 |
|
|
|
|
0 |
25.335 |
0 |
|
|
|
|
|||||
|
|
|
|
|||
|
|
0 |
0 |
72.212 |
|
|
|
|
|
||||
|
|
|
|
(5.6)
Выполним проверку правильности найденных главных напряжений. Для этого необходимо найти инварианты тензора (5.6).
Они |
|
|
должны |
|
совпадать с |
|
найденными ранее значениями |
||||||||||||||||||
I |
1 |
|
1 |
|
2 |
|
3 |
40 |
|
|
|
|
|
|
|
|
|
|
|||||||
I |
|
|
|
|
|
|
|
|
|
|
|
(0) |
|
(0) |
|
(0) |
|
3784 |
|||||||
|
|
|
|
|
|
|
|
2 |
2 |
2 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
2 |
|
|
1 |
|
2 |
|
|
1 |
|
|
|
3 |
|
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
3 |
|
0 |
|
2 |
0 |
105280 |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
0 |
|
0 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
1 |
, |
|
|
I |
2 |
, I |
3 |
|
|
:
Совпадают.
Каждому
значению
|
j |
|
( j
1,2,3)
соответствует
вектор
n |
j |
|
направляющих
косинусов, определяющих направление главных площадок по отношению к исходной системе координат:
n |
j |
cos(x, x |
|
); n |
j |
cos( y, y |
|
); n |
j |
cos(z, z |
|
). |
|
x |
j |
y |
j |
z |
j |
||||||||
|
|
|
|
|
|
|
Компоненты этого вектора направляющих косинусов определяются из системы алгебраических уравнений:
23
( |
|
|
|
) n |
j |
|
|
n |
j |
|
|
n |
j |
0 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
x |
|
|
j |
|
x |
|
xy |
|
y |
|
xz |
|
z |
|
|
yx |
nx |
( y |
j ) ny |
yz |
nz |
|
|
|||||||||
0 |
||||||||||||||||
|
|
j |
|
|
|
|
|
|
|
j |
|
|
|
j |
|
|
zx |
nx |
zy |
ny ( z |
j ) nz |
0 |
|
||||||||||
|
||||||||||||||||
|
|
j |
|
|
|
|
j |
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
которые должны быть дополнены условием нормировки:
n |
j |
|
2 |
n |
j |
|
2 |
n |
j |
|
2 |
1 |
где j 1,2,3. |
|
|
|
|||||||||||
x |
|
y |
|
z |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
(5.7)
(5.8)
NB!
n |
j |
, n |
j |
|
x |
y |
|||
|
|
Определитель, составленный из коэффициентов при неизвестных
j |
системы уравнений (5.7), равен нулю. |
, nz |
Следовательно, три уравнения в (5.7) являются линейно зависимые: одно уравнение (любое) является следствием двух других. Поэтому для определения направляющих косинусов nxj , nyj , nzj любой главной оси нужно
любое из уравнений удалить и к двум оставшимся добавить уравнение (5.8).
Решив полученную систему трех уравнений с тремя неизвестными, найдем
j |
j |
j |
, соответствующие главному напряжению |
направляющие косинусы nx |
, ny |
, nz |
. Положение оставшихся двух осей находят аналогично.
|
j |
|
Нужно иметь в виду, что каждый из направляющих косинусов получается с двумя знаками. Знаки соответствуют повороту осей по часовой стрелке или против часовой стрелки. При этом главные оси занимают одно и то же положение, но направлены в противоположные стороны. При определении положения главных осей нужно оставить одну систему знаков, конкретизировав при этом направления осей.
Например:
1 |
1 |
1 |
составим такую систему уравнений: |
Для nx |
, ny |
, nz |
( x 1 ) n1x xy n1y xz n1zyx n1x ( y 1 ) n1y yz n1z
n1x 2 n1y 2 n1z 2 1
Поделим два первых уравнения на n1 и получим z
00
24
|
|
|
|
|
|
n |
|
|
n |
|
|
|
|
( x |
1 ) |
1 |
xy |
|
1 |
xz |
|||||||
|
1 |
|
|
1 |
|||||||||
|
|
|
|
|
|
|
x |
|
|
|
y |
|
|
|
|
|
|
|
|
nz |
|
nz |
|
|
|||
|
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
yx |
|
n |
|
( y |
1 ) |
n |
1 |
|
|
||||
|
1 |
|
|
yz |
|||||||||
|
|
|
|
x |
|
|
|
|
|
|
y |
|
|
nx |
|
nz |
ny |
|
nz |
nz |
|
|
|||||
|
|
|
|
|
1 |
|
|||||||
1 |
|
2 |
|
|
1 |
|
2 |
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставим найденные ранее значения нормальных, касательных и главных напряжений в эту систему и получим:
|
|
|
|
n |
|
|
n |
|
|
|
|
|
|
|
|
|
1 |
|
|
1 |
|
|
|
|
|
( 70 57.547) |
|
x |
10 |
|
y |
|
10 |
|
||||
1 |
|
1 |
|
|||||||||
|
|
|
|
nz |
|
nz |
|
|
|
|
||
|
1 |
|
|
|
|
|
|
1 |
|
|
||
10 |
n |
x |
( 20 57.547) |
n |
y |
|
|
|||||
|
|
|
22 |
|||||||||
1 |
|
|
1 |
|||||||||
|
nz |
|
|
|
|
|
nz |
|
|
nx |
|
|
ny |
|
nz |
|
1 |
|
1 |
|
2 |
1 |
2 |
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Предоставим программному комплексу Mathcad 15 решить эту систему:
|
1 |
|
|
0.097 |
0.097 |
|
||
n |
x |
|
||||||
|
|
|
|
|
|
|
|
|
1 |
|
0.283 |
0.283 |
|||||
|
n |
y |
|
|
|
|||
|
|
|
|
|||||
|
1 |
|
|
|
0.954 |
0.954 |
|
|
n |
z |
|
||||||
|
|
|
|
|
|
|
|
Возьмем значения из 1-го столбца и найдем значения углов. Углы, которые составляет первая главная ось тензора напряжений с исходными осями координат, находятся вычислением функции arccos от n1x , n1y , n1z :
(x, |
0 |
' |
; ( y, |
0 |
' |
) |
) 95 57 |
|
) 73 55 ; (z, |
||||
1 |
|
|
1 |
|
1 |
|
17 |
0 |
42 |
' |
|
|
.
Аналогично для nx2 , ny2 , nz2
( |
x |
|
2 |
) n2 |
|
xy |
n2 |
|
xz |
n2 |
||||
|
|
|
|
|
x |
|
y |
|
z |
|||||
|
yx |
n2 |
( |
y |
|
2 |
) n2 |
|
yz |
n2 |
||||
|
|
x |
|
|
|
|
|
y |
|
z |
nx2 2 ny2 2 nz2 2 1
0
0 и т.д.
Получим:
25
n |
2 |
|
|
|
0.145 |
0.145 |
|
||
x |
|
||||||||
|
|
|
|
|
|
|
|
||
n |
2 |
|
0.944 |
0.944 |
|||||
|
y |
|
|
||||||
|
|
|
|
||||||
n |
2 |
|
|
|
0.295 |
0.295 |
|
||
z |
|
||||||||
|
|
|
|
|
|
|
|
(x, |
|
0 |
' |
; ( y, |
|
) 160 |
0 |
79 |
' |
; (z, |
|
0 |
' |
2 |
) 81 64 |
|
2 |
|
|
2 |
) 72 84 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
3 |
3 |
|
И для nx |
, ny |
, nz |
|||||
Получим: |
|
|
|||||
n |
3 |
|
|
|
0.984 |
||
x |
|
||||||
|
|
|
|
|
|
||
n |
3 |
|
0.167 |
||||
|
y |
|
|
||||
|
|
|
|||||
|
|
3 |
|
|
|
0.050 |
|
nz |
|
|
|
0.984 |
||
0.167 |
|
|
|
||
|
||
0.050 |
|
|
|
||
|
(x, |
|
) 10 |
0 |
07 |
' |
; ( y, |
) 80 |
0 |
37 |
' |
; (z, |
0 ' |
3 |
|
|
|
|
) 87 11 |
|||||||
|
|
|
|
|
3 |
|
|
|
|
3 |
|
Как проверить правильность найденных nxj , nyj , nzj ?
Проверка правильности вычисления положения главных осей тензора напряжений основана на свойствах матрицы направляющих косинусов. Она относится к ортогональным матрицам и обладает следующими свойствами:
–Определитель ортогональной матрицы равен единице.
–Сумма квадратов элементов, входящих в каждую строку (столбец) равна единице.
–Если рассматривать каждую строку матрицы как вектор-строку, а каждый столбец – как вектор-столбец, то скалярные произведения двух разных векторов-строк (векторов-столбцов) равны нулю.
Воспользуемся первым свойством ортогональных матриц.
Составим определитель из
n |
j |
, n |
j |
, n |
j |
|
x |
y |
z |
||||
|
|
|
и проверим, равен ли он единице:
1 |
1 |
|
1 |
|
0.097 |
0.283 |
0.954 |
|||
n |
x |
n |
y |
n |
z |
|
||||
|
|
|
|
|
|
|
||||
n |
2 |
n |
2 |
n |
2 |
|
0.145 |
0.944 |
0.295 0.9999999986869448238139 1 |
|
x |
y |
z |
||||||||
|
|
|
|
|
|
|
||||
n |
3 |
n |
3 |
n |
3 |
|
0.984 |
0.167 |
0.050 |
|
x |
y |
z |
|
|||||||
|
|
|
|
|
|
|
Даже с учетом всех округлений получилось очень близкое к 1 число.
По заданному критерию прочности Мизеса вычислим эквивалентные напряжения:
экв 1 2 2 2 3 2 1 2 1 3 2 3 113.807 МПа
26
Часть III. Определение усилий в элементах статически неопределимых конструкций
Пример. 3 Решение статически неопределимой задачи с помощью составления уравнения совместности деформаций
Абсолютно жесткий брус закреплен на трех стрежнях и нагружен силой Р. Определить неизвестные усилия в стержнях.
l1 3м; l2 4м; l3 4.5м; F 2 10 4 м2 ; 160 МПа; E 200 10 9 Па
Допустим, что во всех стержнях системы возникают растягивающие усилия. Укажем их на исходной системе, воспользовавшись методом сечений, и рассмотрим равновесие отсеченной нижней части.
Составим три уравнения статики для плоской системы:
X 0
Y 0 : N1 N2 N3 P 0
momO1 0 : N2 a N3 32a P 34a 0
27
Имеем два нетривиальных уравнения статики и три неизвестных усилия. Система 1 раз статически неопределима. Дополним уравнения статики одним уравнением совместности деформаций. Допустим, что брус под действием силы сместится вниз. Отметим отрезки удлинений стержней l1, l2 , l3
Установим геометрическую связь между деформациями стержней.
Проведем горизонтальные линии в трапеции деформаций и запишем уравнения подобия прямоугольных треугольников:
|
|
l |
l |
2 |
|
|
l |
l |
|
|
|
|
1 |
|
= |
2 |
3 |
|
|
||
|
|
|
а |
|
а / 2 |
|
|
|||
|
|
|
|
|
|
|
|
|||
Деформации стержней: |
|
|
|
|
|
|
|
|
|
|
|
|
N l |
|
|
|
N l |
|
|
N l |
|
l |
= |
1 1 |
, l |
= |
|
2 2 |
, l |
= |
3 3 |
|
1 |
|
4EF |
2 |
|
|
5EF |
3 |
|
3EF |
|
|
|
|
|
|
|
|
Подставим в уравнение совместности деформаций:
N l |
|
N l |
|
N l |
|
|
N l |
1 1 |
2 2 |
|
2 2 |
3 3 |
|||
|
|
|
|||||
4EF |
|
5EF |
= |
5EF |
|
3EF |
|
|
а |
|
|
а / 2 |
|||
|
|
|
|
Объединим уравнения статики и совместности деформаций:
N1 N2 N3 P 0 |
|
|
|
|
|
||||||||||
N2 N3 |
|
3 |
P |
3 |
0 |
|
|
|
|
|
|||||
2 |
4 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
3N l |
4N |
l |
2 |
|
|
|
4N |
l |
2 |
|
4.5N l |
|
|||
1 1 |
|
|
2 |
|
=2 |
|
|
2 |
|
|
3 3 |
|
|||
4 |
|
5 |
|
|
5 |
|
|
3 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
Решая эту систему, получим значения неизвестных усилий:
N1 =0.366P, N2 =0.4P, N3 0.231P
28
Пример 4. Применение метода сил для систем, работающих на растяжение-сжатие
Определим усилия, возникающие при нагреве обоих стержней на T считаем абсолютно жестким. Система имеет одну лишнюю связь. основную систему.
Вид канонического уравнения:
|
11 |
X |
1 |
|
0 |
|
|
1T |
|
. Брус BG Начертим
Приложим к основной системе единичную силу
эпюру продольной силы |
N1 . |
X |
1 |
|
1
и построим единичную
Для нахождения единичного перемещения 11 воспользуемся формулой:
|
|
|
|
|
|
|
|
m |
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
ji |
|
|
|
|
|
N j |
l |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
i 1 |
EF |
|
|
|
|
|
i |
|
|
|
|
|
|
|||
i 2 |
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
2 |
|
|
4 a |
|||
|
|
2 |
|
|
|
2 |
|
a |
3 |
|
||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
11 |
|
N j |
l |
|
|
|
|
|
1 |
|
|
a |
|
|
||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
i 1 |
|
EF |
|
|
|
|
i |
|
EF |
|
|
|
|
|
|
|
|
|
|
|
EF |
Для вычисления температурного перемещения воспользуемся формулой:
i 2 |
|
|
|
|
|
||
|
|
|
i Ta 1 Ta |
3 Ta |
3 1 |
||
jT TlN |
1 |
||||||
i 1 |
|
|
|
|
|
Подставляя найденные коэффициенты в каноническое уравнение, определим лишнюю неизвестную:
29
X |
|
|
|
0.18 TEF |
|
1T |
|||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
11 |
|
Усилие во втором стержне можно найти, умножив на единичной продольной силы во втором стержне N1
получим:
X1 значение эпюры от3 . В результате,
С 0.32 TEF
30