Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

m1012

.pdf
Скачиваний:
1
Добавлен:
15.11.2022
Размер:
10.66 Mб
Скачать

Так как для непьезочувствительных материалов n + n0 = n,

где

 

 

 

 

 

 

 

 

 

n n1

n2

C( 1 2 3 ),

(3.31)

то формулу (3.21) можно представить в виде:

 

 

 

C

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( R)

 

 

 

 

 

( 1 2 3 )sds,

(3.32)

x

n

 

x

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

т.е. фиктивные деформации

 

( R)

пропорциональны абсолютному

x

 

пьезооптическому коэффициенту.

Отсюда следует, что при создании идентичных условий нагружения на двух моделях, для которых упругие константы (E,) равны, а пьезооптические коэффициенты разные (С(1), С(2)), можно получить независимую систему уравнений с двумя неизвестными [36]:

 

(1)

 

N

 

2

 

(

 

 

( R)

),

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

0

 

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.33)

 

 

 

N

 

 

 

 

 

 

 

C

(2)

 

 

 

(2)

 

2

2

 

x

 

 

(1)

x

 

 

0

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

( R)

 

 

 

 

 

x

 

 

 

 

 

 

C

 

 

 

 

 

Решая эту систему уравнений, определим

 

x

( R)

). Есте-

 

x

ственно, что устойчивость такой системы к ошибкам измерений существенно зависит от отношения C(2) / C(1). Очевидны также затруднения при реализации этого метода, связанные с получением идентичных моделей и условий их нагружения.

3.6.2. Метод дифракционного соответствия

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

Абсолютные пьезооптические коэффициенты, которые в соответствии с соотношениями Максвелла–Неймана определяют изменение показателя преломления при деформировании, обладают дисперсией, т.е. зависят от длины волны света . Величина фиктивных деформаций (xR) , как было показано выше, пропор-

61

циональна абсолютному пьезооптическому коэффициенту С. Это позволяет, аналогично предыдущему методу, получить систему

уравнений для определения (и разделения) деформаций

 

 

( R)

.

 

x

и x

Для деформированного сечения пространственной модели при двух длинах волн в одном и том же интерферометре получим систему уравнений [36]:

 

(1)

 

N

 

 

2

 

 

 

 

 

( R)

,

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

0

 

x

 

x,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

N

2

 

2

0 x

 

( R)

,

 

 

 

(3.34)

 

x

 

x,2

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

( R)

C

( R)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x,2

 

 

x,1

 

 

 

 

 

 

 

 

 

 

 

 

 

Решая эту систему уравнений,

находим x

и

( R)

.

При этом

x

достаточно одной модели и не возникает проблем с обеспечением идентичности нагружений.

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

Возможности применения этого метода были рассмотрены при определении деформаций однороднодеформируемого растра. При этом интерферограмма растра была получена через деформируемый объем (см. рис. 3.1). Чтобы обеспечить равенство углов дифракции, были использованы He-Ne-лазеры ( = 0,6328 мкм) и He-Cd-лазеры ( = 0,4416 мкм). Установленные в фурье-плоскости 1-й линзы диафрагмы пропускали соответственно плюс-минус вторые дифракционные порядки ( = 0,6328 мкм, m1 = 2) и плюсминус третьи ( = 0,4416 мкм, m2 = 3). Это обеспечивало условие дифракционного соответствия – почти точное равенство углов ди-

фракции, так как

 

 

 

 

 

m1 1

 

m2 2

.

(3.35)

 

 

 

 

n1

 

n2

 

Коэффициент С = C(2)/ C(1) определялся предварительной тарировкой.

62

Ошибки в определении деформированого соcтояния растра при толщине призмы L = 40 мм достигали ~25 %. Впрочем, без учета изменений n они уже на расстояниях, меньших ~10 мм от точки приложения силы, превышали 100 %.

3.6.3.Многоракурсное просвечивание

вприближении малых углов дифракции

При малых углах дифракции лучей в исследуемом сечении модели можно считать, что производная от grad n на линии наблюдения мало отличается от ее значения на линии просвечивания. Если к тому же он слабо меняется по направлению просвечивания, то величина фиктивных деформаций равна (3.22):

 

 

 

L

 

 

 

 

n

 

 

 

 

2

 

 

 

2

 

 

.

( R)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2n

cos

2

 

 

x

2

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

Величина второй производной от n в формуле (3.22) может быть определена при наличии дополнительной информации. Например, при наличии интерферограммы абсолютной разности хода (в общем случае – нескольких таких интерферограмм, полученных при многоракурсном просвечивании) [17, 230, 231]. Действительно, изменение разности хода в этом случае равно:

 

0

 

L

0

 

L

0

 

 

 

 

 

 

 

 

 

(n

n)dz

 

n dz nL,

 

 

 

0

 

 

0

 

(3.36)

где n = n – n0, тогда

( )

 

( n)

L

n

L,

0

 

 

 

 

 

 

 

x

 

x

 

x

 

(3.37)

т.е. дифференцирование поля

 

0

 

позволяет найти x-компоненту

градиента n:

 

 

 

 

 

 

 

n

 

1

( )

,

 

 

 

 

0

 

 

 

 

 

 

 

 

 

x

 

L

x

 

 

 

 

откуда

 

 

 

 

 

 

 

2n

 

1

2 (

0

)

.

x2

L

x2

 

 

 

 

 

 

(3.38)

(3.39)

63

Подставляя последнее выражение в (3.22), получим

 

 

 

L

 

 

 

2

( )

.

 

 

 

 

 

 

 

 

( R)

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

x

 

2n cos

2

 

 

 

x

2

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

(3.40)

Сравнение полученных экспериментальных результатов для задачи Фламана и вычисленных в соответствии с (3.40) показыва-

ет, что при

размерах зоны неоднородного деформирования

L < 40 мм и

регистрируемых пространственных частотах

0 < 100 мм-1

пользование таким приближением не приводит к

появлению существенных ошибок.

3.7. Тестовый пример: полуплоскость под действием плоского штампа

При действии плоского штампа на границу полуплоскости [64] (рис. 3.6, а) изменение показателя преломления для изотропного материала (С1 С2 = С) вдоль оси x равно:

n n n

C(

)

2PC

.

 

 

 

1

2

1

2

 

x

2

a

2

 

 

 

 

 

 

 

 

 

 

 

(3.41)

Рассмотрим, как меняются регистрируемые интерферограммой параметры при просвечивании этой полуплоскости по нормали к ее поверхности.

Если толщина модели после исследуемого сечения L, то изменение оптического пути 0 на этом участке и его производных при нормальном просвечивании равны:

0

0

 

 

 

L

 

 

 

 

 

 

 

L

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

)

 

 

nds

 

 

C(

)ds

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2PC

 

 

 

 

L

 

 

 

 

 

 

2PCL

 

 

 

 

 

 

 

 

 

 

 

 

ds

 

 

,

 

 

 

 

 

x

 

 

a

 

 

 

 

 

x

 

a

 

 

 

 

 

 

 

2

 

2

0

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 0 )

 

 

 

 

 

2PCLx

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x2 a2 )3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 (

)

 

 

2PCL

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

3x2

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

.

x

2

 

 

 

 

(x

2

a

)

3/ 2

 

x

2

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

(3.42)

(3.43)

(3.44)

64

а)

 

P

L

2a

y

 

 

 

Лазер x

б)

Сечение: y = 0

в)

Материал: ОНС

а = 5 мм0 = ±120 мм-1

= 0,633 мк P = 25 Н/мм L = 40 мм

Е = 3,2 103 Н/мм2 С = 3,9 10-5 мм2

г)

д)

Рис. 3.6. Влияние рефракции на точность определения деформаций

взадаче о штампе при плосконапряженном состоянии

Всоответствии с (3.40) величина фиктивных деформаций

(xR):

(xR)

2PCL2

 

 

 

1

 

 

 

 

 

 

3x2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

.

(3.45)

n0

 

(x

2

a

2

)

3/ 2

x

2

a

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

65

Эта зависимость приведена на рис. 3.6, в. На рис. 3.6, б представлены две интерферограммы [54, 176–179]: первая (правая часть фотографии) – интерферограмма абсолютной разности хода при нормальном просвечивании, использована для определения

( R)

;

вторая (левая часть фотографии) – результат интерференции

x

двух ( 2) дифракционных порядков, характеризует распределение

( )

x

( R)

.

Результаты обработки этих интерферограмм для

x

x

сечений y = 0 и x = 2a приведены на рис. 3.6, в, г, д.

Деформации исследуемого сечения x, определяемые как разность между этими двумя зависимостями, хорошо согласуются с аналитическим решением. Погрешность в определении деформаций в зоне экстремума функции x (x) для сечения y = 0 не превышает 7 %.

Следует отметить, что чем меньше деформации сечения и чем больше размеры зоны неоднородного деформирования, тем больше (без учета рефракции) погрешность в определении деформа-

ций. Это

обстоятельство

хорошо иллюстрируют зависимости:

( )

( y) y

( R)

( R)

( y), y

( y) для сечения x = 2a на рис. 3.6, г, д.

y

y

, y

При тех же примерно размерах зоны неоднородного деформирования, но существенно меньших значениях деформаций сеченияy (по сравнению с x) в зоне под штампом суммарные деформа-

ции более чем в три раза превышают деформации исследуемого сечения.

Поскольку чувствительность интерферометра к внутриплоскостным деформациям регулируется изменением пространственной частоты регистрируемой световой волны, то возникает естественный вопрос: каким образом эти изменения отразятся на величине фиктивных деформаций? Получить однозначный ответ на него затруднительно, потому что изменение пространственной частоты волны, связанное с рефракцией в неоднородно деформируемом объекте, зависит от нескольких факторов.

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

66

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

Во-вторых, при определении фиктивных деформаций величина производной от разности хода этих лучей, в соответствии с (3.15), нормируется на величину регистрируемой пространственной частоты. Эта нормировка при увеличении пространственной частоты уменьшает вклад от рефракции. Данная особенность успешно используется для уменьшения влияния рефракции, например, при исследовании температурных задач [196, 197].

