用pipe单元对管道建模中的几点疑问

笑天

小弟初学Ansys，用pipe单元对管道建模有几点疑问，请高手指点，大家都谈点看法，共同提高~~
1.  疑问1： 我对同一管道分别用壳单元和pipe单元建模进行模态分析，结果相差很大，应该是pipe单元分析有误，下面是命令流，请高人指点错误之处？？？
2.  疑问2：用管单元建模后，单元大小Ansys是怎样分配的？建模过程中怎么没有指定单元大小？？
    附：(1)管道形状如图示，参数为：弹性模量为1.139e10 Pa；泊松比取0.35；密度取 1179 kg/m3；内径为25mm；外径为28.5mm；热膨胀系数为2.5×10-5m/m·k。
   (2) 壳单元的模态分析分析结果是可靠的，结果为：
    壳单元分析结果            pipe单元结果为：
    SET   TIME/FREQ           SET   TIME/FREQ   
     1    9.9292               1 0.70668E-03                           
     2    10.865               2 0.97435E-03
                               3 0.10589E-02
     3    17.551               4 0.19878E-02
     4    26.698               5 0.30583E-02
     5    32.381               6 0.41937E-02
     6    36.819               7 0.54694E-02
     7    52.350               8 0.56583E-02
     8    54.917               9 0.70759E-02
     9    79.526              10 0.90870E-02
    10    89.280              11 0.97958E-02
                              12 0.13001E-01

笑天

以下为命令流：
-
fini
/clear  
/TITLE,guandao  
/UNITS,MKS           !MKS system (m, kg, s, °C).
/PREP7               !进入前处理器

MP,EX,1,1.139e10      
MP,PRXY,1,0.35
MP,ALPX,1,2.5e-5         !热膨胀系数
MP,DENS,1,1179


PUNIT,2   !设置为米制单位 (meter A+centimeter B+fraction of cm C/D)
PSPEC,1,,,0.0285,0.0035,    !管子外径28.5CM   壁厚：0.35CM
  
POPT,NC              !管道分析标准: ASME Section III NC, Class 2.
                    
PTEMP,20,20,         !Temperature = 20°

BRANCH,0,0,0,0       !第一根管起点
/PNUM,NODE,1         !显示节点及编号

RUN,0+9+5/10         !meter A+centimeter B+fraction of cm C/D---点2
RUN,1+1+5/10         !点3
RUN,0+10             !转弯1
RUN,0,0,1+31+5/10    !点5
RUN,0,0,0+19         !转弯2
RUN,-1+76            !点7
RUN,-0+30            !转弯3
RUN,0,0+40           !转弯4
RUN,-0+22            !点10
RUN,-0+22            !转弯5
RUN,0,-0+64          !点12
  
/VIEW,1,1,1,1  

BEND,3,4,0.025,2, ,1,
BEND,5,6,0.025,2, ,1,  
BEND,8,7,0.025,2, ,1,
BEND,9,8,0.1,2, ,1,
BEND,10,11,0.1,2, ,1,   

VALVE,10, ,0.695

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

ACEL,0,-9.8,0      !重力加速度
                  
D,1,All

D,2,UY,0   
D,2,UZ,0
D,2,ROTY,0  
D,2,ROTZ,0         !对应K2

D,3,UY,0   
D,3,UZ,0
D,3,ROTY,0  
D,3,ROTZ,0         !对应K3


D,5,UX,0   
D,5,UY,0
D,5,ROTX,0  
D,5,ROTY,0         !对应K5
                  
D,7,UY,0           
D,7,UZ,0           
D,7,ROTY,0         
D,7,ROTZ,0         !对应K7

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOLU
ANTYPE,MODAL
MODOPT,LANB,12
MXPAND,12,,,YES
SOLVE

ferriswanglei

见帮助文档 // Technology Demonstration Guide // 8. Nuclear Piping System Under Seismic Loading // 8.3. Modeling
用pipe单元做试试看

笑天

ferriswanglei 发表于 2011-10-19 10:10
见帮助文档 // Technology Demonstration Guide // 8. Nuclear Piping System Under Seismic Loading // 8. ...

我这个模态分析的命令流就是用pipe单元做的，关键是结果不对，但找不出错误~~~

panqy

圆周大小在SECDATA中设置，长度大小用mesh line设置。
至于分析结果，单元不同，划分就不同，肯定会有些差别。能提供你用shell做分析时的命令流么？

笑天

panqy 发表于 2011-10-19 11:54
圆周大小在SECDATA中设置，长度大小用mesh line设置。
至于分析结果，单元不同，划分就不同，肯定会有些差 ...

SECDATA不可用。见help中，SECDATA, VAL1, VAL2, VAL3, VAL4, VAL5, VAL6, VAL7, VAL8, VAL9, VAL10
Describes the geometry of a section.      Notes
The SECDATA command is divided into these section types: BEAMS, TAPER, SHELLS, PRETENSION, JOINTS, and REINFORCING. The beam types are collectively listed first, followed by the remaining section types.

管单元不能用吧~~  你们应该自己做下，别想当然，pipe单元直接生成的离散模型，不能画网格，应该也不能用你所说的mesh line设置~~    你们觉得呢~~
  

笑天

我找到了一个错误，用RUN,命令时，后面要定义单元数量，不定义的话，就默认为只有一个单元~~  按2cm一个单元，从新定义单元数量后结果提高了一个数量级，但还是不对~~  要再提高3个数量级才正确~~  期待高人指点！！

panqy

笑天 发表于 2011-10-19 17:44
SECDATA不可用。见help中，SECDATA, VAL1, VAL2, VAL3, VAL4, VAL5, VAL6, VAL7, VAL8, VAL9, VAL10
Desc ...

这么说没错，所以能不能用你也应该自己做下

royzhang

纠正一下，应是ansys帮助文档/Mechanical APDL Documentation/Technology Demonstration Guide/9. Nuclear Piping System Under Seismic Loading