Численное моделирование экспериментов по определению типа начальной анизотропии упругих материалов

Автор: Астапов Юрий Владимирович, Христич Дмитрий Викторович

Журнал: Вычислительная механика сплошных сред @journal-icmm

Статья в выпуске: 4 т.8, 2015 года.

Бесплатный доступ

Для материала, в дополнение к известному определению главных осей анизотропии, вводится понятие канонических осей анизотропии как осей декартовой прямоугольной системы координат, в которой тензоры упругих свойств материала имеют наименьшее число независимых ненулевых компонент. Разработана программа экспериментов, дающая возможность установить канонические оси в материалах с различным типом анизотропии. При этом для поиска положения главных осей предлагается выполнить три эксперимента на одноосное сжатие. Выявлено, что у ромбического, моноклинного и триклинного материалов канонические оси анизотропии совпадают с главными осями анизотропии. Для тригонального, тетрагонального и гексагонального материалов однозначно находятся одна главная ось анизотропии - главная поворотная ось, и совпадающая с ней каноническая ось анизотропии. Для отыскания ориентации двух других канонических осей анизотропии требуется еще два эксперимента: на растяжение-сжатие и сдвиг в плоскости, перпендикулярной главной поворотной оси. В кубическом материале главные оси анизотропии выбираются совпадающими с осями лабораторной системы координат. Для конкретизации положения канонических осей анизотропии относительно осей лабораторной системы оказывается достаточным знание результатов двух экспериментов на одноосное сжатие. Разработана программа механических экспериментов для распознания типа начальной упругой анизотропии материала по данным опытов в канонических осях анизотропии в случае, когда эти оси совпадают с осями лабораторной системы координат. Для того чтобы отличить изотропный материал от кубического, необходимо провести эксперимент на двухосное растяжение-сжатие по направлениям любых двух лабораторных осей и эксперимент на сдвиг в содержащей их плоскости. Такие же два эксперимента в плоскости, перпендикулярной главной поворотной оси, позволяют выделить тригональный, тетрагональный и гексагональный материалы. Выполненное компьютерное моделирование задач по идентификации типа начальной упругой анизотропии материала, аналогичных экспериментам, подтвердило работоспособность предложенных программ экспериментов.

Еще

Анизотропные материалы, главные оси анизотропии, программа механических экспериментов, численный эксперимент, моделирование

Короткий адрес: https://sciup.org/14320782

IDR: 14320782   |   DOI: 10.7242/1999-6691/2015.8.4.33

Текст научной статьи Численное моделирование экспериментов по определению типа начальной анизотропии упругих материалов

Деформируемые твёрдые тела при малых деформациях проявляют свойство упругости. При бесконечно малых деформациях и постоянной температуре соотношения, задающие связь между тензорами напряжений и деформаций, асимптотически стремятся к закону Гука:

Е = C ■■ S ,

где ε — линейный тензор деформаций; C — постоянный тензор упругих податливостей четвёртого ранга; S — тензор истинных напряжений. Для материала с каждым типом анизотропии характерна своя структура тензора упругих податливостей C [1–5]. Поэтому тип начальной упругой анизотропии можно выяснить из экспериментов при бесконечно малых деформациях. Проблема состоит лишь в том, что для анизотропного материала необходимо знать, какой точечной группой симметрии обладают его механические свойства, при этом сам материал не обязательно должен быть кристаллическим. Для кристаллов же установление точечной группы симметрии позволяет отнести их к одной из известных кристаллографических систем [2–5]. В качестве объекта исследования рассматривается кубический образец представительных размеров.

В работе [6] задача определения типа анизотропного материала решена в общем виде на основе разложения известного из экспериментов тензора модулей упругости на изотропную (постоянную) часть и части, содержащие два девиатора и нонор. Приведены формулы разложений матриц модулей упругости материалов всех кристаллографических сингоний. Эти разложения для разных типов анизотропии имеют различный вид, поэтому, с опорой на них, проблема идентификации анизотропного материала полностью решается.

Существуют способы нахождения 21 компоненты тензора упругости в лабораторной системе координат с применением акустических волн [10], методов голографической и спекл-интерферометрии [7], экспериментов по статическому однородному деформированию [8]. Однако при этом предварительно должны быть известны значения 21 компоненты тензора упругости или тензора упругих податливостей в некоторой системе координат [6–10]. Как показано в работах [3, 12], тип анизотропии материала распознается, если известно каноническое представление тензора упругости в главных осях анизотропии.

Альтернативный подход к решению проблемы отнесения анизотропных материалов к тому или иному типу состоит в разработке программы механических экспериментов с макрообразцами, по результатам которых можно, не имея конкретных значений упругих констант, провести классификацию материалов по кристаллографическим системам. Такая экспериментальная программа позволит определить тип симметрии свойств анизотропного материала в том случае, когда она априори неизвестна.

2.    Определение положения главных и канонических осей анизотропии материала

