24小时热门版块排行榜

返回列表

岳鑫鑫

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 37.3
帖子: 23
在线: 5.7小时
虫号: 4077555
注册: 2015-09-16
性别: GG
专业: 固体力学

[求助] ANSYS的Do循环十万次输出位移到数组里，从9万多次之后数组里面的位移值全是0

请问各位有没有遇到过这种情况，一个桁架结构，采用ANSYS的Do循环输出最大位移到一个数组里，循环十万次，从9万多次之后数组里面的位移值全是0了，这个怎么回事？谢谢解答！

命令流
! 桁架桥梁结构分析
! 11个输入的分析
finish
/clear

! 文件写入
*Create, dataread, mac
*DIM, truss, array, 100000, 11                   ! 创建样本数据存储的数组
*VREAD, truss(1,1), truss, TXT, , JIK, 11, 100000
(11f7.4)
*End
Dataread
*DIM, weiyi, array, 100000, 1                   ! 创建位移数据存储的数组

/prep7

*get, btim, active, 0, time, CPU

*do, i, 1, 100000, 1 ! 行循环

  ! 进行参数设置
  ! 设桥梁的总长度为L，每段水平桁架长度设为DL，桥梁高度设为H
  L = 32
  DL = L/8
  H = 5.5
  ! 设桥梁钢截面的面积为A，惯性矩为I
  A1 = truss(i, 2)*0.001
  I1 = 3.83E-6
  A2 = truss(i, 3)*0.001
  I2 = 1.87E-6
  A3 = truss(i, 4)*0.001
  I3 = 8.47E-6
  ! 将弹性模量和泊松比设为参数
  e_modu = truss(i, 1)*100000000000
  prxy_Poi = 0.3
  ! 荷载值
  P1 = -10000*truss(i, 5)
  P2 = -10000*truss(i, 6)
  P3 = -10000*truss(i, 7)
  P4 = -10000*truss(i, 8)
  P5 = -10000*truss(i, 9)
  P6 = -10000*truss(i, 10)
  P7 = -10000*truss(i, 11)
  ! 建模分析
  et, 1, link180
  r, 1, A1, I1, , , , ,
  r, 2, A2, I2, 0, 0, 0, 0
  r, 3, A3, I3, 0, 0, 0, 0
  ! 定义钢的弹性模量和泊松比
  mp, ex, 1, e_modu
  mp, prxy, 1, prxy_Poi
  mp, dens, 1, , 7800

  ! 定义关键点
  k, 1, 0, 0, ,
  k, 2, DL, 0, ,
  k, 3, 2*DL, 0, ,
  k, 4, 3*DL, 0, ,
  k, 5, 4*DL, 0, ,
  k, 6, 5*DL, 0, ,
  k, 7, 6*DL, 0, ,
  k, 8, 7*DL, 0, ,
  k, 9, 8*DL, 0, ,
  k, 10, DL, H, ,
  k, 11, 2*DL, H, ,
  k, 12, 3*DL, H, ,
  k, 13, 4*DL, H, ,
  k, 14, 5*DL, H, ,
  k, 15, 6*DL, H, ,
  k, 16, 7*DL, H, ,
  ! 通过几何点生成桥底梁的线
  L, 1, 2 $ L, 2, 3
  L, 3, 4 $ L, 4, 5
  L, 5, 6 $ L, 6, 7
  L, 7, 8 $ L, 8, 9
  ! 生成桥顶梁线和侧梁线
  L, 9, 16 $ L, 15, 16
  L, 14, 15 $ L, 13, 14
  L, 12, 13 $ L, 11, 12
  L, 10, 11 $ L, 1 , 10
  ! 生成桥身弦杆的线
  L, 2 , 10 $ L, 3 , 10
  L, 3 , 11 $ L, 4 , 11
  L, 4 , 12 $ L, 4 , 13
  L, 5 , 13 $ L, 6 , 13
  L, 6 , 14 $ L, 6 , 15
  L, 7 , 15 $ L, 7 , 16
  L, 8 , 16

  ! 选择桥顶梁和侧梁指定单元属性
  lsel, s, , , 9, 16, 1,
  latt, 1, 1, 1, , , ,
  ! 选择桥身弦杆指定单元属性
  lsel, s, , , 17, 29, 1,
  latt, 1, 2, 1, , , ,
  ! 选择桥底梁指定单元属性
  lsel, s, , , 1, 8, 1,
  latt, 1,3,1, , , ,

  ! 划分网格
  allsel, all
  lesize, all, , , 1, , , , , 1
  lmesh, all

  ! 施加位移约束、外力，进行求解
  /solu
  nsel, s, loc, x, 0
  d, all, , , , , , all, , , , ,
  allsel, all
  nsel, s, loc, x, L
  d, all, , , , , , , uy, , , ,
  allsel all
  ! 施加荷载
  fk, 2, fy, P1 $ fk, 3, fy, P2
  fk, 4, fy, P3 $ fk, 5, fy, P4
  fk, 6, fy, P5 $ fk, 7, fy, P6
  fk, 8, fy, P7

  /replot
  allsel, all
  solve

  /post1
  !  weiyi(i) = uy(5)
  *get, weiyi(i), node, 5, U, y ! 返回节点5的Y位置，并将结果作为weiyi存储。

  /prep7
  lclear, all
  ldele, all
  kdele, all

*enddo

*get, etim, active, 0, time, CPU
CPUTIME = etim - btim
*stat, CPUTIME                                  ! 显示总用时长

! 将位移数据导出
*create, datawrite, mac
*cfopen, weiyi, txt
*vwrite, weiyi(1,1)
  (f9.6)
*cfclose
*end
Datawrite@wuming524

回复此楼

» 猜你喜欢

1楼 2019-10-12 16:50:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ansys_apdl

禁虫 (小有名气)

本帖内容被屏蔽

2楼2019-11-23 08:59:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主岳鑫鑫的主题更新

返回列表