[Bezzvyka]

есть 2 пересекающихся цилиндра под прямым углом и развертка выглядит примерно так (извиняюсь за кривизну рисовал в пейнте), нужно написать программу где мы будем вводить радиусы двух цилиндров и нам будет рисоваться эта развертка, уравнение нижней синусойды y=R-(R^2-(2*r*sinx)^2)^(1/2), где х [0; 2*PI*r], верхняя линия равна 2*PI*r, боковые пусть буду r (R-радиус большого цилиндра, r-радиус цилиндра поменьше сообственно которого и производиться развертка. Смог сделать только 3 линии, как написать код для построения этой синусойды не знаю, помогите пожалуйста.

[gomer]

строй полилинию по точкам, так же, как и вручную

[Bezzvyka]

Я тебя не понял, я смог только сделать

(defun C:rzv (/ P1 P2 P3 P4) (setq P1 (getpoint "начальная точка:")) (setq L (getdist P1 "длина квадрата:")) (setq P2 (polar P1 0.0 L)) (setq P3 (polar P2 (/ PI 2) (* 2 PI L))) (setq P4 (polar P3 PI L)) (command "_LINE" P1 P2 P3 P4 "С" ) )

мне бы хотя бы код для построения синусойды y=R-(R^2-(2*r*sinx)^2)^(1/2), где х [0; 2*PI*r]

[VVA]

Добавил в формулу ABS (модуль) [Выделил в коде красным], т.к.
при R=20, r=14 и x=1
формула R^2-(2*r*sinx)=-155,3 и извлечение корня дает ошибку

Код:

[Выделить все]

(defun C:rzv (/ P1 P2 Rbig Rsmall X Y Xmax lst P3 step)
  (initget 3)
  (setq Rbig (getdist "\nРадиус большого цилидра: "))
  (initget 3)
  (while (>= (setq Rsmall (getdist "\nРадиус маленького цилидра: "))
            Rbig)
    (princ "\nРадиус должен быть мешьше чем ")
    (princ Rbig)
    (initget 1)
    )
(initget 3)
(setq Step (getdist "\nШаг по оси X: "))
(initget 1)  
(setq P1 (getpoint "\nНачальная точка: "))
(command "_.PLINE" "_non" P1
         "_non" (setq P2 (polar P1 (* 0.5 pi) Rsmall))
         "_non" (setq P2 (polar P2 0 (* 2 pi Rsmall)))
         "_non" (setq P2 (polar P2 (* 1.5 pi) Rsmall))
         ""
         )
  ;;;Построение синусоиды
  ;;;Формула y=R-(R^2-(2*r*sinx)^2)^(1/2), где x [0; 2*PI*r]
  (setq X (- 0. step)
        Xmax (* 2 pi Rsmall)
        )
  (while (<= (setq X (+ X step)) Xmax)
    (setq P3 (- (* Rbig Rbig)(*(* 2 Rsmall (sin x))(* 2 Rsmall (sin x)))))
    (if (minusp P3)
      (setq P3 (- 0. (sqrt(abs P3))))
      (setq P3 (sqrt P3))
      )
    (setq Y (- Rbig P3))
    (setq lst (cons (list X Y) lst))
    )
  (command "_.PLINE" "_non" P1)
  (foreach pt (reverse lst)
    (command "_non" (mapcar '+ P1 pt))
    )
  (command P2 "")
  (while (> (getvar "CMDACTIVE") 0) (command))
  (princ)
)

[gomer]

VVA,не пойдет программа, развертка строится по количеству опорных точек, а не по расстоянию между ними. Во-вторых, в результате получается черти что, а не развертка.

[VVA]

gomer, количество опорных точек=длина развертки/шаг
Внес изменения в #4. На синусоиду похоже больше.

[Boxa]

Нда... студенту задали задачу, а Вы за него решение пишите. Не стыдно?

[gomer]

Сообщение от Boxa Вы за него решение пишите. Не стыдно?

Не умничайте. VVA развертка не должна пересекать себя, попробуйте диаметры 500, 800 и шаг 156

[VVA]

gomer, Я реализовал формулу. Вспоминать начерталку про развертки мне некогда

[alex8888]

Поскольку я вплотную занимаюсь этим чудом, то сделал для себя несколько функций для разверток. Все это узкоспециализировано под мои задачи, но работает безотказно - не одна сотня деталей сделана по ним. Может быть поможет.
В архив кинул вроде бы все причастные файлы для развертки отвода, ненужное можно выбросить.
Основная функция "at_stutzen", остальные вспомогательные.

[gomer]

где I:AT_LIST->ARRAY?

[alex8888]

Может быть еще какую пропустил, пишите. Так уж повелось - все под разные проекты, не писать же в каждый одно и тоже.

[gomer]

Вроде оно, но в чем разница не могу разобраться пока, вроде формула одна, а результат разный...

[alex8888]

Вроде оно, но в чем разница не могу разобраться пока, вроде формула одна, а результат разный...

А в чем проблема?

[gomer]

Да, похоже что VVA линейные и угловые величины намешал. А общий случай пересечения цилиндров у вас есть?

[alex8888]

Что значит общий случай? Что конкретно?
В моем лиспе можно вместо параметров указать nil или просто () и тогда он будет спрашивать шаг за шагом диаметр трубы, диаметр отвода, высоту от оси трубы до верхней плоскости отвода (это чтобы высоту развертки подсчитать), потом ненужные для вас, но нужные мне текстовые пометки на развертку, можно пропустить нажатием ввода, наконец спросит отступ от продольной оси трубы до оси отвода, в случае если плоскости совпадают, т.е. равно 0, то можно оставить пустым и нажать ввод. Последним запросит точку куда вставить развертку, путем тыка мышкой, от нее потом считается вся развертка. Хотелось бы, чтобы можно было уже потом, после построения указывать где разместить получившуюся полилинию, при этом чтобы она "висела" на курсоре, но при этом не являлась бы блоком - блок потом расчленять надо, а это лишние телодвижения.

[gomer]

Сообщение от alex8888 Что значит общий случай? Что конкретно?

Штуцер может быть и под углом

[alex8888]

Сообщение от gomer Штуцер может быть и под углом

Такая ситуация в моей практике еще не встречалась, поэтому, естественно , я ее не обыгрывал. Единственное, что в расчеты надо принимать угол наклона оси штуцера к оси трубы. В случае, когда цилиндры пересекаются под углом 90 град. синус угла равен 1, в остальных случаях нет и его как-то нужно добавлять к расчетам. Спинным мозгом чувствую, что сложного ничего нет, просто где-то надо на этот синус поделить, возможно в строке 95-96:

Код:

[Выделить все]

 (- Stutzen_hoehe (sqrt (- (expt (* 0.5 Kopf_D) 2.0)(expt (+ (* m (* 0.5 Stutzen_D)) Stutzen_Abstand) 2.0))))

, но не разбирался с этим. Если есть желание разобраться, сообщи результат потом, хорошо?

[VVA]

Сообщение от gomer Да, похоже что VVA линейные и угловые величины намешал.

У меня так же на это есть подозрения.
Бегло просмотрел код alex8888, у него

(setq a (list n) ;создание пустого списка делений длины отвода winkel_list (list (* 2 pi)) ;создание списка углов

а у ТС

для построения синусоиды y=R-(R^2-(2*r*sinx)^2)^(1/2), где х [0; 2*PI*r]

Думаю собака зарыта где-то здесь