Программа моделирования электронной аппаратуры при ударных воздействиях

ровать, а также вернуться на уровень назад.

Автоматизация планов механической обработки поверхностей деталей является дальнейшим шагом развития системы комплексной автоматизации технологической подготовки производства - необходимого инструмента повышения производительности работ машиностроительного предприятия.

Литература

1. Черепашков А.А. Компьютерные технологии. Создание, внедрение и интеграция промышленных автоматизированных систем в машиностроении. Самара: СамГТУ, 2008.

2. Методы обработки резанием круглых отверстий: справочник; [под общ. ред. Б.Н. Бирюкова]. М.: Машиностроение, 1989.

3. Харламов Г.А., Тарапанов А.С. Припуски на механическую обработку: справочник. М.: Машиностроение, 2006.

УДК 621.396.2

ПРОГРАММА МОДЕЛИРОВАНИЯ ЭЛЕКТРОННОЙ АППАРАТУРЫ ПРИ УДАРНЫХ ВОЗДЕЙСТВИЯХ

Д.И. Белов; С.В. Шумарин, к.т.н.; Е.Н. Талицкий, д.т.н.

Рассматривается анализ удароустойчивых конструкций электронной аппаратуры на ЭВМ. Предлагается инструментальная среда для их проектирования, использующая в качестве вычислителя систему конечно-элементного анализа ЛЫБУБ, которая позволяет сократить затраты на натурные испытания и сроки разработки готового изделия.

Ключевые слова: электронная аппаратура, удар, моделирование, испытания, САПР, ЛШУ8, программа, ускорения.

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

Одной из ведущих мировых программных систем, основанных на методе конечных элементов, является система конечно-элементного анализа ЛЫБУБ. Она позволяет проводить анализ прочности и устойчивости конструкций, механики разрушений материалов, учитывать нелинейность задачи (пластичность, деформации, текучесть и др.) [1]. Благодаря этим возможностям ЛШУБ отлично подходит для решения задач определения динамических характеристик ЭА, подвергаемой ударному воздействию.

Анализ удароустойчивости конструкций ЭА в системах конечно-элементного анализа включает три основных этапа: ввод данных, конечно-элементный анализ и вывод результатов. Исходными данными являются трехмерная геометрическая модель блока или шкафа ЭА (ГОСТ Р 516762000), значения механических свойств используемых материалов и параметры граничных условий, то есть схема крепления. Конечно-элементный анализ заключается в создании и последующем

расчете сеточной модели конструкции. На основании сравнения действующих ускорений и напряжений конструкции ЭА с допустимыми значениями принимается решение об ее удароустойчиво-сти. При необходимости конструкция модифицируется, процесс моделирования повторяется.

На каждом из этапов конструктор сталкивается с рядом существенных трудностей. Во-первых, при подготовке адекватной модели исследуемой ЭА важно правильно учесть демпфирующие свойства материалов и конструкций. В ЛЫБУБ демпфирование задается с помощью нескольких коэффициентов, выбор которых зависит от типа решаемой задачи и воздействующей на модель нагрузки. К тому же эти величины связаны с частотой колебаний, поэтому требуются дополнительные вычисления для определения их численных значений. Во-вторых, при задании граничных условий необходимо преобразовать единицы измерения амплитуды ударного импульса. В техническом задании на ЭА удар задается как зависимость ускорения от времени, тогда как в ЛШУБ предусмотрен ввод амплитуды только в единицах перемещения. В-третьих, выбор необходимой длины массива исходных данных, вида и размера сетки конечных элементов, определяющих точность и длительность расчетов, зависит от конкретных условий. Также требуется преобразование результатов расчета ЛШУБ: из табличных данных «перемещение - время» в совмещенные графики ускорений воздействия и отклика конструкции. Кроме

того, конструктору ЭА необходимо значительное время на освоение системы ANSYS.

Решить подобные проблемы возможно, создав специализированные инструментальные среды для математических решателей систем конечно-элементного анализа. Это позволяет, в частности, рассчитывать конструкции ЭА с помощью мощных систем конечно-элементного анализа, но с использованием стороннего графического интерфейса с радиоконструкторской терминологией. Примером такой инструментальной среды для системы ANSYS является разработанная авторами программа ShockCover (№ 2010612028; зарегистр. 17.03.2010 в Федеральной службе по интеллектуальной собственности, патентам и товарным знакам РФ), которая обеспечивает создание, анализ и модификацию ударозащищенных конструкций ЭА при значительно меньших сроках моделирования и подготовки пользователя по сравнению с использованием стандартного интерфейса.

Основная функция ShockCover заключается в подготовке полноценной модели ЭА c заданным ударным воздействием, в проведении ее расчета в ANSYS в фоновом режиме и обработке результатов. Программа включает в себя графический интерфейс, модули численного дифференцирования (Differ) и интегрирования (Integral), а также БД моделей типовых конструкций ЭА и материалов. Общая структура ShockCover представлена на рисунке 1.