В работах [3, 11, 12] была предложена программа экспериментов для идентификации типа упругой симметрии анизотропного материала, структура которого в общем случае заранее неизвестна. Базовым в программе является эксперимент по выяснению положения в материале главных осей анизотропии. Согласно В.В. Новожилову [13], главными осями анизотропии называются главные оси тензора напряжений, возникающих в анизотропном материале в ответ на чисто объёмную деформацию. В работе [14] показано, что главные оси анизотропии отождествимы с главными осями тензора деформаций, возникающих в материале при гидростатическом сжатии, когда тензор напряжений S = S 0 , где S 0 — шаровой тензор.

В статье [15] описана возможность установления главных осей анизотропии из трёх экспериментов на сжатие. Тензор напряжений S 0 представляется в виде суммы

S 0 = S 1 + S 2 + S 3 , (2)

где S 1 = - 1 1 e 1 e 1, S 2 = - 1 2 e 2 e 2, S 3 = - 1 3 e 3 e 3 — тензоры, каждый из которых описывает сжатие одинаковым напряжением t, = t вдоль одной из осей некоторой фиксированной декартовой прямоугольной системы координат O x'y'z' с базисом e 1, e 2, e 3, которую в дальнейшем будем называть лабораторной. Полагаем, что оси лабораторной системы проходят по рёбрам рассматриваемого кубического макрообразца. Каждому тензору напряжений S 1 , S 2 , S 3 соответствует тензор деформаций ε 1 , ε 2 , ε 3 , связанный с ним законом Гука (1):

  • £ 1 = C ■■ S 1 , £ 2 = C ■■ S 2, £ 3 = C ■■ S 3 . (3)

В силу линейности закона Гука, с учётом формул (2) и (3), получаем:

  • £ о = C ■ ■ S о = C ■ ■ ( S i + S 2 + S з ) = £ 1 + £ 2 + £ 3 . (4)

Аддитивность представления (4) разрешает отождествить опыт на всестороннее сжатие с тремя экспериментами на одноосное сжатие, в которых требуется измерить все компоненты тензора деформации образца. Выполнение трех таких одноосных экспериментов позволяет идентифицировать главные оси анизотропии материала как главные оси тензора деформаций (4), который описывает отклик образца на всестороннее сжатие. Если не представляется возможным измерить в экспериментах все компоненты тензоров деформаций ε 1 , ε 2 , ε 3 , то недостающие следует выразить через геометрические параметры деформированных образцов [15].

Главные оси a 1 , a 2 , a 3 тензора деформаций ε 0 (4) являются главными осями анизотропии материала и определяются методами линейной алгебры. В результате нахождения главных осей этого тензора становится известной ориентация главных осей анизотропии материала в пространстве относительно лабораторной системы координат.

В случае кратных собственных значений тензора ε 0 главные оси анизотропии отыскиваются в общем случае с точностью до трёх параметров. Однако тензор упругих податливостей C имеет наименьшее число независимых констант исключительно в базисе, который присущ элементам симметрии свойств материала. Для кристаллов такой базис связывают с кристаллофизическими осями [4, 13, 16], а для композитных материалов или древесины — с тремя осями анизотропии по преимущественным структурным направлениям [17].

В работе [18] показано, что в случае, когда внутреннее строение некоторого анизотропного материала заранее неизвестно, базис декартовой прямоугольной системы координат для элементов симметрии свойств этого материала может быть установлен из механических экспериментов. Оси этой системы координат названы каноническими осями анизотропии материала. Их базисные векторы обозначим как k 1 , k 2 , k 3 . Канонические оси анизотропии — это оси такой декартовой прямоугольной системы координат, в которой тензоры, описывающие свойства материала, имеют наименьшее число независимых ненулевых констант. Канонические оси анизотропии с отвечающими им элементами симметрии кристаллической структуры всегда совпадают с главными осями анизотропии (см. [13]). Однако главные оси анизотропии, найденные в соответствии с понятием В.В. Новожилова, допускают произвол в выборе их ориентации и в общем случае могут не совпадать с каноническими осями.

Рассмотрим три возможных случая кратности главных значений тензора деформаций ε 0 (4) и для каждого из них покажем способы конкретизации положения канонических осей анизотропии.

Итак, если по результатам экспериментов (3) оказывается, что все собственные значения тензора ε 0 различны ( S 1 * s 2 * s 3 ), то его собственные векторы распознаются однозначно и задают ориентацию главных осей анизотропии a 1 , a 2 , a 3 . Канонические оси анизотропии k 1 , k 2 , k 3 в этом случае совпадают с главными осями анизотропии, а материал по своим свойствам может быть отнесён к ромбическому, моноклинному или триклинному. Методика идентификации этих трёх материалов разработана в статье [19].

В случае, когда два из трёх собственных значений тензора деформаций (4) одинаковы ( S 1 = s 2 * S 3 ), материал является одноосным кристаллом — тригональным, тетрагональным или гексагональным. При этом однозначно определяется только одна главная ось анизотропии — главная поворотная ось a 3 . Базисные векторы a 1 и a 2 могут быть выбраны произвольно, а базисные векторы k 1 , k 2 необходимо согласовывать с боковой поворотной осью; для вектора к 3 справедливо равенство: к 3 = a 3. Базисы a 1, a 2, a 3 и k 1 , k 2 , k 3 связаны ортогональным тензором поворота

