论坛 › C03：ANSYS--结构软件（Mechanical\Stuctural\Rigid Dynamic › 发一个变厚度复合材料壳的例子

amani 发表于 2008-2-10 17:12:38

发一个变厚度复合材料壳的例子

/prep7
ET,1,SHELL91,32,1
KEYOPT,1,11,0,
KEYOPT,1,8,1
ET,2,SHELL91,30,1
KEYOPT,2,11,0,
KEYOPT,2,8,1
MP,EX,1,135000
MP,EY,1,8800
MP,EZ,1,8800
MP,PRXY,1,0.33
MP,PRYZ,1,0.33
MP,PRXZ,1,0.33
MP,GXY,1,4470
MP,GYZ,1,4470
MP,GXZ,1,4470
RECTNG,0,10,0,10,
RECTNG,10,20,0,10,
AGLUE,all
LESIZE,all, , ,20, , , , ,1
!*
ASEL,S,LOC,X,0,10
MSHAPE,0,2D
MSHKEY,0
TYPE, 1
AMESH,all
allsel,all
ASEL,S,LOC,X,10,20
MSHAPE,0,2D
MSHKEY,0
TYPE, 2
AMESH,all
allsel,all
CSWPLA,11,0,1,1,
esys,11
*GET,MXNODE,NODE,,NUM,MAX,
*GET,MXELE,ELEM,,NUM,MAX,
*DIM,THICK,,MXNODE
ESEL,S,TYPE,,1
NSLE,S
*DO,ELE,1,MXELE
*IF,ESEL(ELE),EQ,1,THEN
    R,ELE
    RMODIF,ELE,1,32,0, , , ,0
    RMODIF,ELE,13,1,45,0.125,0,0,0, !1
    RMODIF,ELE,19,1,135,0.125,0,0,0, !2
    RMODIF,ELE,25,1,90,0.125,0,0,0, !3
    RMODIF,ELE,31,1,135,0.125,0,0,0, !4
    RMODIF,ELE,37,1,45,0.125,0,0,0, !5
    RMODIF,ELE,43,1,90,0.125,0,0,0, !6
    RMODIF,ELE,49,1,135,0.125,0,0,0, !7
    RMODIF,ELE,55,1,45,0.125,0,0,0, !8
    RMODIF,ELE,61,1,0,0.125,0,0,0, !9
    RMODIF,ELE,67,1,45,0.125,0,0,0, !10
    RMODIF,ELE,73,1,135,0.125,0,0,0, !11
    RMODIF,ELE,79,1,90,0.125,0,0,0, !12
    RMODIF,ELE,85,1,135,0.125,0,0,0, !13
    RMODIF,ELE,91,1,45,0.125,0,0,0, !14
    RMODIF,ELE,97,1,0,0.125,0,0,0, !15
    RMODIF,ELE,103,1,45,0.125,0,0,0, !16
    RMODIF,ELE,109,1,135,0.125,0,0,0, !17
    RMODIF,ELE,115,1,90,0.125,0,0,0, !18
    RMODIF,ELE,121,1,45,0.125,0,0,0, !19
    RMODIF,ELE,127,1,135,0.125,0,0,0, !20
    RMODIF,ELE,133,1,90,0.125,0,0,0, !21
    RMODIF,ELE,139,1,135,0.125,0,0,0, !22
    RMODIF,ELE,145,1,45,0.125,0,0,0, !23
    RMODIF,ELE,151,1,0,0.125,0,0,0, !24
    RMODIF,ELE,157,1,45,0.125,0,0,0, !25
    RMODIF,ELE,163,1,135,0.125,0,0,0, !26
    RMODIF,ELE,169,1,90,0.125,0,0,0, !27
    RMODIF,ELE,175,1,45,0.125,0,0,0, !28
    RMODIF,ELE,181,1,135,0.125,0,0,0, !29
    RMODIF,ELE,187,1,90,0.125,0,0,0, !30
    RMODIF,ELE,193,1,135,0.125,0,0,0, !31
    RMODIF,ELE,199,1,45,0.125,0,0,0, !32
*ENDIF
*ENDDO
*DO,NODE,1,MXNODE
*IF,NSEL(NODE),EQ,1,THEN
    THICK(node) = 0.125-0.0125*NX(NODE)
*ENDIF
*ENDDO
RTHICK,THICK(1),105,106,107,108
RTHICK,THICK(1),111,112,113,114
*DO,NODE,1,MXNODE
    THICK(node) =0

