Оружие будущего

Форум про оружие будущего
Текущее время: Вт май 21, 2019 5:41 pm

Часовой пояс: UTC + 3 часа




Начать новую тему Ответить на тему  [ Сообщений: 64 ]  На страницу Пред.  1, 2, 3, 4
Автор Сообщение
 Заголовок сообщения:
СообщениеДобавлено: Ср окт 08, 2008 9:53 pm 
Админ
Аватар пользователя

Зарегистрирован: Пн июл 23, 2007 12:02 pm
Сообщений: 10860
Откуда: Пермь
sqr - квадрат
sqrt - квадратный корень


Вернуться наверх
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: Ср окт 08, 2008 9:58 pm 
СМЕРТЬ!
Аватар пользователя

Зарегистрирован: Пн апр 28, 2008 7:24 pm
Сообщений: 233
Откуда: Москва, Вокзал у Припяти
вот блин. Для меня sqr и sqrt было одно и тоже...

_________________
Адский Гауссган!


Вернуться наверх
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: Пн окт 20, 2008 6:13 pm 
СМЕРТЬ!
Аватар пользователя

Зарегистрирован: Пн апр 28, 2008 7:24 pm
Сообщений: 233
Откуда: Москва, Вокзал у Припяти
Короче скрипт для максимы я написал, но почему-то конечный результат получается не таким каким должен быть. Видимо где-то вкралась ошибка.
Код:
U:450;    /*В напряжение на конденсаторах*/
C:500;    /*мкФ ёмкость конденсаторов*/
R:100;    /*мОм сопротивление цепи(без катушки)*/
D:500;    /*мкм диаметр провода*/
d1:8;     /*мм внутренний диаметр катушки*/
d2:24;    /*мм внешний диаметр катушки*/
l:16;     /*мм длинна катушки*/
Lcoil:0;/*Гн индуктивность катушки*/
Rcoil:0;/*Ом сопротивление катушки*/
N:0;   /*шт количество витков*/
/*------------- CODE START ---------*/
C:C*10^(-6);/*Перевод величин*/
R:R*10^(-3);
D:D*10^(-6);
d1:d1*10^(-3);
d2:d2*10^(-3);
l:l*10^(-3);
po:0.0172*10^(-6);/*Удельное сопротивление меди*/
/*----------------------------------*/
dm:(d1+d2)/2;
S:%pi*D^2/4;
if N=0 then N:(d2-d1)*l/D^2;
Len:%pi*N*dm;
if Rcoil=0 then Rcoil:po*Len/S;
if Lcoil=0 then Lcoil:2*(dm*N)^2/(dm+3*(l+(d2-d1)));
L:Lcoil;
R:Rcoil+R;
/*--------- ГРАФИКИ ---------------*/
if 1/(L*C)>(R/(2*L))^2 then flag:0;
if flag=0 then I(t):=U/(sqrt(1/(L*C))*L)*sin(sqrt(1/(L*C))*t)*exp(-R/(2*L)*t);
if flag=0 then plot2d(I(t),[t,0,6*10^(-3)]);
if 1/(L*C)<(R/(2*L))^2 then flag:1;
if flag=1 then I(t):=U/(2*sqrt((R/(2*L))^2-1/(L*C))*L)*(exp((sqrt((R/(2*L))^2-1/(L*C))-R/(2*L))*t)-exp((-sqrt((R/(2*L))^2-1/(L*C))-R/(2*L))*t));
if flag=1 then plot2d(I(t),[t,0,6*10^(-3)]);

Если найдёте ошибку в формулах, сообщите. Ещё осталось дописать оптимальный диаметр, повидимому просто найти максимум.

ДОБАВЛЕНО: Кажется понял где. Количество витков получилось слишком маленьким, пока исправляю.

Уже получился график......

_________________
Адский Гауссган!


Вернуться наверх
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: Пн окт 20, 2008 9:07 pm 
In Trance
Аватар пользователя

Зарегистрирован: Сб дек 22, 2007 2:45 pm
Сообщений: 4880
Откуда: Украина, Днепропетровск
Llevellyn
Тут насколько мне не изменяет мой моск, в твоей проге нужно создать цикл, как в фемме т.е. пересчитывать индуктивность катушки и силу втягивания снаряда через некоторый промежуток времени. Иначе у тебя получится что-то на подобии CoilGun Simulator.

_________________
ZADROT MODE COMPLETE
MISSION PASSED


Вернуться наверх
 Профиль  
 
Показать сообщения за:  Сортировать по:  
Начать новую тему Ответить на тему  [ Сообщений: 64 ]  На страницу Пред.  1, 2, 3, 4

Часовой пояс: UTC + 3 часа


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
cron

 

Яндекс цитирования
Русская поддержка phpBB