[LISP程序]：求助空间钢束计算的LSP程序?

zhouyxcs

在预应力混凝土结构的设计中常遇到空间钢束的长度和引伸量的计算.一般是分别绘出立面和平面图.这时如何选取了这两个图就能计算出长度.另外还有引伸量的计算.求助高手.

maplesu

其实求出空间的钢束长度，要是能准确画出就知道长度，但是难点就在于怎么画出，同时有平弯竖弯的是后这个曲线就是两个拄面的交线，是很难画出的。在施工中估计也很难精确作出，其实在施工中一般应该都是采取球面的做法。其实我一般是把平弯和竖弯尽量不同时发生，也就是说可以在竖弯完成以后在进行平弯，这样分别计算长度就容易计算，而且需要有平弯的时候几乎都是要弯到腹板中的吧，所以在到竖弯点前完成平弯也是可以很好的不和腹板打架的。以上都是个人的意见。

至于引伸量我倒是有一个程序了，不过做一端张拉和两端张拉（需要是对称的）还可以，两段不对称的就不行了，因为每边的损失是不同的，很难找到向两边力相等的那个点，加之两边的引伸量也不同。所以我也没有做到。

你看看可以不？？这个也是只有单纯的一个方向的弯。


[PHP]
(defun c:scl ()
  (setvar "cmdecho" 0)  
  (setq Ay 1668)
  (setq Eg 195000)
  (setq u 0.25)
  (setq k 0.0015)
  (setq MPa 1860)
  (setq pv (* 0.75 (* MPa Ay)))
  (setq en (entsel "objet"))
  (setq sum 0)
  (while (/= en nil)
    (command "lengthen" (car en) "")
    (setq l (getvar "perimeter"))
    (setq en_data (entget (car en)))
    (setq kind (cdr (assoc 0 en_data)))
    (if (= kind "ARC")
      (progn
        (setq qd (CDR (ASSOC 50 en_data)))
        (setq zd (CDR (ASSOC 51 en_data)))
        (setq a (abs (- zd qd)))
      )
    )
    (if (= kind "LINE")
      (setq a 0)
    )
    (if (= kind "SPLINE")
      (progn
        (setq n 0)
        (setq i 0)
        (setq az 0)
        (setq aa 0)
        (setq aaa 0)
        (repeat (length en_data)
          (setq pp (nth n en_data))
          (setq key (car pp))          
          (if (= key 10)
            (progn
              (setq x (cadr pp))
              (setq y (caddr pp))
              (if (/= i 0)
                (progn
                   (setq zx (- x xx))
                  (setq zy (- y yy))
                  (setq a1 (atan (abs zy) (abs zx)))
                  (if (> i 1) (setq aaa (abs (- a1 aa))))                    
                  (setq aa a1)                  
                )
                (setq a1 0)
              )
              (setq xx x)
              (setq yy y)
              (setq i (1+ i))
              (setq az (+ az aaa))
            )
          )
          (setq n (1+ n))
        )
        (setq a az)
      )
    )   
    (setq pv (* pv (/ (+ 1 (exp (- 0 (+ (* u a) (* k (/ l 100)))))) 2)))
    (setq dl (/ (* pv l) (* Ay Eg)))
    (setq sum (+ sum dl))
    (setq en (entsel "objet:"))
  )
  (alert (strcat "伸长量=" (rtos sum 2 2)))
  (prin1)
)
(prin1)
[/PHP]


还有我把需要输入的参数值直接复值的，需要改的话在程序里面改就可以了，钢束图要准确的画在cad中，然后一个一个选取每段钢束就可以了。（对称的只取一半，然后乘2就可以了）