|
|
|
|
|
Название проекта: |
|
Мне нужно смоделировать движение стального стержня через магнит
|
|
Кто разместил: |
|
Внешний проект с 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 |
|
Категория проекта: |
|
|
|
Бюджет проекта: |
|
|
|
|
|
|
|
|
|
| Проект |
Открыт |
Добавление страниц товаров Elementor WP
Категория: Сайты, WordPress |
19 Apr 2026 09:57 GMT |
Телеграмм бот по поиску Бригад строительных
Категория: Программирование, Python
Бюджет: 100 000 руб |
19 Apr 2026 09:45 GMT |
Авитолог
Категория: Маркетплейс менеджмент, Авито Молл |
19 Apr 2026 09:42 GMT |
Локальный сайт с мини рассказами.
Категория: Сайты, Сайт «под ключ» |
19 Apr 2026 09:41 GMT |
Разработка ПО для Квадракоптера
Категория: Программирование, Google Android |
19 Apr 2026 09:27 GMT |
Написать короткий текст на картинке
Категория: Дизайн, Дизайн для соцсетей |
19 Apr 2026 08:31 GMT |
Иллюстрация водяной юлы
Категория: Рисунки и иллюстрации
Бюджет: 5 000 руб |
19 Apr 2026 08:28 GMT |
Сборка Windows-приложения из готового Flutter-проекта
Категория: Программирование, Интерактивные приложения |
19 Apr 2026 08:15 GMT |
Web MVP (Joomla)
Категория: Сайты, Создание MVP |
19 Apr 2026 08:03 GMT |
Парсинг данных с приложения
Категория: Программирование, Парсинг данных
Бюджет: 5 000 руб |
19 Apr 2026 07:55 GMT |
Консультация менеджера WB с опытом работы в категории Духи
Категория: Маркетплейс менеджмент, Wildberries |
19 Apr 2026 07:15 GMT |
Составление ведомости маткриалов пр сметам (для всех)
Категория: Тексты, Технический перевод |
19 Apr 2026 07:14 GMT |
Провести техническую seo оптимизацию сайта на битрикс
Категория: Сайты, SEO-оптимизация
Бюджет: 8 000 руб |
19 Apr 2026 07:08 GMT |
Внедрить чат бот (для всех)
Категория: AI — искусственный интеллект, YandexGPT |
19 Apr 2026 06:45 GMT |
Нужен таргетолог со своим рекламным кабинетом
Категория: Социальные сети, Реклама в соцсетях |
19 Apr 2026 06:45 GMT |
|
|
Все проекты
|
|
|
|
|
|
