24小时热门版块排行榜    

查看: 1161  |  回复: 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

金虫 (正式写手)

★ ★
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的回帖

coldwind042

金虫 (正式写手)

hehe,哥们也是岩土工程的阿,多交流阿!我以前用dp准则是采用42单元
5楼2008-12-05 10:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

好啊,学岩土的,不过用ansys用的不是很多,呵呵。
多交流,常联系!!

[ Last edited by sunxiao on 2008-12-5 at 23:08 ]
西南交大隧道方向招硕、博士
6楼2008-12-05 23:06:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luohui5151

铁虫 (小有名气)

我是岩土的山东的啊

你好你好啊
一辈子做好一件事情
7楼2008-12-12 23:53:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

楼上的也好,多联系哦
西南交大隧道方向招硕、博士
8楼2008-12-13 01:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sunxiao 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +5 雨打竹帘 2026-08-30 9/450 2026-08-30 22:23 by 余韵清
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
信息提示
请填处理意见