doi: 10.5862/MCE.53.2
Определение остаточной несущей способности железобетонных балок на стадии эксплуатации по критерию
прочности арматуры и бетона
Д.т.н., профеесор В.С. Уткин,
Вологодский государственный университет
Аннотация. Экспериментально рассмотрен теоретический метод оценки остаточной несущей способности индивидуальной железобетонной балки на стадии эксплуатации по критериям прочности рабочей арматуры и прочности бетона сжатой зоны балки. Использованы интегральные методы испытаний балки и вероятностные методы описания случайных величин.
В качестве меры несущей способности принята предельная нагрузка в интервальном виде на стадии эксплуатации, приведен анализ выбора ее значения из полученных границ интервала. Рассмотрен переход от экспериментально полученного значения предельной нагрузки к предельной эксплуатационной нагрузке.
Приведены примеры расчетных схем балок с эксплуатационной нагрузкой и указания по определению их предельных значений.
Ключевые слова: несущая способность; железобетонная балка; прочность арматуры; прочность бетона; предельная нагрузка; расчетная схема
Введение
Железобетонные балки (и балочные плиты) входят в состав большинства зданий и сооружений. Одной из мер их качества и конструкционной (механической) безопасности эксплуатации служит остаточная несущая способность по различным критериям: прочности арматуры и бетона, жесткости, трещинообразованию, ширине раскрытия трещин. Под несущей способностью элементов конструкций понимается их способность выдерживать наибольшие нагрузки, обеспечивая безопасное функционирование конструкции на стадии эксплуатации. По стандарту ГОСТ Р 54257-2010 «Надежность строительных конструкций и оснований» под несущей способностью понимается «максимальный эффект воздействия, реализуемый в строительном объекте без превышения предельных состояний». Мерой несущей способности является значение наибольшей нагрузки (воздействия), которую способна выдержать конструкция по критериям работоспособности, не доходя до состояния отказа. По межгосударственному стандарту ГОСТ 31937-2011 «Здания и сооружения. Правила обследования и мониторинга технического состояния», вступившего в силу с 01.01.2014 г., мерой несущей способности конструкции является наибольшая нагрузка, при которой отсутствует недопустимый риск разрушения.
Из приведенной информации следует, что основной задачей при оценке несущей способности конструкций является определение значения наибольшей, предельной нагрузки или воздействия. В процессе эксплуатации конструкций несущая способность несущих элементов понижается [1] вследствие старения (деградации) материалов, образования неисправностей и повреждений и т. д. Соответственно, уменьшается значение предельной нагрузки на нее по условиям безопасности эксплуатации. Это приводит к снижению надежности, сокращению остаточного временного ресурса и других показателей, характеризующих безопасность зданий и сооружений, что свидетельствует об актуальности обсуждаемой проблемы - оценки остаточной (оперативной) несущей способности железобетонных конструкций, в частности балок, на стадии их эксплуатации.
Постановка задачи
Для определения несущей способности железобетонной балки необходимо выявить наименьшую предельную нагрузку по всем критериям ее работоспособности (по всем видам предельных состояний). Критериями для железобетонных балок по действующим нормативным документам являются прочность арматуры и бетона, жесткость балки, трещинообразование в растянутой зоне и сжатых областях бетона и ширина раскрытия трещин. В настоящей статье предлагается рассмотреть проблему определения остаточной несущей способности только по критерию прочности арматуры и бетона.
Обзор источников по проблеме
В СП 13-102-2003 «Правила обследования несущих строительных конструкций зданий и сооружений», а также в стандартах ГОСТ 31337-2011 и ивО 13882:2010 (Е) при оценке несущей способности конструкций по тому или иному критерию работоспособности предписывается проводить теоретические поверочные расчеты по результатам неразрушающих испытаний механических и физических характеристик материалов конструкций на стадии эксплуатации. В ГОСТ 31937-2011 приведена информация о видах дефектов и повреждений в железобетонных конструкциях и о возможном снижении их несущих способностей. В работе [2] приведен теоретический анализ несущей способности железобетонных балок после их частичного разрушения по критерию прочности бетона, а также приведена методика экспериментальных исследований балок. В работе [3] исследуется несущая способность железобетонных балок при изгибе по трещинообразованию с использованием теории механики разрушения.
Расчетам несущей способности эксплуатируемых железобетонных конструкций с дефектами и повреждениями посвящена работа А.Н. Бедова и В.Ф. Сапрыкина [4], где также приводятся методики испытаний бетона и арматуры в условиях эксплуатации железобетонных конструкций и статистической обработки результатов измерений. В [5] рассматриваются методы испытаний железобетонных конструкций. В работе [6] предложен способ неразрушающего определения несущей способности металлических конструкций, а в [7] - железобетонных конструкций. В качестве контролируемых параметров в [7] используются деформации арматуры и бетона
£ь. Однако в качестве предельных значений деформаций приняты средние значения без учета
изменчивости этих случайных величин. Так, принято еъ д = 0,002 и „р = К, /Е3 , где -
расчетное сопротивление арматуры, - среднее значение модуля упругости стали арматуры.
Кроме этого, зависимости нагрузки от деформаций для арматуры приняты линейными, в то время как они криволинейные для арматуры при ее совместной работе с бетоном.
Полезная информация о повреждениях и дефектах в железобетонных конструкциях и об их влиянии на несущую способность приведена в работе [8], однако в ней отсутствует информация о методах оценки несущей способности железобетонных элементов.
Наиболее эффективным методом определения остаточной несущей способности железобетонных балок является метод интегральных испытаний конструкций с измерением контролируемых параметров в совокупности с измерением механических и физических характеристик стали арматуры и бетона балки. Методика проведения интегральных испытаний конструкций рассмотрена в общем виде в работах [5] и др. В работе [9] приводятся методы определения несущей способности металлических и железобетонных конструкций при различных видах деформирования при ограниченных (малых) объемах информации о контролируемых параметрах с помощью методов теории возможностей [10]. При ограниченной информации описание неопределенностей этими методами оправдано [10, 11]. Однако из-за более осторожного подхода и использования целого класса или множества распределений вероятностей при описании неопределенностей в этом методе не удается выявить возможный резерв несущей способности конструкций, что снижает его экономическую эффективность в случаях, когда такой резерв имеется.
Описание исследования
Предлагается новый метод определения остаточной несущей способности железобетонной балки и балочных плит по критериям прочности арматуры и бетона на стадии эксплуатации. По этому методу балка нагружается пробной (испытательной) нагрузкой, направленной для безопасности испытаний противоположно эксплуатационной нагрузке, которую удалить по тем или иным причинам невозможно. Нагрузка (разгрузка) с помощью гидравлического домкрата с манометром прикладывается в месте наибольшего изгибающего момента от эксплуатационной нагрузки.
Наибольшее (безопасное) значение пробной нагрузки Етах определяется из равенства
Мр <МП, где М - наибольший изгибающий момент от эксплуатационной нагрузки (с учетом
" ч я.
собственного веса балки) и Мр - изгибающий момент от нагрузки Ешах в том же сечении балки, определяемые методами строительной механики.
Испытание балки нагрузкой Ц <¥т;хх осуществляется многократным (>10) воздействием на балку не менее чем при трех разных значениях ¥ для возможности проведения статистического анализа и применения известного «правила трех сигм» [12]. Предварительно на участке приложения пробной нагрузки ¥ удаляется защитный слой бетона на участке длиной 100-200 мм. Для измерения деформаций на крайних стержнях арматуры и в сжатой зоне бетона наклеивают тензорезисторы в соответствии с указаниями [7]. Для стабилизации деформаций железобетонная балка выдерживается под нагрузкой до четырех часов. Деформации по [5]
определяются по формуле £ = {^-Яо^кЦ), где Ц),^ - омическое (электрическое)
сопротивление тензорезисторов до и после приложения нагрузки и ее выдержки, к -коэффициент тензочувствительности тензорезисторов.
По результатам измерений деформаций ев и еь от каждого значения нагрузки ¥{ при числе измерений т >10 находят средние значения и средние квадратические отклонения Sj деформаций на каждом уровне пробной нагрузки . Используя «правило трех сигм» [12, 13 и др.], находят наибольшие возможные значения деформаций £•. тах = г. + с вероятностью 0,997 при каждом уровне от нагрузки . По методу наименьших квадратов [13 и др.], используя для описания случайных величин е8 и еъ по [14] нормальный закон распределения, находят с помощью компьютерных программ функции ¥в= ^ ^ и /2( + 33£ .
На рисунке 1 условно показаны эти функции в виде кривых, например, для арматуры балки. Также на рисунке 1 условно изображены графики функций плотностей распределения
вероятностей ^ и доверительные интервалы .
Зависимость между нагрузкой и деформациями в арматуре представлена на рисунке 1 криволинейной, т. к. при совместной работе арматуры с бетоном в растянутой зоне балки нарушается линейная зависимость между сг и е, наблюдаемая при испытаниях арматуры. Также
учтено возрастание среднеквадратического отклонения измеряемой величины [15].
с ростом среднего значения
35, 35
Рисунок 1. Графики функциональной зависимости нагрузки ¥ от деформаций в арматуре ¥1е ^ + и условные функции плотностей вероятностей р
Значения предельных деформаций £зпр,£Ьпр по условиям прочности стали арматуры и бетона находят по априорной информации. Так, по СП 52-101-2003 при двухлинейной диаграмме <т-£ для стальной арматуры недопустимо превышение среднего значения деформации,
равного £ =0,025. Этому значению деформации будет соответствовать верхнее значение предельной нагрузки р* как меры несущей способности балки по критерию прочности арматуры.
Значение ее находят из уравнения ^ = ^ = £
У,.Щ> _
Так как предельная деформация по СП 52-101-2003 определяется по результатам испытаний серии образцов арматуры, то ее приходится оценивать по средним значениям и доверительному интервалу, равному 65^, если использовать «правило трех сигм». Значение наименьшей предельной деформации (худший вариант) определяется как -3¿V. . В этом
случае предельная нагрузка (нижнее значение) определяется из уравнения К" — £,„„ — З^ ,
пр ^ ^ь.пр ь}1р ^
как показано на рисунке 1, т. е. по пессимистическому (наихудшему) варианту.
Известно [14], что коэффициент вариации Уа для предела текучести стали арматуры для одного и того же профиля с одного и того же завода равен Уа = 0,01 -г- 0,04 . Примем Уа = 0,025 . Среднее квадратическое отклонение прочности арматуры находим из Ба=у(7-ат или, при отсутствии информации о значениях предела текучести ат для арматуры по условию Бсг = усгКяп, где нормативное значение прочности стали арматуры, которое принимается
по СП 52-101-2003 или определяется испытаниями арматуры балки по [4]. Коэффициент вариации модуля упругости стали арматуры по [14] колеблется от 0,02 до 0,06. По СП 52-1-1-2003
принято значение = 2-105 МПа.
Примем для модуля упругости стали арматуры Ея коэффициент вариации по [16] равным
--С -3
0,04. Тогда 3Е =Е3у = 2-10 -0,04 = 8-10 МПа. При линейной зависимости между £3 и сг^
имеем ея=ая/Ея Например, при
Используя метод линеаризации [14], найдем 3Е =
2
Я
(У_
Ж
Я
а3 = ^ п = 400 МПа
1
Л2
2-10
•10
400
>л
>
имеем
= 0,5-10"
^ = уа(73 = 0,025 • 400 = 10 МПа и
3£„ = 1,5-10"
Тогда
имеем
£„^-3^=0,025 - 0,0015 = 0,0235, и эта деформация будет соответствовать нижнему значению предельной нагрузки ^, значение которой находится графически по рисунку 1 или
аналитически из функции = /9 ^ — . Верхнее значение предельной нагрузки р* находят также графически по рисунку 1 или аналитически из функции Е= ^
Косвенным образом в рассмотренном методе определения остаточной несущей способности железобетонной балки учитывается работа растягиваемого бетона балки (в отличие от СНиП 52-01-2003, введенного с 2013 г., где несущая способность балки рассчитывается без учета работы растянутого бетона). Тем самым достигается реализация существующего резерва несущей способности балок на стадии эксплуатации.
В ответственных случаях для определения параметров необходимо использовать апостериорную информацию, для чего проводят испытания арматуры балки по рекомендациям [4]. В результате мера несущей способности балки по критерию прочности арматуры характеризуется
интервалом значений [17, 18] предельной нагрузки . При более осторожном подходе
принимают Епр =ЕЦр.
2
3
3
5
Если рассматривать Fnp в интервалeF^p<Fnp<F°p
как случайную величину ¥ , то для
ее вероятностного описания можно использовать равномерное распределение [12] с функцией плотности вероятности (при обозначении Р = Х)\
/х:
i,
если x < F^p >
1/ Fв - Fн еслиFн < x < Fв
11 Fnp 1 np если 1 np ^ x ^ 1 np
если x > F^p
Если считать ^ и р^ точными крайними значениями для ¥ и средние значения принять
Т7® , то для вероятностного описания К можно использовать усеченное
как /'„,=о,5(.;;
интервальное распределение [11, 19] в виде границ нижних и верхних функций множества возможных распределений вероятностей отказа балки по прочности арматуры, представленных условно на рисунке 2 и имеющих следующие аналитические выражения:
о,
x-m
x
PX
x-a 1,
0,
b-m„
b — x
если х<т„
если mx<x<b ;
если х>Ъ
если x< a
если a<x<mx
если x m
(1)
где а и Ъ - нижняя и верхняя границы случайной величины X ; т - среднее значение X .
Значение предельной нагрузки р , принятое из интервала
F
np '
a, ь.
, или для нашей задачи
н f в
np> np .
можно характеризовать значением границ вероятностей отказа из интервала вероятностей \_,Р .
Р(х)1
0 Х = Кр
Рисунок 2. Границы множества распределений вероятностей Р_4(~2Р случайной величины X, усеченного интервального распределения а — Р^ Ъ — Р^р
0
1
По такой же методике устанавливается значение предельной нагрузки на железобетонную балку по критерию прочности бетона балки. Для этого одновременно с измерением деформаций в арматуре измеряются деформации в сжатой зоне бетона с учетом требований в [7]. По результатам испытаний балки методом наименьших квадратов [13] (по компьютерной
программе) находят функции ^ = ^Сь и > графический вид которых
аналогичен функциям на рисунке 1. Значения предельных деформаций бетона устанавливаются по априорной информации. Известно, что по СП 52-101-2003 значение предельной деформации бетона при сжатии и непродолжительном действии нагрузки составляет £ьо = 0,002. По [20]
коэффициент вариации прочности бетона при сжатии можно принять равным 0,04. Немецкие исследователи Х. Рюш и др. нашли, что прочность бетона на сжатие хорошо описывается нормальным распределением. Модуль упругости бетона при сжатии по СП 52-101-2003 зависит от
о
класса бетона. Так, например, при классе бетона ВЗО имеем Еь =34,5-10 МПа.
Методика определения класса бетона в изделии на стадии эксплуатации хорошо описана
п
в [4]. Там же описана оценка средней прочности Кт='^Кл/п и среднеквадратического
¿=1
отклонения прочности бетона по известной [13] формуле $ =
п
П- 1„ или по
2 = 1
размаху [13] в зависимости от объема информации. Перспективным является определение модуля упругости бетона методом динамического индектирования [21] применительно к эксплуатируемым железобетонным конструкциям. Использование эмпирической формулы для определения модуля упругости бетона в этом методе позволяет определить среднее значение модуля упругости и его среднее квадратическое отклонение для статистического анализа.
Таким образом, имеем среднее значение предельной деформации бетона ^ь
применяя метод линеаризации, его среднеквадратическое отклонение:
3 Р —
°пр
'142 \Еъ)
л2
г»
\ЕЪ )
К 2 КЕ
Из двух значений предельной нагрузки по критериям прочности арматуры и бетона в интервальной форме для оценки несущей способности балки на практике используется
(при осторожном подходе). Однако
наименьшая предельная нагрузка
Р ■ =Р"
пр 1П111 пр
действительная предельная нагрузка находится в интервале и остается неизвестной. В ряде обоснованных случаев, как было описано выше, возможно идти на некоторый риск [22, 23] и
принимать Рпр >Р"1р.
Для повышения точности определения предельной нагрузки (несущей способности) для железобетонной балки необходимо провести несколько независимых испытаний балки и получить
некоторое множество интервалов значений
Г рв
Гпрг, гпрг
Обрабатывая это множество
с использованием теории свидетельств Демпстера - Шейфера [24, 25], можно получить интервал предельной нагрузки и оценить его функцией доверия и функцией правдоподобия
Упр .
р1рпр ^ [26]. Если при этом будут известны в интервальном виде и эксплуатационные нагрузки,
то можно рассчитать надежность (вероятность безотказной работы) балки и любой другой конструкции, используя методику из [26]. Эти замечания показывают перспективность приведенного метода определения предельной нагрузки (несущей способности) несущих элементов любых конструкций.
Полученное значение предельной нагрузки Г не будет иметь практического значения, если нельзя выразить через нее эксплуатационную нагрузку на балку в стадии эксплуатации.
ГОСТ 31937-2011 (от 01.01.2014 г.) предписывает на стадии эксплуатации в результате обследования выявлять «реальную расчетную схему здания или сооружения и его отдельных конструкций для проверочных расчетов несущей способности».
Таким образом, возникает проблема в выявлении уточненной расчетной схемы балки на стадии эксплуатации и перехода от Г к предельной эксплуатационной нагрузке, например,
np
равномерно распределенной ^ и с учетом собственного веса самой конструкции. Рассмотрим этот вопрос на примерах.
Пример 1. Балка свободно (шарнирно) оперта и загружена равномерно распределенной нагрузкой, показанной на рисунке 3.
qnp
cL -rb п
// . l ~77. 7
г {
Рисунок 3. Расчетная схема балки
Предельную нагрузку а найдем из равенства Мг =М^ , где М^ - наибольший
11гР Чир Рпр Рпр
изгибающий момент от предельной нагрузки Р , полученной по вышеописанному методу; М„ -
Пр дпр
изгибающий момент от нагрузки дп + дсе (дсв - собственный вес балки). В качестве д}
np
принимается дП (нижнее значение, определяемое по нижнему значению ) или большее значение, но с некоторым риском.
Из расчетной схемы по рисунку 3 имеем М„ = \11„ + qr,a /, /8 и = ——. Отсюда
Чпр ^нр — пр ^
2 F" 2 FH
н _ пр и _ пр
Чпр """ Чс.в. ~ J и Чпр ~ , Чс.в.
l
2 К
Верхнее значение предельной нагрузки д°р=—у^-~Чсе ■ Несущая способность балки характеризуется интервалом нагрузки §,"„,¿7■ Действующая нагрузка на балку должна
отвечать условию д < д"пр при осторожном (пессимистическом) варианте несущей способности.
Если <7 в интервале ее значений принять как случайную величину дпр, то при среднем
np
н
пр
дпр ее можно описывать усеченным интервальным распределением по
значении дпр рисунку 2.
Пример 2. На рисунке 4 показана уточненная расчетная схема балки с учетом упругой податливости опор при наличии угловой жесткости опорных закреплений.
в
qnp
Рисунок 4. Уточненная расчетная схема балки
В такой схеме место с наибольшим изгибающим моментом неизвестно. Поэтому предварительно перед испытаниями балки нагрузкой ¥ приходится выявлять угловую жесткость опорных закреплений и только после этого решать статически неопределимую задачу и находить в балке Мтах и его место. В этом месте и прикладывается испытательная (пробная) нагрузка.
Описание метода выявления угловой жесткости опор А и В и дальнейший расчет балки по рисунку 4 можно найти в работе [9].
Заключение
1. Предложен новый метод определения несущей способности железобетонных балок и балочных плит по критериям прочности арматуры и бетона балки на стадии эксплуатации.
2. Способ в исполнении отличается безопасностью, т. к. вместо нагружения балки при интегральных испытаниях используется ее разгружение.
3. В качестве меры несущей способности использована предельная нагрузка в интервальной форме, которую выдерживает балка с определенным уровнем вероятности.
4. Предложенный метод может быть использован для других видов предельных состояний железобетонной балки, для других конструктивных форм и различных видов воздействий.
Литература
1. Рибицки Р. Повреждения и дефекты в строительных конструкциях. Пер. с нем. М.: Стройиздат, 1982. 432 c.
2. Байда Д.Н. Остаточная несущая способность железобетонных балок после их частичного разрушения: Автореф. дисс.... канд. техн. наук. Киев, 2005. 21 с.
3. Лотыш В.В. Определение несущей способности железобетонных балок при изгибе с учетом трещинообразования. Анализ развития трещин: Автореф. дисс.. канд. техн. наук. Львов, 2001. 17 с.
4. Бедов А.И., Сапрыкин В.Ф. Обследование и реконструкция железобетонных и каменных конструкций эксплуатируемых зданий и сооружений. М.: Изд-во АСВ. 1995. 192 с.
5. Землянский А.А. Обследование и испытание зданий и сооружений. М.: Изд-во АСВ, 2001. 240 с.
6. Патент на изобретение № 2460057. Способ неразрушающего определения несущей способности строительных конструкций / Уткин В.С., Редькин А.Н. М., 2012, бюл. № 24.
7. Патент на изобретение № 2275613. Способ неразрушающего контроля несущей способности железобетонных конструкций / Уткин В.С., Кошелева Ж.В. М., 2006, бюл. №12.
8. Mohebimoghaddan B., Dianat S.H. Evolution of the corrosion and strength of concrete exposed to sulfate solution // International Journal of Civil Engineering and Technology. 2012. Vol. 3. Issue 2. Pp. 158-206.
9. Уткин В.С., Уткин Л.В. Несущая способность и надежность строительных конструкций. Вологда: ВоГУ, 2000. 152 с.
10. Utkin L.V., Kozine I. On new cautious structural reliability models in the framework of imprecise probabilities // Structural Satefy. 2010. Vol. 32. Issue 6. Pp. 411-416.
11. Utkin V.S. Safety analysis of the soil beds of foundations based on bearing-capacity criterion // Soil Mechanics and Foundation Engineering. 2014. Vol. 51, No. 1 (Russian Original No.1, January-February, 2014). Pp. 9-16.
12. Шишкин Н.Ф. Основы метрологии, стандартизации и контроля качества: Учеб. пособие. М.: Изд-во стандартов, 1987. 320 с.
13. Джонсон Н., Лион Ф. Статистика и планирование эксперимента в технике и науке. Пер. с англ. М.: Мир, 1980. 510 с.
14. Шпете Г. Надежность несущих строительных конструкций. Пер. с нем. М.: Стройиздат, 1994. 288 с.
15. Ермолаев Н.Н., Михеев В.В. Надежность оснований сооружений. Л.: Стройиздат, 1976. 152 с.
16. Аугусти Г., Баратта А., Кашмати Ф. Вероятностные методы в строительном проектировании. Пер. с англ. М.: Стройиздат, 1988. 580 с.
17. Yan Wang. Imprecise probabilities based on generalized intervals for system reliability assessment // International Journal of Reliability and Safety. 2010. Vol. 4. No. 4. Pp. 319-342.
18. Ning-Cong Xiao, Hong-Zhong Huang, Yan-Feng Li, Zhonglai Wang, Xiao-Ling Zhang. Non-probabilistic reliability sensitivity analysis of the model of structural systems with s interval variables whose state of dependence is determined by constrains // Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability. 2013. Vol. 227. Issue 5. Pp. 491-498.
19. Ярыгина О.В. Методы расчета надежности железобетонных конструкций в составе зданий и сооружений при ограниченной статистической информации: Автореф. дисс.... канд.техн. наук. Санкт-Петербург, 2013. 23 c.
20. RQsch H., Sell R., Rackwity R. Statistische Analyse der Betonfestigkeit. Berlin: Dentschee Ausschup fQ Stanleton, 1969. 206 p.
21. Мацулевич О.В., Рудницкий В.А., Литвиновский Д.А. Определение модуля упругости бетона методом динамического индектирования // Неразрушающий контроль и диагностика. 2012. №4. С. 46-56.
22. Aven T. Interpretations of alternative uncertainty representations in a reliability and risk analysis context // Reliability Engineering & System Safety. 2011.Vol. 96. Issue 3. Pp. 353-360.
23. Aven T., Zio E. Some considerations on the treatment of uncertainties in risk assessment for practical decision making // Reliability Engineering & System Safety. 2011. Vol. 96. Issue 1. Pp. 64-74.
24. Dempster A.P. Upper and lower probabilities induced by a multivalved mapping // Annals of Mathematical Statistics. 1967. Vol. 38. No.2. Pp. 325-339.
25. Shafer G.A. Mathematical Theory of Evidence. Princeton: Princeton University Press, 1976. 297 p.
26. Zhang Z., Jiang C., Han X., Dean Hu, Yu S. A response surface approach for Structure reliability analysis using evidence theory // Advanced in Engineering Software. 2014. Pp. 37-45.
Владимир Сергеевич Уткин, г. Вологда, Россия Тел. раб.: +7(8172)518396; эл. почта: [email protected]
© Уткин В.С., 2015