[問題] ANSYS Solid46建纏繞角度model的問題
小弟最近在搞一個複合材料的應力分析。
但分析出來的結果，感覺好像錯很大。
我在想是因為solid46這元素我還不太會用。
問題描述:
一個三維模型的圓筒只去模擬三十度，
而圓筒的材質 內壁=鋁 外壁=十層不同纏繞角度的金屬纖維(繞桶身綑綁)
內壁給壓力 40MPa，外壁給0.1MPa。
其他四個璧面給對稱
但模擬結果的趨勢卻是錯誤的，
爬了很多文，但還是不知道錯在哪邊?!
我覺得因該是纏繞角度，但我實在不知道怎麼搞，局部座標?!
煩請有使用SOLID46的高高高高手神神神人大大賜教!!!!
以下為完整的apdl，能否請板上大大。幫小弟我解答一下 問題是出在哪邊?!
! IL=內半徑
! AL=鋁厚度
! ONELAYER=1層碳纖維的厚度
! LAYERNU=碳纖維層數
! NL=碳纖維總厚度
! PVHL=桶身一半總長
! L1AN=第一層纏繞的角度
finish
/clear
IL=44
AL=1.8
ONELAYER=0.42
LAYERNU=10
NL=ONELAYER*LAYERNU
PVHL=80
!Inner layer to outer layer
L1AN=90
L2AN=-90
L3AN=45
L4AN=-45
L5AN=90
L6AN=-90
L7AN=45
L8AN=-45
L9AN=90
L10AN=-90
/prep7
!=======================鋁========================
et,1,solid95
MP,EX,1,70E3 !Al
MP,NUXY,1,0.3 !Al
TB,BKIN,1 !
TBDATA,,246,0 !yielding
!=======================================================
!=======================複合材料=========================
ET,2,solid46,,,,,2,4
! 各層纏繞角度
R,2,LAYERNU,1 !et=2 Layer=10 symmetric stacking
RMORE
RMORE,2, L1AN, ONELAYER, 2, L2AN, ONELAYER
RMORE,2, L3AN, ONELAYER, 2, L4AN, ONELAYER
RMORE,2, L5AN, ONELAYER, 2, L6AN, ONELAYER
RMORE,2, L7AN, ONELAYER, 2, L8AN, ONELAYER
RMORE,2, L9AN, ONELAYER, 2, L10AN, ONELAYER
! ET2 材質
MP,EX,2,154100
MP,EY,2,11410
MP,EZ,2,11410
MP,PRXY,2,0.33
MP,PRYZ,2,0.49
MP,PRXZ,2,0.49
MP,GXY,2,7092
MP,GYZ,2,3792
MP,GXZ,2,7092
!========================================================
k,1, PVHL, IL,0
k,2, -PVHL, IL,0
k,3, -PVHL, IL+AL,0
k,4, PVHL, IL+AL,0
k,5, PVHL, IL+AL+NL,0
k,6, -PVHL, IL+AL+NL,0
A,1,2,3,4
A,3,4,5,6
!==============中心軸==================
k,7, 0, 0,0
k,8,10, 0,0
L,7,8
!===========面拉成體 30度===================
APLOT
FLST,5,2,5,ORDE,2
FITEM,5,1
FITEM,5,-2
ASEL,S, , ,P51X
FLST,2,2,5,ORDE,2
FITEM,2,1
FITEM,2,-2
FLST,8,2,3
FITEM,8,7
FITEM,8,8
VROTAT,P51X, , , , , ,P51X, ,30, ,
!==================VGLUE======================================
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VGLUE,P51X
!==========================================
/VIEW,1,1,2,3
!========內壁面 a3 40MPa====================================
FLST,2,1,5,ORDE,1
FITEM,2,3
SFA,P51X,1,PRES,40
!========外壁面 a9 0.1MPa=================================
FLST,2,1,5,ORDE,1
FITEM,2,9
SFA,P51X,1,PRES,0.1
!========旁邊壁面對稱 a 11 7 8 6 1 2 4 10====================================
FLST,2,8,5,ORDE,7
FITEM,2,1
FITEM,2,-2
FITEM,2,4
FITEM,2,6
FITEM,2,-8
FITEM,2,10
FITEM,2,-11
DA,P51X,SYMM
!==========體積指定材料===================================
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 2, 2, 2, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 1, 2, 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
ALLSEL,ALL
!=============mesh====================================
SMRT,1
MSHAPE,1,3D
MSHKEY,0
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!==================================
ALLSEL,ALL
!=============================================
FINISH
/SOL
ALLSEL,ALL
/REPLOT
/STATUS,SOLU
SOLVE
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 58.114.87.156
推
11/16 16:04, , 1F
11/16 16:04, 1F
→
11/16 16:05, , 2F
11/16 16:05, 2F
→
11/16 16:06, , 3F
11/16 16:06, 3F
→
11/18 23:58, , 4F
11/18 23:58, 4F