Исходными данными, необходимыми для начала работы с ShockCover, являются модель исследуемой конструкции ЭА, механические параметры используемых материалов, а также математическое описание ударного воздействия. ANSYS запускается в пакетном, фоновом режиме с командным файлом (макрос расчета), формируемым ShockCover [1]. Полученные в результате расчета динамические характеристики конструкции (ускорения и напряжения) выводятся пользователю в собственном графическом интерфейсе.

Как видно из схемы, ввод исходных данных ShockCover начинается с задания характеристик ударного воздействия, представленных в ТЗ, а именно формы, амплитуды и длительности импульса ускорения (в модуле графического интерфейса «Определение ударного импульса»). Их автоматическое преобразование к едини-

цам перемещения, необходимое для ANSYS, происходит путем двойного интегрирования данных в модуле Integral. Далее создается модель конструкции ЭА в одноименном модуле. Это может выполняться двумя способами: либо с помощью импорта готовой модели, выполненной в САПР трехмерного черчения (в формате IGES или Parasolid), либо с помощью собственной БД типовых конструкций. Кроме того, в этом модуле определяются механические свойства используемых в конструкции материалов: плотность, модуль упругости, коэффициент Пуассона и коэффициент механических потерь. Значения этих величин могут быть взяты из встроенной БД материалов или введены пользователем вручную.

После подготовки исходных данных Shock-Cover автоматически формирует макрос динамического расчета и запускает ANSYS в пакетном режиме, по окончании которого программа получает файл результатов. Его обработка происходит в модуле «Вывод результатов». Преобразование полученных действующих перемещений конструкции к единицам ускорения осуществляется путем двойного дифференцирования в модуле Differ. Окончательно действующие ускорения представляются в виде графика зависимости от времени, совмещенного с графиком исходного ударного импульса, что позволяет оценить степень их возрастания и сравнить с допустимыми значениями.

ТЗ на ударное воздействие

Форма, амплитуда и длительность импульса ускорений

ShockCover

Integral

Конструкция

ЭА

Файл в формате IGES или Parasolid

Данные о материалах конструкции

Механические свойства материалов

\z

Графический интерфейс

»

БД типовых конструкций

Модуль «Определение ударного импульса»

БД материалов

Differ

Модуль «Создание модели»

Ж

Модуль «Вывод результатов»

Макрос расчета

А V

Файл результатов

ANSYS

Динамические характеристики ЭА

Microsoft Office Excel

Рис. 1. Структурная схема системы ShockCover

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

Рассмотрим алгоритм работы программы ShockCover, представленный на рисунке 2. Условно его можно разделить на этапы формирования модели ЭА, задания ударной нагрузки и получения динамических характеристик.

Формирование модели исследуемой ЭА начинается с обработки информации, введенной конструктором. Программа определяет, загружен ли

19

Задание демпфирования

Ввод количества конечных элементов

12

Сетка конечных элементов

13

Ввод мест крепления детали

16

Модальный расчет

15

Задание крепления

3

да J Импорт CAD-модели

6

БД

материалов

9

Присвоение механических свойств

_

14

Ввод типа крепления детали

20

Сборка деталей конструкции

21 Ввод характеристик ударного импульса

22

Создание массива ускорений

23

Двойное интегрирование массива ускорений

24

Создание массива перемещений

25

Выбор места удара

26

Динамический расчет

27

Вывод результатов

Конец

Рис. 2. Алгоритм формирования макроса расчета

файл, содержащий модель аппаратуры, или необходимо использовать БД для построения ее геометрии (блоки 1-3). В блоках 4-6 задаются характеристики конструкционных материалов. Далее вычисляются демпфирующие свойства конструкционных элементов ЭА. Для этого модель разбивается на отдельные детали (блок 7), для каждой из которых определяется собственная частота колебаний (СЧК) (блок 17), с помощью модального анализа ANSYS (блок 16). Выполнению каждого такого расчета предшествуют процедуры присвоения механических свойств модели (блок 9), создание сетки конечных элементов (блоки 10-12) и задание крепления детали, соответствующее ее креплению в общей конструкции ЭА (блоки 1315). После получения количественных значений СЧК вычисляется и задается коэффициент, определяющий демпфирование (блоки 18-19). После определения демпфирующих свойств всех деталей программа переходит к заключительной части формирования модели ЭА - сборке всех ее конструктивных элементов (блок 20).

Следующим этапом является приложение ударного воздействия к исследуемой конструкции. Для этого после считывания исходных данных об исходном импульсе (блок 21) создается массив ускорений (блок 22). После двойного интегрирования данного массива (блок 23) формируется таблица «перемещение - время» (блок 24), которая определяет значения ударного импульса во времени для ANSYS. В блоке 25 выбирается точка или поверхность приложения ударного воздействия, что завершает формирование макроса динамического расчета.

Заключительный этап работы программы ShockCover - получение динамических характеристик конструкции ЭА. Для этого проводится динамический расчет созданной модели в системе ANSYS (блок 26). После его окончания результаты вычислений действующих амплитуд и напряжений элементов конструкции выводятся в виде графиков (блок 27).

