Sapunov_Prikladnaya_teoriya_uprugosti_2008
.pdfгде N1 и N2 − значения главных напряжений на бесконечности; α − угол, который главная ось 1 составляет с осью x . Постоянная
Cможет быть принята равной нулю. Введем в рассмотрение новую функцию
Ω(z)=Φ(z)+ z Φ′(z)+Ψ(z) ,
которая на основании имеющихся соотношений будет иметь вид:
Ω(z)= |
κ(X +iY ) 1 |
|
|
|
|
|
1 |
|
|||
|
|
|
|||||||||
|
|
|
|
+ Γ + Γ′+Ο |
|
. |
|||||
2π(1 |
+ κ) z |
|
|||||||||
|
|
|
|
|
z2 |
|
Комплексная функция Ψ(z) может быть определена через новую функцию Ω(z) следующим образом:
Ψ(z)=Ω(z)−Φ(z)+ z Φ′(z) ,
что позволяет представить напряжения через комплексные функции Φ(z) и Ω(z). В частности, соотношение (1.8, а)
σy −iτxy = Φ(z)+Φ(z)+ z Φ′(z)+ Ψ(z) |
|
будет иметь вид: |
|
σy −iτxy = Φ(z)+Ω(z )+(z − z) Φ′(z) . |
(1.8, б) |
Аналогично можно преобразовать и соотношение, определяющее перемещения (1.2), если вместо функции ψ(z) ввести новую функцию
ω(z)= ∫Ω(z)dz = z Φ(z)+ ψ(z)+const .
Будем иметь:
71
|
|
|
|
|
|
|
|
|
|
′ |
(z)−ψ(z) |
|
|||||
2μ(u +iv) = κϕ(z)− zϕ |
|
|||||||
2μ(u +iv) = κϕ(z)−ω(z )−(z − z) |
|
|
(1.9, а) |
|||||
Φ(z)+ const . |
Перейдем теперь к решению первой основной задачи. Будем
считать, что заданы внешние нагрузки X + , Y + , X − , Y − , приложенных ко всей границе L ; верхними индексами " + " и " −" отме-
чены значения нагрузок, принимаемые соответственно на верхних и нижних краях разрезов. Кроме того, будем считать заданными
постоянные Γ = Γ = B (постоянная C принята равной нулю) и Γ′ = B′+ iC′ , т.е. считаем заданными напряжения на бесконечности.
Согласно соотношению (1.8, б) граничные условия в поставленной задаче принимает вид:
Φ+(t)+Ω−(t)=Y +−i X + , Φ−(t)+Ω+(t)=Y −−i X − .
Складывая и вычитая полученные уравнения, имеем:
[Φ(t) +Ω(t)]+ +[Φ(t) +Ω(t)]− = 2p(t) ,
[Φ(t) −Ω(t)]+ − [Φ(t) −Ω(t)]− = 2q(t) ,
где 2p (t), 2q (t) − заданные на L функции:
2p (t)=[Y ++Y −]−i [X + + X −] , 2q (t)=[Y +−Y −]−i [X + − X −] .
(1.10)
(1.11)
Так как Φ(∞)− Ω(∞)= −Γ′, то общее решение граничной задачи (1.11) определяется формулой:
Φ(z)−Ω(z)= |
1 |
∫ |
q (t) |
dt − |
|
′ . |
|
||
Γ |
(1.12) |
||||||||
πi |
t − z |
|
|||||||
|
L |
|
|
|
|
72
Далее, вводя в рассмотрение функцию
Χ(z)= |
n (z − a |
k |
)1/ 2(z −b |
|
)1/ 2 |
, |
|
|
|
|
|||||
|
|
|
∏ |
|
|
|
k |
|
|
|
|
|
|
||
|
|
|
k =1 |
|
|
|
|
|
|
|
|
|
|
|
|
запишем общее решение граничной задачи (1.10) 1: |
|
||||||||||||||
Φ(z)+ Ω(z)= |
|
1 |
|
|
∫ |
Χ+(t)p (t) |
dt + |
2Pn |
(z) |
. |
|||||
πi Χ(z) |
|
t − z |
|
Χ(z) |
|
||||||||||
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
Здесь Pn (z) обозначает полином степени не выше n : |
|
||||||||||||||
P |
n |
(z)= C |
0 |
zn +C zn−1 |
+ ... +C |
n |
, |
|
|
||||||
|
|
|
|
1 |
|
|
|
|
|
|
(1.13)
(1.14)
апод Χ+(t) подразумевается значение, принимаемое функцией
Χ(z) на левой стороне L .
Из формул (1.12) и (1.14) следует:
Φ(z)= Φ0 (z)+ PΧn((zz))− 12 Γ′ , Ω(z)= Ω0 (z)+ PΧn((zz))+ 12 Γ′ ,
где
Φ0 (z)= |
1 |
|
∫ |
Χ+(t)p (t) |
dt + |
1 |
|
2πi Χ(z) |
t − z |
2πi |
|||||
|
L |
|
|||||
|
|
|
Χ+(t)p (t) |
|
|
||
Ω0 (z)= |
1 |
|
∫ |
dt − |
1 |
||
2πi Χ(z) |
t − z |
2πi |
|||||
|
L |
|
|||||
|
|
|
|
|
|
(1.15)
∫ |
q (t) |
|
dt , |
|
t − z |
||||
L |
q (t) |
|
||
∫ |
dt . |
|||
t − z |
||||
L |
|
|
|
1 Решение дается без вывода. Полное решение задачи приводится в работе [7 ].
73
Остается определить полином Pn (z). Для определенности будем считать, что под функцией Χ(z) подразумевается ветвь, имеющая при больших значениях z вид:
Χ(z)= zn + an−1zn−1 + ...
Коэффициент C0 сразу определяется по первой из формул (1.15) и по условию Φ(∞)=Γ:
C0 = Γ + 12 Γ′ .
Остальные коэффициенты должны быть определены из условия однозначности перемещений, комплексное представление которых дано формулой:
2μ(u +iv) = κϕ(z)−ω(z )−(z − z)Φ(z)+ const .
Приведенное соотношение показывает, что выдвинутое условие будет выполняться, если функция κϕ(z)−ω(z ) будет возвращаться к своим первоначальным значениям, когда точка z описывает замкнутые контуры Λk , охватывающие отрезки akbk = Lk ( k =1, 2, ..., n ). Стягивая контуры Λk к отрезкам Lk , условие од-
нозначности перемещений можно привести к следующим равенствам:
2(κ +1)∫ |
Pn (t) |
dt + κ ∫[Φ0+ (t)−Φ0− (t)]dt + ∫[Ω0+ (t)−Ω0− (t)]dt = 0 , |
||
Χ(t) |
||||
L |
L |
L |
||
k |
|
k |
k |
которые представляют собой систему линейных уравнений относительно постоянных C1 , C2 , . . . , Cn .
74
В частном случае, когда края разрезов (щелей) свободны от на-
пряжений (растяжение пластины, ослабленной разрезами), функ-
ции Φ 0 (z) и Ω0(z) равны нулю и решение принимает достаточно простой вид:
Φ(z)= |
Pn (z) |
1 |
|
|
Ω(z)= |
Pn (z) |
1 |
|
|
|
||
Γ′ , |
Γ′ . |
|
||||||||||
|
− |
2 |
|
+ |
2 |
(1.16) |
||||||
Χ(z) |
Χ(z) |
В этом случае коэффициенты полинома Pn (z) определяются соотношениями:
C0 = Γ + |
1 Γ′ |
, ∫ |
Pn (t) |
dt = 0 |
( k =1, 2, ..., n ) . |
|
|||||||||||||||||
|
|
||||||||||||||||||||||
|
|
|
|
|
2 |
|
|
L |
Χ+(t) |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
При n =1 |
(растяжение пластины, ослабленной одним разре- |
||||||||||||||||||||||
зом), полагая a1 = −l , |
b1 = l |
( l − полудлина разреза), получим со- |
|||||||||||||||||||||
всем простые формулы: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
(2Γ + |
|
′)z |
− 1 |
|
|
Ω(z)= |
(2Γ + |
|
′)z |
+ |
1 |
|
′ . |
|
||||||||
Φ(z)= |
Γ |
|
′ , |
Γ |
|
||||||||||||||||||
Γ |
Γ |
(1.17) |
|||||||||||||||||||||
|
2 z |
2 |
−l |
2 |
|
2 |
|
|
2 |
z |
2 |
−l |
2 |
|
2 |
|
|
||||||
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Представленные функции определяют напряженное состояние неограниченной пластины с одиночным разрезом.
1.7. Напряженное состояние в вершине разреза (трещины) при растяжении пластины
В расчетах на прочность элементов конструкций и сооружений с трещинами отправной точкой является исследование напряжений и деформаций, возникающих в них под действием приложенных нагрузок. При этом особый интерес представляет область в непосредственной близости к вершине (кончику) трещины, размеры которой
75
малы по сравнению с размерами тела и трещины, поскольку именно в этой области начинается дальнейшее разрушение. В рамках линейной механики разрушения, исходящей из модели идеально упругого тела и представляющей трещину в виде разреза нулевой толщины, поверхности которого свободны от напряжений, задача определения напряженного состояния в вершине трещины сводится к граничной задаче теории упругости для плоскости с полубесконечным разрезом при соответствующем нагружении.
Рассмотрим задачу о растяжении плоскости с одиночным разрезом (трещиной), ориентированным перпендикулярно растягивающим напряжениям σ, приложенным на бесконечности, а нагрузка на верхнем и нижнем краях разреза равна нулю. Разрез занимает
область y = 0 , | x | ≤ l . |
|
Поскольку в рассматриваемом случае N1 = σ , |
N2 = 0 , α = π/ 2 , |
будем иметь: B = σ/ 4 , Γ = σ/ 4 (постоянная C |
равна нулю) и |
Γ′ = B′+ iC′ = σ/ 2 . |
|
Соответственно, комплексные функции Φ(z) и Ω(z), определяемые соотношениями (1.17), принимают вид:
Φ(z)= |
|
σ z |
|
− |
σ |
, |
Ω(z)= |
|
σ z |
|
+ |
σ . |
||||
|
|
2 |
|
2 |
4 |
|
|
2 |
|
2 |
||||||
2 |
z |
−l |
|
|
|
2 |
z |
−l |
|
4 |
||||||
|
|
|
|
|
|
|
|
|
|
Напомним, что полученные соотношения определяют напряженное состояние неограниченной пластины с одиночным разре-
зом при больших значениях z . Для перехода к задаче определе-
ния напряженного состояния в вершине трещины (к задаче для плоскости с полубесконечным разрезом) введем полярные координаты r , θ, поместив начало координат в вершине разреза, и новую
комплексную переменную ζ = z −l = re iθ |
( z = ζ +l ). Комплексные |
|||||
функции Φ(z) и Ω(z) будут иметь вид: |
|
|
|
|||
Φ(z)= |
|
σ(ζ +l) |
|
− |
σ |
, |
|
ζ(ζ + 2l) |
4 |
||||
2 |
|
|
76
Ω(z)= |
|
σ(ζ +l) |
+ |
σ . |
|
ζ(ζ + 2l) |
|||
2 |
|
4 |
При малых (сравнительно с длиной трещины) значениях ζ асимптотическое представление комплексных функций Φ (z) и Ω(z) в окрестности вершины разреза позволяет получить:
Φ(z) |
|
ζ |
|
→0 |
= Ω(z) |
|
ζ |
|
→0 |
= |
σ πl |
= |
KI e −iθ/ 2 |
, |
|
|
|
|
|||||||||||
|
|
|
|
2 2πζ |
2 2πr |
|||||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
где введено обозначение: KI = σ πl − коэффициент интенсивно-
сти напряжений для трещины нормального отрыва (для трещины типа I ).
Напряженное состояние в вершине трещины определим, используя первое из соотношений (1.8) и соотношение (1.8, а):
σx +σy = 2 [Φ(z)+Φ(z)] ,
σy −iτxy = Φ(z)+Ω(z )+(z − z) Φ′(z) .
Используя второе из представленных зависимостей, для нормального напряжения σy можем получить:
|
|
σy = Re[Φ(z)+Ω(z)+(z − z) |
Φ′(z)] = |
|
|
|||||||||
|
KI |
|
|
|
iθ |
iθ |
|
|
3iθ |
|
||||
= |
|
Re exp |
− |
|
|
+ exp |
|
|
−i sin θ exp |
2 |
|
|
||
2 2πr |
|
|
||||||||||||
|
|
|
|
2 |
|
2 |
|
|
|
|
|
|
KI |
|
θ |
|
θ |
|
3θ |
|||
σy = |
|
cos |
2 |
1 |
+sin |
2 |
sin |
2 |
. |
||
2πr |
|||||||||||
|
|
|
|
|
|
|
Аналогичным образом находим соотношение, определяющее касательное напряжение τxy :
77
−iτxy = |
K |
I |
|
−i sin θcos |
3θ |
|
τxy = |
|
K |
I |
cos |
θ |
sin |
θ |
cos |
3θ |
. |
||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
2πr |
2 |
2 |
2 |
|
||||||||||||||
|
2 2πr |
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||||
Соответственно, для нормального напряжения σx будем иметь: |
|
||||||||||||||||||||
|
|
|
|
σx = 2 [Φ(z)+ |
Φ(z)] −σy |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
KI |
|
|
|
θ |
|
θ |
|
3θ |
|
|
|
|
|
|
|
||
|
|
σx = |
|
|
cos |
2 |
1−sin |
2 |
sin |
. |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
2πr |
|
|
|
2 |
|
|
|
|
|
|
|
Обратим внимание, что введенный параметр KI = σ πl являет-
ся мерой сингулярности напряжений в области, прилегающей к вершине трещины, и играет для этой области такую же роль, как коэффициенты концентрации напряжений для макрообъемов тела. В отличие от коэффициентов концентрации коэффициенты интен-
сивности напряжений имеют размерность( МПа м или кг/мм3/2 ). Очевидно, что можно говорить о трещинах поперечного (тип трещины II ) и продольного (тип III ) сдвига.
Особенностью поля напряжений в вершине трещины является то, что асимптотика напряженного состояния не зависит от длины трещины, формы тела и схемы действующих нагрузок. С другой стороны, интенсивность асимптотического распределения напряжений определяется только коэффициентом интенсивности напряжений, являющимся функцией проложенной нагрузки, геометрии тела и трещины, но не зависящим от координат r , θ точки вблизи вершины разреза.
Таким образом, коэффициенты интенсивности напряжений являются базовыми характеристиками напряженно-деформирован- ного состояния в вершине трещины. Определение этих величин для тел с трещинами на основе решения соответствующих краевых задач является обширной самостоятельной областью математической теории упругости.
78
2. Изгиб тонких плит (пластин)
Пластиной (плитой) будем называть тело призматической или цилиндрической формы, у которого высота h мала по сравнению с размерами основания.
Пластины, толщина которых не превышает 15 наименьшего
размера основания, относятся к тонким. Расчеты таких пластин ведут на основе классической теории изгиба, базирующейся на гипотезах Кирхгофа. Расчеты пластин, толщина которых превышает 15 наименьшего размера основания, ведут на основе теории тол-
стых плит.
Плоскость, делящая толщину пластины пополам, называется срединной. Перемещения, которые получают точки срединной плоскости при нагружении пластины нагрузкой, перпендикулярной к срединной плоскости (при изгибе пластины), называют прогибами. Срединная плоскость после деформирования пластины перехо-
дит в срединную поверхность.
Пластину считают жесткой, если величина максимального прогиба (стрелы прогиба) не превышает 15 толщины. При построе-
нии теории изгиба жесткой пластины можно пренебречь напряжениями растяжения (сжатия), равномерно распределенными по толщине пластины (мембранными напряжениями). Для гибких пластин, наряду с изгибными напряжениями, учет мембранных напряжений необходим. Пластину принято считать абсолютно гибкой или мембраной, если стрела прогиба превышает толщину в 5 и более раз. При расчете мембраны можно пренебречь изгибными напряжениями по сравнению с мембранными напряжениями в срединной плоскости.
Построение инженерной теории изгиба тонких пластин базируется на следующих основных допущениях (гипотезах):
− прямая, нормальная к срединной плоскости пластины до изгиба, переходит в нормаль к срединной поверхности после изгиба (гипотеза Кирхгофа − Лява);
−напряжением, действующим в направлении, перпендикулярном срединной поверхности, можно пренебречь;
−срединная плоскость остается после изгиба нейтральной.
79
2.1. Изгиб прямоугольных пластин
Выберем систему координат так, чтобы плоскость xy являлась
срединной плоскостью пластины, а ось z располагалась в направлении ее толщины (рис. 2.1). Размеры пластины показаны на рисунке.
Гипотезу о прямых нормалях в рассматриваемом случае можно
трактовать как требование неизменности прямых углов в плоскостях xz и yz , что означает отсутствие сдвигов γxz и γyz в этих
плоскостях ( γxz = γyz = 0 ) и равенство нулю соответствующих ка-
сательных напряжений τxz и τyz ( τxz = τyz = 0 ).
Вторая гипотеза позволяет утверждать, что на площадках, параллельных срединной плоскости, можно пренебречь нормальным напряжением σz ( σz = 0 ). Соответственно, можно пренебречь взаимным нажатием горизонтальных слоев пластины и принять, что линейная деформация εz равна нулю ( εz = 0 ). Как следствие,
можно утверждать, что вертикальные перемещения точек (прогибы w ) в пределах нормали (по толщине пластины) остаются постоянными и равными вертикальным перемещениям точек срединной поверхности w = w(x, y). Поскольку принято, что срединная плоскость остается после изгиба нейтральной, для точек этой плоскости имеем u0 = 0 , v0 = 0 , w0 = w(x , y).
Принятые допущения существенно упрощают анализ напряжен- но-деформированного состояния изгибаемой пластины и позволяют в качестве расчетной схемы использовать элемент срединной плоскости пластины вместо ее элемента объема.
80