Серия «Физико-химические проблемы и высокие технологии строительного материаловедения» Выпуск № 1(10), 2015
16. Тарасевич Ю.И. Адсорбция на глинистых минералах. Киев: Наукова думка, 1975.
350 с.
17.Бондаренко, А. В. Адсорбционные свойства активированного монтмориллонита: Дисс. канд. хим. наук : 02.00.04 / А. В. Бондаренко. Липецк, 2002. 174 с.
18.Бельчинская, Л. И. Комплексная природозащитная технология использования лаковых композиций и смол в деревообрабатывающей промышленности: Дис. д-ра техн. наук: 11.00.11 / Л. И. Бельчинская. Екатеринбург, 1996. 530 с.
19.Арипов Э.А. Природные минеральные сорбенты, их активирование и
модифицирование, Ташкент: Фан, 1970. С. 56-62.
20.Смирнов А.Д. Сорбционная очистка воды. - Л., Химия, 1982, 168 с.
21.Тарасевич Ю.И. Природные сорбенты в процессах очистки воды. Киев: Наукова думка, 1981. 207 с.
22.Стрельникова, О. Ю. Адсорбция формальдегида из водного раствора на природных и модифицированных минеральных сорбентах / О. Ю. Стрельникова, Л. И. Бельчинская, Н. А. Ходосова // Научные ведомости БелГУ. Серия Естественные науки. 2011.
№15 (110). Вып. 16. С. 103-108.
23.Myriam M.et al. Structural and textural modifications of palygorskite and sepiolite under acid treatment. - Myriam M., Sua´rez M. & Martı´n-Pozas J.M. (1998) Clays and Clay Minerals, 46, P. 225-231.
24.Gala´n E. Properties and applications of palygorskite-sepiolite clays. (1996) Clay Minerals, 31.
25.Oliveira, R. et al. The increase of surface area of a Brazilian palygorskite clay activated with sulfuric acid solutions using a factorial design / Oliveira, R. N.I, Acchar, W.; Soares, G. D. A.; Barreto, L. S // Mat. Res. 2013. Vol.16 № 4.
________________________________________________________________________________
Ходосова Наталия Анатольевна - к.х.н, доцент кафедры химии Воронежского государственного архитектурно-строительного университета. E-mail: nhodosova@mail.ru Тел. 89155846434
Бельчинская Лариса Ивановна – д.т.н., проф., заведующая кафедрой химии Воронежского государственного лесотехнического университета им. Г.Ф. Морозова. E-mail: chem@vglta.vrn.ru Тел. 8(473)2-53-76-59
Фролов Вадим Юрьевич – студент 2 курса строительно-технологического факультета Воронежского государственного архитектурно-строительного университета. E-mail: frol656@ya.ru
71
Научный Вестник Воронежского государственного архитектурно-строительного университета
МЕТОДЫ КОНТРОЛЯ И ДИАГНОСТИКИ МАТЕРИАЛОВ И ИЗДЕЛИЙ
УДК 539.319: 539.219.2
А.А. Трещѐв, А.В. Морозова
РЕШЕНИЕ ЗАДАЧИ ТИПА КИРША ДЛЯ ПЛАСТИНЫ ИЗ АНИЗОТРОПНОГО РАЗНОСОПРОТИВЛЯЮЩЕГОСЯ МАТЕРИАЛА
Произведена постановка задачи типа Кирша для одноосно сжатой пластины конечных размеров с отверстием из ортотропного разносопротивляющегося материала, адаптированы определяющие соотношения, представленные в монографии А.А. Трещева [1] для данного класса задач, приведѐн вывод разрешающих уравнений. На основе построенной модели с использованием метода конечных элементов решена конкретная задача.
Ключевые слова: задача Кирша, одноосное сжатие, плоское напряжение, ортотропный разносопротивляющийся материал, концентрация напряжений, перфорированная пластина
A.A. Treshchev, A.V. Morozova
THE SOLUTION OF THE TASK LIKE THE KIRSCH'S PROBLEM FOR A PLATE MADE OF ANISOTROPIC MATERIAL WITH MULTIMODULUS BEHAVIOR
The solution of the task like the Kirsch's problem for the plate of eventual sizes with hole made of orthotropic material with multimodulus behavior is raising, the defining relations submitted in the monograph of A.A. Treshchov [1] are adapted for this class of problems. The article contains the conclusion of resolving equations. The concrete problem is decided on the basis of the built model with the use of method of finite elements.
Keywords: Kirsch's problem, monoaxial stretching, flat tension, ortotropic material with multimodulus behavior, concentration of tension, perforated plate.
Введение. Концентрация напряжений является одним из основных факторов, влияющих на прочность перфорированных элементов, поэтому исследование влияния различных концентраторов на напряженно-деформированное состояние конструктивных элементов из анизотропных разносопротивляющихся материалов является одной из актуальных задач механики деформированного твѐрдого тела.
Теоретическая часть. Рассматривается равновесие тонкой физически нелинейной пластины с отверстием радиусом r и толщиной h в направлении оси X3. Она имеет конечные размеры 2a и 2b по направлению основных осей декартовой системы координат координат X1, X2 и загружена в своей плоскости сжимающей равномерно распределѐнной по толщине нагрузкой px, которая приложена к внешнему контуру Г2 в соответствии с рисунком 1. Нагружение считается простым при активной деформации, поэтому для описания свойств материала принимаются нелинейные уравнения состояния, представленные через параметры нормированного пространства.
Здесь все функции, описывающие напряженно-деформированное состояние, зависят только от координат X1 и X2. Вектор перемещений имеет при этом две основные ненулевые составляющие u1 = u и u2 = υ в направлении осей X1 и X2, соответственно. Перемещения
u3 w в направлении оси Х3 есть, но они имеют второстепенное значение и определяются после решения краевой задачи тривиально: w e33х3 . При этом деформации e33 не зависят
от координаты Х3 и настолько малы, что в дальнейшем рассматриваться не будут.
________________________________________________________________________________
© Трещѐв А.А., 2015
72
Серия «Физико-химические проблемы и высокие технологии строительного материаловедения» Выпуск № 1(10), 2015
px=const
y |
Г3 |
||
|
|||
|
|
|
|
a |
|
|
|
|
|
|
=const |
|
|
|
x |
|
|
|
p |
0 |
|
|
x |
r |
|
|
|
Г1 |
|
|
|
Г2 |
Г2 |
|
|
|
b |
|
|
Г3 |
|
|
|
|
|
|
|
Рис. 1. Расчетная схема пластины
В рамках теории механики сплошных сред принимаются уравнения связи между основными компонентами тензора малой деформации и перемещениями в условиях плоского напряженного состояния для пластинки в виде:
e11 = u,1 ; e22 = u,2; 12 = u,2 + u,1. |
(1) |
По всей толщине пластинки и на обеих еѐ плоских поверхностях имеем равенство:
33 13 23 0 |
(2) |
Остальные компоненты тензора напряжений ij (i, j 1,2) не зависят от координаты
X3 и имеют постоянное по толщине h ненулевое значение.
Уравнения связи между компонентами тензора малых деформаций и напряжениями в условиях плоского напряженного состояния с учетом равенства (2) в форме метода переменных «параметров упругости» Биргера приводятся к виду [2]:
{e} [C]{ }, |
(3) |
где {е} {е11е22 12}Т – вектор-столбец ненулевых деформаций;{ } { 11 22 12}Т – векторстолбец ненулевых напряжений; [C] – матрица податливостей ортотропного
разносопротивляющегося материала, зависящих от вида напряженного состояния и от достигнутого уровня напряжений, компоненты которой определяются через тензорные параметры:
Ckkkk Akkkk Bkkkk kk (здесь и далее по индексам не суммировать);
Ciijj Aiijj Biijj ( ii jj) при i j ;Cijij Aijij Bijij 
2 ij при i j ,
где Aijkm, Bijkm – константы, определяемые из экспериментов по деформированию образцов
материала [1]; ij – нормированные напряжения, ij ij / S ij / |
|
|
ij ij [1]. |
||
Независимых констант будет 8 и они запишутся так [1]: |
|
|
73 |
|
|
Научный Вестник Воронежского государственного архитектурно-строительного университета
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1/Ε |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Α1111 = 0,5 1/Ε |
|
|
|
1/Ε |
; Β1111 = 0,5 1/Ε |
|
|
|
|
|
; Α2222 |
= 0,5 1/Ε |
2 |
1/Ε |
2 |
; |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1/E |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
B2222 = 0,5 1/E |
2 |
2 |
|
; A1122 0,5 |
ν |
|
|
/E |
|
|
ν |
/E |
; |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12 |
|
|
1 |
|
|
12 |
|
|
2 |
|
|
|
|
|
|
||||||||
B |
0,5 |
ν /E ν |
|
/E |
; A |
|
|
|
1/E |
|
1/E |
|
|
1 |
1/E 1/E |
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
1122 |
|
|
|
|
|
|
|
12 1 |
12 |
|
|
2 |
|
|
|
|
1212 |
|
|
|
|
|
12 |
|
|
|
|
|
12 |
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1/ E |
|
|
2 |
|
|
/ E |
|
|
/ E |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1/ E |
|
2 |
|
|
21 |
|
; |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
21 |
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,125 |
|
|
|
|
|
|
|
|
|
1/E |
|
|
|
|
|
|
1/E |
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
2 1/E 1/E |
|
|
|
2 |
1/E |
1/E |
|
2 |
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
1212 |
|
|
|
|
|
|
12 |
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
2 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
/E |
|
|
ν |
|
|
|
|
|
/E |
|
ν |
|
/E |
|
; ν |
|
/E |
|
|
ν |
|
|
/E |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
4 ν |
|
|
|
/E |
; ν |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
здесь (i, j 1,2), |
|
|||||||||||||||||||||||||||||||||||
|
21 |
|
1 |
|
|
|
|
|
21 |
1 |
|
|
ij |
|
|
|
j |
|
|
ji |
|
|
i |
|
|
|
|
ij |
|
j |
|
|
|
|
ji |
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
где E j – модули упругости при растяжении и сжатии в направлениях соответствующих главных осей анизотропии; ij – коэффициенты поперечной деформации при растяжении и
сжатии в направлениях соответствующих главных осей анизотропии; Eij –модули упругости
при сжатии и растяжении в направлениях под углом 45◦ к соответствующим главным осям анизотропии.
Для плоского напряженного состояния в ортотропном теле через каждую точку проходят две взаимно перпендикулярные плоскости упругой симметрии, в силу чего на основании закона упругости для однородного упругого разносопротивляющегося материала запишем:
e11 (А1111 В1111 α11) σ11 А1122 В1122 (α11 α22 ) σ22 ; e 22 А1122 В1122 (α11 α 22 ) σ11 (А2222 В2222 α22 ) σ22 ; e12 (А1212 В1212 
2 α12 ) 12.
Обращая уравнения (3) с учетом геометрических соотношений (1) получим зависимость между напряжениями и деформациями при плоском напряженном состоянии:
{ } [D( )]{е}, |
(4) |
где [D( )] [C] 1 – матрица жесткостей материала, |
зависящая от вида напряженного |
состояния и степени нагружения Принимая за основу те или иные физические соотношения, мы не изменяем
уравнения статико-геометрической природы [1]. Основные положения и зависимости классической теории упругости остаются справедливыми, поэтому уравнения равновесия для данной задачи при условном отсутствии объемных сил представим следующим образом
[3]:
11,1 12,2 0 ; |
12,1 22,2 0. |
(5) |
|
При статических граничных условиях на контуре пластинки имеем: |
|
||
р1 11 cos( x1) 12 cos( x2 ) ; |
(6) |
||
р2 22 cos( x2 ) 12 cos( x1), |
|||
|
|||
где cos( xk ) - косинусы углов между осями |
xk и внешней нормалью к площадке, |
на |
|
которую действуют усилия с проекциями рk , (k 1,2). |
|
||
74 |
|
|
|
Серия «Физико-химические проблемы и высокие технологии строительного материаловедения» Выпуск № 1(10), 2015
Сохраняется одно из общих уравнений неразрывности деформаций Сен-Венана:
е11,22 е22,11 12,12 |
(7) |
Выражения (1), (5), (3) содержат восемь уравнений с восьмью неизвестными, которые можно решить по методу сил или по методу перемещений. Для метода сил заменив в зависимостях (7) деформации через напряжения (3) в сочетании с уравнениями равновесия
(5) получим систему трех уравнений, имеющую три неизвестных напряжения. Таким образом придѐм к третьему уравнению из разрешающей системы:
С1212 12,12 С1212,12 12 С1212,2 12,1 С1212,1 12,2 |
С1111,22 С1122,11 11 (С1122,22 |
С2222,11 22 С1111 11,22 С1122 11,11 22,22 С2222 22,11 2 С1111,2 11,2 |
|
С1122,2 22,2 С1122,1 11,1 С2222,1 22,1 . |
(8) |
Для метода перемещений в два уравнения равновесия (5) подставляем выражения напряжений (4), в которых сразу заменяем деформации согласно зависимостям (1). В результате этих преобразования получаем два уравнения с двумя неизвестными перемещениями [4]:
D11u,11 D11,1u,1 D12υ,12 D12,1υ,2 D33(u,22 υ,12 ) D33,2 (u,2 υ,1 ) 0 |
|
D22υ,22 D22,2υ,2 D12u,12 D12,2u,1 D33(u,12 υ,22 ) D33,1(u,2 υ,1 ) 0 , |
(9) |
где Dij – компоненты матрицы жесткостей материала.
Для решения задачи принимается плоская конечно-элементная модель с двумя степенями свободы в узле. Вся область, ограниченная контуром пластинки разбивается треугольными симплекс конечными элементами [4]. В месте ожидаемой концентрации напряжений, т.е. около отверстия сетка конечных элементов сгущается. Перемещения в произвольной точке элемента u(x1, x2 ) и υ(x1, x2 ) представляются через перемещения
узлов элемента {U} следующим образом:
|
|
|
|
|
|
|
u |
|
υ T [N]{U}, |
|
|
(10) |
|||
где {U} u |
i |
υ |
u |
j |
υ |
j |
u |
k |
υ |
k |
Т |
; i, j, k – номера узлов конечного элемента. |
|||
|
i |
|
|
|
|
|
|
|
|
|
|||||
Поля перемещений в пределах выбранного в задаче треугольного элемента |
|||||||||||||||
определяются |
с |
|
помощью |
интерполяционного |
полинома |
первой степени [5]: |
|||||||||
u 1 2x1 3x2; |
υ 4 5x1 6x2 . Функции формы [N] принимаются в следующем |
||||||||||||||
виде: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
[N] |
Ni |
0 N j 0 |
Nk |
0 |
(11) |
|||
|
|
|
|
|
|
|
|
|
|
, |
|||||
|
|
|
|
|
|
|
|
|
|
0 Ni 0 N j |
0 Nk |
|
|||
где Ni ai bix1 cix2 /(2 ), |
|
(i j k), |
площадь |
треугольного элемента; |
|||||||||||
ai x1jx2k x1kx2j ; |
bi x1i x2k ; ci x1k x2j . |
|
|
|
|||||||||||
Пользуясь уравнениями (1) и продифференцировав матрицу [N] , получим зависимости между деформациями и перемещениями узлов элемента:
75