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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 【求助】关于DP材料的一个简单例子，计算不收敛的原因探讨
51/1返回列表
查看: 1127  |  回复: 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的回帖

coldwind042

金虫 (正式写手)
hehe，哥们也是岩土工程的阿，多交流阿！我以前用dp准则是采用42单元
赞一下 回复此楼
5楼2008-12-05 10:00:49
已阅   回复此楼   关注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的回帖

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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +7 Tide man 2026-08-10 8/400 2026-08-11 15:49 by 铁帽子农民
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +5 布布和一二 2026-08-10 6/300 2026-08-11 10:15 by 医学老男孩
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目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 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
信息提示
关闭
请填处理意见
关闭