明经CAD社区

 找回密码
 注册

QQ登录

只需一步,快速开始

搜索
查看: 1089|回复: 7

转子抽速公式编程

[复制链接]
发表于 2016-5-6 11:31:26 | 显示全部楼层 |阅读模式
求大神指教: 下面是螺杆泵的抽速公式,要怎样调用程序才能求出最后的sth
(defun qiu_acos (num);;定义arccos
  (if (= num 0) (atan 1 0)
                (atan (/ (sqrt (- 1 (* num num))) num))))

(defun tan(mm) ;;定义tan
  (/ (sin mm) (cos mm)))

(defun chousu (r r0 r1 R n p / s2 s4 s1 s3 c)
  (setq aaa (/ r0 R))
  (setq bbb (/ r1 R))
  (setq cc (qiu_acos aaa))
  (setq dd (qiu_acos bbb))
  (setq ee (tan cc))
  (setq s2 (* (/ (* R R) 2) (- pi (- (/ pi 3) (sqrt 3)) (- ee cc))))
  (setq s4 (* (/ (* r r) 2) (+ pi (- (/ pi 3) (sqrt 3)) (- ee cc))))
  (setq s1 (* (/ 1 6) (* r0 r0) (expt ee 3)))
  (setq s3 (- (+ (* (+ (* 2 r1 r1) (* R R)) f) (* (* r1 r1) (sin (* 2 f))) (* (/ 1 4) (* R R) (sin (* 4 f)))) (+ (* 3 r1 R (sin f)) (* r1 R (sin (* 3 f))))))
  (setq f dd)
  (setq c (- 1 (/ (- (+ s1 s2 s4) s3) (* pi R R))))
  (setq sth (* 2 n c p pi R R))
  )

(setq qq 30 ww 30 tt1 60 yy 90 uu 3000 ii 35)
(chousu qq ww tt1 yy uu ii)
(princ sth)


 楼主| 发表于 2016-5-6 11:52:46 | 显示全部楼层
求大神啊求大神
发表于 2016-5-6 18:45:35 | 显示全部楼层
1、函数chousu定义中不允许有两个R参数(r和R是相同意义的),可将R改为R3等
2、求s3算式中,f参数无值,当然不能求出S3值,
 楼主| 发表于 2016-5-7 12:13:51 | 显示全部楼层
USER2128 发表于 2016-5-6 18:45
1、函数chousu定义中不允许有两个R参数(r和R是相同意义的),可将R改为R3等
2、求s3算式中,f参数无值,当然 ...

你好,谢谢你的回复。
程序里面是给f赋值了的
(setq dd (qiu_acos bbb))
(setq f dd)
将R改为R3之后显示参数太少
 楼主| 发表于 2016-5-7 12:15:45 | 显示全部楼层
USER2128 发表于 2016-5-6 18:45
1、函数chousu定义中不允许有两个R参数(r和R是相同意义的),可将R改为R3等
2、求s3算式中,f参数无值,当然 ...

您好,感谢您的回复
我是给f赋值了的
(setq f dd)
我将R改为R3之后结果还是显示参数太少
发表于 2016-5-7 17:34:01 | 显示全部楼层
  1. ;;; 求大神指教: 下面是螺杆泵的抽速公式,要怎样调用程序才能求出最后的sth
  2. (defun qiu_acos (num);;定义arccos
  3.   (if (= num 0) (atan 1 0)
  4.                 (atan (/ (sqrt (- 1 (* num num))) num))))

  5. (defun tan(mm) ;;定义tan
  6.   (/ (sin mm) (cos mm)))

  7. (defun chousu (r r0 r1 R3 n p / s2 s4 s1 s3 c)
  8.   (setq aaa (/ r0 R3))
  9.   (setq bbb (/ r1 R3))
  10.   (setq cc (qiu_acos aaa))
  11.   (setq dd (qiu_acos bbb))
  12.   (setq ee (tan cc))
  13.   (setq s2 (* (/ (* R3 R3) 2) (- pi (- (/ pi 3) (sqrt 3)) (- ee cc))))
  14.   (setq s4 (* (/ (*  r  r) 2) (+ pi (- (/ pi 3) (sqrt 3)) (- ee cc))))
  15.   (setq s1 (* (/ 1 6) (* r0 r0) (expt ee 3)))
  16.   (setq f dd)
  17.   (setq s3 (- (+ (* (+ (* 2 r1 r1) (* R3 R3)) f) (* (* r1 r1) (sin (* 2 f))) (* (/ 1 4) (* R3 R3) (sin (* 4 f)))) (+ (* 3 r1 R3 (sin f)) (* r1 R3 (sin (* 3 f))))))
  18.   (setq c (- 1 (/ (- (+ s1 s2 s4) s3) (* pi R3 R3))))
  19.   (setq sth (* 2 n c p pi R3 R3))
  20.   )

  21. (setq qq 30 ww 30 tt1 60 yy 90 uu 3000 ii 35)
  22. (chousu qq ww tt1 yy uu ii)
  23. (princ sth)

  24. ;;; (setq f dd)要放到(setq s3 ... 之前
  25. ;;; 现可以计算其结果= 1.23464e+025 你自己核算一下公式.
 楼主| 发表于 2016-5-8 16:10:42 | 显示全部楼层
感谢贵人!!!
 楼主| 发表于 2016-5-8 16:18:29 | 显示全部楼层
USER2128 发表于 2016-5-7 17:34

感谢贵人!!!
另外还有个问题,如果已知了sth的值,怎样调用程序求R3呢?
您需要登录后才可以回帖 登录 | 注册

本版积分规则

小黑屋|手机版|CAD论坛|CAD教程|CAD下载|联系我们|关于明经|明经通道 ( 粤ICP备05003914号 )  
©2000-2023 明经通道 版权所有 本站代码,在未取得本站及作者授权的情况下,不得用于商业用途

GMT+8, 2024-11-24 00:17 , Processed in 0.177658 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表