*ENDDO ALLSEL,ALL
ESEL,S,TYPE,,2
NSLE,S
*DO,ELE,1,MXELE
*IF,ESEL(ELE),EQ,1,THEN
    R,ELE
    RMODIF,ELE,1,30,0, , , ,0
    RMODIF,ELE,13,1,45,0.125,0,0,0, !1
    RMODIF,ELE,19,1,135,0.125,0,0,0, !2
    RMODIF,ELE,25,1,90,0.125,0,0,0, !3
    RMODIF,ELE,31,1,135,0.125,0,0,0, !4
    RMODIF,ELE,37,1,45,0.125,0,0,0, !5
    RMODIF,ELE,43,1,90,0.125,0,0,0, !6
    RMODIF,ELE,49,1,135,0.125,0,0,0, !7
    RMODIF,ELE,55,1,45,0.125,0,0,0, !8
    RMODIF,ELE,61,1,0,0.125,0,0,0, !9
    RMODIF,ELE,67,1,45,0.125,0,0,0, !10
    RMODIF,ELE,73,1,135,0.125,0,0,0, !11
    RMODIF,ELE,79,1,90,0.125,0,0,0, !12
    RMODIF,ELE,85,1,135,0.125,0,0,0, !13
    RMODIF,ELE,91,1,45,0.125,0,0,0, !14
    RMODIF,ELE,97,1,0,0.125,0,0,0, !15
    RMODIF,ELE,103,1,45,0.125,0,0,0, !16
    RMODIF,ELE,109,1,135,0.125,0,0,0, !17
    RMODIF,ELE,115,1,90,0.125,0,0,0, !18
    RMODIF,ELE,121,1,45,0.125,0,0,0, !19
    RMODIF,ELE,127,1,135,0.125,0,0,0, !20
    RMODIF,ELE,133,1,90,0.125,0,0,0, !21
    RMODIF,ELE,139,1,135,0.125,0,0,0, !22
    RMODIF,ELE,145,1,45,0.125,0,0,0, !23
    RMODIF,ELE,151,1,0,0.125,0,0,0, !24
    RMODIF,ELE,157,1,45,0.125,0,0,0, !25
    RMODIF,ELE,163,1,135,0.125,0,0,0, !26
    RMODIF,ELE,169,1,90,0.125,0,0,0, !27
    RMODIF,ELE,175,1,45,0.125,0,0,0, !28
    RMODIF,ELE,181,1,135,0.125,0,0,0, !29
    RMODIF,ELE,187,1,90,0.125,0,0,0, !30
*ENDIF
*ENDDO
*DO,NODE,1,MXNODE
*IF,NSEL(NODE),EQ,1,THEN
    THICK(node) = 0.125-0.0125*(NX(NODE)-10)
*ENDIF
*ENDDO
RTHICK,THICK(1),99,100,101,102
RTHICK,THICK(1),105,106,107,108
*DO,NODE,1,MXNODE
THICK(node) =0
*ENDDO
ALLSEL,ALL
/ESHAPE,1.0
EPLOT

wood257 发表于 2008-2-25 09:22:07

不错，感谢搂主，不知道搂主有没有做过带倒角的复合材料结构件，网格是怎么划分的

wood257 发表于 2008-2-25 11:04:51

RTHICK,THICK(1),105,106,107,108
RTHICK,THICK(1),111,112,113,114？？？？？
什么意思，应该是15，16，17，18
表示节点在单元中的实常数位置，并且为什么两个rthick
15 TK(I) Layer thickness at node I for layer 1
16 TK(J) Layer thickness at node J for layer 1
17 TK(K) Layer thickness at node K for layer 1
18 TK(L) Layer thickness at node L for layer 1

rocketma 发表于 2008-2-26 09:32:40

学习一下

meky121 发表于 2008-10-19 12:46:30

一个薄壁圆筒，复合材料共三层，外层厚度函数变化；

gyangpen 发表于 2008-10-20 19:50:31

好，雪中送碳

真是高人

cxh@ 发表于 2009-3-15 10:15:15

做得太好了！

sh_lin30 发表于 2009-3-15 11:15:45

shell91/99在新版12.0ansys中已经不列入单元文档了
估计以后得人，不见的会知道这两种单元

12.0ansys采用281替代这样两种单元
问题，281的单元截面定义和实常数赋值，好像都不能完全替代91/99在构建变厚度复合材料壳体方面的功能，不知道ansys废止91/99的用意何在？
还是，281的建模功能没用被如我这样的人开发出来？

liuhongliang 发表于 2009-5-24 21:33:25

偶是刚学ansys的新手，先保存下楼主的命令流，希望以后能看懂

weijenfary 发表于 2009-10-13 23:24:32

楼主说的也太专业了吧，俺是新学者，什么都看不懂:'(

cugxiong 发表于 2009-11-2 22:29:30

资源共享共同提高

yanzifu 发表于 2010-3-11 16:35:21

请问楼主，为什么要用两个RTHICK

getgoing 发表于 2010-7-27 21:59:17

靠，不得不顶啊

zhouad 发表于 2010-11-29 12:03:04

顶一个先，很是不错的，实际例子

yaozhenhua 发表于 2011-8-12 23:10:18

顶一下，最近就是做这方面例子，感谢楼主分享

zpj551 发表于 2011-9-4 00:58:11

偶是刚学ansys的新手，先保存下楼主的命令流，希望以后能看懂

doonyli 发表于 2011-9-4 11:07:58

*DIM,THICK,TABLE,2,1,1,Z,,,11
*SET,THICK(1,0),0
*SET,THICK(2,0),80
*SET,THICK(1,1),0.2
*SET,THICK(2,1),0.6
SECTYPE,1,SHELL
SECOFFSET,TOP
SECDATA,1,1,0,3,1
SECFUNCTION,%THICK%
用function也可以 不过显示厚度后看起来不那么美观

6l-d6 发表于 2011-12-23 11:06:42

谢谢楼主！学习下！

sl_hzx 发表于 2012-4-2 16:39:33

好深奥，看不懂，可是还是顶一下

amare 发表于 2012-4-24 18:47:40

好东西！！！！！！！！1

查看完整版本: 发一个变厚度复合材料壳的例子