应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 【求助】关于DP材料的一个简单例子，计算不收敛的原因探讨
51/1返回列表
查看: 954  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
西南交大隧道方向招硕、博士
1楼 2008-12-04 01:06:22
已阅   回复此楼   关注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的回帖
查看全部 8 个回答

coldwind042

金虫 (正式写手)
★ ★
kuhailangyu(金币+1,VIP+0):thanks~~!
lei0736(金币+1,VIP+0):谢谢 加一个太少了
你这个是非线性问题，必须先要考虑问题的真实性，否则是很难收敛的。你先可以试算一下线弹性静力，看看问题的变形。你这个问题主应力已经达到0.4e9pa，严重超过了你的材料强度。
建议依据实际问题来设计算例，确定荷载及材料参数。欢迎有结果再拿出来讨论！
赞一下(20人) 回复此楼
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的回帖

coldwind042

金虫 (正式写手)
hehe，哥们也是岩土工程的阿，多交流阿！我以前用dp准则是采用42单元
赞一下 回复此楼
5楼2008-12-05 10:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