home *** CD-ROM | disk | FTP | other *** search
/ Garbo / Garbo.cdr / mac / util / tkcmpsts.sit / mechloads / mechloads
Encoding:
Text File  |  1990-07-26  |  6.5 KB  |  792 lines  |  [TKSM/TKSR]
  1. TK+1
  2. =v
  3. #2
  4. :n
  5.  
  6. :s
  7. b
  8. :c
  9. --------
  10. #3
  11. :n
  12.  
  13. :s
  14. b
  15. :c
  16. Loadings
  17. #4
  18. :n
  19.  
  20. :s
  21. b
  22. :c
  23. --------
  24. #5
  25. :n
  26. N11
  27. :s
  28. b
  29. :u
  30. N/m
  31. :d
  32. N/m
  33. :c
  34. stress resultants
  35. #6
  36. :n
  37. N22
  38. :s
  39. b
  40. :u
  41. N/m
  42. :d
  43. N/m
  44. #7
  45. :n
  46. N12
  47. :s
  48. b
  49. :u
  50. N/m
  51. :d
  52. N/m
  53. #8
  54. :n
  55. M11
  56. :s
  57. b
  58. :u
  59. N
  60. :d
  61. N
  62. :c
  63. moment resultants
  64. #9
  65. :n
  66. M22
  67. :s
  68. b
  69. :u
  70. N
  71. :d
  72. N
  73. #10
  74. :n
  75. M12
  76. :s
  77. b
  78. :u
  79. N
  80. :d
  81. N
  82. #12
  83. :n
  84.  
  85. :s
  86. b
  87. :c
  88. ----------------
  89. #13
  90. :n
  91.  
  92. :s
  93. b
  94. :c
  95. Laminate Strains
  96. #14
  97. :n
  98.  
  99. :s
  100. b
  101. :c
  102. ----------------
  103. #15
  104. :n
  105. eps11
  106. :s
  107. b
  108. :u
  109. e
  110. :d
  111. ╡e
  112. :c
  113. in-plane strains
  114. #16
  115. :n
  116. eps22
  117. :s
  118. b
  119. :u
  120. e
  121. :d
  122. ╡e
  123. #17
  124. :n
  125. eps12
  126. :s
  127. b
  128. :u
  129. e
  130. :d
  131. ╡e
  132. #18
  133. :n
  134. kap11
  135. :s
  136. b
  137. :u
  138. 1/m
  139. :d
  140. 1/m
  141. :c
  142. curvatures
  143. #19
  144. :n
  145. kap22
  146. :s
  147. b
  148. :u
  149. 1/m
  150. :d
  151. 1/m
  152. #20
  153. :n
  154. kap12
  155. :s
  156. b
  157. :u
  158. 1/m
  159. :d
  160. 1/m
  161. =u
  162. =l
  163. #1
  164. :n
  165. NM
  166. :c
  167. matrix of stress/moment resultants
  168. :v
  169. E
  170. #2
  171. :n
  172. ek
  173. :c
  174. matrix of laminate strains/curvatures
  175. :v
  176. E
  177. #3
  178. :n
  179. spl11up
  180. :d
  181. MPa
  182. :u
  183. Pa
  184. :v
  185. E
  186. #4
  187. :n
  188. spl11dn
  189. :d
  190. MPa
  191. :u
  192. Pa
  193. :v
  194. E
  195. #5
  196. :n
  197. spl22up
  198. :d
  199. MPa
  200. :u
  201. Pa
  202. :v
  203. E
  204. #6
  205. :n
  206. spl22dn
  207. :d
  208. MPa
  209. :u
  210. Pa
  211. :v
  212. E
  213. #7
  214. :n
  215. spl12up
  216. :d
  217. MPa
  218. :u
  219. Pa
  220. :v
  221. E
  222. #8
  223. :n
  224. spl12dn
  225. :d
  226. MPa
  227. :u
  228. Pa
  229. :v
  230. E
  231. #9
  232. :n
  233. epl11up
  234. :d
  235. ╡e
  236. :u
  237. e
  238. :v
  239. E
  240. #10
  241. :n
  242. epl11dn
  243. :d
  244. ╡e
  245. :u
  246. e
  247. :v
  248. E
  249. #11
  250. :n
  251. epl22up
  252. :d
  253. ╡e
  254. :u
  255. e
  256. :v
  257. E
  258. #12
  259. :n
  260. epl22dn
  261. :d
  262. ╡e
  263. :u
  264. e
  265. :v
  266. E
  267. #13
  268. :n
  269. epl12up
  270. :d
  271. ╡e
  272. :u
  273. e
  274. :v
  275. E
  276. #14
  277. :n
  278. epl12dn
  279. :d
  280. ╡e
  281. :u
  282. e
  283. :v
  284. E
  285. #15
  286. :n
  287. spp11up
  288. :d
  289. MPa
  290. :u
  291. Pa
  292. :v
  293. E
  294. #16
  295. :n
  296. spp11dn
  297. :d
  298. MPa
  299. :u
  300. Pa
  301. :v
  302. E
  303. #17
  304. :n
  305. spp22up
  306. :d
  307. MPa
  308. :u
  309. Pa
  310. :v
  311. E
  312. #18
  313. :n
  314. spp22dn
  315. :d
  316. MPa
  317. :u
  318. Pa
  319. :v
  320. E
  321. #19
  322. :n
  323. spp12up
  324. :d
  325. MPa
  326. :u
  327. Pa
  328. :v
  329. E
  330. #20
  331. :n
  332. spp12dn
  333. :d
  334. MPa
  335. :u
  336. Pa
  337. :v
  338. E
  339. #21
  340. :n
  341. epp11up
  342. :d
  343. ╡e
  344. :u
  345. e
  346. :v
  347. E
  348. #22
  349. :n
  350. epp11dn
  351. :d
  352. ╡e
  353. :u
  354. e
  355. :v
  356. E
  357. #23
  358. :n
  359. epp22up
  360. :d
  361. ╡e
  362. :u
  363. e
  364. :v
  365. E
  366. #24
  367. :n
  368. epp22dn
  369. :d
  370. ╡e
  371. :u
  372. e
  373. :v
  374. E
  375. #25
  376. :n
  377. epp12up
  378. :d
  379. ╡e
  380. :u
  381. e
  382. :v
  383. E
  384. #26
  385. :n
  386. epp12dn
  387. :d
  388. ╡e
  389. :u
  390. e
  391. :v
  392. E
  393. =f
  394. #1
  395. :n
  396. MECH_DOC
  397. :c
  398. version 1.0
  399. :t
  400. Rule
  401. #2
  402. :n
  403. pelc
  404. :c
  405. Ply strains, laminate coordinates
  406. :t
  407. Procedure
  408. :i
  409. nplies
  410. #3
  411. :n
  412. pslc
  413. :c
  414. Ply stresses, laminate coordinates
  415. :t
  416. Procedure
  417. :i
  418. nplies
  419. #4
  420. :n
  421. tformse
  422. :c
  423. Transform stress/strain
  424. :t
  425. Procedure
  426. :i
  427. nplies
  428. #5
  429. :n
  430. lameps
  431. :c
  432. Laminate strains
  433. :t
  434. Procedure
  435. #1
  436. :n
  437. MECH_DOC
  438. :z
  439. "This is the Mechanical Loadings model of the TK!Composites package.
  440. "Copyright 1990 by Barry J. Berenberg.
  441. "See Documents model for distribution information.
  442.  
  443. "This model allows you to calculate stresses and strains in a
  444. "laminated plate by entering loadings as stress and moment resultants.
  445. "Outputs are laminate strains and curvatures, and ply strains and
  446. "stresses in both laminate and ply coordinates.
  447.  
  448. "To use this model, first load the laminate model, then "Merge..."
  449. "this model to it.  Set up the laminate in the laminate interactive 
  450. "table, enter the loadings on the variable sheet, and solve.  List
  451. "solving for families of laminates is not available yet, but will be
  452. "soon.
  453.  
  454. "Only laminate strains will appear on the variable sheet.  Ply stresses
  455. "and strains are shown in the tables selam and seply.  Laminated plate
  456. "theory assumes ply stresses and strains vary linearly through the
  457. "ply thickness.  The tables show stress and strain at the upper and
  458. "lower surface of each ply
  459. /E
  460. E
  461. #2
  462. :n
  463. pelc
  464. :z
  465. for i=1 to nplies
  466.    zup := 'h[i]
  467.    zdn := zup - 'tply[i]
  468.    'epl11up[i] := 'ek[1] - zup*'ek[4]
  469.    'epl11dn[i] := 'ek[1] - zdn*'ek[4]
  470.    'epl22up[i] := 'ek[2] - zup*'ek[5]
  471.    'epl22dn[i] := 'ek[2] - zdn*'ek[5]
  472.    'epl12up[i] := 'ek[3] - zup*'ek[6]
  473.    'epl12dn[i] := 'ek[3] - zdn*'ek[6]
  474. next i
  475. /E
  476. E
  477. #3
  478. :n
  479. pslc
  480. :z
  481. for i=1 to nplies
  482.    'spl11up[i] := 'Et1111[i]*'epl11up[i] + 'Et1122[i]*'epl22up[i] + 2*'Et1112[i]*'epl12up[i]
  483.    'spl11dn[i] := 'Et1111[i]*'epl11dn[i] + 'Et1122[i]*'epl22dn[i] + 2*'Et1112[i]*'epl12dn[i]
  484.    'spl22up[i] := 'Et1122[i]*'epl11up[i] + 'Et2222[i]*'epl22up[i] + 2*'Et2212[i]*'epl12up[i]
  485.    'spl22dn[i] := 'Et1122[i]*'epl11dn[i] + 'Et2222[i]*'epl22dn[i] + 2*'Et2212[i]*'epl12dn[i]
  486.    'spl12up[i] := 'Et1112[i]*'epl11up[i] + 'Et2212[i]*'epl22up[i] + 2*'Et1212[i]*'epl12up[i]
  487.    'spl12dn[i] := 'Et1112[i]*'epl11dn[i] + 'Et2212[i]*'epl22dn[i] + 2*'Et1212[i]*'epl12dn[i]
  488. next i
  489. /E
  490. E
  491. #4
  492. :n
  493. tformse
  494. :z
  495. for i=1 to nplies
  496.    theta := 'angle[i]
  497.    'spp11up[i] := 'spl11up[i]*cosd(theta)^2 + 'spl22up[i]*sind(theta)^2 + 2*'spl12up[i]*sind(theta)*cosd(theta)
  498.    'spp11dn[i] := 'spl11dn[i]*cosd(theta)^2 + 'spl22dn[i]*sind(theta)^2 + 2*'spl12dn[i]*sind(theta)*cosd(theta)
  499.    'spp22up[i] := 'spl11up[i]*sind(theta)^2 + 'spl22up[i]*cosd(theta)^2 - 2*'spl12up[i]*sind(theta)*cosd(theta)
  500.    'spp22dn[i] := 'spl11dn[i]*sind(theta)^2 + 'spl22dn[i]*cosd(theta)^2 - 2*'spl12dn[i]*sind(theta)*cosd(theta)
  501.    'spp12up[i] := ('spl22up[i]-'spl11up[i])*sind(theta)*cosd(theta) + 'spl12up[i]*(cosd(theta)^2-sind(theta)^2)
  502.    'spp12dn[i] := ('spl22dn[i]-'spl11dn[i])*sind(theta)*cosd(theta) + 'spl12dn[i]*(cosd(theta)^2-sind(theta)^2)
  503.    'epp11up[i] := 'epl11up[i]*cosd(theta)^2 + 'epl22up[i]*sind(theta)^2 + 2*'epl12up[i]*sind(theta)*cosd(theta)
  504.    'epp11dn[i] := 'epl11dn[i]*cosd(theta)^2 + 'epl22dn[i]*sind(theta)^2 + 2*'epl12dn[i]*sind(theta)*cosd(theta)
  505.    'epp22up[i] := 'epl11up[i]*sind(theta)^2 + 'epl22up[i]*cosd(theta)^2 - 2*'epl12up[i]*sind(theta)*cosd(theta)
  506.    'epp22dn[i] := 'epl11dn[i]*sind(theta)^2 + 'epl22dn[i]*cosd(theta)^2 - 2*'epl12dn[i]*sind(theta)*cosd(theta)
  507.    'epp12up[i] := ('epl22up[i]-'epl11up[i])*sind(theta)*cosd(theta) + 'epl12up[i]*(cosd(theta)^2-sind(theta)^2)
  508.    'epp12dn[i] := ('epl22dn[i]-'epl11dn[i])*sind(theta)*cosd(theta) + 'epl12dn[i]*(cosd(theta)^2-sind(theta)^2)
  509. next i
  510. /E
  511. E
  512. #5
  513. :n
  514. lameps
  515. :z
  516. "Multiplies ABDn matrix by NM matrix by dotting columns of
  517. "ABDn with NM.  Not true matrix multiplication--only works
  518. "because ABDn is symmetric.
  519.  
  520. for i=1 to 6
  521.    'ek[i] := dot('ABDn[i],'NM)
  522. next i
  523. "Curvatures should be -D*'NM
  524. for i=4 to 6
  525.    'ek[i] := -'ek[i]
  526. next i
  527. /E
  528. E
  529. =r
  530. #1
  531. :r
  532. place('NM,1)=N11
  533. #2
  534. :r
  535. place('NM,2)=N22
  536. #3
  537. :r
  538. place('NM,3)=2*N12
  539. #4
  540. :r
  541. place('NM,4)=M11
  542. #5
  543. :r
  544. place('NM,5)=M22
  545. #6
  546. :r
  547. place('NM,6)=2*M12
  548. #7
  549. :r
  550. call lameps()
  551. #8
  552. :r
  553. eps11='ek[1]
  554. #9
  555. :r
  556. eps22='ek[2]
  557. #10
  558. :r
  559. eps12='ek[3]
  560. #11
  561. :r
  562. kap11='ek[4]
  563. #12
  564. :r
  565. kap22='ek[5]
  566. #13
  567. :r
  568. kap12='ek[6]
  569. #14
  570. :r
  571. call pelc(nplies)
  572. #15
  573. :r
  574. call pslc(nplies)
  575. #16
  576. :r
  577. call tformse(nplies)
  578. =p
  579. =t
  580. #1
  581. :n
  582. selam
  583. :v
  584. Vertical
  585. :s
  586. Show
  587. :f
  588. 1
  589. :t
  590. "Ply stresses and strains, laminate coordinates"
  591. :u
  592. #1
  593. :l
  594. spl11up
  595. :w
  596. 10
  597. :h
  598. ""
  599. #2
  600. :l
  601. spl11dn
  602. :w
  603. 10
  604. :h
  605. ""
  606. #3
  607. :l
  608. spl22up
  609. :w
  610. 10
  611. :h
  612. ""
  613. #4
  614. :l
  615. spl22dn
  616. :w
  617. 10
  618. :h
  619. ""
  620. #5
  621. :l
  622. spl12up
  623. :w
  624. 10
  625. :h
  626. ""
  627. #6
  628. :l
  629. spl12dn
  630. :w
  631. 10
  632. :h
  633. ""
  634. #7
  635. :l
  636. epl11up
  637. :w
  638. 10
  639. :h
  640. ""
  641. #8
  642. :l
  643. epl11dn
  644. :w
  645. 10
  646. :h
  647. ""
  648. #9
  649. :l
  650. epl22up
  651. :w
  652. 10
  653. :h
  654. ""
  655. #10
  656. :l
  657. epl22dn
  658. :w
  659. 10
  660. :h
  661. ""
  662. #11
  663. :l
  664. epl12up
  665. :w
  666. 10
  667. :h
  668. ""
  669. #12
  670. :l
  671. epl12dn
  672. :w
  673. 10
  674. :h
  675. ""
  676. E
  677. #2
  678. :n
  679. seply
  680. :v
  681. Vertical
  682. :s
  683. Show
  684. :f
  685. 1
  686. :t
  687. "Ply stresses and strains, ply coordinates"
  688. :u
  689. #1
  690. :l
  691. spp11up
  692. :w
  693. 10
  694. :h
  695. ""
  696. #2
  697. :l
  698. spp11dn
  699. :w
  700. 10
  701. :h
  702. ""
  703. #3
  704. :l
  705. spp22up
  706. :w
  707. 10
  708. :h
  709. ""
  710. #4
  711. :l
  712. spp22dn
  713. :w
  714. 10
  715. :h
  716. ""
  717. #5
  718. :l
  719. spp12up
  720. :w
  721. 10
  722. :h
  723. ""
  724. #6
  725. :l
  726. spp12dn
  727. :w
  728. 10
  729. :h
  730. ""
  731. #7
  732. :l
  733. epp11up
  734. :w
  735. 10
  736. :h
  737. ""
  738. #8
  739. :l
  740. epp11dn
  741. :w
  742. 10
  743. :h
  744. ""
  745. #9
  746. :l
  747. epp22up
  748. :w
  749. 10
  750. :h
  751. ""
  752. #10
  753. :l
  754. epp22dn
  755. :w
  756. 10
  757. :h
  758. ""
  759. #11
  760. :l
  761. epp12up
  762. :w
  763. 10
  764. :h
  765. ""
  766. #12
  767. :l
  768. epp12dn
  769. :w
  770. 10
  771. :h
  772. ""
  773. E
  774. =n
  775. =g
  776. #1
  777. :c
  778. .0000009999999999999999999
  779. :i
  780. 10
  781. :v
  782. Yes
  783. :a
  784. Yes
  785. :r
  786. Yes
  787. :h
  788. No
  789. :t
  790. 1
  791. %Tv,f,10,4,0,1,0
  792.