Full Text

Введение Вблизи поверхностей твердых тел на наноуровне проявляются особые физико-механические свойства, которые могут значительно отличаться от свойств материала вдали от поверхности - начиная отличием в химических связях для приповерхностных атомов и наличия свободных радикалов и заканчивая присутствием оксидных пленок и повышенной концентрацией примесных атомов [1-3]. Один из основных способов для учета поверхностных явлений - использование, помимо уравнений теории упругости, уравнения поверхностной упругости [1-4]. Достаточно общий метод решения задач теории упругости, учитывающий поверхностные эффекты, заключающийся в использовании характерных нестандартных граничных условий, вытекающих из уравнений поверхностной упругости, был развит в работах [1-6]. Однако к настоящему времени число полученных аналитических решений ограничено, что, возможно, связано с весьма «неудобным» видом получаемых граничных условий. Значительная часть полученных решений для задач, учитывающих поверхностные эффекты, ограничено довольно простыми геометриями, такими как сферическая пора [7-9], плоскость с одиночным круговым отверстием [10], одиночная пластина [11], однако относительный вклад поверхностных эффектов для круговых и сферических поверхностей минимален. В работах [12, 13] рассматривались задачи о плоскости с двумя круговыми отверстиями и в неконцентрической трубе, где влияние поверхностных эффектов было более существенным. Из недавних работ, посвященных исследованию вклада поверхностных эффектов, отметим [14-18]. В настоящей работе рассмотрена задача о круговом отверстии в подвергнутой растяжению полуплоскости; в данной задаче, несмотря на довольно простую геометрию, из-за достаточно малого расстояния между границами (т.е. между отверстием и границей полуплоскости) можно ожидать более весомого вклада поверхностных эффектов. Отметим, что теория поверхностной упругости, используемая в настоящей работе, является наиболее распространенной, но не единственной. В данной теории, ставшей уже традиционной, в дополнение к объемным напряжениям и деформациям, предполагается существование дополнительных поверхностных напряжений, действующих подобно напряжениям в мембранах, совпадающих с внешними поверхностями тел либо внутренними поверхностями раздела. Данные поверхностные напряжения связываются с объемными деформациями посредством аналога закона Гука, а дополнительные поверхностные деформации при этом не вводятся. Альтернативные теории (например, [19]) рассматривают дополнительные поверхностные (интерфейсные) деформации, связываемые с объемными напряжениями. Недавно предложены также варианты теории, объединяющие данные подходы [20-22]. Отметим также, что полученные решения представляют интерес не только сами по себе, но и в качестве вспомогательных решений при решении более сложных задач, например, задач о взаимодействии пор с растущими трещинами [23-24]. 1. Геометрия, биполярные координаты, общие соотношения Рассмотрим полуплоскость, содержащую круговое отверстие и растягиваемую вдоль ее границы постоянными напряжениями, приложенными на бесконечности (рис. 1). В данной задаче выполнены условия плоской деформации, а на контуре отверстия и границе полуплоскости вводятся дополнительные соотношения, характеризующие поверхностные эффекты (см. п. 2). В случае классической упругости задача была решена в [25-27]. В настоящей работе получено обобщение данного решения с учетом поверхностных эффектов на контуре кругового отверстия и границы полуплоскости. Следуя работе [25], введем биполярную систему координат (α, β), связанную с декартовой системой координат (х, y) следующим образом: (1) где Полюс O1 расположен на расстоянии a от полуплоскости. Масштабный коэффициент имеет вид (2) Рис. 1. Геометрия задачи Fig. 1. The problem geometry Координатные линии представляют собой набор окружностей с центрами на оси Oy. Окружности, у которых значение соответствуют положительным значениям, лежат выше оси Ox, а отрицательным - ниже оси Ox, координатная линия - соответствует оси Ox. Координатные линии представляют собой дуги окружностей, проходящих через точку O1, ортогональные окружностям Справа от оси Oy значения а с левой стороны - значение соответствует оси Oy, кроме отрезка OO1, где координата β терпит разрыв, а именно: слева от оси Oy а справа Непосредственно в точке O1 - а β не определено. Будем считать, что известное постоянное значение соответствует контуру отверстия, тогда, для радиуса отверстия R и расстояния между границей полуплоскости и центром отверстия d имеют место следующие соотношения: (3) Компоненты тензора напряжений определяются через бигармоническую функцию напряжений Φ(α,β) следующим образом [26]: (4) (5) Бигармоническая функция напряжения Φ(α, β) для плоской задачи, в которой содержится круговое отверстие, представляется в виде суперпозиции двух функций [27]: (6) где функция удовлетворяет заданным напряжениям на бесконечности, а функция F подбирается так, чтобы бигармоническая функция удовлетворяла известным условиям на контуре отверстия и равнялась нулю на бесконечности, т.е. Для рассматриваемой задачи функция Φ(α, β) имеет вид [27] (7) Для удовлетворения граничных условий на контуре отверстия и на границе полуплоскости представим напряжения в виде рядов Фурье: (8) где коэффициенты ck, bk, dk выражаются через параметры An, Bn, Cn, DnJ. 2. Постановка задачи Задача рассматривается в постановке плоской деформации. Всюду внутри тела предполагаются выполненными основные уравнения теории упругости (уравнения равновесия, совместности и закона Гука), что автоматически обеспечивается введением бигармонической функции (4). В случае классической упругости (т.е. без учета поверхностных эффектов) граничные условия имеют вид Целью задачи является нахождение неизвестных коэффициентов An, Bn, Cn, DnJ и вычисление значения на контуре отверстия. Однако если учитывать поверхностные эффекты на границах, то данные граничные условия могут не равняться нулю. Для решения задачи необходимо определить граничные условия с учетом поверхностной упругости. В рассматриваемой теории изменения любых физических свойств, в частности упругих постоянных, с уменьшением размеров тела моделируются введением дополнительных свойств, приписываемых бесконечно-тонкой поверхности, а объемные свойства предполагаются не зависящими от размера тела. Таким образом, с уменьшением размеров тела относительный вклад его поверхности возрастает, что приводит к изменению эффективных свойств с уменьшением размера. В общем случае, с использованием обобщенного закона Лапласа - Юнга [28], граничные условия на поверхностях, редуцирующихся вследствие выполнения условий плоской деформации в линии границы полуплоскости и контура отверстия, представимы в виде (9) где - тензор объемных напряжений; - вектор нормали к границе; - тензор поверхностных напряжений; - поверхностный градиент [28]. Для криволинейной границы поверхностный градиент поверхностных напряжений представляется следующим образом: (10) где R1, R2 - радиусы кривизн; h1, h2 - метрические коэффициенты; - орты криволинейной ортогональной системы координат (α1, α2). Для рассматриваемой геометрии расчетной области (см. рис. 1), когда цилиндрическая поверхность образована окружностью контура отверстия, для в биполярной (бицилиндрической) системе координат (9), параметры имеют вид [29] (11) (12) Аналогично для границы полуплоскости, для параметры имеют вид (13) (14) Следовательно, компоненты на границе полуплоскости не испытывают влияния поверхностных эффектов и равны нулю. Учитывая это, можно упростить выражение для бигармонической функции (7). Подставив (7) в (5) при учитывая (14), получаем (15) тогда выражение для бигармонической функции можно переписать в виде (16) где An, En, J - коэффициенты, подлежащие определению. Объединив два уравнения в (12), получим (17) Здесь поверхностное напряжение по-прежнему остается неизвестным. Для получения граничного условия выразим поверхностное напряжение через известные соотношения (законы Гука и Шаттлворса): (18) где - тензор модулей поверхностной упругости; - окружная поверхностная деформация; E, ν - соответственно модуль Юнга и коэффициент Пуассона (в данной формуле нет суммирования по повторяющемуся индексу). Из уравнений (18) следует (19) Для характерных значений модуля Юнга и констант поверхностной упругости величину R всегда можно выбрать такой, чтобы параметр ε был малым. С учетом выражения (19) граничные условия (12), (17) выражаются следующим образом: (20) (21) Решение уравнений (20), (21) получаем путем разложения напряжений по малому параметру ε: (22) Приравнивая члены с равными степенями малого параметра при ε, получаем рекуррентные системы для определения коэффициентов An, En, J: (23) 3. Нулевое приближение (классическая упругость) Нулевое приближение поставленной задачи (23) является решением задачи классической упругости [26] и имеет вид (24) В данных выражениях коэффициенты, являющиеся аналогами коэффициентов в уравнении (16), отмечены дополнительным верхним индексом 0, указывающим на соответствие классическому известному решению [26] т.е. решению без учета поверхностных эффектов. Выражения (24) являются нулевым приближением в разложении по параметру решения задачи, с учетом поверхностной упругости. 3. Применение рекуррентных соотношений С помощью соотношений (23) задача сводится к решению набора последовательных задач для каждого члена разложения по малому параметру ε. Если известно (m - 1)-е решение (т.е. известны для компонент напряжения, то аналогично (8), с помощью (5), (6) получаем (25) где имеют вид (26) (27) где - символ Кронекера. Граничные условия (последние два уравнения в (23)) могут быть представлены следующим образом: (28) где (29) Последние уравнения в (29) образуют систему для нахождения коэффициентов через (30) где (31) Приравнивая коэффициенты при синусах и косинусах в (30), аргументом функции которой является вторая координата β в биполярной системе координат (в дальнейшем для упрощения записи опуская верхний индекс (m)), получаем (32) где (33) Предположив J известным и применяя процедуру [26], а именно - умножая третье уравнение из (32) на и суммируя для всех n, а также учитывая условие сходимости рядов получаем (34) Из первого уравнения (32) непосредственно следует равенство (35) При значения вычисляются непосредственно (36) где (37) Таким образом, используя уравнения (31), (34), находим оставшиеся неизвестные коэффициенты (38) где (39) Кроме того, используя уравнение (35) и первое уравнение в (38), определяем коэффициенты J, A1: (40) Подставляя найденные коэффициенты (38), (40) в уравнение (27), получим необходимую компоненту напряжения Полученное решение (38), (40) можно рассматривать как обобщенное решение в [26]. 5. Результаты Для решения задачи упругие постоянные материала принимались соответствующими константам алюминия: ГПа, а константы, характеризующие свойства поверхности, - упругими постоянными, полученными в [5]: откуда модуль поверхностной упругости рассчитывался как Концентрации напряжений в точках A, B (см. рис. 1) вычислялись для различных значений относительного расстояния между отверстием и границей полуплоскости (рис. 2). Радиус отверстия для расчетов полагался равным нм. Рис. 2. Расположение отверстия в полуплоскости Fig. 2. Location of the hole in the half-plane На рис. 3, а, б представлены зависимости концентрации напряжений с учетом и без учета поверхностного напряжения в точках А и В соответственно. Зависимости b1, a1 соответствуют значениям концентрации напряжения с учетом и без учета поверхностного напряжения в точке A; b2, a2, - в точке B (см. рис. 3, б). а б Рис. 3. Концентрации напряжений с учетом (a1,2) и без учета (b1,2) поверхностного напряжения: а - в точке А; б - в точке В Fig. 3. Concentration of stresses with account of (a1,2) and without account of (b1,2) of the surface stress, (a) is in the point A, (b) is in point B На рис. 4 представлены зависимости с1, с2, соответствующие значению поверхностного напряжения на контуре отверстия при относительно малом и относительно большом расстояниях между отверстием и границей полуплоскости. Из полученных графиков видна разница между концентрациями напряжений с учетом и без учета поверхностного напряжения. При отдалении кругового отверстия от границы полуплоскости наблюдается уменьшение влияния границы полуплоскости на концентрацию напряжения на круговом отверстии. Рис. 4. Поверхностное напряжение в зависимости от координаты β на контуре отверстия Fig. 4. Surface stress depending on coordinate β at the contour of the hole Заключение Получено и исследовано решение задачи теории упругости о полуплоскости, содержащей круговое отверстие с учетом поверхностных эффектов. При помощи разложения в ряды Фурье переменных, представленных в биполярной системе координат, с использованием уравнений теории упругости, дополненных уравнениями поверхностной упругости, рассчитаны величины концентрации напряжений на контуре отверстия. Полученные выражения можно рассматривать как обобщенное решение задачи в традиционной постановке, не учитывающей вклад поверхностных эффектов. Проведены сравнения концентраций напряжений с учетом и без учета вклада поверхностных эффектов в различных точках контура отверстия. Рассмотрен вклад, вносимый поверхностными эффектами в зависимости от относительного расстояния между отверстием и границей полуплоскости. Показано, что, несмотря на довольно простую геометрию, благодаря достаточно малому расстоянию между отверстием и границей полуплоскости значения концентрации напряжений с учетом и без учета поверхностных эффектов существенно отличаются друг от друга.

About the authors

D. V Gandilyan

Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences

References

Statistics

Views

Abstract - 323

PDF (Russian) - 215