| 查看: 1042 | 回复: 1 | ||
[求助]
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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
基金不中，共勉
已经有4人回复
-
国自然面上复盘~欢迎讨论
已经有5人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
-
2026国自然函评费到账
已经有25人回复
ansys_apdl
禁虫 (小有名气)
|
本帖内容被屏蔽 |
2楼2019-11-23 08:59:29
投票:
