|
|
|
|
|
Название проекта: |
Мне нужно смоделировать движение стального стержня через магнит
|
|
Кто разместил: |
Внешний проект с 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 |
|
Категория проекта: |
|
|
Бюджет проекта: |
|
|
|
|
|
|
|
|
| Проект |
Открыт |
Футуристичный дизайн дисков для машин
Категория: Дизайн, Фирменный стиль
Бюджет: 8 000 руб |
05-Apr-2026
14:48 GMT |
Работы в еду 1с фреш
Категория: Аутсорсинг и консалтинг, Ввод и обработка данных/текста
Бюджет: 2 000 руб |
05-Apr-2026
14:32 GMT |
Сделать макет логотипа
Категория: Дизайн, Дизайн / Логотипы
Бюджет: 500 руб |
05-Apr-2026
14:28 GMT |
Специалист Колл-центра
Категория: Реклама и Маркетинг, Продажи и генерация лидов |
05-Apr-2026
14:24 GMT |
Жудожник чтобы продолжить Мультфильм Gravity Falls
Категория: Дизайн, AR дополненная реальность |
05-Apr-2026
14:10 GMT |
Консультация маркетолога
Категория: Реклама и Маркетинг, Маркетинговая стратегия |
05-Apr-2026
14:06 GMT |
Восстановить сайт из Резервной копии на хостинге Джино
Категория: Программирование, Веб-программирование |
05-Apr-2026
13:57 GMT |
Разработка проектной документации с дизайном и сметами
Категория: Дизайн, Архитектура |
05-Apr-2026
13:47 GMT |
Создать реалистичную ИИ-девушку
Категория: AI — искусственный интеллект, AI-аватары
Бюджет: 120 000 руб |
05-Apr-2026
13:32 GMT |
Рисунок для плаката в виде комикса
Категория: Рисунки и иллюстрации
Бюджет: 2 500 руб |
05-Apr-2026
13:16 GMT |
Нужен дизайн сайта
Категория: Дизайн, Веб-дизайн |
05-Apr-2026
13:16 GMT |
Дизайны ковров под мышки
Категория: Дизайн, Трикотажный и текстильный дизайн
Бюджет: 7 500 руб |
05-Apr-2026
13:15 GMT |
Анимации графиков криптовалюты
Категория: Анимация, 2D Анимация
Бюджет: 6 000 руб |
05-Apr-2026
13:14 GMT |
Ищу монтажера видео MInecraft
Категория: Аудио/Видео/Фото, Видеомонтаж
Бюджет: 150 руб |
05-Apr-2026
13:14 GMT |
DevOps–аудит и стабилизация AI–инфраструктуры (n8n, Docker, LLM)
Категория: Автоматизация бизнеса, N8n |
05-Apr-2026
13:13 GMT |
|
|
Все проекты
|
|
|
|
|
|