На основании изложенных причин указать заранее, каков будет итог влияния вышеперечисленных факторов при изменении пространственной частоты регистрируемой волны, в общем случае, невозможно. Однако, как по-

казывают численные

экспери-

 

менты, зависимость эта довольно

 

слабая.

В качестве

иллюстрации

 

на рис. 3.7 приведена зависи-

 

мость

фиктивных

деформаций

 

( R)

(x)

в сечении y = 0

от про-

 

x

 

странственной частоты волн, ди-

Рис. 3.7. Зависимость фиктивных

фрагировавших в

вертикальной

деформаций от пространственной ча-

плоскости xoz.

 

 

стоты регистрируемой волны

 

При малых частотах дифрагировавших волн ( < 100 мм-1)

изменение пространственной частоты регистрируемой световой волны почти не отражается ни на величине, ни на характере распределения фиктивных деформаций.

На частотах ~300 мм-1 и более уже заметно уменьшение величины фиктивных деформаций. Кроме того, с ростом пространственной частоты происходит смещение экстремума и постепенное выглаживание функции.

Таким образом, использование стандартного приема – уменьшения чувствительности к влиянию рефракции за счет повышения чувствительности к деформациям [196, 197] – не дает

67

желаемого эффекта. Изменения происходят почти синхронно. Кроме того, при увеличении пространственной частоты возникает и проблема «мертвых зон», связанная с полным внутренним отражением света от поверхности объекта. Поэтому в дальнейшем для регистрации, как правило, используются волны с частотами, не превышающими ~300 мм-1.

3.8. Численные эксперименты при исследовании пространственных объектов

Как было ранее отмечено, фиктивные деформации в общем случае определяются соотношением (3.32):

где

 

 

 

( R)

 

 

 

 

 

 

 

 

x

 

I (x, y, z)

2

 

1

1

 

3

C

s

 

2

I (x, y, z)

 

 

 

 

 

 

 

 

 

1

2

sds,

(3.46)

 

 

 

 

 

n

 

 

 

x

 

0

 

 

 

 

 

 

 

 

 

 

 

 

– первый инвариант тензора напряже-

ний.

В условиях сложного объемного деформирования величина фиктивных деформаций будет определяться положением исследуемого сечения в модели, геометрией схемы регистрации и может носить довольно сложный характер.

Проследим зависимость фиктивных деформаций от положения исследуемого сечения в случае действия плоского, круглого в плане штампа на границу полупространства [50, 177].

На рис. 3.8 представлена зависимость первого инварианта тензора напряжений и его вторых производных для нескольких вертикальных сечений (y = const), лежащих в плоскости zoy (P = 4 кН, a = 10 мм). Характерной, уже встречавшейся ранее при рассмотрении деформированого состояния двумерного штампа, особенностью является изменение знака вторых производных I1(x, y, z) для некоторых сечений. Однако в отличие от плоской задачи их величина меняется еще и по направлению просвечивания.

Рассмотрим динамику изменений локальных значений фиктивных деформаций, предположив, что направление просвечивания перпендикулярно указанным сечениям (т.е. параллельно oy). Для этого, считая z = const параметром, будем переходить «от удаленных сечений со значениями y >> a к последующим (y → 0) по направлению просвечивания. Так, для z = a/2 вторые производные, а следовательно, и локальные значения фиктивных деформаций сначала малы и уменьшаются, достигая нулевого значения

68

при приближении к значениям y = a. Далее они меняют знак и резко возрастают по модулю, затем уменьшаются и, повторно изменив знак, возрастают до экстремального значения при у = 0. Поведение фиктивных деформаций для последующих значений y определяется с учетом их симметрии относительно значения y = 0. Конечное значение величины фиктивных деформаций определяется как сумма их локальных значений по направлению просвечивания в зоне только после исследуемого сечения. Поэтому эта величина будет существенно зависеть от положения сечения внутри объекта.

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

I1(x, y, z), МПа

Плоскость x = 0 Сечения: y = const

d 2I1 /dz2, МПа/мм2

Рис. 3.8. Распределение первого инварианта тензора напряжений I1(x, y, z) и его второй производной d 2I1 /dz 2 в задаче действия плоского,

круглого в плане штампа на границу полупространства

69

На рис. 3.9 представлены эпюры деформаций для нескольких сечений, а также величины фиктивных деформаций для них. Для этого разность хода между дифрагировавшими в этих сечениях лучами определялась численным интегрированием по направлению их распространения; пределы интегрирования: от исследуемого сечения до 10а (а радиус штампа). Величины фиктивных деформаций определялись в соответствии с (3.15) дифференцированием этой разности.

Рис. 3.9. Зависимость величины фиктивных деформаций от положения исследуемого сечения

Анализ полученных результатов показывает, что величина фиктивных деформаций существенно зависит от положения исследуемого сечения как относительно зоны резких неоднородностей, так и относительно зоны регистрации. В частности, для симметричных относительно штампа сечений (например, y = а,1,25а) фиктивные деформации могут не только отличаться в несколько раз, но и иметь разные знаки. Для сечений y = 1,25а в зоне z > 1,5a фиктивные деформации составляют ~30 % от z, но

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]