|
|
|
|
|
Название проекта: |
|
Мне нужно смоделировать движение стального стержня через магнит
|
|
Кто разместил: |
|
Внешний проект с weblancer.net
|
|
Открыт: |
|
29-Nov-2025 08:03 GMT |
|
Описание: |
|
<p># Сохраните как simulate_<a href="http://magnetization.py">magnetization.py</a> и запустите: python simulate_<a href="http://magnetization.py">magnetization.py</a></p><p>import numpy as np</p><p>import pandas as pd</p><p>import matplotlib.pyplot as plt</p><p>mu0 = 4*np.pi*1e-7</p><p>params = {'Ms':1.6e6, 'a':1.0e5, 'k':5.0e4, 'c':0.2, 'alpha':1e-3}</p><p>sigma = 1.0e6</p><p>rod_thickness = 0.01</p><p>magnet_peak_B = 0.5</p><p>magnet_center = 0.0</p><p>magnet_sigma = 0.02</p><p>def B_external(x):</p><p> return magnet_peak_B <i> np.exp(-0.5 </i> ((x - magnet_center)/magnet_sigma)**2)</p><p>def M_anh(H_eff, Ms, a):</p><p> return Ms*np.tanh(H_eff/a)</p><p>tau = mu0 <i> sigma </i> (rod_thickness/2)**2</p><p>print("tau (s) =", tau)</p><p>def integrate_JA_time(H_ts, dt, params, tau):</p><p> Ms = params; a = params; k = params; c = params; alpha = params</p><p> M = 0.0</p><p> M_hist = np.zeros_like(H_ts)</p><p> for i, H_ext in enumerate(H_ts):</p><p> dH_dt = 0.0 if i==0 else (H_ts-H_ts)/dt</p><p> H_eddy = -tau * dH_dt</p><p> H_eff = H_ext + H_eddy</p><p> Man = M_anh(H_eff + alpha*M, Ms, a)</p><p> delta = np.sign(dH_dt) if dH_dt != 0 else 1.0</p><p> denom = (k <i> delta - alpha </i> (Man - M))</p><p> dM_dH = 0.0 if abs(denom) < 1e-12 else (Man - M) / denom</p><p> dM_dt = dM_dH * dH_dt</p><p> M = M + dM_dt * dt</p><p> M = c*Man + (1-c)*M</p><p> M_hist = M</p><p> return M_hist</p><p>velocities = np.array()</p><p>results = </p><p>for v in velocities:</p><p> x_start, x_end = -0.08, 0.08</p><p> total_time = abs(x_end-x_start)/v</p><p> nsteps = max(2000, int(total_time/1e-5))</p><p> t = np.linspace(0, total_time, nsteps)</p><p> dt = t-t</p><p> x_t = x_start + v*t</p><p> B_t = B_external(x_t)</p><p> H_t = B_t / mu0</p><p> M_t = integrate_JA_time(H_t, dt, params, tau)</p><p> results.append({'velocity_m_s':v, 'mean_abs_M_A_m':np.mean(np.abs(M_t)), 'peak_abs_M_A_m':np.max(np.abs(M_t))})</p><p> plt.figure()</p><p> plt.plot(t, M_t/params)</p><p> plt.xlabel("time (s)")</p><p> plt.ylabel("M(t)/Ms")</p><p> plt.title(f"v={v} m/s")</p><p> plt.grid(True)</p><p><a href="http://plt.show">plt.show</a>()</p><p>df = pd.DataFrame(results)</p><p>print(df)</p><p><a href="http://df.to">df.to</a>_csv("magnetization_vs_velocity_improved.csv", index=False)</p><p>print("Saved: magnetization_vs_velocity_improved.csv")</p> |
|
Project ID:
|
|
4816636 |
|
Категория проекта: |
|
|
|
Бюджет проекта: |
|
|
|
|
|
|
|
|
|
| Проект |
Открыт |
Бухгалтерская услуга
Категория: Аутсорсинг и консалтинг, Бухгалтерия |
30 Nov 2025 11:05 GMT |
|
Мне надо бот сигнал
|
30 Nov 2025 11:02 GMT |
Накрутить ботов в группу ВК
Категория: Социальные сети, VK |
30 Nov 2025 11:00 GMT |
Edycja zdjecia
Категория: Обработка фото |
30 Nov 2025 10:59 GMT |
Seo продвинуть сайт в топ 3 google и яндекс
Категория: Продвижение сайтов (SEO, GEO), Арбитраж трафика
Бюджет: 10000 руб |
30 Nov 2025 10:46 GMT |
Разработчик на Битрикс
Категория: Программирование, Веб-программирование |
30 Nov 2025 10:44 GMT |
Рерайт на антиплагиат.
Категория: Тексты, Рерайтинг
Бюджет: 1500 руб |
30 Nov 2025 10:44 GMT |
|
Реклама, лиды в телеграм канал по криптообмену
|
30 Nov 2025 10:38 GMT |
Нарезчик на стримы (для всех)
Категория: Аудио/Видео/Фото, Видеомонтаж |
30 Nov 2025 10:38 GMT |
ЮРИСТ в области наружной рекламы
Категория: Реклама и Маркетинг, Медиапланирование |
30 Nov 2025 10:27 GMT |
Проектування одноповерхового будинку 8х9 м з плоским дахом, виконаний з газоблоку, з терасою. - 5000UAH
Категория: Дизайн интерьеров, Чертежи и схемы |
30 Nov 2025 10:24 GMT |
Парсер новых товаров для доски объявлений на запросах
Категория: Программирование, Создание скриптов
Бюджет: 45000 руб, для всех |
30 Nov 2025 10:20 GMT |
|
Создание презентации для автосалона на китайском и украинском
|
30 Nov 2025 10:19 GMT |
Картинки из нейросетей
Категория: Дизайн, Баннеры
Бюджет: 50000 руб |
30 Nov 2025 10:06 GMT |
Разработка лого для ювелирного бренда
Категория: Дизайн, Дизайн / Логотипы |
30 Nov 2025 09:58 GMT |
|
|
Все проекты
|
|
|
|
|
|
