АКТУАТОР

Искусственная мышца на основе композита СВМПЭ–крейзинг с атакамитом и аусетической структурой: концепция резонансного управления

Аннотация

В работе предложена концепция искусственной мышцы, сочетающей три физических эффекта: (1) обратимый межкристаллитный крейзинг сверхвысокомолекулярного полиэтилена (СВМПЭ) с высокой обратимостью деформации (50–85%); (2) магнитокалорический эффект атакамита (Cu₂Cl(OH)₃), обеспечивающий управляемое охлаждение композита в импульсном магнитном поле; (3) аусетическую re-entrant сотовую структуру с отрицательным коэффициентом Пуассона для усиления механического отклика. Управление предлагается осуществлять резонансным электромагнитным полем с частотой 57,142857 Гц, что, согласно гипотезе α-Триады, соответствует фундаментальному резонансу тетраэдральной решётки. Ожидаемые характеристики прототипа: относительное сокращение ≥30%, усилие ≥50 кг/см², ресурс ≥100 000 циклов, КПД ≥85%.

Ключевые слова: СВМПЭ, крейзинг, атакамит, магнитокалорический эффект, аусетическая структура, отрицательный коэффициент Пуассона, искусственная мышца, резонансное управление, α-Триада, 57 Гц

1. Введение

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

В настоящей работе мы объединяем три эффекта: (1) обратимый межкристаллитный крейзинг СВМПЭ; (2) магнитокалорический эффект атакамита; (3) аусетическую геометрию. Управление предлагается осуществлять резонансным полем с частотой 57,142857 Гц.

2. Физические основы

2.1. Межкристаллитный крейзинг СВМПЭ

Межкристаллитный крейзинг (МКК) СВМПЭ — это процесс деформации плёнок в присутствии физически активных жидких сред (ФАЖС), приводящий к образованию фибриллярно-пористой структуры с наноразмерными порами. При последующей релаксации напряжений происходит низкотемпературная спонтанная обратная деформация.

Ключевые параметры процесса:

• Обратимость деформации: 50–85%
• Восстановление пористой структуры при вторичном нагружении: до ~20 об.%
• Размер пор: нанометровый диапазон (до 10 нм)
• Эффект открытия/закрытия пор при циклическом нагружении

Это делает СВМПЭ-крейзинг идеальной матрицей для управляемого актуатора: при охлаждении матрица усаживается, при нагреве — возвращается в исходное состояние.

2.2. Магнитокалорический эффект атакамита

Атакамит (Cu₂Cl(OH)₃) — минерал с уникальной магнитной структурой: ионы меди образуют длинные цепочки из соединённых треугольников («пилообразные цепочки»). Из-за треугольной геометрии спины ионов не могут выстроиться полностью антипараллельно — возникает магнитное фрустрированное состояние.

При воздействии магнитного поля этот фрустрированный порядок разрушается, что вызывает значительное охлаждение материала — температура падает почти вдвое. Исследователи из TU Braunschweig и HZDR подтвердили этот эффект в импульсных магнитных полях.

Для нашего применения это означает: подавая импульсное магнитное поле на композит СВМПЭ + атакамит, мы можем управлять температурой материала, а следовательно — его механической деформацией.

2.3. Аусетическая структура (re-entrant honeycomb)

Ауксетики — материалы с отрицательным коэффициентом Пуассона: при сжатии они становятся тоньше в поперечном направлении (в отличие от обычных материалов, которые расширяются). Для re-entrant сотовых структур (вывернутые пчелиные соты) отрицательный коэффициент Пуассона может достигать –5.

В нашей конструкции это означает:

• При осевом сжатии ячейки втягиваются внутрь, а не распирают
• Плотность материала в рабочей зоне увеличивается
• Передача усилия становится более эффективной

3. Гипотеза резонансного управления на частоте 57,142857 Гц

Согласно α-Триаде (резонанс–связь–материя), фундаментальная частота

