|
|
|
|
|
Название проекта: |
Мне нужно смоделировать движение стального стержня через магнит
|
|
Кто разместил: |
Внешний проект с 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 |
|
Категория проекта: |
|
|
Бюджет проекта: |
|
|
|
|
|
|
|
|
| Проект |
Открыт |
АВТОНОМНАЯ ЭКОСИСТЕМА БЕЗ ПОДПИСОК СВОЯ
Категория: Программирование, Fullstack |
31-Jan-2026
11:02 GMT |
Перекрашивание техники, размещение на цветном фоне
Категория: Дизайн, Веб-дизайн |
31-Jan-2026
11:00 GMT |
Финансовое задание с расчетом
Категория: Аутсорсинг и консалтинг, Финансовый консультант
Бюджет: 1 500 руб |
31-Jan-2026
10:43 GMT |
|
Перевод статьи с английского на русский
|
31-Jan-2026
10:30 GMT |
SEO продвижение сайта по ГДЗ (для всех)
Категория: Продвижение сайтов (SEO, GEO), Поисковые системы |
31-Jan-2026
10:18 GMT |
Внести правки в сайт на ModX
Категория: Сайты, MODx |
31-Jan-2026
10:14 GMT |
|
Перенос презентации из Figma в PowerPoint
|
31-Jan-2026
10:11 GMT |
Размещение отзывов
Категория: Реклама и Маркетинг, SERM (управление репутацией) |
31-Jan-2026
10:10 GMT |
Написать скрипт на BAS
Категория: Программирование, Создание скриптов |
31-Jan-2026
10:09 GMT |
|
Написать скрипт на BAS
|
31-Jan-2026
10:07 GMT |
битрикс правки на сайте
Категория: Сайты, 1С Битрикс |
31-Jan-2026
09:58 GMT |
|
Нужен ручной тестер
|
31-Jan-2026
09:53 GMT |
N8N автоматизация (для всех)
Категория: AI — искусственный интеллект, N8n |
31-Jan-2026
09:09 GMT |
Создание 3D модели для фрезеровки на ЧПУ (барельеф) декоративной накладки по фото с доработками
Категория: 3D Графика, 3D Моделирование
Бюджет: 2 000 руб |
31-Jan-2026
08:36 GMT |
Создание MVP
Категория: Mobile, Создание MVP |
31-Jan-2026
08:31 GMT |
|
|
Все проекты
|
|
|
|
|
|