Q 3 = cos ф ( к 1 к 1 + к 2 к 2 ) + sin ф ( к 1 к 2 - к 2 к 1 ) + к 3 к 3

так, что а ' = к ' Q 3 , где i = 1,2,3.

Чтобы выяснить, каков угол ф в случае тригонального материала, необходимо провести один эксперимент: двухосное растяжение–сжатие вдоль векторов a 1 , a 2 . В этом опыте тензор напряжений имеет вид S 4 = t 4 ( а 1 а 1 - а 2 а 2 ) , а измеряемые компоненты тензора деформаций вычисляются по формулам:

S 11 S 22 t 4 ( C 1111 C 1122 ), S 33 S 12 0, S 1 3 = 2 1 4 C 1123 sin3 ф, s 23 = 2 1 4 C 1123 cos3 ф .

Из выражений (6) следует формула для угла ф в тригональном материале:

ф ( тр ) = ' arctg ^.

3         S 23

Определение угла ф, характерного для тетрагонального материала, требует проведения двух экспериментов: на растяжение–сжатие вдоль векторов a1 , a2 и на сдвиг в плоскости этих векторов, который осуществляется при выполнении растяжения-сжатия под углом 45° к направлениям а1, а2. В первом из этих экспериментов тензор напряжений составляет S4 = 14 (а1 а1 - а2а2), а измеряемые компоненты тензора деформаций равняются

= -Б 22 = 1 4 ( ( С - C H22 )COS 2 2 ф + 2 C 1212 Sin 2 2 ф ) , ^ 12 = t 4 ( 2 C 1212 - ( C 1111 - C 1122 ) ) sin 2 ф cos 2 ф, Е 13 = Б 23 = ^ 33 = 0.

Во втором опыте тензор напряжений имеет вид S5 = t5 (a1 a2 + a2a1), а измеряемые деформации находятся как sn =—е22 = 15(2C1212 -(C1111 -C1122))sin2фcos2ф,

1 12 = 1 5 ( ( C 1111 - C 1122 )sin2 2 ф + 2 C 1212 cos2 2 ф ) ,                                  (9)

^ 13 = 1 23 = 1 33 = 0.

Из (8), (9) получаются следующие формулы для угла ф в случае тетрагонального материала:

ф ( т ) =- 'arctg 2 1 12 t 5 — или ф ( т ) =- 'arctg    2 1" t 4— .                    (10)

  • 4         S 11 t 5 - S 12 t 4                         4         S 11 t 5 - S 12 t 4

Ориентация канонических осей анизотропии относительно лабораторной системы координат определяется по углу ф (см. формулы (7) и (10)) для тригонального и тетрагонального материалов, соответственно, из соотношений:

к 1 = a 1 cos ф- a 2 sin ф ,    к 2 = a 'sin ф + a 2cos ф ,    к 3 = a 3.

Упругие свойства гексагонального материала инвариантны относительно любых поворотов вокруг вектора к 3 = a 3 [2]. В гексагональном материале измеряемые компоненты тензора деформаций составляют: 1 11 = -1 22 = t 4( C 1111 - C 1122), 1 12 = 1 1 3 = 1 23 = 1 33 = 0 (так как для такого материала C 1212 = ( C am - C 1122)/2) при действии напряжений S 4 = t 4 ( a1 a 1 - a 2 a 2 ) [3]; 1 12 = 1 5( C 1111 - C 1122), 1 11 = 1 22 = 1 13 = 1 23 = 1 3 3 = 0 — при S 5 = 1 5 ( a 1 a 2 + a 2 a 1).

Если в результате эксперимента по установлению главных осей анизотропии по В.В. Новожилову получается, что главные значения тензора деформаций (4) одинаковы ( е 1 = е 2 = е 3), то в качестве главных осей анизотропии a 1 , a 2 , a 3 следует рассматривать оси лабораторной системы координат, в которой проводился эксперимент. Такой выбор векторов a 1 , a 2 , a 3 позволяет использовать данные экспериментов на одноосное сжатие для нахождения канонических осей анизотропии k 1 , k 2 , k 3 в кубическом материале.

Взаимная ориентация векторных базисов a i и k i задается ортогональным тензором поворота Q = qj к‘кj : a ' = к‘ Q ( i = 1,2,3 ). Для компонент тензора Q выполняются тождества:

q 11 + q 21 + q 31    1,              q 12 + q 22 + q 32    1,              q 13 + q 23 + q 33    1,

  • q 11 q 12 + q 21 q 22 + q 31 q 32 = 0,      q 11 q 13 + q 21 q 23 + q 31 q 33 = 0,      q 12 q 13 + q 22 q 23 + q 32 q 33 = 0.

Покажем, что для определения положения канонических осей анизотропии k i относительно осей лабораторной системы a i достаточно использовать результаты двух экспериментов на одноосное сжатие (3) в лабораторной системе координат.

В первом эксперименте S 1 = - t a 1 a 1. Измеряемые деформации выражаются через константы податливости и компоненты тензора Q следующим образом:

1 11 =- 1 [ Q 1 ( C 1111 - ( C 1122 + 2 C 1212 ) ) + C 1122 + 2 C 1212 ] ,

