Упругие волны, излучаемые при торможении тела, движущегося по границе полупространства Текст научной статьи по специальности «Физика»

Электронный журнал «Техническая акустика» http://www.ejta.org

2019, 1

Ю. М. Заславский

Институт прикладной физики Российской академии наук.

Россия 603950, Н. Новгород, ул. Ульянова, 46, тел.(831)4164764, факс: (831)4365976. E-mail: zaslav@appl.sci-nnov.ru

Упругие волны, излучаемые при торможении тела, движущегося по границе полупространства

Получена 15.01.2019, опубликована 22.02.2019

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

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

ВВЕДЕНИЕ

Анализ упругих объемных и поверхностных рэлеевских волн, возбуждаемых при торможении твердого тела, падающего на земную поверхность, проводился в ряде исследований и связан, в частности, с проблемой интерпретации данных регистрации сейсмических откликов на метеоритное воздействие [1-3]. В работе [3] выполнены оценочные расчеты параметров сейсмического отклика при соударении падающего тела с поверхностью полупространства и последующем резком замедлении. При оценках сейсмического эффекта предполагалось вертикально направленное движение тела с высокой начальной скоростью, неупругое соударение с границей полупространства. Расчеты имели целью получить предварительные оценки параметров сейсмического отклика, в связи с чем использовалась огрубленная модель. Не принималось во внимание прохождение участка пути сквозь слой земной атмосферы, на котором принципиален учет «черенковского» механизма излучения акустических и сейсмических волн (при разогреве и взрыве паров в атмосфере). Не учитывались деформационно-пластические и тепловые явления (разрушение среды, внедрение тела вглубь и образование кратеров) при вхождении тела в толщу осадочных пород. Поэтому представленные данные дают завышенные оценки уровня сейсмического толчка, который прогнозируется по результатам количественных расчетов.

В настоящей работе в продолжение указанных исследований анализируются характерные особенности упругих волн, излучаемых в твердое полупространство за счет реакции на воздействие, касательное к свободной плоской границе со стороны тела, движущегося с резким замедлением прямолинейно вдоль нее. Моделируется ситуация, аналогичная предыдущей [3], когда не сгоревший в атмосфере метеорит, представляя твердую фазу, приземляется в некоторой точке сухопутной территории по касательной к поверхности (в чем и состоит отличие от ранее рассмотренного случая). Как и в предыдущем рассмотрении, не затрагиваются упомянутые выше эффекты, связанные с движением по атмосфере, со сложными физико-механическими явлениями при столкновении с твердой границей (рыхление и нагрев поверхностного слоя грунта и т. д.), требующие специального анализа. Реальные траектории тел, падающих на земную поверхность, имеют обе компоненты, т.е. их скорость направлена под произвольным углом к границе. Поэтому подход, в котором рассматриваются по отдельности случаи с нормальной (вертикальной) и с горизонтальной (касательной) ориентацией вектора начальной скорости при соприкосновении с земной поверхностью, представляет собой идеализацию, удобную для анализа. Для сравнения с предыдущими результатами важно получить оценки уровня вибросейсмического отклика (толчка), сопровождающего эффект горизонтально ориентированного приземления метеорита на ровную земную поверхность, моделируемую в простейшем случае плоской границей однородного упругого полупространства. Очевидно, что при замедленном движении тела, «скользящего» по твердой границе, кроме указанных потерь, доля кинетической энергии затрачивается на генерацию упругих поверхностных волн, локализованных в приповерхностной области грунта, и на излучение объемных волн в земную толщу.

1. СЕЙСМИЧЕСКИЙ ОТКЛИК, ОБУСЛОВЛЕННЫЙ ПРИХОДОМ ПОВЕРХНОСТНОЙ ВОЛНЫ РЭЛЕЯ

Кратко изложим решение задачи возбуждения упругих волн при воздействии на границу в течение короткого промежутка времени (в пределе - стремящегося к нулю) силы, противоположной силе трения (торможения), действующей на тело массы M, движущегося прямолинейно по плоской границе однородного твердого полупространства с начальной скоростью V < cR, где сд - скорость рэлеевской волны, ct - скорость поперечной волны, с - скорость продольной волны, и с равномерным замедлением, равным a (рис.1). В принципе, можно принять начальную скорость, равной скорости звука в воздухе, т. е. полагать, что V = сзвуш. Вместе с тем, к числу

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