Для оценки достоверности результатов, получаемых в системе Shock-Cover, проведено натурное испытание блока ЭА, схема которого представлена на рисунке 3. Эксперимент проводился на ударном стенде свободного падения [2]. Погрешность поддержания испытательного режима составляла ±20 %. Измерения проводились в центре верхней печатной платы, на которой устанавли-

вался вибропреобразователь.

Основными элементами анализируемого блока ЭА являются корпус и крышка, выполненные из алюминиевого сплава Д16, а также три стекло-текстолитовые печатные платы без электрорадиоэлементов. Механические свойства материалов данных деталей приведены в таблице. Воздействию механического удара - одиночного импульса ускорения в форме полуволны синусоиды с амплитудой ускорения 500 g и длительностью 1 мс - подвергалось основание блока. Общее время исследования составляло 4 мс для возможности анализа состояния конструкции после прекращения действия ударного импульса.

Материал Модуль упругости, Н/м2 Коэффициент Пуассона Плотность, кг/м3 Коэффициент механической прочности Предел текучести, Н/м2

Стеклотекстолит 3*1010 0,22 2000 0,02 3,5*108

Д16 6,9*1010 0,33 2800 0,02 3,6*108

В программе ^НоскСоует задавались исходные данные для расчета блока ЭА (рис. 4). Для построения его геометрического образа использовалась БД типовых конструкций. На созданную модель накладывалась сетка конечных элементов размерностью 30 разбиений по каждому направлению оси, после чего был произведен расчет.

На рисунке 5 представлены результаты испытания и динамический расчет блока ЭА в программе ЗкоскСоувг и при эксперименте.

Определение ударного воздействия БлскЭА №1 Блок ЭА NF2

Создание модели Вывод результатов

Геометрические параметры

Импорт конструкции ЭА IGES PARASOLID

ф Толщина

71уро1

Геометрические размеры Ширина платы|0,053 ,м Длина платы |0.0вЭ м Высота корпуса,м Толщина корпуса|0.002 ,м

Количество печатных плат I Параметры печатных плат

Материал корпуса Модуль Юнга |б,9е10 ,Н/м*м К-нт Пуассона 0,33 Плотность [2700 КМП |0.02

Толщина платы

Ширина платы

Длина платы

N- Платы Модуль Юнга.Н/м "м^оф-нт Пуассона Плотность.кг/м*м |кмп Ь"олщина.м Уровень.м

1 ЗеЮ 0.22 2000 0.02 0,0015 0.01

2 ЗеЮ 0.22 2000 0.02 0.0015 0.03

3 ЗеЮ 0.22 2000 |0.02 00015 005

Сетка конечных элементов

Рис. 4. Исходные данные для анализа блока ЭА в программе ShockCover

Т

Эксперимент — — — 8ЬоскСоуег

Рис. 5. Результаты исследования блока ЭА при воздействии ударного импульса с амплитудой ускорения 500 g

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

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

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

Литература

1. Морозов Е.М. [и др.]. ANSYS в руках инженера: Механика разрушения. М.: ЛЕНАНД, 2010. 456 с.

2. Белов Д.И. [и др.]. Моделирование ударных процессов в электронной аппаратуре с помощью системы конечно-элементного анализа «АК8У8». Проектирование и технология ЭС. 2009. № 1. С. 32-34.

УДК 519.673

ПРОГРАММА ДЛЯ ПОДДЕРЖКИ ПРИНЯТИЯ РЕШЕНИЙ О СТРУКТУРЕ МОДЕЛЕЙ ФИЗИКО-ХИМИЧЕСКИХ СИСТЕМ

С.Б. Коныгин, к.т.н.

Описана программа для моделирования физико-химических систем методом вероятностного клеточного автомата. Приведен пример использования программы для построения модели процессов окисления металла.

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

Совершенствование существующих и создание новых технологий тесно связаны с разработкой их математических моделей. Особую роль при этом играют модели физико-химических процессов (ФХП), описывающих изменения химического состава и структуры вещества. Они широко используются при разработке процессов в химической и нефтегазовой промышленности, микроэлектронике, в материаловедении и т.д. Отличительной особенностью данных технологий является наличие широкого круга химических компонентов и параллельно протекающих ФХП как в объеме, так и на поверхности раздела фаз [1].

В связи с этим одной из ключевых задач, возникающих при построении указанных моделей, является выбор их структуры, то есть совокупности учитываемых химических компонентов, ФХП и взаимосвязей между ними.

Для принятия решения о целесообразности включения того или иного ФХП или компонента в модель ее разработчику необходимо иметь возможность быстро оценить степень влияния данных факторов на целевые показатели описываемого процесса. Удачным решением указанной про-

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

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

Согласно базовому принципу предложенного метода [1-3], микроструктуре моделируемой системы ставится в соответствие сетка ячеек ВКА, каждая из которых представляет собой вероятностный автомат. Типам частиц ставятся в соответствие состояния ячеек S, а случайным процессам