【原创】有预应力的盘压杆结构模态分析

liuyiliang823

我想做这个盘压杆机构的带预应力的模态分析，该结构由横向的三角架和纵梁还有斜拉杆组成，三角架和纵梁之间耦合了在重合节点处的沿x，y，z三个方向的平移自由度，和一个z方向的转动自由度，这是我第一次写命令，不知道什么地方出了错，麻烦各位高手帮我看看，谢谢，命令流如下：
/prep7
ET,1,44      !定义三角架以及纵梁属性
MP,EX,1,5.7e10
MP,EY,1,1.24e10
MP,EZ,1,1.24e10
MP,PRXY,1,0.208
MP,PRYZ,1,0.058
MP,PRXZ,1,0.208
MP,GXY,1,6.19e9
MP,GYZ,1,6.19e9
MP,GXZ,1,6.19e9
MP,DENS,1,2000
SECTYPE,1,BEAM,RECT,,
SECDATA,0.003,0.003,4,4
  
ET,2,10      !定义斜拉杆属性
R,2,25E-6,0.01!定义斜拉杆的横截面积以及盘压杆上的初应力
MP,EX,2,1e10
MP,DENS,2,100
  
ET,3,21      !定义质量快属性，
R,3,1/3,1/3,1/3,,,
  
!建立盘压杆结构模型
/VIEW,1,,
CSYS,1
*AFUN, DEG
*DO,I,1,71
N,(I-1)*3+1,0.125,90,(I-1)*0.145
N,(I-1)*3+2,0.125,210,(I-1)*0.145
N,(I-1)*3+3,0.125,330,(I-1)*0.145
*ENDDO
TYPE,1
MAT,1
SECNUM,1
  
*DO,I,1,71
E,(I-1)*3+1,(I-1)*3+2
E,(I-1)*3+2,(I-1)*3+3
E,(I-1)*3+3,(I-1)*3+1
*ENDDO
  
TYPE,2
MAT,2
R,2
*DO,I,1,70
E,(I-1)*3+1,(I-1)*3+5
E,(I-1)*3+2,(I-1)*3+6
E,(I-1)*3+3,(I-1)*3+4
*ENDDO
  
K,1,0.125,90,0
K,2,0.125,210,0
K,3,0.125,330,0
K,4,0.125,90,70*0.145
K,5,0.125,210,70*0.145
K,6,0.125,330,70*0.145
LSTR,1,4
LSTR,2,5
LSTR,3,6
  
TYPE,1
MAT,1
SECNUM,1
LESIZE,ALL,,,700,,,,
LMESH,ALL
  
TYPE,3
R,3
E,211
E,212
E,213
  
CPINTF,UX
CPINTF,UY
CPINTF,UZ
CPINTF,ROTZ
  
/SOLU      !先做静态分析，载荷就是由于斜拉杆产生的力
D,1,ALL    !约束盘压杆底端三角架的三个节点的全部自由度
D,2,ALL
D,3,ALL
  
SOLVE
/POST1
plnsol,s,eqv,0,1
  
还有几个问题向大家请教，
1，在该模型中，三角架我用的是beam44单元，在每个节点有六个自由度，斜拉杆用的是link10单元，每个节点有3个自由度，是不是我把斜拉杆和三角架通过共用节点直接相连接就可以让斜拉杆相对于三角架自由转动？希望高手能指点一下，
2，我如果想把44号单元用189号单元代替，我应该怎样用命令设置keypoint
3，我调试这个命令流时出现错误，不知道为什么，
  
希望大家能帮帮我，谢谢你们！！！

foransys

这位老兄用画笔写字的功夫堪称一流，呵呵！

liuyiliang823

：）不好意思，是不是我没有说明白自己得问题啊？我又换了纵梁以及盘压杆得单元类型。
具体请见命令流文件！
/prep7
ET,1,189      !定义三角架以及纵梁属性
MP,EX,1,5.7e10
MP,EY,1,1.24e10
MP,EZ,1,1.24e10
MP,PRXY,1,0.208
MP,PRYZ,1,0.058
MP,PRXZ,1,0.208
MP,GXY,1,6.19e9
MP,GYZ,1,6.19e9
MP,GXZ,1,6.19e9
MP,DENS,1,2000
SECTYPE,1,BEAM,RECT,,
SECDATA,0.003,0.003,4,4
R,1,,
  
ET,2,8      !定义斜拉杆属性
R,2,25E-6,0.0001
MP,EX,2,1e10
MP,DENS,2,100
  
ET,3,21      !定义质量快属性
R,3,1/3,1/3,1/3,,,
  
/VIEW,1,1,1,1
CSYS,1
*AFUN,DEG
K,1,0.125,90,0
K,2,0.125,210,0
K,3,0.125,330,0
K,1000,0,0,5
LSTR,1,2
LSTR,2,3
LSTR,3,1
  
LATT,1,1,1,,1000,,1
LESIZE,ALL,,,4
LMESH,ALL
  
LGEN,71,ALL,,,0,0,0.145,,0
  
TYPE,2
MAT,2
REAL,2
*DO,I,1,70
E,(I-1)*36+1,(I-1)*36+38
E,(I-1)*36+2,(I-1)*36+50
E,(I-1)*36+14,(I-1)*36+37
*ENDDO
  
*GET,NKP,KP,,NUM,MAX
K,NKP+1,0.125,90,0
K,NKP+2,0.125,210,0
K,NKP+3,0.125,330,0
K,NKP+4,0.125,90,0.145*70
K,NKP+5,0.125,210,0.145*70
K,NKP+6,0.125,330,0.145*70
LSTR,NKP+1,NKP+4
LSTR,NKP+2,NKP+5
LSTR,NKP+3,NKP+6
  
LSEL,S,LINE,,214,216,1
LATT,1,1,1,,1000,,1
LESIZE,ALL,,,1400
LMESH,ALL
  
TYPE,3
REAL,3
LSEL,S,LOC,Z,0.145*70
NSEL,R,LOC,X,0.125
E,2521
E,2522
E,2534
  
ALLSEL
NROTAT,ALL
CPINTF,UX,
CPINTF,UY,
CPINTF,UZ,
CPINTF,ROTY,
CPINTF,ROTZ,
  
/SOLU
DK,1,ALL
DK,2,ALL
DK,3,ALL
SOLVE
  
有个问题不大明白，希望高手能帮帮我啊！！！！！！！！！！！！！！
就是在斜拉杆上加了预应力（拉力）
我觉得ansys求解得盘压杆得变形和实际不相符合！不知道为什么啊？

liuyiliang823

求解结果如图所示！我不明白为什么在斜拉杆上加上预应力（拉力）之后，为什么盘压杆得顶端出现了扩张啊，我觉得应该从低端到顶端，盘压杆在截面方向，即在三角架平面内，是逐渐缩小得啊，各位大侠帮我看看这个问题到底怎么回事，谢谢你们！！！！！！！！！！！！！

liuyiliang823

还有一个问题，就是从应力云图结果上看，无论从求解的节点平动自由度，或者转动自由度看，节点得求解结果都很小，，为什么从变形图上还有动画显示看到产生得变形那么大啊？？

harbinliying

可能是显示结果的比例放大了，可以设置为按真实比例来显示就可以了