движется равнозамедленно. Анализ сейсмического отклика проводится на основе рассмотрения модели однородного полупространства при упрощающем предположении о преобладающей роли поверхностных рэлеевских и поперечных объемных SV - (вертикальной поляризации) волн. В процессе торможения тела также должны возбуждаться SH - (горизонтальной поляризации) волны, вместо которых в реальных условиях, как правило, наблюдаются поверхностные волны Лява. Однако для их рассмотрения требуется модель среды со структурой, содержащей хотя бы единственный слой поверх полупространства с пониженной плотностью или скоростью упругих волн (для адекватного описания реальной слоистой структуры), что выходит за рамки настоящего подхода, базирующегося на упрощенной модели однородного упругого полупространства. В этой связи их анализ может быть выполнен отдельно. Роль продольных Р -волн при касательном воздействии пренебрежимо мала [4, 5].

SV - м>аув SV - м>аув

Рис. 1. Расчетная схема анализа

Запись интегральных выражений, описывающих всю совокупность волновых возмущений на границе полупространства, возбуждаемых силовым воздействием, является отправной при вычислении и анализе отклика. Она предполагает свертку спектра воздействия с функцией Грина задачи. Вычисление интегрального выражения (контурный интеграл) дает явную запись в виде суммы слагаемых, содержащих временную зависимость, характерную для каждого волнового типа, который определяется вкладом в интеграл, вносимым «окрестностью соответствующей особой точки», лежащей на контуре интегрирования [6-9].

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

промежутка времени постоянной силы, касательной к границе, при этом в задаче отсутствует фиксированная, навязанная источником частота осцилляций. Будем рассматривать случай Уа/а ^ 0, т.е. промежуток времени от начала торможения тела

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

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