f0=400/7≈57,142857 Гц связана с резонансной структурой тетраэдральной решётки. Мы предполагаем, что воздействие переменного магнитного поля на этой частоте на композит СВМПЭ + атакамит позволит:

1. Эффективно возбуждать магнитные моды атакамита (резонансное усиление магнитокалорического эффекта)
2. Синхронизировать механический отклик СВМПЭ-матрицы с частотой поля
3. Достичь максимальной энергоэффективности преобразования

Данная гипотеза требует экспериментальной проверки и является предметом настоящего исследования.

4. Конструкция и материалы

4.1. Матрица

Пористый СВМПЭ-крейзинг (плёнка или волокно) с открытыми порами 10–50 нм, полученный деформацией в присутствии ФАЖС с последующей стабилизацией.

4.2. Активный наполнитель

Синтезированный атакамит (Cu₂Cl(OH)₃) с подтверждённой кристаллической структурой (РФА), диспергированный в порах СВМПЭ.

4.3. Геометрия

Re-entrant сотовая структура (вывернутые пчелиные соты) с отрицательным коэффициентом Пуассона, вырезанная лазером из композитного листа.

4.4. Управление

Двухканальный генератор с частотой 57,142857 Гц, регулируемой скважностью импульсов (ШИМ) и возможностью сдвига фазы 0°/180° для двух актуаторов (основная мышца и антимышца).

4.5. Конструкция «мышца + антимышца»

Система состоит из двух идентичных актуаторов, работающих в противофазе:

• Основная мышца — сокращается при подаче поля (охлаждение → усадка СВМПЭ)
• Антимышца — сокращается в противоположном направлении при переключении фазы, обеспечивая активный возврат за миллисекунды без ожидания пассивного нагрева

5. Ожидаемые характеристики

Предполагается, что разработанный образец будет демонстрировать следующие параметры: относительное сокращение активного участка составит не менее 30% от исходной длины; удельное усилие на единицу поперечного сечения достигнет не менее 50 кг/см²; время сокращения и время активного возврата (при переключении на антимышцу) не превысит 20 мс в каждом случае; циклический ресурс выдерживания составит не менее 100 000 циклов «сокращение – возврат»; эффективность преобразования электрической энергии в механическую работу будет не ниже 85%; работоспособность сохраняется в диапазоне температур от –50 до +60 °C.

6. Заключение

Предложенная концепция объединяет три экспериментально подтверждённых эффекта:

1. Обратимый крейзинг СВМПЭ (обратимость 50–85%, нанопоры до 10 нм)
2. Магнитокалорический эффект атакамита (охлаждение почти вдвое в импульсном поле)
3. Аусетическую re-entrant структуру (отрицательный коэффициент Пуассона до –5)

Управление на частоте 57,142857 Гц, согласно гипотезе α-Триады, может обеспечить резонансное усиление эффекта и высокую энергоэффективность. Система из двух противофазных актуаторов (мышца + антимышца) позволяет достичь активного возврата за миллисекунды, что критически важно для робототехнических применений.

Симуляция:

python

import numpy as np

import matplotlib.pyplot as plt

# ============================================================

# ПАРАМЕТРЫ СИСТЕМЫ

# ============================================================

# Частота и период

f0 = 400/7 # 57.142857 Гц

T_period = 1/f0 # 17.5 мс

# Скважность импульсов (доля времени, когда поле включено)

D = 0.3 # рабочий режим

# Длительность импульса и паузы

t_imp = D * T_period

t_pause = (1 - D) * T_period

# Параметры тепловой модели

C_eff = 2.0 # эффективная теплоёмкость композита, Дж/(К·см³)

lambda_heat = 0.5 # коэффициент теплообмена с окружающей средой, Вт/(К·см³)

T_env = 20.0 # температура окружающей среды, °C

T_initial = 20.0 # начальная температура

# Магнитокалорический эффект атакамита (охлаждение при включении поля)

DeltaT_cool = -60.0 # °C, пиковое охлаждение при импульсе (зависит от амплитуды поля)

