| 查看: 6751 | 回复: 4 | |||
[求助]
求职！ANSYS瞬态响应分析：完全法和模态叠加法的问题 已有1人参与
|
分别用完全法和模态叠加法做了一根悬臂梁的瞬态响应分析 ，但是结果不同，请帮忙看看，我的问题出在哪里 。谢谢~~http://muchongimg.xmcimg.com/data/b ... _1463561510_685.png 如图 ，载荷是在右端施加一个时间极短的脉冲载荷（一个周期的高频率正弦载荷） 建模部分命令流： FINISH /CLEAR /PREP7 /eshape,1 ET,1,BEAM4 !用BEAM4梁单元 R,1,0.0001,0.00000000083,0.0021,0.01,0.01 !定义1号实常数（BEAM4单元） MP,EX,1,2.1E11 !弹性模量 MP,PRXY,1,0.3 !泊松比 MP,DENS,1,7800 !密度 !******几何模型建立 K,1,0,0,0 K,2,0.5,0,0 K,3,1,0,0 LSTR,1,2 !1号线 LSTR,2,3 !2号线 LSEL,S,,,1,2 !选中1、2号线 LATT,1,1,1 !为已选中的线分配材料模型、实常数、单元类型 LESIZE,ALL,,,100 !将已选中的线划分为100份 LMESH,ALL !线单元划分 !******施加约束 D,1,,,,,,ALL,,,,,！约束左端为固支 采用完全法分析的命令流： /PREP7 *SET,N_CIRCLE,20!计算的周期数 *SET,NSTP_C,30!每周期的载荷步数 *SET,FRE,10000!载荷频率 *SET,CURR_T,0 *SET,DT,1.0/FRE/NSTP_C *SET,PI,3.141592653 /SOLU ANTYPE,4 !指定分析类型：动力学分析 TRNOPT,FULL !完全法 TIMINT,ON !时间积分效应 OUTRES,ALL,LAST !所有子步，所有记录写入DATABASE *DO,I,0,N_CIRCLE-1,1!外部循环：共计算N_CIRCLE个周期 *DO,J,1,NSTP_C,1!内部循环：每个周期NSTP_C个载荷步 CURR_T=(I*NSTP_C+J)*DT!计算当前载荷步的时间点 *IF,I,EQ,0,THEN ALLSEL,ALL NSEL,S,LOC,X,1 !选择加载点 F,ALL,FY,-100*SIN(2*PI*FRE*CURR_T)!施加当前时间点的载荷 ALLSEL,ALL *ENDIF ALLSEL,ALL TIME,CURR_T !设置当前载荷步的时间点 AUTOTS,1 !自动时间步长 KBC,0 !载荷步之间线性插值：递增方式 DELTIM,DT/10,DT/12,DT/8 !本载荷步之中指定时间步长 LSWRITE,I*NSTP_C+J !写载荷步 *ENDDO *ENDDO LSSOLVE,1,N_CIRCLE*NSTP_C!求解所有载荷步 FINISH 采用模态叠加法分析瞬态响应的命令流： !模态分析 !********************** /SOLU ANTYPE,2!指定为模态分析 MODOPT,LANB,10 !指定提取特征值数及所用方法 MXPAND,10,,,YES!指定扩展的模态阶数 SOLVE FINISH /PREP7 *SET,N_CIRCLE,20!计算的周期数 *SET,NSTP_C,30!每周期的载荷步数 *SET,FRE,10000!载荷频率 *SET,CURR_T,0 *SET,DT,1.0/FRE/NSTP_C *SET,PI,3.141592653 /SOLU ANTYPE,4 !指定分析类型：动力学分析 TRNOPT,MSUP,10 !模态叠加法 TIMINT,ON !时间积分效应 OUTRES,ALL,LAST !所有子步，所有记录写入DATABASE *DO,I,0,N_CIRCLE-1,1!外部循环：共计算N_CIRCLE个周期 *DO,J,1,NSTP_C,1!内部循环：每个周期NSTP_C个载荷步 CURR_T=(I*NSTP_C+J)*DT!计算当前载荷步的时间点 *IF,I,EQ,0,THEN ALLSEL,ALL NSEL,S,LOC,X,1 !选择加载点 F,ALL,FY,-100*SIN(2*PI*FRE*CURR_T)!施加当前时间点的载荷 ALLSEL,ALL *ENDIF ALLSEL,ALL TIME,CURR_T !设置当前载荷步的时间点 AUTOTS,1 !自动时间步长 KBC,0 !载荷步之间线性插值：递增方式 DELTIM,DT/10,DT/12,DT/8 !本载荷步之中指定时间步长 LSWRITE,I*NSTP_C+J !写载荷步 *ENDDO *ENDDO LSSOLVE,1,N_CIRCLE*NSTP_C!求解所有载荷步 FINISH !!!!模态扩展 /SOLU EXPASS,ON NUMEXP,all OUTRES,all,all SOLVE FINISH Main Menu-TimeHist Posrproc打开绘制 2号节点（梁中位置）的时程曲线： 完全法： http://muchongimg.xmcimg.com/data/b ... _1463562564_853.png 模态叠加法： http://muchongimg.xmcimg.com/data/b ... _1463562564_383.png 可以看出这差太多了。。。怎么回事能有大神帮忙看看么  QQ图片20160518164958.png  f.png  m.png |
回复此楼» 收录本帖的淘帖专辑推荐
 ANSYS命令流 |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有21人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有8人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
dongmiaol2楼2016-05-19 08:16:09
3楼2016-05-19 14:40:33
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
4楼2016-05-20 07:19:26
5楼2018-09-12 11:19:45
