Калькулятор тепловых процессов в металлокерамике

Пиковая мощность, Вт

Длит. импульса t_pulse, с

Период повторения t_rep, с

Темп. окруж., T_amb, °C

R_jc кристалл→плата, °C/W

R_cp плата→радиатор, °C/W

R_r∞ радиатор→окр., °C/W

C_pl плата, J/°C

C_rad радиатор, J/°C

Отклик Параметры

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

Плата рассматривается как тонкая пластина, для которой собственные частоты зависят от жёсткости изгиба, геометрии, граничных условий и плотности материала. При импульсном воздействии в временной области рассчитывается переходная (временная) характеристика, а затем её спектр FFT используется для выделения резонансных частот.

Основные величины и обозначения

• a, b — размеры платы по осям X и Y (м);
• h — толщина пластины (м);
• E — модуль Юнга (Па);
• ν — коэффициент Пуассона;
• ρ — плотность материала (кг/м³);
• D — жёсткость изгиба платы (N·m);
• m — масса пластинки (кг) или поверхностная плотность ρ·h (кг/м²).

Формулы

1. Жёсткость изгиба пластины:
$$D = \frac{E\,h^{3}}{12\,\bigl(1-\nu^{2}\bigr)}$$
2. Приближённая формула собственной частоты для мода (m,n) прямоугольной пластины

Зависит от граничных условий — пример для просто поддерживаемой пластины; в более точных задачах используют численные методы / FEM:
$$\omega_{mn} \;\approx\; \pi^{2}\,\sqrt{\frac{D}{\rho\,h}}\;\Bigl(\frac{m^{2}}{a^{2}} + \frac{n^{2}}{b^{2}}\Bigr)
$$
$$
f_{mn}=\frac{\omega_{mn}}{2\pi}$$
Здесь \(m,n\) — целые числа (номера мод), а точность аналитической формулы зависит от реальных граничных условий (закрепление/запирание/свободные края).

Для точного определения собственных частот и форм колебаний рекомендуется использовать FEM.

Как работает калькулятор

1. Пользователь вводит геометрию (a, b, h) и свойства материала (E, ν, ρ) — или выбирает готовый материал (например, AlN, Al2O3 и т.п.).
2. Калькулятор вычисляет жёсткость \(D\) и приближённые частоты \(f_{mn}\) для наборов (m,n) (по умолчанию 1≤m,n≤4 — это можно расширить).
3. Если включён режим «Импульсный отклик», модель формирует временной сигнал уравнения движения (линейная система с демпфированием) и интегрирует во времени (метод Newmark/прямой интегратор).
4. Полученный временной отклик преобразуется в частотную шкалу через FFT — пики спектра соответствуют модальным частотам \(f_{mn}\).
5. Результат — таблица мод, список резонансных частот и графики: временная форма отклика и амплитудный спектр. Дополнительно — режим «анимация» формы мод (синусоидальная визуализация) при выбранной частоте.

Тип входных данных

• Геометрия: указывайте реальные размеры в метрах (например 0.05 для 50 мм).
• Толщина: при малых h (менее 0.5 мм) учтите, что реальная упаковка металлизации и слои могут влиять на эффективную жёсткость.
• Материал: для металлокерамики чаще всего используются подложки на основе AlN, Al₂O₃ или DBC (керамика + медные шины). Есть опция «металлизированный слой» — указывайте суммарную толщину металла и его свойства, если хотите учитывать слоистую модель.
• Демпфирование: в реальности платы имеют внутреннее демпфирование — задайте коэффициент затухания ξ (например 0.001…0.05) для получения реалистичных пиков в спектре.
• Возбуждение: одноимпульс (широко­полосный) — лучший вариант для выделения всех мод; узкополосный импульс (короткий гармонический пакет) — для локализации близких мод.

Пример расчёта — жёсткость изгиба

Рассмотрим алюминий-нитридную плату AlN толщиной \(h=1\) мм, \(E\approx 320\) ГПа, \(\nu\approx 0.24\).
Тогда жёсткость изгиба:
$$D = \frac{320\cdot10^{9}\cdot(0.001)^{3}}{12\,(1-0.24^{2})} \approx
$$
$$
\approx 28.3\ \text{N·m}$$
Далее по формуле для \(f_{mn}\) можно получить приближённые значения собственных частот — при необходимости калькулятор выдаст численные значения для выбранных (m,n).

Интерпретация результатов

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

Типичные ошибки и советы проектировщику

• Неправильные граничные условия. В реальной плате края не всегда строго «просто поддерживаемые» — крепления, пайка и корпуса меняют частоты.
• Игнорирование металлизации и слоистой структуры — тонкие медные дорожки и контактные площадки могут существенно изменить локальную жёсткость и массу.
• Недооценка демпфирования — без учёта затухания пики в спектре будут нереалистично острыми.
• Смена размеров: удвоение толщины повышает жёсткость ~в 8 раз (через \(h^{3}\)), что сильно сдвигает частоты вверх.

Какие режимы расчёта доступны

1. Аналитический — приближённые формулы для первых мод (быстро, грубая оценка).
2. Временной — численное интегрирование отклика на импульс + FFT (практичный режим для инженерной работы).
3. Слоистая модель — учет металлизации (двухслойная плата) — более точное прогнозирование частот.
4. FEM экспорт — экспорт параметров для расчёта в ANSYS/COMSOL при необходимости детального исследования.

Таблица свойств металлокерамических подложек

Программа предназначена для быстрого инженерного анализа: выявления потенциальных резонансов, сравнения вариантов толщин и материалов, оценки влияния демпфирования и металлизации. Для критичных проектов всегда подтверждайте расчёты численным FEM-анализом и лабораторными измерениями.