те (г, г) - Ма8(х - уаг + аг 2/2Ж , (1)

где М — масса тела.

Спектральное разложение тензора напряжений по частоте (частотно-временное преобразование), вычисляемое как Фурье-компонента силового воздействия на поверхность полупространства при резком торможении тела, координата которого описывается как х = У0г - аг2/2, может быть записано:

ад ад ад а

Тх2(кх,ку,с) = Ма\ 11вш-*хх-*'у8{х-У0г + аг2 /2)5(у)йхйуйг = Ма}2/2л . (2)

0 -ад-ад 0

Представив интеграл (2) в виде, удобном для приближенных оценок,

а

Я(кх, ку, с) = |ег(ш-кху°) (1 + 1акхг2 / 2>Й,

0

для спектра воздействия можно получить выражение:

ТХ2 (кх, ку ,с) = МаЯ (кх, ку ,с),

Я (кх, ку ,с) = 1

kV (3)

\ a 2a a

1 - cos("- 'УоГо + kZLsn^-kloV a 2a a

где kx = k cos p, ky = k sin p, k - проекция волнового вектора на границу z = 0, р - азимутальный угол.

Проекция волнового вектора k является переменной интегрирования в интегральных выражениях, с помощью которых далее вычисляются волновые смещения.

В формуле для спектра силового воздействия фигурируют слагаемые с точностью до квадратичных значений относительно оУ^/а, а при последующих вычислениях в

предположении кУ^ /а << 1 будут удержаны члены, содержащие первую степень, и

отброшены слагаемые с более высокими степенями.

Для описания поля упругих смещений в поверхностной волне далее необходимо использовать функцию Грина (для поверхностной Я -волны), с помощью которой может быть описано решение задачи волнового возбуждения полупространства переменной силой, касательной к границе, воздействующей в течение короткого промежутка времени. Явный вид функции Грина представлен в [10-13]. При свертке спектра силового воздействия по частоте и по волновым числам с указанной функцией, можно получить интегральное выражение, описывающее волновые смещения на свободной границе.

Вертикальная компонента упругих волновых смещений в общем случае может быть представлена как:

ад ад ад

и* (X, 7, * = 0, г) = 111 (кх, ку, со)ТХ2 (кх, ку, с)е+1кхХ+,куйкхйкуйо, (4)

— ад—ад—ад

где X = гсовщ, У = гвтщ - координаты точки наблюдения на плоскости * = 0,

(кх, ку ,о) - функция Грина. Указанные координаты много больше координат, соответствующих местоположению тела X >> х, У >> у, кроме того, предполагается дальняя зона кхХ >> 1, к У >> 1, где

применимы асимптотические оценки интегралов.

С учетом явного вида упомянутой функции Грина формула (4) может быть преобразована:

- iMV 777 k (2k2 - к? - 2v,vt )e +'kyJS (к , к ,é)

uz (X,Y, z = 0, t) = -=MVL f f f kí^-t----( x' y' )dkxdkyda, (5)

zV 7 pc? JJJ (2k2 - k2)2 - 4k 2v,vt x y

~ t -7 -7 -7 V t ' l t

гда v, k 2 - kj , vt k 2 - kt2 .

При вычислении интегралов в (5) предполагается kr >> 1 и используются следующие элементарные соотношения:

dkx = dk cos р- k sin (pdp dky = dk sinp + k cospdp

Это позволяет представить формулу (5) в развернутой записи:

dkxdky = kdkdp.

-iMV0 772¡k2(2k2 -k2 -2v¡vt)e-'^+ikrcos(-v)S(k,p,ffl)cospj;j ,

uz(X,Y,z = 0,t) =-I I I —-t-' tJ , ,---v ' '-t-dkdpdrn . (6)

z ( ) pc] -J7J0 f (2k2 - k2)2 - 4k\vt ( ()

Интегрирование по азимутальному углу р выполняется методом стационарной фазы:

2— 2— , кг 2 I , —

Г '1 / \ Г ¡кг-1—р Л 1кг-1—

|соБре с°8(р-^)ёр = соб^!е 2 ( 2А— соБЩв 4 . (7)

0 0 ' кг

С учетом проведенного интегрирования формула (6) преобразуется следующим образом:

и = -. 24—МУо с°8^ Ц к2(2к2 -к? -2у1у1 (к¿кг-шдёп. (8)

2 ре] Ц ((2к2 - к,2)2 - 4^,к2)4кг

На рис. 2 представлен замкнутый контур на плоскости комплексной переменной к, используемый для интегрирования в (8), о чем упоминалось выше, а также нанесены особые точки (полюс кК и две точки ветвления к1, к), а также выходящие из них

линии разрезов, обеспечивающих однозначность подынтегрального выражения на каждом Римановом листе, что строго обосновано в [6-9].

■ .................

1т к

V

^ ! Ь,;.................

к, к. кК Яе к

0

Рис. 2. Контур (пунктир) для вычисления интеграла по волновому числу к

При вычислении интеграла в (8) по переменной к необходимо взять вычет в точке полюса кК = а/еЕ и получить выражение, соответствующее вкладу, обусловленному поверхностными рэлеевскими волнами:

у—МГ0 с°в^} М(2к2К - к2 - 2у1у, ук*г-а8К (кК ,щ,а)

-1 -г ёа .

4ре,-л/2ге^ 1( 2 ,2 „,,4к'й - ЗкЖ + к;) + 2к;к; (9)

I 5- " (2*5-кг)2 J

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

и2 =

acR cosy

-т--)a

cR V0

V„

)(1—0cosy)

ie—

-ад

(sin % + i(1 - cos^))

w, = ■

M

f

8pc4V2

V j (1 —0cosy)2

2

2 - % - 2-c t

22 4 - 3| CR + %

ci c

t у

+ 2

22 С ct

f 2 Л 2

2 - С2

v ct у

(10)

В результате вычисления получается формула, описывающая осциллограмму отклика как функцию времени, в виде «интегрируемой особенности», имеющую знакопеременный характер:

c

R

v

у

3

inMa2c\ cosy

w, =

С c 2^

2 - С1,

v ct у

16 pc4 ^ cRV0r (1 - Vo cos y)

2

2 - - 22 cR 2

4 - 3

- +

V c

2 "R -2 t у

c

+ 2 cR

22 ct ct

f ri\2 v ct у

(11)

if

t -

Vv

-R у

'V V Vo(1 - cos y)

i

Y2

a

cn

if

Vv

V

a

V

(1 — cosy)

VV Vo(1 - ^cosy)

1

Y1

a

c

R

R

v

r

r

t

>

c

c

R

R

причем в двух корневых слагаемых в (11) предполагаются ненулевые значения при t > г/ск в первом слагаемом (в фигурных скобках) и при t > г/ск + (V /а)(1 - (Уо/ск )ео8у) - во втором, а при задержках, меньших указанных, отклик является нулевым (см. рис. 3 а). Интервал задержек г = (К0/а)(1 - (¥0 /ск )ео8у)

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

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

а б

Рис. 3. Осциллограммы волновых откликов (качественный вид) а - рэлеевская .R-волна, б - поперечная ^F-волна

По формуле (11) выполнены количественные оценки уровня сейсмического отклика, выраженного в единицах колебательной скорости смещений грунта vz — (a/V0 )uz, что

обычно принято в расчетах. Однако при выводе формулы (11) не учтена диссипация энергии при распространении волны, что ведет к огрублению оценки в сторону ее завышения. Так, при массе тела M —10т, начальной скорости V0 - 340 м/с,

замедлении a —100 м/с2 на дистанции 1 км от начальной точки касания на поверхности грунта z — 0 (р — 2г/смъ, е{ — 2000 м/с, ct — 577.35 м/с ) возникнет импульсный волновой отклик с уровнем ~ 4- 5 мкм/с, имеющий характерную (эффективную) длительность т = VQ/a — 5c. Наблюдение такого сигнала в условиях надежного превышения над уровнем микросейсм можно обеспечить с использованием одиночного низкочастотного сейсмоприемника без использования пространственного накопления на дистанциях до нескольких километров (без применения сейсмокосы).

2. СЕЙСМИЧЕСКИЙ ОТКЛИК, ОБУСЛОВЛЕННЫЙ ПРИХОДОМ ПОПЕРЕЧНОЙ SV-ВОЛНЫ

Расчет и оценка колебательного отклика в предположении доминирования поверхностной рэлеевской волны в пренебрежении диссипативного затухания могут рассматриваться только как предварительные и для незначительных удалений от места остановки тела. Для более удаленных дистанций может представлять интерес анализ возможности приема упругих колебаний в виде ^F-волны (объемная волна вертикальной поляризации), которая по сравнению с поверхностной проникает в более глубокие слои земной толщи с пониженными поглощающими характеристиками и далее распространяется в градиентной среде по лучевым траекториям. При этом не исключается возвращение волнового луча на нулевой горизонт, т. е. на свободную земную поверхность, при котором открывается перспектива ее регистрации приемниками, установленными более удаленно, чем в расчете на прием рэлеевской

волны. В пренебрежении искривлением луча можно оценить уровень отклика по £К-волне, вновь используя модель однородного полупространства. В этом случае следует учитывать поглощение волны по мере распространения, что нетрудно сделать на последнем этапе расчета.

Упругие смещения ив в сдвиговой ЗУ -волне (в вертикальном сечении и с

ориентацией по в, перпендикулярной радиусу-вектору) можно вычислить по схеме, рассмотренной применительно к Я -волне. Необходимо воспользоваться функцией Грина и спектром воздействия, пересчитанными на случай возбуждения волны указанного типа [10-13]:

/Ма

ив = — РС

Ш +^ —-к.

-ю-ю-ю к \\2к кг

к\\2к2 - к/)2 + 4к,К,к2)

\а - к(У0 в соъщ)У{) ( \а - к(У0 в соъщ)У{)

--—-+ /11 - СОБ

а - - 1 а V а

где ЗЗУ к ,щ,в,т) =

\а - к У Бт в СОБ^)

12

Интегрирование по кх. к выполняется подобно уже рассмотренному, причем ключевую роль, определяющую вклад в интеграл, вносит окрестность вблизи точки ветвления к,. Это позволяет записать:

(

тМа сощсоъв

2 БШ2 в - 1 + 2сО8вБШ в

иа =■

рс2 Я

\2мп2в-1)2 + 4 соБв бш2 в

V

\ - мп2 в

У

1 —-БтвсоБ^ V С

Л

Я . а | а ©Я

-1%\г—)--£--

С У0 У0 2 с У0

\1--соб^ Бт в

ю е с'

(бШ^ + /\1 - СОБ^))

£

(13)

где Я - радиус-вектор в точку наблюдения в сферической системе координат.

В экспоненциальный множитель подынтегрального выражения введен декремент © для учета поглощения при распространении £К-волны.

Интегрирование по частоте в (13) проводится так же как представлено в (10):

2

-ю

лМа cosy cos в

2sin в-1 + 2cosesine

Ид =-

pc2 R

(2sin2 в -1)2 + 4cosesin2 в

~ R 2ct (t--) ct

arctg > - arctg'

0R

cL c1

R Vl

ct a

- sin2 в

V

V

1 —0 sin в cosy

V c

(14)

0R

причем ненулевые значения в первом слагаемом в фигурных скобках в (14) предполагаются при t > R/ct, а во втором - при t > R/ct + (V0 /а)(1 - (V0 ¡ct) cos y sin в).

По формуле (14) построена осциллограмма сейсмического отклика в поперечной SV -волне (см. рис. 3б). Сравнение осциллограмм на рис. 3а, б показывает, что принципиальным отличием откликов в поверхностной и в объемной волнах является знакопеременный вид, при котором хвостовая часть как бы «отделяется» от основного тела импульса в поверхностной волне в противоположность случаю импульсного отклика в объемной волне, который характеризуется однополярной зависимостью. Последнее обстоятельство также указывает на своеобразие, заключающееся в знакопостоянном характере отклика по объемным сейсмическим волнам, отличающем их от двуполярного акустического отклика, в который преобразуется любое исходное возмущение с первоначально различной импульсной формой по мере распространения с удалением на бесконечность в однородной безграничной атмосфере.

Заметим, что если начальная скорость V0 > cñ, то при торможении тела,

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

c

ЗАКЛЮЧЕНИЕ

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

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

Работа выполнена в рамках государственного задания ИПФ РАН (проект 0035-2019-0019).

ЛИТЕРАТУРА

1. Добрынина А. А., Саньков В. А., Чечельницкий В. В., Черных Е. Н. Сейсмические эффекты метеорита Челябинск. В кн.: Метеорит Челябинск - год на Земле. Материалы Всероссийской научной конференции / Составитель А.Н. Антипин. Челябинск: Каменный пояс, 2014. С. 668-677.

2. Селезнев В. С., Лисейкин А. В., Еманов А. А., Белинская А. Ю. Метеороид Челябинск (взгляд сейсмолога). В кн.: Метеорит Челябинск - год на Земле. Материалы Всероссийской научной конференции / Составитель А.Н. Антипин. Челябинск: Каменный пояс, 2014. С. 678-681.

3. Заславский Ю. М., Заславский В. Ю. К оценке сейсмического эффекта при падении на Землю метеорита. Техническая акустика. Электронный журнал http://www.ejta.org 2011, 4.

4. Гущин В. В., Докучаев В. П., Заславский Ю. М., Конюхова И. Д. О распределении мощности между различными типами излучаемых волн в полубезграничной упругой среде // Исследование Земли невзрывными сейсмическими источниками. М.: Наука. 1981. С. 113-118.

5. Заславский Ю. М. Нагрузка с циклическим перемещением точки воздействия как источник сейсмических волн-гармоник / Физика земли 1994. №2. С. 52.

6. Ewing E., Jardetsky W., Press F. Elastic waves in layered media. N-Y.: McGrow-Hill, 1957. 454 p.

7. Аки К., Ричардс П. Количественная сейсмология. Т.1, 2. М.: Мир, 1991. 880 с.

8. Петрашень Г. И. Распространение волн в упругих средах. Л.: Наука, 1980. 324 с.

9. Miller G. F., Pursey H. The field and radiation impedance of mechanical radiators on the free surface of semi-infinite isotropic solid // Proc. Royal Soc., Ser.A, 1954, V. 223, №1155, P. 521-541.

10. Шемякин Е. И., Файншмидт В. М. Распространение волн в упругом полупространстве, возбужденном поверхностной касательной силой // Ученые записки

ЛГУ, №177, Сер. Математические науки. Вып. 28. Динамические задачи теории упругости, IV, 1954. С. 148-179.

11. Krylov V. V. Generation of ground vibrations by accelerating and braking road vehicles. ACUSTICA - acta acustica vol.82 (1996) 642-649.

12. Заславский Ю. М. Волны Лява, возбуждаемые движущимся источником / Акуст. ж. 2016. т. 62, № 1, с. 87-94.

13. Заславский Ю. М. Излучение сейсмических волн вибрационными источниками. Изд-во ИПФ РАН Н.Новгород, 2007. - 200 с.