I 33 =- 1 [ Q 5 ( C 1 13 =- 1 [ Q 8 ( C i

-

-

■ ( C 1122 + 2 C 1212 ) ) + C 1122 ] , ( C 1122 + 2 C 1212 ) ) ] ,

1 22 = - 1 [ Q 4 ( C 1111 - ( C 1122 + 2 C 1212 ) ) + C 1122 ] ,

S 12

S 23

= - 1 Q 7 ( C

-

( C 1122 + 2 C 1212 ) ) ; ,

=- 1 [ Q 9 ( C 1111 - ( C 1122 + 2 C 1212 )£,

где обозначено        Q 1    q n + q 12 + q 13 ,        Q 4    q n q 21 + q 12 q 22 + q 13 q 23 ,        Q 5    q 11 q 31 + q 12 q 32 + q 13 q 33 ,

Q 7    q 21 q 11 + q 22 q 12 + q 23 q 13 ,

Q 8    q 31 q 11 + q 32 q 12 + q 33 q 13 ,

Q 9    q 21 q 31 q 11 + q 22 q 32 q 12 + q 23 q 33 q 13 .

Во втором эксперименте S 2 = - t a 2 a 2, и измеряемые компоненты тензора деформаций е 2 связываются с компонентами тензора Q через константы податливости формулами вида:

SH =— t [ Q 4 ( С цц ( C 1122 + 2 ^212 ) ) + C 1122 + 2 C 1212 ] ,

3 = — t [ Q 6 ( C 1111 ( C 1122 + 2 C 1212 ) ) + C 1122 ] ,

Ё 23 =— t [ Q 11 ( C 1111 ( C 1122

_£ =— t [ Q2 ( C 312 =— t [ Q10 ( C

( C 1122 + 2 C 1212 ) ) + C 1122 ] ,

( C 1122 + 2 C 1212 ) ) ] ,

+ 2 C 1212 ) ) ] ,    S 13 =— t [ Q 12 ( C 1111 ( C 1122 + 2 C 1212 ) ) ] .

Здесь Q 2    q 21 + q 22 + q 23 , Q 6    q 21 q 31 + q 22 q 32 + q 23 q 33 , Qw   q 11 q 21 + q 12 q 22 + q 13 q 23 , Q u    q 31 q 21 + q 32 q 22 + q 33 q 23 ,

Q 12 = q 11 q 31 q 21 + q 12 q 32 q 22 + q 13 q 33 q 23 .

Для отыскания девяти компонент тензора Q независимых соотношения из (12), (13):

используем шесть соотношений (11) и четыре

t [ Q 7 ( C 1111 ( C 1122 + 2 C 1212 ) ) ] = S 12 ,

t [ Q 9 ( C 1111 ( C 1122 + 2 C 1212 ) ) ]=Ё 23 ,

t [ Q 8 ( C 1111 ( C 1122 + 2 C 1212 ) ) ] = S 13 ,

t [ Q 10 ( C 1111 ( C 1122 + 2 C 1212 ) ) ]=V

Исключая из них множитель (Cm1 — (C1122 + 2C1212)), получаем три уравнения относительно компонент тензора Q :

л3       л3   л л3 — £13 („ „3 , „ „33 \ q31 q11 + q32 q12 + q33q13 = — ( q21 q11 + q22q12 + q23q13 ), £12

л л л2 , „ „ „2         „2 _ S23 I „ 3      33 \

q21 q31 q11 + q22q32q12 + q23 q33q13 = — ( q21 q11 + q22q12 + q23q13 ), _

3         3          3 _Л2/",, „3         33 \ q11 q21 + q12q22 + q13q23 = — (q21 q11 + q22q12 + q23q13 ).

S 12

Численное решение системы уравнений (11), (14) позволяет найти компоненты тензора Q , то есть установить ориентацию канонической системы координат в кубическом материале относительно лабораторной системы координат по измеряемым в опытах деформациям.

Таким образом, разработана система механических экспериментов, с помощью которых можно в однородном кубическом образце определять направления канонических осей анизотропии, связанных с элементами симметрии свойств материала. Положение этих осей относительно лабораторной системы координат представляется для одноосных кристаллов соотношениями (5), (7), (10), для кубического материала — соотношениями (11), (14). При этом сначала распознаются тригональный, тетрагональный и кубический материалы. Осуществляется это по характерным для каждого типа материала ненулевым компонентам тензора деформаций, измеряемым в экспериментах, и выражениям (6), (8), (9), (12), (13). В случае совпадения канонических осей анизотропии с осями лабораторной системы координат для конкретизации типа анизотропии в материале требуется проведение дополнительных экспериментов.

3.    Программа экспериментальной идентификации изотропного и кубического материалов

Пусть в опытах на одноосное сжатие в лабораторной системе координат достоверно находятся сдвиговые деформации, например, £12, e12, тогда материал является кубическим, а ориентация его канонических осей анизотропии определяется соотношениями (11), (14). Если же измеряемые сдвиговые деформации не превышают погрешности измерений и можно считать, что 3^ = £j = 0 (i * j), то, как показано выше, канонические оси анизотропии k1 , k2 , k3 совпадают с главными осями a1 , a2 , a3 и осями лабораторной системы координат. В этом случае, для того чтобы отличить изотропный материал от кубического, требуется провести эксперимент на двухосное растяжение–сжатие по направлениям векторов a 1, a2 с тензором напряжений S4 = 14 (а1 а1 — а2а2) и эксперимент на сдвиг с тензором напряжений S5 = 15 (a1 a2 + a2a1) в плоскости этих векторов.

В качестве реакции на напряжения S 4 и в изотропном, и в кубическом материалах, в соответствии с представлениями их тензоров податливостей C [3], возникают деформации:

^ .) = 8 < к ) = 1 4 ( C nn - C ^) ( a 1 a 1 - a 2 a 2 ) = ( C nn - C ^) 5 4 ,                      (15)

где Cijkl — компоненты тензора упругих податливостей.

В опыте на сдвиг при напряжениях S 5 в изотропном материале деформации составляют:

e 5 из ) = 2 1 5 C 1212 ( a 1 a 2 + a 2 a 1 ) = 1 5 ( C im C 1122 ) ( a 1 a 2 + a 2 a 1 ) = ( C nn C 1122 ) 5 5.

