[LISP程序]：一个矩形截面受弯计算程序

yangleisd

2. [php];--------------------------------------+
   
3. ;文件名：M-AS.lsp                      +
   
4. ；计算梁板受弯配筋计算　　　　　　　 　+
   
5. ;--------------------------------------+
   
6. (defun c:M-AS( / c lx b h gj fy ksib m mode fc h0 tt w as)
   
7. (inivar) ;初始化系统变量
   
8. (textscr)
   
9. (setq c (if (setq temp (getreal "\n砼强度等级:<20>"))  temp 20)
   
10.        lx (if (= (setq temp (getstring "计算板(1)/梁(2):<1>")) "") "1" temp)
    
11.        b (if (= lx "2") (getreal "梁宽(mm)：") 1000)
    
12.        h (if (= lx "1") (getreal "板厚(mm):")
    
13.                         (getreal "梁高(mm):"))
    
14.        gj (if (setq temp (getreal "钢筋等级:<2>"))  temp 2)
    
15.        fy (if (= gj 1) 210 310)
    
16.        ksib (if (= gj 1) 0.614 0.544) ;极限受压区高度
    
17.        m (getreal "弯矩(KN-m):")
    
18.        mode (getint "\n按旧规范GB10-89(0)/按新规范GB50010-2002(1)计算:<1>")
    
19.        mode (if mode mode 1)
    
20. ) ;setq end
    
21. (if (null (type concrete)) (load "concrete.lsp"))
    
22. (setq fc (concrete "fc" c mode)) ;查询混凝土轴心受压强度
    
23. (if (= mode 1) (setq fc (* fc (alfa1 c)))) ;如果按新规范
    
24. (setq h0 (if (= lx "1") (- h 20) (- h 35)) ;计算有效受压区高度
    
25.        tt (/ (* m 1e6) (* fc b h0 h0))      
    
26.        w (- 1 (sqrt (- 1 (* 2.0 tt))))
    
27.        as (/ (* fc b h0 w) fy )
    
28. ) ;setq end
    
29. (if (<= w ksib)
    
30.    (progn
    
31.       (princ (strcat "\n相对受压区高度= " (rtos w 2 4)))
    
32.       (princ (strcat "\n配筋面积= " (rtos as 2 2) "(mm2)" ))
    
33.    ) ;progn end
    
34.    (progn
    
35.       (princ (strcat "\n相对受压区高度=" (rtos w 2 4) " >极限受压区高度=" (rtos ksib 2 4)))
    
36.       (princ "\n截面不够,请加大截面!")
    
37.    ) ;progn
    
38. ) ;if end
    
39. (resvar) ;还原系统变量
    
40. (princ)
    
41. ) ;defun end[/php]

[php];--------------------------------------+<br /> ;文件名：M-AS.lsp                      +<br /> ；计算梁板受弯配筋计算　　　　　　　 　+<br /> ;--------------------------------------+<br /> (defun c:M-AS( / c lx b h gj fy ksib m mode fc h0 tt w as) <br /> (inivar) ;初始化系统变量<br /> (textscr)<br /> (setq c (if (setq temp (getreal "\n砼强度等级:<20>"))  temp 20) <br />        lx (if (= (setq temp (getstring "计算板(1)/梁(2):<1>")) "") "1" temp)<br />        b (if (= lx "2") (getreal "梁宽(mm)：") 1000)<br />        h (if (= lx "1") (getreal "板厚(mm):") <br />                         (getreal "梁高(mm):"))<br />        gj (if (setq temp (getreal "钢筋等级:<2>"))  temp 2)<br />        fy (if (= gj 1) 210 310)<br />        ksib (if (= gj 1) 0.614 0.544) ;极限受压区高度<br />        m (getreal "弯矩(KN-m):")<br />        mode (getint "\n按旧规范GB10-89(0)/按新规范GB50010-2002(1)计算:<1>")<br />        mode (if mode mode 1)<br /> ) ;setq end<br /> (if (null (type concrete)) (load "concrete.lsp")) <br /> (setq fc (concrete "fc" c mode)) ;查询混凝土轴心受压强度<br /> (if (= mode 1) (setq fc (* fc (alfa1 c)))) ;如果按新规范<br /> (setq h0 (if (= lx "1") (- h 20) (- h 35)) ;计算有效受压区高度<br />        tt (/ (* m 1e6) (* fc b h0 h0))      <br />        w (- 1 (sqrt (- 1 (* 2.0 tt))))<br />        as (/ (* fc b h0 w) fy ) <br /> ) ;setq end<br /> (if (<= w ksib) <br />    (progn <br />       (princ (strcat "\n相对受压区高度= " (rtos w 2 4)))<br />       (princ (strcat "\n配筋面积= " (rtos as 2 2) "(mm2)" ))<br />    ) ;progn end<br />    (progn <br />       (princ (strcat "\n相对受压区高度=" (rtos w 2 4) " >极限受压区高度=" (rtos ksib 2 4)))<br />       (princ "\n截面不够,请加大截面!")<br />    ) ;progn<br /> ) ;if end<br /> (resvar) ;还原系统变量 <br /> (princ)<br /> ) ;defun end[/php]

njsabc

非常值得学习的源程序