CC BY
56
Измерение. Мониторинг. Управление. Контроль
УДК 531.2.084.2:004.94
М. А. Фролов, А. В. Самошин
МЕТОДИКА РАСЧЕТА ЧУВСТВИТЕЛЬНОГО ЭЛЕМЕНТА ТИПА «МЕМБРАНА-ПО ДУШКА-СТАКАН» МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ
M. A. Frolov, A. V. Samoshin
DESIGN PROCEDURE OF THE SENSITIVE ELEMENT OF TYPE «MEMBRANE-PILLOW-GLASS» THE FINITE ELEMENT METHOD
Аннотация. Рассмотрены методики расчета датчика давления методом конечных элементов в осесимметричной постановке и аналитическим способом, их преимущества и недостатки. Произведено сравнение результатов расчета с экспериментальными данными.
Abstract. Design procedures of the gauge of pressure by a method of final elements in ax-isymmetric statements and are considered analytical way. Their advantages and lacks are considered. Comparison of results of calculation with experimental data is made.
Ключевые слова: датчик давления, метод конечных элементов, выходной сигнал, чувствительный элемент, мембрана, подушка, эффективная площадь.
Key words: the pressure sensor, finite element method, the target signal, sensitive element, membrane, pillow, the effective area.
Решение задач механики на ранних этапах разработки и проектирования датчикопреобразующей аппаратуры (ДНА) с целью получения оптимальных геометрических размеров и стабильных метрологических характеристик является важным этапом современного проектирования.
Все измерения начинаются с восприятия измеряемых величин и формирования измерительного сигнала, который далее подвергается необходимым преобразованиям. Функцию восприятия входной величины выполняет чувствительный элемент.
Существует множество вариантов конструктивных исполнений ЧЭ: мембрана, мембрана-балка, «мембрана-колпачок», «кольцо», «стакан» и т.д. В данной статье рассмотрены методики, позволяющие решать аналитическим методом и методом конечных элементов (МКЭ) с помощью программного модуля (НМ) Simulation программы SolidWorks задачу определения требуемого выходного сигнала.
В качестве примера для апробации методики расчета выходного сигнала был взят ЧЭ датчика давления тензорезистивного типа ЛХ-412, представляющий собой узел, включающий в себя три детали: мембрана-подушка-стакан, - с наклеенными на внешней цилиндрической поверхности тензорезисторами.
Нростейшая схема ЧЭ в виде стакана изображена на рис. 1.
2013,№4(6)
Рис. 1. Схема ЧЭ в форме стакана
Как видно из рис. 1, основными геометрическими размерами ЧЭ, непосредственно влияющими на значение выходного сигнала, являются толщина мембраны (И), внутренний диаметр стакана (!), внешний диаметр стакана (О). Стрелкой указано направление воздействия нагрузки, в нашем случае - это давление измеряемой среды.
Аналитический метод решения данной задачи подразумевает следующий порядок действий:
1. На первом этапе производится расчет воспринимающей рабочее давление мембраны с целью определения ее тягового усилия. Тяговое усилие характеризует способность мембраны передавать силу объекту, находящемуся на пути ее перемещения. Тяговое усилие выражают как [1]
Р = ЯэфР ,
где Язф - эффективная площадь мембраны; Р - давление измеряемой среды.
(1)
Эффективная площадь для подобного типа мембран [1]
$эф = Ц • (D2 + D • d + d 2),
(2)
где О - диаметр внешней части мембраны (возле жесткой заделки); ! - диаметр жесткого центра (диаметр подушки).
2. На втором этапе производится расчет относительной радиальной и тангенциальной деформаций [2]:
Єг =■
Є =
F ,
SCE;
^F
SIE ,
(3)
где el - относительная радиальная деформация; F - сила (получена на первом этапе расчета); <SCT - площадь поперечного сечения стакана; E - модуль упругости материала стакана; ed - относительная тангенциальная деформация; ц - коэффициент Пуассона.
3. На третьем этапе в зависимости от ранее определенной схемы преобразователя производится преобразование величин деформации в выходной сигнал датчика (мВ, мА и т.д.). В случае несоответствия полученного значения выходного сигнала заданным требованиям технического условия производится изменение геометрических размеров узла «мембрана-подушка-стакан», и вновь повторяется расчет, начиная с первого этапа. Такой аналитический метод решения принято называть прямым.
На практике наиболее удобен обратный путь решения, т.е., изначально задавшись требуемым значением выходного сигнала (относительной деформацией), определить геометрические размеры стакана и мембраны.
Приведенные выше методики расчета имеют ряд недостатков и связаны со спецификой настройки выбранного в качестве примера датчика ЛХ-412, где после сборки узлов ЧЭ производится обдавливание мембраны избыточным давлением, превышающим номинальное значение (обычно 1,25^1,5-^Рном), вследствие чего мембрана подвергается пластической деформации. Обдавливание необходимо для выборки зазоров, образовавшихся в процессе сборки и из-
43
готовления узлов (в зависимости от конструктивного исполнения датчика компенсация всех зазоров может происходить с использованием различных компенсаторов, которые могут устанавливаться между подушкой и стаканом либо стаканом и основанием). После процесса об-давливания мембраны происходит изменение ее эффективной площади, и как следствие сила, определенная на первом этапе расчета, не будет соответствовать реальной. Затруднен выбор оптимальной толщины мембраны (h), поскольку большинство аналитических формул сводят расчет мембраны к определению ее максимального прогиба или относительных деформаций, исходя из ее основных геометрических размеров и физико-механических свойств (толщины, диаметра, модуля упругости, коэффициента Пуассона).
Решая аналогичную задачу методом конечных элементов, мы исключаем перечисленные выше недостатки, поскольку имеется возможность внесения любого изменения в геометрию модели (тем самым копируя геометрию мембраны после воздействия на нее избыточного давления).
Рассматриваемая конструкция является осесимметричной, следовательно, наиболее рационально использовать 2Б-упрощение расчетной модели. Модуль расчета Simulation программы SolidWorks позволяет решать данную задачу в 2Б-интерпретации. Суть 2Б-упрощения заключается в следующем: программа относительно выбранной осевой линии и плоскости строит оболочечную модель по контуру 3Б-модели (рис. 2). Дальнейший расчет ведется только в сечении, построенном программой, исходя из условия равенства в любом другом сечении.
Данное упрощение многократно сокращает время решения задачи за счет уменьшения количества степеней свободы конечных элементов, что позволяет использовать оптимальную по размеру сетку конечных элементов для улучшения сходимости результатов.
Решение данной задачи МКЭ в современных САПР представляет собой работу с пре- и постпроцессорами. В препроцессоре происходит задание граничных условий, имитирующих работу датчика. Для нашей модели, исходя из конструктивных особенностей датчика, достаточно наложить всего четыре граничных условия (ГУ), а именно:
1) задание связей между деталями расчетной модели. По умолчанию в ПМ Simulation используется глобальный контакт - «Связанные». Рекомендуется добавить дополнительный контакт компонентов - «Нет проникновения» между воспринимающей мембраной и подушкой;
2) задание воздействующей силы. На внешнюю сторону мембраны (со стороны воздействия измеряемой среды) задается граничное условие - «Давление». В нашем случае величина данного ГУ равна верхнему пределу измерения датчика;
3) задание ограничений. Согласно схеме установки датчика на объекте выбрано ГУ -«Жесткая заделка». Данное ГУ накладывается на грань расчетной модели, которая в реальной конструкции имеет резьбу и ввинчивается при монтаже на объекте;
4) для улучшения сходимости и повышения точности расчета создается управление для сетки конечных элементов, локально меняющее размер. В нашем случае необходимо уменьшить величину конечного элемента на гранях мембраны и стакана.
Работа с препроцессором закончена. Решение задачи происходит в автоматическом режиме. После решения задачи пользователь начинает работу с постпроцессором программы, в котором только необходимо выбрать интересующие величины (в нашем случае относительные радиальные и тангенциальные деформации).
Рис. 2. Пример 2Б-упрощения
2013,№4(6)
45
Преобразование величины относительных деформаций в выходной сигнал идентично аналитическому решению.
Для апробации предложенного подхода к расчету ЧЭ данной конструкции воспользуемся действующей конструкцией датчика ЛХ 412.
Проведем расчет аналитическим методом и МКЭ апробируемого датчика с пределом измерений 15 МПа. Решим задачу аналитическим способом, для этого подставим соответствующие значения О и ! в формулу (2), получим 5эф = 1,1-Ю-4 м2. При данном значении эффективной площади тяговое усилие Р по формуле (1) составляет 1617 Н. По формулам (3) определим значения относительных деформаций: = 12,3-Ю-4; = 3,7-Ю-4. Выходной сигнал
для данных значений относительных деформаций составляет 19 мВ [2]. По результатам моделирования при данных геометрических размерах ЧЭ значение выходного сигнала датчика составляет 21 мВ. По результатам экспериментальных замеров выходной сигнал датчика составляет ~ 22,4 мВ. Погрешность моделирования в данном случае составляет 6,3 %, погрешность аналитического расчета - 15,2 %.
Предложенный для решения осесимметричных задач подход доказал свою эффективность, показав наименьшее расхождение с экспериментальными данными, а также высокую производительность. Высокая точность расчета позволила сократить время настройки ДПА за счет минимального разброса выходных характеристик, что в конечном счете привело к сокращению времени изготовления. Дополнительным преимуществом расчета МКЭ является тот факт, что расчетные модели являются основой для создания 3Б-моделей датчиков, которые на последующих этапах проектирования путем создания ассоциативных видов преобразуются в чертежи, которые являются основой конструкторской документации (КД), сокращая тем самым время ее разработки.
Использование современных САПР на стадии проектирования ответственных деталей и узлов современных датчиков давления позволяет избежать ошибок проектирования, которые впоследствии могут являться причиной отказа ракетно-космической техники.
Список литературы
1. Водяник, В. И. Эластичные мембраны / В. И. Водяник. - М. : Машиностроение, 1974. -131 с.
2. Разработка теории и методов расчета упругих элементов сложной конструкции /
Е. П. Осадчий, А. И. Тихонов [и др.]. - Пенза, 1978. - 871 с.
Фролов Михаил Алексеевич
начальник лаборатории, Научно-исследовательский институт физических измерений E-mail: [email protected]
Frolov Mihail Alekseevich
chief of laboratory, Scientific research Institute of physical measurements
Самошин Алексей Владимирович
ведущий инженер, Научно-исследовательский институт физических измерений E-mail: [email protected]
Samoshin Alexey Vladimirovich
leading engineer,
Scientific research Institute of physical measurements
УДК 531.2.084.2:004.94 Фролов, М. А.
Методика расчета чувствительного элемента типа «мембрана-подушка-стакан» методом конечных элементов / М. А. Фролов, А. В. Самошин // Измерение. Мониторинг. Управление. Контроль. - 2013. - № 4 (6). - С. 42-45.