Тогда

(4)       (5)

£ 11      £ 12

, t 4 t 5

так как для изотропного материала C 1212 = C 2323 = C 1313 = ( C m1 C 1122)/2 [3].

В кубическом материале при сдвиге с S 5 в плоскости векторов a 1 , a 2 появляются деформации, описываемые тензором:

E < к ) = 2 1 5 C 1212 ( a 1 a 2 + a 2 a 1 ) = 2 C 1212 5 5.                                       (17)

Для компонент тензора упругих податливостей кубического материала C1212 = C2323 = C1313 * * (C1111 - C1122)/2 [3], п0ЭТ0мУ е5к) = 2C121255 * (C1111 — C1122)55.

Из выражений (15), (17), (18) следует, что для кубического материала имеет место неравенство:

(4)       (5)

£ 11    . £12

.

t 4 t 5

Итак, если в опытах на двухосное растяжение–сжатие и сдвиг в плоскости векторов a 1 , a 2 выполняется равенство (16), то материал изотропный, если справедливо неравенство (19), то материал кубический. В двух указанных экспериментах компоненты тензоров деформаций ε 4 , ε 5 выражаются через измеряемые геометрические параметры деформированных образцов [15].

4.    Программа экспериментальной идентификации одноосных кристаллов

Под одноосными кристаллами [2] будем понимать кристаллы, имеющие главную поворотную ось третьего, четвёртого или шестого порядков, то есть относящиеся к тригональной, тетрагональной или гексагональной кристаллографическим системам. Для материалов с такими типами анизотропии канонические оси определены в разделе 2. При канонических осях анизотропии k 1 , k 2 , не совпадающих с главными осями анизотропии a 1 , a 2 , тригональный материал можно отличить от тетрагонального и гексагонального по наличию в измерениях при двухосном растяжении–сжатии ненулевых компонент деформаций £ 13 и £ 23. У гексагонального материала в эксперименте на растяжение-сжатие вдоль векторов a 1 , a 2 с тензором напряжений 5 4 = t 4 ( a 1 a 1 a 2 a 2 ) измеряемые компоненты тензора деформаций имеют вид:

£ 11       £ 22     t 4 ( C 1111    C 1122 ),        £ 12     £ 13     £ 23     £ 33     0.

Гексагональный материал отличается от тетрагонального тем, что в указанном опыте ему соответствует £ 12 = 0, в то время как у тетрагонального деформация £ 12 * 0 и вычисляется согласно выражению (8).

Рассмотрим случай, когда главные оси анизотропии, установленные экспериментально, совпадают с каноническими осями анизотропии одноосных кристаллов, то есть ф = 0 . В этом случае в эксперименте на растяжение-сжатие вдоль векторов a 1 , a 2 с тензором напряжений 5 4 = t 4 ( a 1 a 1 a 2 a 2 ) тензоры деформаций, в соответствии с представлениями их тензоров податливостей C [3], составляют:

– для тригонального материала

^) = 1 4 ( C

a 1 a 1 a 2 a 2 ) + 2 1 4 C 112 3 ( a 2 a 3 + a 3 a 2 ) ;

– для тетрагонального и гексагонального материалов

= . т) = = 4 2 ) = 1 4 ( C

a a a a

) = (C

C 1122 ) S 4 .

В опыте на сдвиг при нагружении с S 5 = t 5 ( a 1 a 2 + a 2 a 1 ) возникают деформации:

– в тригональном материале

( тр )

= 5

t 5 ( C 1

C 1122 ) ( a 1 a 2 + a 2 a 1 ) + 2 1 5 C 112 3 ( a 1 a 3 + a 3 a 1 ) ,

– в тетрагональном материале е5т) = 215 C1212 (a1 a2 + a2a1) = 2C1212S5.

В гексагональном материале при напряжениях S 5 в плоскости векторов a 1 , a 2 деформациям отвечает тензор

