24小时热门版块排行榜    

查看: 1124  |  回复: 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的回帖

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,严重超过了你的材料强度。
建议依据实际问题来设计算例,确定荷载及材料参数。欢迎有结果再拿出来讨论!
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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 确定了,国自然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
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目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
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +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
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见