小弟最近在搞一个复合材料的应力分析。 但分析出来的结果，感觉好像错很大。 我在想是因为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