|
|
|
|
|
Название проекта: |
Мне нужно смоделировать движение стального стержня через магнит
|
|
Кто разместил: |
Внешний проект с 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 |
|
Категория проекта: |
|
|
Бюджет проекта: |
|
|
|
|
|
|
|
|
| Проект |
Открыт |
|
Менеджер по привлечению клиентов
|
15-Feb-2026
21:07 GMT |
Нужен конверсионный сайт на Tilda под Яндекс Директ (бухгалтерские услуги) (для всех)
Категория: Сайты, Веб-дизайн |
15-Feb-2026
20:11 GMT |
|
Создание персонажа
|
15-Feb-2026
20:00 GMT |
Сделать лендинг + настроить все формы и интеграции (заявка → оплата → календарь → CRM)
Категория: Сайты, Zero-coding |
15-Feb-2026
19:15 GMT |
Дизайн интерьера кухни-гостиной
Категория: Дизайн, Интерьеры |
15-Feb-2026
19:15 GMT |
Создать ИИ агента для программы Superslicer
Категория: Программирование, Прикладное программирование
Бюджет: 10 000 руб |
15-Feb-2026
19:15 GMT |
Продвижение в строительной отрасли
Категория: Реклама и Маркетинг, Контекстная реклама |
15-Feb-2026
19:14 GMT |
Умная поддержка в боте через n8n
Категория: AI — искусственный интеллект, N8n |
15-Feb-2026
19:14 GMT |
Поиск контактов элитных сервисов (клининг, капельницы, салоны массажа
Категория: Аутсорсинг и консалтинг, Обслуживание клиентов и поддержка |
15-Feb-2026
19:14 GMT |
Автоматизация через Make.com
Категория: Автоматизация бизнеса, Make.com |
15-Feb-2026
19:13 GMT |
Срочный перенос сервера с VPS на VDS
Категория: Сайты, Техническая поддержка
Бюджет: 10 000 руб |
15-Feb-2026
19:13 GMT |
установить brainycp на centos7
Категория: Программирование, Системное программирование |
15-Feb-2026
19:12 GMT |
Разработка ПОД (проектов сноса) 5 зданий
Категория: Инжиниринг, ArchiCAD
Бюджет: 50 000 руб |
15-Feb-2026
19:11 GMT |
создание сайта
Категория: Сайты, Интернет-магазины |
15-Feb-2026
19:11 GMT |
Экспрессивно зачитать текст на камеру
Категория: Аудио/Видео/Фото, Диктор |
15-Feb-2026
19:10 GMT |
|
|
Все проекты
|
|
|
|
|
|