# Параметры механической деформации

# Коэффициент термической усадки СВМПЭ-крейзинга (отрицательное расширение при охлаждении)

alpha_T = -0.008 # 1/°C, ~0.8% усадки на градус (ориентировочно)

# Максимальная обратимая деформация крейзинга

eps_max = 0.35 # 35% при охлаждении от 20 до -40°C

# Усиление от аусетической структуры (отрицательный Пуассон)

ausetic_gain = 1.5 # коэффициент усиления деформации

# Параметры нагрузки

load_kg_per_cm2 = 50.0 # целевое усилие

# Время симуляции

t_total = 0.5 # секунд (500 мс)

dt = 0.0001 # шаг по времени, с (0.1 мс)

N_steps = int(t_total / dt)

# ============================================================

# ИНИЦИАЛИЗАЦИЯ

# ============================================================

time = np.linspace(0, t_total, N_steps)

# Для двух каналов: 0 - основная мышца, 1 - антимышца

T = np.zeros((2, N_steps))

eps = np.zeros((2, N_steps))

force = np.zeros((2, N_steps))

# Начальные условия

T[0, 0] = T_initial

T[1, 0] = T_initial

eps[0, 0] = 0.0

eps[1, 0] = 0.0

# Фазовый сдвиг для антимышцы (180° = половина периода)

phase_shift = 0.5 * T_period

# Переменные для хранения истории температуры для гистерезиса

T_prev = T_initial

# ============================================================

# ЦИКЛ СИМУЛЯЦИИ

# ============================================================

for i in range(1, N_steps):

t = time[i]

# Определяем фазу для основной мышцы (канал 0)

# Поле включено, если t mod T_period < t_imp

mod_time0 = t % T_period

field_on0 = (mod_time0 < t_imp)

# Для антимышцы (канал 1) - сдвиг фазы

mod_time1 = (t + phase_shift) % T_period

field_on1 = (mod_time1 < t_imp)

# Тепловой баланс для каждого канала

for ch in [0, 1]:

field_on = field_on0 if ch == 0 else field_on1

# Мощность охлаждения/нагрева

if field_on:

# Поле включено: атакамит охлаждается (отрицательный источник тепла)

P_field = -C_eff * abs(DeltaT_cool) / t_imp # средняя мощность охлаждения за импульс

else:

P_field = 0.0

# Теплообмен с окружающей средой

P_env = lambda_heat * (T_env - T[ch, i-1])

# Изменение температуры

dT_dt = (P_field + P_env) / C_eff

T[ch, i] = T[ch, i-1] + dT_dt * dt

# Ограничение температуры (не ниже -136°C, не выше +60°C)

T[ch, i] = np.clip(T[ch, i], -136.0, 60.0)

# Механическая деформация (зависит от температуры)

# Используем модель: eps = eps_max * (T_ref - T) / (T_ref - T_min)

# где T_ref = 20°C (начальная), T_min = -40°C (макс усадка)

T_ref = 20.0

T_min = -40.0

if T[ch, i] <= T_ref:

eps_thermal = eps_max * (T_ref - T[ch, i]) / (T_ref - T_min)

else:

eps_thermal = 0.0 # при нагреве выше начальной температуры - расширение не учитываем

# Применяем аусетическое усиление (для сжатия)

eps_ausetic = eps_thermal * ausetic_gain

# Ограничение по максимальной обратимой деформации крейзинга (85%)

eps_ausetic = np.clip(eps_ausetic, 0.0, 0.85)

# Деформация с учетом знака (отрицательная - сокращение)

eps[ch, i] = -eps_ausetic # отрицательная деформация = сокращение

# Усилие (пропорционально модулю упругости и деформации)

# Для простоты считаем, что усилие линейно зависит от деформации

# Коэффициент подогнан под целевое усилие 50 кг при деформации 30%

E_mod = 500.0 # кг/см² (условный модуль упругости)

