小弟最近在搞一個複合材料的應力分析。 但分析出來的結果，感覺好像錯很大。 我在想是因為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