24小时热门版块排行榜    

查看: 1122  |  回复: 7
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

[交流] 【求助】关于DP材料的一个简单例子,计算不收敛的原因探讨

/PREP7  
!*  
!*****************************************赋DP材料参数
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2600
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,3e9
MPDATA,PRXY,1,,0.32
TB,DP,1,,,  
TBMODIF,1,1,3e7
TBMODIF,1,2,20  
TBMODIF,1,3,
!**********************************************建立单元,平面应变分析**************  
ET,1,PLANE183   
!*  
KEYOPT,1,1,0
KEYOPT,1,3,2
KEYOPT,1,6,0
KEYOPT,1,10,0   
!****************************************建模并划分网格
RECTNG,0,50,0,50,   
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,50, , , , ,1
!*  
MSHAPE,0,2D
MSHKEY,0
!*  
MSHKEY,1
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH
!*****************************************施加荷载并求解  
/SOL
FLST,2,1,4,ORDE,1   
FITEM,2,1   
!*  
/GO
DL,P51X, ,UY,0  
FLST,2,2,4,ORDE,2   
FITEM,2,2   
FITEM,2,4   
!*  
/GO
DL,P51X, ,UX,0  
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO
!*  
SFL,P51X,PRES,9e8,  
acel,,10
SOLVE   

问题描述:一个正方形,左右两边及下边界约束住相应方向的位移,上边界施加9e8pa的均布荷载,进行求解。
只是为了学习DP准则,所以材料参数没有实际参照。

为什么总是不收敛?请高手指教!!!不胜感激!!

[ Last edited by woshilsh on 2008-12-4 at 08:27 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

西南交大隧道方向招硕、博士
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

楼上的也好,多联系哦
西南交大隧道方向招硕、博士
8楼2008-12-13 01:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

coldwind042

金虫 (正式写手)

★ ★
kuhailangyu(金币+1,VIP+0):thanks~~!
lei0736(金币+1,VIP+0):谢谢 加一个太少了
你这个是非线性问题,必须先要考虑问题的真实性,否则是很难收敛的。你先可以试算一下线弹性静力,看看问题的变形。你这个问题主应力已经达到0.4e9pa,严重超过了你的材料强度。
建议依据实际问题来设计算例,确定荷载及材料参数。欢迎有结果再拿出来讨论!
2楼2008-12-04 19:26:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

楼上说的很有道理,我施加的荷载的确太大,在线弹性材料计算中,变形已经完全没有意义。所以我将荷载去除,并不再对重力进行计算,而是简单的在上表面施加了5000帕的压力,但是仍然不能收敛,请指教。Y向变形只有0.128*10*-5。不胜感激!
命令流如下:
/PREP7  
!*  
!*****************************************赋DP材料参数
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2600
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,3e9
MPDATA,PRXY,1,,0.32
TB,DP,1,,,  
TBMODIF,1,1,3e7
TBMODIF,1,2,20  
TBMODIF,1,3,
!**********************************************建立单元,平面应变分析**************  
ET,1,PLANE183   
!*  
KEYOPT,1,1,0
KEYOPT,1,3,2
KEYOPT,1,6,0
KEYOPT,1,10,0   
!****************************************建模并划分网格
RECTNG,0,50,0,50,   
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,50, , , , ,1
!*  
MSHAPE,0,2D
MSHKEY,0
!*  
MSHKEY,1
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH
!*****************************************施加荷载并求解  
/SOL
FLST,2,1,4,ORDE,1   
FITEM,2,1   
!*  
/GO
DL,P51X, ,UY,0  
FLST,2,2,4,ORDE,2   
FITEM,2,2   
FITEM,2,4   
!*  
/GO
DL,P51X, ,UX,0  
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO
!*  
SFL,P51X,PRES,5e3,  
SOLVE
西南交大隧道方向招硕、博士
3楼2008-12-04 23:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

引用回帖:
Originally posted by coldwind042 at 2008-12-4 19:26:
你这个是非线性问题,必须先要考虑问题的真实性,否则是很难收敛的。你先可以试算一下线弹性静力,看看问题的变形。你这个问题主应力已经达到0.4e9pa,严重超过了你的材料强度。
建议依据实际问题来设计算例,确 ...

问题基本解决,经验证不是DP材料的问题,主要原因在于单元问题,谢谢指教!
西南交大隧道方向招硕、博士
4楼2008-12-05 06:54:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么网上很多人说本周 12号出结果 +4 瞬息宇宙 2026-08-10 4/200 2026-08-11 02:02 by 旁观2020
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 8/400 2026-08-10 22:07 by demonstreets
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 15/750 2026-08-10 19:29 by 要乐观耀哥
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见