force[ch, i] = E_mod * abs(eps[ch, i]) # кг/см²

# ============================================================

# ВИЗУАЛИЗАЦИЯ РЕЗУЛЬТАТОВ

# ============================================================

fig, axes = plt.subplots(3, 1, figsize=(12, 10))

# 1. Температура

ax1 = axes[0]

ax1.plot(time, T[0, :], 'b-', label='Основная мышца', linewidth=1.5)

ax1.plot(time, T[1, :], 'r-', label='Антимышца', linewidth=1.5)

ax1.axhline(y=-40, color='k', linestyle='--', alpha=0.5, label='Рабочий нижний предел')

ax1.axhline(y=20, color='k', linestyle='--', alpha=0.5, label='T среды')

ax1.set_ylabel('Температура, °C')

ax1.set_xlabel('Время, с')

ax1.set_title(f'Температура актуаторов (D={D})')

ax1.legend()

ax1.grid(True, alpha=0.3)

# 2. Деформация (относительное сокращение)

ax2 = axes[1]

ax2.plot(time, eps[0, :], 'b-', label='Основная мышца', linewidth=1.5)

ax2.plot(time, eps[1, :], 'r-', label='Антимышца', linewidth=1.5)

ax2.axhline(y=-0.3, color='k', linestyle='--', alpha=0.5, label='Целевое сокращение 30%')

ax2.set_ylabel('Относительная деформация (отрицательная = сокращение)')

ax2.set_xlabel('Время, с')

ax2.set_title('Деформация актуаторов')

ax2.legend()

ax2.grid(True, alpha=0.3)

# 3. Усилие

ax3 = axes[2]

ax3.plot(time, force[0, :], 'b-', label='Основная мышца', linewidth=1.5)

ax3.plot(time, force[1, :], 'r-', label='Антимышца', linewidth=1.5)

ax3.axhline(y=50, color='k', linestyle='--', alpha=0.5, label='Целевое усилие 50 кг/см²')

ax3.set_ylabel('Усилие, кг/см²')

ax3.set_xlabel('Время, с')

ax3.set_title('Развиваемое усилие')

ax3.legend()

ax3.grid(True, alpha=0.3)

plt.tight_layout()

plt.show()

# ============================================================

# ВЫВОД СТАТИСТИКИ

# ============================================================

# Определим установившийся режим (последние 0.1 с)

steady_start = int(0.9 * N_steps)

eps0_steady = eps[0, steady_start:]

eps1_steady = eps[1, steady_start:]

T0_steady = T[0, steady_start:]

T1_steady = T[1, steady_start:]

print("="*50)

print("РЕЗУЛЬТАТЫ СИМУЛЯЦИИ")

print("="*50)

print(f"Скважность D = {D}")

print(f"Средняя температура основной мышцы: {np.mean(T0_steady):.1f} °C")

print(f"Средняя температура антимышцы: {np.mean(T1_steady):.1f} °C")

print(f"Среднее сокращение основной мышцы: {-np.mean(eps0_steady)*100:.1f}%")

print(f"Среднее сокращение антимышцы: {-np.mean(eps1_steady)*100:.1f}%")

print(f"Максимальное усилие основной мышцы: {np.max(force[0, :]):.1f} кг/см²")

print(f"Максимальное усилие антимышцы: {np.max(force[1, :]):.1f} кг/см²")

print(f"Время сокращения (10-90%): ~{np.argmax(eps[0, :] < -0.27)*dt*1000:.1f} мс")

print("="*50)

Что можно менять

• D – скважность (0.05–0.8). Чем больше D, тем холоднее.
• DeltaT_cool – сила магнитокалорического эффекта (зависит от амплитуды поля).
• ausetic_gain – усиление от аусетической структуры.
• t_total – длительность симуляции.

Дальнейшие работы направлены на экспериментальную верификацию предложенной концепции и оптимизацию параметров композита.

Авторы: Марков М.В., Приходько С.М.

Проект: Resonance Group