= 5 2 ) = 2 1 5 C 1212 ( a 1 a 2 + a 2 a 1 ) = 1 5 ( C im C 1122 ) ( a 1 a 2 + a 2 a 1 ) = ( C im C 1122 ) S 5 ,                (24)

так как для гексагонального материала C 1212 = ( C 1111 C 1122)/2 [3].

Если канонические оси анизотропии k i совпадают с главными осями анизотропии a i , то в эксперименте на двухосное растяжение–сжатие тригональный материал отличим от тетрагонального и гексагонального по наличию ненулевых компонент деформаций б 23 (см. выражения (21), (22)).

В соответствии с (22)–(24) распознание тетрагонального и гексагонального материалов в канонических осях анизотропии можно выполнить по аналогии с кубическим и изотропным материалами, то есть использовать соотношения, похожие на (16), (19):

– если

(4)       (5)

Ь 11      ь 12

, t 4 t 5

то C 1212 = ( С цц C 1122)/2 , и материал гексагональный;

– если

(4)       (5)

_ 1*_ * Ь 12

t 4 t 5

то C 1212 * ( С цц C 1122)/2, и материал тетрагональный.

Таким образом, если два из трёх собственных значений тензора деформаций, возникающих при всестороннем сжатии, совпадают и при двухосном растяжении–сжатии в плоскости собственных векторов, соответствующих равным собственным значениям, возникают сдвиговые деформации в перпендикулярной плоскости векторов a 2 , a 3 , то материал тригональный. Полученные условия (21), (25) и (26) по данным двух экспериментов на двухосное растяжение–сжатие позволяют узнать, каким является исследуемый материал — тригональным, тетрагональным или гексагональным.

Анализ влияния погрешности измерений на точность выполнения критериев идентификации проведён в статьях [20, 21]. Показано, что при значениях погрешностей экспериментальных измерений, характерных для современной испытательной техники, отклонения от точных равенств в предложенных критериях (16), (25) составляют не более 5%.

5.    Моделирование задач, аналогичных экспериментам по определению типа анизотропного материала

Для компьютерного моделирования задач, соответствующих описанным экспериментам, разработан программный продукт [22], который осуществляет расчёт и визуализацию деформаций материального образца при различных видах нагружений, предусмотренных предложенными программами экспериментов. Ориентация канонических осей анизотропии относительно лабораторной системы координат и компоненты тензора упругих податливостей для анизотропных материалов различных типов в программе задаются, но считаются «неизвестными» для экспериментатора (пользователя). Результатом работы программы является тензор деформаций ε , который определяется из закона Гука (1) по заданным напряжениям на гранях куба. Программа предусматривает, что приложенные нагрузки и вызванные ими деформации измеряются с некоторой погрешностью, максимальная величина которой оговаривается изначально.

Приведём примеры моделирования экспериментов по распознанию типа начальной упругой анизотропии материала.

  • 5.1.    Моделирование базового эксперимента

Рассматривается кубический образец, рёбра которого направлены вдоль лабораторной системы координат с базисом e 1 , e 2 , e 3 . Проводятся три опыта на одноосное сжатие образца по программе базового эксперимента (2). Деформации измеряются с погрешностью Аб = 0,01 - 10 4 . Результаты представлены в Таблице 1.

Таблица 1. Опыты на одноосное сжатие – 1-й уровень номинальных напряжений

Опыт

Направление сжатия

Вид тензора напряжений

Ненулевая компонента напряжений

Заданные напряжения с учетом погрешности измерений, Па

Измеренные деформации

Обозначение

Номинальное значение, Па

1

Вдоль оси e 1

S 1 = S ,‘1) e' e'

(1)

S 11

- 10 7

S Л) =- 0,9918 - 107

б j =

Г- 1,59    0     0 '

0    0,58   0

v 0      0   0,58 ;

- 10 - 4

2

Вдоль оси e 2

S 2 = S 22) e2 e2

(2)

S22

S 22) =- 1,0056 - 107

(2)

Б У

г 0,59    0      0

0    - 1,60   0

0     0    0,59

- 10 - 4

3

Вдоль оси e 3

S 3 = S 3 3 3) e3 e3

(3)

S 33

S 33) =- 0,9957 - 107

Г 0,58   0     0 1

б j 3) =| 0    0,58     0

( 0     0    - 1,59 у

- 10 - 4

Компоненты тензора деформаций £0 = £1 + £2 + £3 являются откликом материала на всестороннее сжатие:

(0) б у

(1)      (2)      (3)

бу + Бij + бу

( - 0,42

I 0

- 0,44

0 1

0      - 10 - 4 .

- 0,42 J

Согласно (27) главные значения тензора деформаций составляют: б, =-0,42-10-4, б2 =-0,44 -10-4, б3 = -0,42 -10-4. С точностью до 5% их можно считать равными друг другу. Тогда исследуемый материал следует отнести к кубическому или изотропному. При этом главные оси анизотропии с достаточной степенью точности совпадают с лабораторными осями координат. При сжатии одинаковыми напряжениями вдоль различных лабораторных осей осевые, а также поперечные деформации одинаковы с точностью до погрешности измерений, а сдвиговые деформации меньше Аб . В этом случае, как показано в Разделе 2, канонические оси анизотропии материала совпадают с главными осями анизотропии и с лабораторными осями координат.

Для определения типа материала в соответствии с критериями (16), (19) проведём эксперименты: в направлении канонических осей анизотропии к 1 = e 1, к 2 = e 2 — на двухосное растяжение-сжатие; в проходящей через них плоскости — на сдвиг, который реализуем как растяжение–сжатие под углом 45 ° к ней.

Опыт 4 растяжение–сжатие вдоль осей e 1 , e 2 . Номинальный тензор напряжений имеет вид:

S 4 = t 4 ( e 1 e 1

e 2 e 2 ) , где 1 4 = 10 7 Па. Измеренные напряжения и деформации, соответственно, составляют:

s '1” = 1,0061 - 10 7 Па,

S 24) =- 0,9996 - 10 7 Па,

(4) б у

2,20

I 0

0    0 1

-2,18 0 - 10 - 4. 0    0 ,

Опыт 5 сдвиг в плоскости осей e 1 , e 2 . Номинальный тензор напряжений представляется как

S 5 = t 5 ( e 1 e 2 + e 2 e 1 ) , где 1 5 = 10 7 Па. Измеренные напряжения и деформации равняются:

S 1(25 = 0,9997 - 10 7 Па,

(5) б у

Г 0 0,70

0,70 0 A

0 - 10 - 4 .

0,

лоа 2,20 - 10 4    0,70 . 10 -4    _                             _

В соответствии с критерием (19) ---------^---------- поэтому исследуемый материал является

1,0061 - 10 7    0,9997 - 107

кубическим.

  • 5.2.    Моделирование базового эксперимента с другими номинальными значениями напряжений

Как и в предыдущем примере, проводятся три опыта на одноосное сжатие образца по программе базового эксперимента (2) с номинальными значениями напряжений S ® = S 22) = S 3(3) = - 10 6 Па. Погрешность измерения деформаций та же: Аб = 0,01 - 10 - 4 .

Таблица 2. Опыты на одноосное сжатие – 2-й уровень номинальных напряжений

Опыт

Направление сжатия

Вид тензора напряжений

Ненулевая компонента напряжений

Заданные напряжения с учетом погрешности измерений, Па

Измеренные деформации

Обозначение

Номинальное значение, Па

1

Вдоль оси e 1

S 1 = S (1) e 1 e 1

(1)

S 11

- 106

S 1(1) =- 1,0058 - 106

Г- 1,42 0,09   0 A

б j ° =1 0,09 0,40    0 I- 10 - 4

( 0      0    0,91 J

2

Вдоль оси e 2

S 2 = S ( 2 2) e 2 e 2

(2)

S22

S 22) =- 1,0063 - 106

Г 0,41 - 0,09   0 A

б ? =1- 0,09 - 1,41    0 I- 10 - 4

( 0      0    0,91 J

3

Вдоль оси e 3

S 3 = S 3 ( 3 3) e3 e3

(3)

S33

s 3 3 3 ) =- 0,9974 - 106

Г 0,91   0      0 A

Б (j 3) =1 0     0,91      0     I- 10^

( 0     0    - 1,87 J

Компоненты тензора деформаций £0 = £1 + £ 2 + £3, которые являются откликом материала на всестороннее сжатие, в этом случае составляют:

(1)      (2)      (3)

Бу + Бу + бу

Г- 0,1

- 0,1

0 A

0       - 10 4 .

- 0,05 J

Главные значения тензора деформаций найдем из (28): Б 1 = б2 =- 10 5, б3 =- 0,5 - 10 5. При этом главные оси тензора деформаций совпадают с лабораторными осями координат. Такой материал может быть тригональным, тетрагональным или гексагональным, поскольку все они имеют главную ось анизотропии, соответствующую однократному главному значению б3 , совпадающую с поворотной осью симметрии и, следовательно, с канонической осью k 3 .

Для определения угла поворота ф канонических осей анизотропии к 1, к 2 относительно лабораторных

осей e 1 , e 2 выполним эксперименты на двухосное растяжение–сжатие вдоль осей e 1 , e 2 и на сдвиг в плоскости этих осей.

Опыт 4 двухосное растяжение–сжатие вдоль осей e 1 , e 2 . Номинальный тензор напряжений

записывается как S 4 = t 4 ( e 1 e 1

e 2 e 2 ) , где t 4 = 106 Па. Измеренные напряжения и деформации равняются:

S 1 ( 1 4) = 1,0058 - 106 Па, S 24) = - 1,0063 - 106 Па,

б

Г 1,83

- 0,18

- 0,18 0 '

- 1,81 0 0    0 ?

- 10 - 4 .

Опыт 5 сдвиг в плоскости осей e 1 , e 2 . Номинальный тензор напряжений представляется как

S 5 = t 5 ( e 1 e 2 + e 2 e 1 ) , где 1 5 = 106 Па. Измеренные напряжения и деформации, соответственно, имеют вид:

S (5) = 1,0074 - 10 6 Па,

(5) ъ м

Г 0,18

1,83 ч 0

1,83

- 0,18 0

0 J

- 10 - 4.

Так как в опыте 4 считаем, что не возникают сдвиговые деформации е 13 и б 23 (их значения меньше погрешности измерения деформаций и поэтому не фиксируются), исследуемый материал не является тригональным.

Далее, определим угол ф по формуле (10): ф®- 0,391 рад. Погрешность найденного значения угла ф по сравнению со значением, заданным изначально в программе, составляет 2,3%. Угол ф отличен от нуля, поэтому исследуемый материал является тетрагональным.

Таким образом, численное моделирование экспериментов по идентификации типа начальной упругой анизотропии материала показывает, что разработанные программы экспериментов, при доступной точности измерений, способны правильно устанавливать тип анизотропного материала.

Работа выполнена при финансовой поддержке РФФИ (проекты № 13-01-97501-р_центр_а, 14-01-31138-мол_а, 15-01-01875-а) и Министерства образования и науки РФ (госзадание № 467).

Список литературы Численное моделирование экспериментов по определению типа начальной анизотропии упругих материалов

  • Грин А.Е., Адкинс Дж. Большие упругие деформации и нелинейная механика сплошной среды. -М.: Мир, 1965. -456 с.
  • Ландау Л.Д., Лифшиц Е.М. Теоретическая физика: Теория упругости. -М.: Наука, 1987. -Т. 7. -248 с.
  • Маркин А.А., Соколова М.Ю., Христич Д.В. Процессы упругопластического конечного деформирования. -Тула: Изд-во ТулГУ, 2011. -374 с.
  • Сиротин Ю.И., Шаскольская М.П. Основы кристаллофизики: учеб. пособие. -М.: Наука, 1979. -640 с.
  • Черных К.Ф. Введение в анизотропную упругость. -М.: Наука, 1988. -192 с.
  • Остросаблин Н.И. Линейные инвариантные неприводимые разложения тензора четвертого ранга модулей упругости//Динамика сплошной среды. -2002. -№ 120. -С. 149-160.
  • Цвелодуб И.Ю. К определению упругих характеристик однородных анизотропных тел//ПМТФ. -1994. -Т. 35, № 3. -С. 145-149.
  • Hayes M. A simple statical approach to the measurement of the elastic constants in anisotropic media//J. Mater. Sci. -1969. -Vol. 4, no. 1. -P. 10-14.
  • Jarić J.P. On the conditions for the existence of a plane of symmetry for anisotropic elastic material//Mech. Res. Commun. -1994. -Vol. 21, no. 2. -P. 153-174.
  • Norris A.N. On the acoustic determination of the elastic moduli of anisotropic solids and acoustic conditions for the existence of symmetry planes//Q. J. Mechanics Appl. Math. -1989. -Vol. 42, no. 3. -P. 413-426.
  • Маркин А.А., Соколова М.Ю. Определение типа исходной анизотропии и распространение частного постулата Ильюшина на начально анизотропные материалы//Устойчивость, пластичность, ползучесть при сложном нагружении: сб. науч. тр. -Тверь: Изд-во ТвГТУ, 2000. -Вып. 2. -С. 66-71.
  • Маркин А.А., Соколова М.Ю. Термомеханические модели обратимого конечного деформирования. -Тула: Изд-во ТулГУ, 2010. -268 с.
  • Новожилов В.В. Теория упругости. -Л.: Судпромгиз, 1958. -370 с.
  • Соколова М.Ю., Христич Д.В., Генералова Е.М. К вопросу об определении главных осей анизотропии//Современные проблемы математики, механики, информатики: материалы Международной научной конференции. -Тула: Изд-во ТулГУ, 2008. -С. 300-301.
  • Христич Д.В, Каюмов Р.А., Мухамедова И.З. Программа экспериментов по определению главных осей анизотропии материала//Известия КазГАСУ. -2012. -№ 3. -С. 216-224.
  • Ньюнхем Р.Э. Свойства материалов. Анизотропия, симметрия, структура. -М.-Ижевск: НИЦ «Регулярная и хаотическая динамика», 2007. 652 с.
  • Победря Б.Е. Механика композиционных материалов. -М.: Изд-во МГУ, 1984. -336 с.
  • Христич Д.В. К вопросу об определении главных осей анизотропии материала//Известия ТулГУ. Естественные науки. -2014. -№ 2. -С. 203-213.
  • Христич Д.В. Критерий экспериментальной идентификации ромбического, моноклинного и триклинного материалов//Известия ТулГУ. Естественные науки. -2013. -№ 3. -С. 166-178.
  • Христич Д.В. Критерий экспериментальной идентификации гексагонального, тригонального и тетрагонального материалов//Вестник КГТУ им. А.Н. Туполева. -2013. -№ 2-1. -С. 67-72.
  • Христич Д.В. Критерий экспериментальной идентификации изотропного и кубического материалов//Известия ТулГУ. Естественные науки. -2012. -№ 3. -С. 110-118.
  • Христич Д.В. Компьютерное моделирование экспериментов по определению типа начальной анизотропии упругих материалов//Известия ТулГУ. Естественные науки. -2014. -№ 4. -С. 110-119.
Еще
Статья научная