Обеспечение безопасности оборудования и обслуживающего персонала объектов нефтегазовой отрасли от воздействия ударных волн
На правах рукописи
ТРОПКИН СЕРГЕЙ НИКОЛАЕВИЧ
ОБЕСПЕЧЕНИЕ БЕЗОПАСНОСТИ ОБОРУДОВАНИЯ И
ОБСЛУЖИВАЮЩЕГО ПЕРСОНАЛА ОБЪЕКТОВ НЕФТЕГАЗОВОЙ
ОТРАСЛИ ОТ ВОЗДЕЙСТВИЯ УДАРНЫХ ВОЛН
Специальность 05.26.03 – Пожарная и промышленная безопасность
(нефтегазовая отрасль)
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата технических наук
Уфа – 2013 2
Работа выполнена в ФГБОУ ВПО «Уфимский государственный нефтяной технический университет» на кафедре «Технологические машины и оборудование».
Научный руководитель доктор технических наук Тляшева Резеда Рафисовна
Официальные оппоненты: Идрисов Роберт Хабибович доктор технических наук, профессор, Государственное унитарное предприятие «Институт проблем транспорта энергоресурсов»
РБ / отдел «Безопасность эксплуатации трубопроводных систем», заведующий отделением Самигуллин Гафур Халафович кандидат технических наук, доцент, Национальный минерально-сырьевой университет «Горный» / кафедра транспорта и хранения нефти и газа, заведующий кафедрой
Ведущая организация ООО «Центр исследований экстремальных ситуаций», г. Москва
Защита состоится «24» января 2014 г. в 14-30 на заседании диссертационного совета Д 212.289.05 при ФГБОУ ВПО «Уфимский государственный нефтяной технический университет» по адресу: 450062, Республика Башкортостан, г. Уфа, ул. Космонавтов, 1.
С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «Уфимский государственный нефтяной технический университет».
Автореферат разослан «» декабря 2013 года
Ученый секретарь диссертационного совета Ризванов Риф Гарифович
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность В настоящее время к объектам нефтегазовой отрасли проявляют все больший интерес. Это связано с тем, что количество аварий, заканчивающихся взрывами и приводящих к поражению персонала и оборудования увеличивается, а также с постоянным ужесточением требований по обеспечению безопасности промышленных объектов. Это привело к тому, что резко возросла потребность в создании новых и усовершенствовании существующих методов и способов повышения безопасности эксплуатации оборудования, а также обеспечения защищенности обслуживающего персонала.
На предприятиях нефтегазовой отрасли и нефтехимии обращается значительное количество пожаро- и взрывоопасных веществ, которые способны образовывать взрывоопасные облака, несущие потенциальную опасность.
Необходимо оценивать создаваемую ими реальную угрозу в случае возникновения аварийной ситуации, а также иметь возможность прогнозировать воздействие ударной волны на оборудование и разработать ряд мер для обеспечения взрывобезопасности.
Помимо вышеуказанного, многие объекты нефтегазовой отрасли проектировались и возводились в 60-70х годах ХХ века, когда требования к обеспечению безопасности и защищенности персонала были существенно ниже существующих. После проведения экспертизы промышленной безопасности было установлено, что для обеспечения защищенности данных объектов на современном уровне требуется возведение специальных защитных устройств в условиях существующей плотной застройки предприятий нефтегазовой отрасли.
В этой связи актуальной становится проблема повышения устойчивости оборудования, зданий и сооружений путем разработки конструкции устройства защиты от действия ударной волны.
Область исследования соответствует требованиям паспорта специальности 05.26.03 – Пожарная и промышленная безопасность (нефтегазовая отрасль): 6 Исследование и разработка средств и методов, обеспечивающих снижение пожарной и промышленной опасности технологических процессов, предупреждения пожаров и аварий, тушения пожаров;
7 Разработка технических средств защиты людей от пожаров и производственного травматизма;
9 Исследование процессов протекания аварий, условий их каскадного и катастрофического развития, разработка методов оценки различных воздействий, проявляющихся в процессе развития аварий на нефтегазовых объектах.
Целью диссертационной работы является повышение защищенности объектов высокой важности (зон размещения персонала, оборудования) от воздействия ударных волн применением нового типа защитного устройства на предприятиях нефтегазовой отрасли в условиях существующей плотной застройки пространства территории предприятия, а также обеспечением метода определения его эффективности.
Для достижения цели решаются следующие задачи:
1 Анализ потенциально опасных зон на предприятии;
определение предельных значений параметров ударных волн, которые могут появиться в случае возникновения чрезвычайных ситуаций на объектах нефтегазовой отрасли.
2 Создание метода построения численной модели поведения защитного устройства и его конструктивных элементов под действием ударной волны.
3 Создание конечно-элементной модели взаимодействия ударной волны с защитным устройством.
4 Решение полностью связанной задачи «прочность-газодинамика»
воздействия ударной волны на защитное устройство с помощью программного комплекса SIMULIA Abaqus.
5 Проведение численных экспериментов, позволяющих определить зависимость конструктивных параметров защитного устройства от различных параметров ударной волны, с целью минимизации избыточного давления, приходящегося на защищаемый объект.
Объект исследования:
Техническое средство, обеспечивающее защиту персонала и оборудования на объектах нефтегазовой отрасли от воздействия ударной волны.
Методы исследования:
В качестве основного теоретического метода был выбран метод численного моделирования связанных задач газодинамики и прочности, позволяющий учитывать методы и уравнения механики твёрдого деформируемого тела, механики сплошных сред и теории ударных волн.
Имитационное моделирование взаимодействия ударной волны с защитным устройством производилось на основе метода связанных прочностных и газодинамических расчётов Coupled Euler Lagrange.
Научная новизна 1 Создан метод построения численной модели анализа поведения защитных конструкций под действием ударной волны и определения их эффективности.
Решена полностью связанная задача воздействия ударной волны на защитное устройство, учитывающая изменения параметров обтекания устройства в зависимости от степени его деформирования или разрушения.
2 Впервые на основе численного моделирования предложен и обоснован тип защитной конструкции, представляющей из себя стальную сварную ячеистую перегородку, в каждую ячейку которой установлен элемент-гаситель, состоящий из набора особым образом расположенных листов, образующих канал типа «гребенка». Данная защитная конструкция обеспечивает безопасность персонала и оборудования объектов нефтегазовой отрасли от воздействия ударных волн.
3 Определены конструктивные параметры гасителя защитного устройства в зависимости от значения избыточного давления и профиля ударной волны для случая дефлаграции. Определены конструктивные параметры и номограмма расположения защитного устройства относительно защищаемого объекта.
На защиту выносятся 1 Конструкция защитного устройства, обеспечивающего безопасность оборудования и персонала нефтегазовой отрасли от воздействия ударных волн.
2 Метод построения численной модели анализа поведения защитных устройств под действием ударной волны и определения их эффективности.
3 Численные модели воздействия ударной волны на защитное устройство.
4 Номограммы зависимости эффективности защитного устройства от его конструктивных параметров.
Практическая ценность Разработано техническое устройство для защиты обслуживающего персонала от производственного травматизма.
Создан метод построения численной модели воздействия ударной волны на сооружения, находящиеся на территории объектов эксплуатации нефтегазовой отрасли, который можно использовать при проектировании новых защитных устройств, при реконструкции и оценке состояния защищенности уже существующих объектов.
Создана параметрическая конечно-элементная модель взаимодействия ударной волны и защитного устройства.
Достоверность проведенных исследований Достоверность исследований обеспечивается используемой в них нормативной базой;
соответствием результатов расчета изначально наложенным ограничениям;
обоснованными современными расчетными методами, а также высокой степенью соответствия результатов математического моделирования и практических экспериментов, полученных путем решения задач верификации.
Апробация результатов работы Основные положения и результаты диссертационной работы докладывались на международной научно-практической конференции «Инженерные системы» (г. Москва, 2010, 2011 г.), 64-ой научно-технической конференции молодых ученых и аспирантов (Уфа, 2013 г.), VII-ой научно практической конференции «Промышленная безопасность на взрывопожароопасных и химически опасных производственных объектах» (г.
Уфа, 2013 г.).
Публикации Основные результаты диссертационных исследований опубликованы в 8 статьях и 1 тезисе, в том числе 2 статьи в ведущих рецензируемых научных журналах из перечня ВАК РФ.
Объем и структура диссертации Диссертация состоит из введения, 4 глав, заключения, списка использованных источников из 105 наименований и приложения, содержит 162 страницы, включая 29 рисунков и 35 таблиц.
КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы диссертации, определена цель и задачи исследования, сформулированы научная новизна, практическая ценность и основные положения, выносимые на защиту, дана общая характеристика работы.
В первой главе производится анализ опасностей, возникающих на объектах эксплуатации нефтегазовой отрасли. Исследуются статистические данные о взрывах и пожарах на предприятиях нефтегазовой отрасли, причины их возникновения. Исследуются типы наиболее опасных объектов, их энергетическая насыщенность и возможные причины возникновения аварийных ситуаций.
Определён один из наиболее характерных типов взрывного превращения для углеводородных смесей – дефлаграция.
Во второй главе производится оценка параметров взрывных волн, возникающих при аварийных ситуациях на предприятиях нефтегазовой отрасли. Производится анализ существующих методов оценки параметров взрывных волн, особенностей их применения.
На основе патентного поиска и обзора литературных источников анализируются существующие типы конструкций, обеспечивающих защиту объектов от воздействия взрывных волн. Исследуется область их применения, особенности установки и эксплуатации. Рассматриваются существующие методы анализа защитных конструкций, оценки их эффективности и поведения при воздействии взрывной волны.
На основании произведенного анализа делается вывод об отсутствии разработанного защитного устройства, удовлетворяющего требованиям размещения в условиях существующей плотной застройки объектов нефтегазовой отрасли.
Предлагается конструкция нового типа защитного устройства, описываются принципы его работы и основные конструктивные особенности.
Конструкция защитного устройства представляет собой сварную стальную ячеистую перегородку, состоящую из балочных профилей, стальных листов и разрушаемых элементов – гасителей (рисунок 1). Защитное устройство рассеивает энергию ударной волны за счёт ее частичного отражения, а также за счёт поглощения энергии ударной волны в специальных ячейках – гасителях.
Гаситель представляет собой набор особым образом расположенных стальных листов, формирующих канал типа «гребенка», обеспечивающих при его прохождении снижение давления на неотраженной части фронта ударной волны. При прохождении ударной волны гасители могут упруго-пластически деформироваться и за счёт этого также обеспечивать поглощение энергии ударной волны. Таким образом, после прохождения взрывной волной защитной конструкции значительно снижается давление на её фронте. За счёт этого обеспечивается безопасность оборудования и персонала объектов эксплуатации нефтегазовой отрасли.
Для применения защитной конструкции в условиях существующих объектов необходимо определение его конструктивных параметров в зависимости от опасностей, возможных на данном объекте, а также конфигурации защищаемого объекта.
В третьей главе приводится описание метода исследования на основе численного моделирования поведения защитного устройства при действии ударной волны.
1 - рама защитного устройства, 2 - гаситель Рисунок 1 – Конструкция защитного устройства Для исследования поведения защитной конструкции необходимо решение связанной задачи газодинамики и прочности. Решение подобной задачи аналитическими и эмпирическими методами является крайне сложным, поскольку необходимо решение связанных нелинейных задач, либо проведение масштабного и крайне дорогостоящего экспериментального исследования.
Обосновывается эффективность решения задачи с помощью технологий численного моделирования.
Совместное решение задач газодинамики и прочности реализуется за счёт проведения связанных расчётов гидрогазодинамики и прочности типа Fluid Structure Interaction (FSI) – реализованного между такими программными комплексами как SIMULIA Abaqus и FlowVision, Ansys Mechanical и Ansys CFX и др., а также непосредственно в программном комплексе SIMULIA Abaqus.
Достижением последних лет в области численного моделирования является внедрение в универсальные программные комплексы численного моделирования, такие как SIMULIA Abaqus и LS-DYNA, технологии совместного решения задач прочности и гидрогазодинамики – Coupled Euler Lagrange – связанный метод Эйлера-Лагранжа. Данная технология позволяет производить моделирование задач гидрогазодинамики на Эйлеровых сетках и задачи динамической прочности на Лагранжевой сетке в рамках одного пространства решения, обеспечивая передачу параметров между сетками по типу «давление-перемещение». Данная технология позволяет решает класс задач, в которых обтекание тела напрямую зависит от его текущей деформированной конфигурации.
Приводится описание метода исследования – для решения задачи анализа воздействия взрывной волны на защитное устройство и определения его эффективности, а также для оценки степени защищенности объектов. В программном комплексе Abaqus с помощью явного решателя Abaqus/Explicit создаются и исследуются следующие модели:
1) газодинамическая модель распространения ударной волны;
2) динамическая прочностная модель поведения гасителя;
3) совмещенная газодинамическая-прочностная модель воздействия ударной волны на защитное устройство и расположенный за ним объект.
Моделирование газодинамики: решается газодинамическая задача прохождения ударной волны через расчётную область.
Воздух в расчётной области задается с помощью модели идеального газа.
Взрывная волна моделируется как скачок давления на боковой грани расчётной области по определенной амплитуде.
В результате расчёта в исследуемой области становится известен закон распределения давлений, а также скоростных и температурных характеристик при прохождении ударной волны.
Моделирование задачи динамической прочности явным методом:
исследуется воздействие взрывной волны на конструкцию с помощью приложения на поверхности скачка давления, изменяющегося по времени в соответствии с профилем ударной волны. Создается трехмерная или квази двухмерная модель.
В области решения могут быть использованы любые модели материала, применяемые в механике твердого деформируемого тела, различные критерии пластичности, вязкопластичности, хрупкого разрушения и т.п.
Выбор соответствующей модели поведения материала определяется исходя из типа материала, из которого изготавливается конструкция, и наличия необходимых данных для интересующей модели материала.
Конечно-элементная модель закрепляется в пространстве, и на выбранные поверхности прикладывается изменяющееся по времени давление.
В результате расчёта определяется динамическое напряженно деформированное состояние исследуемого объекта, его перемещения, скорости и ускорения.
Решение полностью связанной задачи взаимодействия ударной волны с конструкцией: после создания и проверки адекватности газодинамической и прочностной модели создается единая модель для анализа взаимодействия «газ конструкция». Для этого сетки моделей совмещаются в единой геометрической области, как показано на рисунке 2.
Для созданной модели учитываются все граничные условия, характерные для газодинамической и прочностной модели.
Рисунок 2 – Единая расчётная область задачи взаимодействия «газ-конструкция», включающая интересующий объем газа и исследуемую конструкцию Для передачи взаимодействия от газа к конструкции используется технология «Основного контакта». Данная технология позволяет автоматически отслеживать контактные взаимодействия в выбранной расчётной области, как между газовой средой и твердым телом, так и между твердыми телами.
Результатом решения задачи являются поля распределения давлений в газовой среде по времени на границах исследуемой конструкции, с учетом изменения формы деформируемых поверхностей, скоростные и температурные характеристики, а также напряженно-деформированное состояние конструкции от воздействия давления ударной волны.
Верификация моделей является необходимым шагом при проведении исследований с помощью численных методов. Сравнение результатов численного моделирования производится с известными экспериментальными или аналитическими решениями.
Оценка достоверности и точности результатов расчета при использовании рассмотренной выше связанной модели взаимодействия ударной волны в газовой среде с деформируемой преградой проведена на основе аналитического решения приведенного в работе «Убежища гражданской обороны:
Конструкции и расчёт» авторов В.А. Котляревский, В.И. Ганушкин, А.А.
Костин и др. Исследовалось время разрушения стены и сравнивалось с имеющимися экспериментальными данными и результатами аналитических расчётов.
Внешний вид модели приведен на рисунке 3.
Высота 1,7 м Ширина 1,32 м Длина 4 м Кладка:
Высота 26 кирпичей Длина 16 кирпичей Рисунок 3 – Общий вид модели Исследовалось поведение ограждающих конструкций из кирпича толщиной 12 см при значениях избыточного давления на фронте ударной волны равных 80 кПа и 107 кПа.
Расчёт разбивается на два этапа. На первом этапе с помощью решателя Abaqus/Standard выполняется решение статической контактной задачи по определению напряженно-деформированного состояния кирпичной кладки под действием собственного веса.
На втором этапе с помощью решателя Abaqus/Explicit явным методом решалась динамическая задача воздействия взрывной волны на кирпичный заполнитель.
Результаты численного моделирования в сравнении с экспериментальными данными и данными расчёта аналитическими методами представлены в таблице 1.
С помощью описанной выше методики была верифицирована задача динамической прочности. Для верификации полной газодинамической прочностной модели была создана модель стены с учетом газодинамики. Форма разрушения стены представлена на рисунке 4.
Таблица 1 – Сравнение результатов расчета методом конечных элементов с аналитическим решением Избыточное Численное Погрешность Натурный Аналитическое давление P, решение в численного эксперимент решение кПа решения, % ABAQUS 17,0 мс 17,2 мс 17,5 мс 80 2, 14,0 мс 14,2 мс 15,0 мс 107 7, Полученное решение для связанной задачи «газ-конструкция»
сравнивалось с ранее решенной задачей динамической прочности. Сравнение перемещений и скоростей кирпичей в центральной части стены представлено на рисунке 5, а и б. По представленным графикам хорошо видно совпадение решения, полученного в динамической и связанной постановке.
1 мс 80 мс 50 мс 100 мс Рисунок 4 – Обрушение частей кирпичной стены в интервале времени 100 мс В четвертой главе представлены результаты применения метода для исследования поведения защитного устройства при воздействии ударной волны.
б) а) Перемещение, м Скорость, м/с 2 Время, с Время, с (секунды) 1 - задача динамической прочности, 2 - задача «газ-конструкция»
Рисунок 5 – Сравнение результатов расчётов а) перемещения;
б) скорости Для проведения исследования на основе расчётной методики РД-03-409 01 и РД 09-540-03 были рассчитаны параметры взрывной волны для объектов III класса опасности, содержащих в себе более 5000 кг взрывоопасных веществ (бензин). Были рассчитаны параметры ударной волны для случая дефлаграции.
Параметры ударной волны представлены в таблице 2.
Таблица 2 – Параметры ударной волны Радиус от Эффективный Амплитуда Амплитуда Длительность Длительность эпицентра энергозапас фазы фазы фазы сжатия, фазы взрыва, м горючей сжатия, разрежения, с разряжения, с смеси, Дж кПа кПа 50 2,2E+11 198 30 0,0616 0, На рисунках 6 и 7 приведено распределение полей давлений и скоростей при прохождении ударной волны через гаситель защитного устройства. На рисунке 8 представлены графики изменения давления в различных точках защитного устройства.
В таблице 3 приводятся максимальная и минимальная величины давления в выбранных для исследования точках. Таким образом подтверждается эффективность снижения давления на фронте ударной волны при прохождении гасителя.
Рисунок 6 – Деформированная конфигурация гасителя, поле давлений Рисунок 7 – Деформированная конфигурация гасителя, поле скоростей Был выполнен анализ напряженно-деформированное состояние (НДС) элементов конструкции защитного устройства. НДС и зона пластических деформаций представлены на рисунке 9. При воздействии взрывной волны максимальные напряжения достигают 250 МПа, что соответствует переходу к пластическому деформированию конструкции для использованной в расчётах стали Ст3. Зона пластической деформации находится в точке крепления первой пластины гасителя, выделенной на рисунке 9 овалом.
Таблица 3 – Давления в выбранных точках гасителя, кПа Номер точки 1 2 3 4 5 6 Амп. фазы 185 185 132 101 60 8 сжатия Амп. фазы -29 -29 -22 -17 -11 -6 - разряжения Давление, кПа Время, с (секунды) Рисунок 8 – Давление на фронте ударной волны в различных точках гасителя На основе проведенного расчёта вычисляется коэффициент локального поглощения для гасителя kloc =Pin/Pout, где Pin – давление на фронте ударной волны до гасителя, Pout – давление на фронте ударной волны после прохождения гасителя. Коэффициент локального поглощения для выбранной конфигурации гасителя составил 18.
На основании проведенного исследования гасителя, можно сделать предположение о непроницаемости защитного устройства для ударных волн в рамках выбранного значения. Таким образом, при исследовании обтекания ударной волной всего защитного устройства можно сделать допущение о его непроницаемости, то есть можно рассматривать защитное устройство в качестве сплошной стены.
Напряжения Деформации Рисунок 9 – Распределение напряжений по Мизесу и эквивалентных пластических деформаций в гасителе На основе предложенного метода создана модель взаимодействия ударной волны с незащищенным объектом – операторной. Габаритные размеры операторной: высота - 4 м, длина - 15 м. Фронтальная стена операторной расположена в 65 метрах от места взрыва. Картина воздействия волны на операторную представлена на рисунке 10.
Распределение давлений имеет несколько характерных зон. Зона срыва потока – точка, в которой развивается максимальное давление от действия ударной волны. Центральная часть поверхности – основная часть поверхности, в которой распределение давлений равномерно.
Оценивались максимальные значения величины давления ударной волны, действующие на переднюю стену, крышу и заднюю стену здания, и зоны срыва потоков. Величины давлений представлены в таблице 4.
Таблица 4 – Давление ударной волны на различные зоны конструкции Свободно стоящая операторная Максимальное давление 178 кПа Фронтальная стена 125 кПа Крыша 62 кПа Задняя стена 36 кПа T=0,024 с T=0,075 с Фронт ударной волны Операторная Операторная T=0,051 с T=0,099 с Операторная Операторная Рисунок 10 – Взаимодействие ударной волны с операторной по времени Согласно РД-03-409-01 давления превышающие 100 кПа приводят к полному разрушению здания. Таким образом, в текущей конфигурации операторная будет разрушена. Для обеспечения безопасности операторной необходимо снизить приходящее на нее давление от ударной волны.
На основе предложенного метода производится расчёт взаимодействия ударной волны с защитным устройством и операторной.
Исследовались различные варианты высоты защитного устройства: 4 м;
5,5 м;
7 м;
8,5 м и варианты расположения защитного устройства относительно защищаемого объекта – на расстоянии 2 м, 4 м, 6 м, 8 м, 10 м, 12 м.
На основе результатов моделирования были определены значения коэффициента глобального поглощения давления ударной волны защитным устройством в зависимости от конструктивных параметров устройства и его расположения, представленные на рисунке 11. В результате исследования были определены зависимости давлений на защитное устройство и операторную от высоты защитного устройства и расстояния между защитным устройством и фронтальной стеной операторной, представленные на рисунке 12.
Высота устройства, м Рисунок 11 – Зависимость глобального коэффициента поглощения от конструктивных параметров и расположения защитного устройства Высота устройства, м Рисунок 12 - Зависимость давления фронта ударной волны на операторную в центральной точке фронтальной плоскости от высоты и расположения защитного устройства ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ 1 Достигнута цель исследования – создано и обосновано защитное устройство нового типа, используемое для повышения защищенности объектов высокой важности от воздействия ударных волн в условиях плотной застройки территории предприятия.
2 Проведен анализ статистики аварий на объектах эксплуатации нефтегазовой отрасли. Выявлены наиболее опасные производственные объекты и их энергонасыщенность.
3 Создан и верифицирован метод построения модели численного анализа поведения защитного устройства и его конструктивных элементов под действием ударной волны.
4 Созданы конечно-элементные модели прохождения ударной волны через выделенное пространство, динамическая прочностная модель гасителя защитного устройства, связанная модель «газодинамики-прочности»
воздействия ударной волны на защитное устройство.
5 Проведённый анализ параметров взрывных волн в зависимости от ситуационного плана позволил прогнозировать характеристики ударной волны для использования в предложенном методе.
6 Предложена конструкция защитного устройства, позволяющая эффективно снижать давление на фронте ударной волны. Одной из главных особенностей представленной конструкции является ее компактность и возможность использовать в случае плотной застройки.
7 Предложен метод исследования поведения защитного устройства от действия взрывной волны. Решены газодинамическая, динамическая прочностная и связанная динамическая задача «газ-конструкция», позволяющие оценить прочность предложенного защитного устройства и его эффективность.
8 Получены коэффициенты поглощения защитного устройства, позволяющие количественно оценить эффективность работы как составляющих конструкцию защитного устройства гасителей, так и общую эффективность поглощения энергии взрывной волны защитным устройством.
Основные результаты диссертационной работы опубликованы в следующих научных трудах:
в ведущих рецензируемых научных журналах:
1 Тропкин, С.Н. Оценка и обеспечение взрывоустойчивости оборудования нефтегазовой отрасли [Текст] / С.Н. Тропкин, Р.Р. Тляшева, Е.А.
Гостёнова, М.И. Кузеев // Башкирский химический журнал, 2011. – Т. 18. – № 1. – С. 118-124.
Тропкин, С.Н. Верификация математической модели разрушения части здания от воздействия ударной волны [Электронный ресурс] / С.Н. Тропкин, Р.Р. Тляшева, М.И. Баязитов, З.Р. Рафикова, И.Р. Кузеев // Электронный научный журнал «Нефтегазовое дело». - 2013. – №1. – С. 476-486. Режим доступа: http://www.ogbus.ru/authors/TropkinSN/TropkinSN_1.pdf в других изданиях Тропкин, С.Н. Исследование динамического поведения аппаратов колонного типа при сейсмических воздействиях [Текст] / С.Н. Тропкин, Р.Р.
Тляшева, И.Р. Кузеев // Мировое сообщество: проблемы и пути решения. 2010. – Т. 88. – С. 52-56.
Гостёнова, Е.А. Исследование динамического поведения ректификационной колонны с трубопроводной обвязкой при воздействии воздушной взрывной волны с помощью программного комплекса SIMULIA ABAQUS [Текст] / Е.А. Гостёнова, И.Р. Кузеев, Р.Р. Тляшева, С.Н. Тропкин // Инженерные системы–2010: Материалы международной научно-практической конференции. – М.: Изд-во РУДН, 2010. – С. 88 - 90.
5 Тропкин, С.Н. Разработка защитного устройства операторной станции при воздействии воздушной взрывной волны с помощью программного комплекса ABAQUS [Текст] / С.Н. Тропкин, Р.Р. Тляшева, М.И. Баязитов, О. Л. Смольников // Инженерные системы–2011: Материалы 4-ой международной научно-практической конференции. – М.: Изд-во РУДН, 2011.
– С. 107- Тропкин, С.Н. Обеспечение защищенности операторной нефтеперерабатывающего предприятия от воздействия ударной волны [Текст] / С.Н. Тропкин, Р.Р. Тляшева, М.И. Баязитов, О.Л. Смольников // Мировое сообщество: проблемы и пути решения. - 2011. – Т. 30. – С. 71-73.
7 Тляшева Р.Р. Обеспечение защищенности обслуживающего персонала предприятий нефтегазопереработки [Текст] / Р.Р.Тляшева, И.Р. Кузеев, С.Н.
Тропкин, М.И. Баязитов, И.Б. Сиражетдинов // Развитие инновационной инфраструктуры университета: Материалы III международного научного семинара - Уфа: Изд-во «Реактив», 2012. – С. 3 - 5.
8 Кузеев И.Р. Метод оценки взрывоустойчивости оборудования нефтегазовой отрасли [Текст] / И. Р. Кузеев, Р.Р. Тляшева, М.И. Баязитов, С.Н.
Тропкин // Химическая техника. – 2013. - №6. – С. 30- 9 Кузеев, И.Р. Оценка взрывоустойчивости оборудования нефтегазовой отрасли [Текст] / Р.Р.Тляшева, М.И.Баязитов, С.Н.Тропкин // Обеспечение промышленной безопасности на взрывопожароопасных и химически опасных производственных объектах: Сборник материалов VII научно-практической конференции. –Уфа: Изд-во УГНТУ, 2013. – С. 92-99.