版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

应助: 0 (幼儿园)
贵宾: 6.013
金币: 16317
红花: 10
帖子: 2050
在线: 82.4小时
虫号: 45802
注册: 2004-05-15
性别: GG
专业: 岩土工程
管辖: 仿真模拟

[交流] 【求助】关于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

荣誉版主 (著名写手)

小木虫钓鱼协会会长

应助: 0 (幼儿园)
贵宾: 6.013
金币: 16317
红花: 10
帖子: 2050
在线: 82.4小时
虫号: 45802
注册: 2004-05-15
性别: GG
专业: 岩土工程
管辖: 仿真模拟

楼上的也好，多联系哦

赞一下

回复此楼

西南交大隧道方向招硕、博士

8楼2008-12-13 01:27:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 8 个回答

coldwind042

金虫 (正式写手)

应助: 0 (幼儿园)
贵宾: 0.03
金币: 543.5
红花: 3
帖子: 630
在线: 175.2小时
虫号: 154293
注册: 2006-01-03
性别: GG
专业: 水工结构和材料及施工

★ ★
kuhailangyu(金币+1,VIP+0):thanks~~!
lei0736(金币+1,VIP+0):谢谢加一个太少了

你这个是非线性问题，必须先要考虑问题的真实性，否则是很难收敛的。你先可以试算一下线弹性静力，看看问题的变形。你这个问题主应力已经达到0.4e9pa，严重超过了你的材料强度。
建议依据实际问题来设计算例，确定荷载及材料参数。欢迎有结果再拿出来讨论！

赞一下(20人)

回复此楼

2楼2008-12-04 19:26:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

应助: 0 (幼儿园)
贵宾: 6.013
金币: 16317
红花: 10
帖子: 2050
在线: 82.4小时
虫号: 45802
注册: 2004-05-15
性别: GG
专业: 岩土工程
管辖: 仿真模拟

楼上说的很有道理，我施加的荷载的确太大，在线弹性材料计算中，变形已经完全没有意义。所以我将荷载去除，并不再对重力进行计算，而是简单的在上表面施加了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

荣誉版主 (著名写手)

小木虫钓鱼协会会长

应助: 0 (幼儿园)
贵宾: 6.013
金币: 16317
红花: 10
帖子: 2050
在线: 82.4小时
虫号: 45802
注册: 2004-05-15
性别: GG
专业: 岩土工程
管辖: 仿真模拟

引用回帖:

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 个回答