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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 实体模型分析
191/2返回列表12下一页
查看: 1166  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

1988yyb

新虫 (初入文坛)

[求助] 实体模型分析

新学ansys，最近做了个模型，有以下假设：
1、本模型未考虑挡土墙结构与土的接触面；
2、挡墙开挖施工范围外的原始土层、墙前后填土层、墙基夯实填土设不同参数；
3、挡土墙结构变形缝间距12.60m，挡墙此处自由端、但土层为单向约束。
4、衬砌结构和土体均采用实体单元SOLID95；土、结构均为弹性。
出现以下问题：
1、选单元类型，当土体和挡墙都选择solid45，或solid95时，结果能够收敛，可是后处理查看竖直方向的位移时，与实际并不相符合。（未定义D-P）
2、土体和挡墙分别选择solid45和solid65时，它们都是8节点6面体单元，结果提示警告错误。(未定义D-P)
3、土的本构模型一般选用D-P模型，可是用TB，DP,2 $TBDATA,1,0,10,20，加上去后就不收敛了。
4、扶肋做网格划分时，出现过警告。
模型的边界条件是：看第四张图片和第五张图片
第一张图片：按理来说，我只加了重力荷载，挡墙的填土应该有沉降才对呀，咋全是往上呢？
5、各位大神们，我的金币就这么多了，请各位帮个忙。

Y轴方向的位移图



模型介绍图



挡墙



整个模型介绍图



扶壁式挡墙



网格划分警告
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-03-24 21:47:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1988yyb

新虫 (初入文坛)
大神们，帮帮忙呀。。。最少“吱一声”呀
赞一下 回复此楼
2楼2013-03-25 20:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小七-vip

新虫 (小有名气)
★
xiegangmai: 金币+1, 谢谢参与！ 2013-03-26 16:54:09
网格加密试试
回复此楼
3楼2013-03-26 15:08:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
1988yyb: 金币+12, ★★★很有帮助 2013-03-27 11:35:58
xiegangmai: 金币+2, 谢谢建议！ 2013-03-27 21:49:57
你用的三种 SOLID 单元，两种是 8 节点六面体，一种是 20 节点六面体，它们的界面上是不协调的。最好采用有相同节点数的实体单元。
赞一下 回复此楼
4楼2013-03-27 09:48:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1988yyb

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by htbbzzg at 2013-03-27 09:48:11
你用的三种 SOLID 单元，两种是 8 节点六面体，一种是 20 节点六面体，它们的界面上是不协调的。最好采用有相同节点数的实体单元。

老师，您也在这呀，已经按建议改了，就是左侧与实际情况有点不相符合，土体受重力作用应该有沉降呀，可现在整体向上位移，至少有个分界面吧，还有个就是土的本构模型加不上去，
赞一下 回复此楼
5楼2013-03-27 11:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1988yyb

新虫 (初入文坛)
1988yyb: 回帖置顶 2013-03-27 14:15:57
引用回帖:
3楼: Originally posted by 小七-vip at 2013-03-26 15:08:03
网格加密试试

已经加密了。。。还是有问题的。。。
赞一下 回复此楼
6楼2013-03-27 11:36:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小七-vip

新虫 (小有名气)
★
xiegangmai: 金币+1, 谢谢参与！ 2013-03-27 21:50:03
本构方程有问题吧
赞一下 回复此楼
7楼2013-03-27 17:01:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1988yyb

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by 小七-vip at 2013-03-27 17:01:59
本构方程有问题吧

我也觉得有可能，可是用的本构方程，就是D-P 模型，ansys里本身就具有的呀，不是自己定义的。。。只需要按要求输入参数就可以了。。。晚些时候，我把命令流附上，望您能指导一下。。。
赞一下 回复此楼
8楼2013-03-27 18:04:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)
忘了一点，在 ANSYS 中，重力是通过 重力加速度来施加的，而重力与重力加速度方向相反。估计你定义的重力加速度是向下的，因而重力是向上的，造成位移也向上。如是这种情况，将重力加速度方向颠倒一下即可。
赞一下 回复此楼
9楼2013-03-28 08:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1988yyb

新虫 (初入文坛)
引用回帖:
9楼: Originally posted by htbbzzg at 2013-03-28 08:25:04
忘了一点，在 ANSYS 中，重力是通过 重力加速度来施加的，而重力与重力加速度方向相反。估计你定义的重力加速度是向下的，因而重力是向上的，造成位移也向上。如是这种情况，将重力加速度方向颠倒一下即可。

重力加速度，我是定义向上的，这个没有问题，现在还有个问题是，加上了本构模型d-p，就不收敛了。现在换种方法建模型。按老师，您说的:
finish
/clear
/filname,truss,1
/title,the anlysis of structual
/replot
/prep7              !进入前处理
et,1,solid65         !钢筋混凝土单元  8节点6面体单元
et,2,solid45        !8节点6面体单元   土体单元
!et,1,solid95        !20节点6面体单元  土体单元
!et,2,solid186       !20节点6面体单元   混凝土单元

mp,ex,1,2.6e10        !钢筋混凝土材料
mp,prxy,1,0.167
mp,dens,1,2500

      
mp,ex,2,7e7            !挡墙前、后填土
mp,prxy,2,0.32
mp,dens,2,1800
!mp,mu,3,0.4
!TB,DP,2
!TBDATA,2,20,15
!r,2

mp,ex,3,9e7           !挡墙基础回填土（夯实）
mp,prxy,3,0.28
mp,dens,3,2000
!mp,mu,4,0.4
!r,3
!TB,DP,3
!TBDATA,3,20,15




MP,EX,4,4e7          !挡墙开挖范围外的原始土层      
mp,prxy,4,0.32
mp,dens,4,1900
!mp,mu,5,0.5
!TB,DP,4
!TBDATA,4,0,25,25

!创建关键点、线、面

k,4                     !坐标原点（墙角内边缘点）   
k,3,-0.35
k,2,0,7.4
k,1,-0.35,7.4
k,5,-1.05
k,6,-1.05,-0.6
k,8,3.2
k,7,3.2,-0.6
k,9,-3.05,-0.6
k,10,-3.05,-2.1
k,11,5.2,-2.1
k,12,5.2,-0.6
k,13,25,-2.1
k,14,-10,-2.1
k,15,-10,0.9
k,17,25,7.4
k,16,-0.35,0.9
k,18,-3.05,0.9
k,19,3.2,7.4
k,21,5.2,7.4
k,20,-1.05,0.9
k,22,-10,-0.6
k,23,-3.05,0
k,24,0,0.9
k,25,-1.05,-2.1
k,26,5.2
k,27,25
k,28,25,-0.6
k,29,3.2,-2.1
k,30,-10
k,31,-0.35,-0.6
k,32,0,-0.6
k,33,-0.35,-2.1
k,34,0,-2.1
k,35,1.6,3.7,0


l,1,2      
l,1,16
l,2,24
l,24,4
l,4,32
l,4,8
l,8,7
l,7,12
l,7,32
l,32,31
l,31,6
l,6,5
l,5,3
l,3,16
l,3,4
l,3,31
l,16,20
l,20,5
l,20,18
l,18,23
l,23,9
l,9,10
l,18,15
l,15,30
l,30,22
l,22,14
l,14,10
l,10,25
l,6,9
l,25,33
l,33,34
l,34,29
l,29,11
l,11,12
l,12,26
l,11,13
l,13,28
l,28,27
l,27,17
l,17,21
l,21,26
l,21,19
l,19,8
l,19,2
l,2,35
l,35,8
l,8,26
l,26,27
l,12,28
l,7,29
l,32,34
l,31,33
l,6,25
l,5,23
l,23,30
l,9,22

l,16,24
l,24,35
l,35,19


al,1,2,3,57     
al,57,14,15,4
al,15,16,10,5
al,5,6,7,9
al,16,13,12,11
al,11,52,30,53
al,10,52,31,51
al,9,51,32,50
al,8,34,33,50
al,13,14,17,18
al,18,19,20,54
al,12,54,21,29
al,22,29,53,28
al,44,45,59
al,43,59,46
al,41,42,43,47
al,7,8,35,47

al,20,23,24,55
al,21,55,25,56
al,22,56,26,27

al,39,40,41,48
al,38,48,35,49
al,37,49,34,36
al,3,45,58
al,4,6,46,58


!形成3D几何模型、划分体单元
allsel
vext,24,25,1,,,3         !扶壁至变形缝段回填土（3.0m、MAT2）
vatt,2,,2
VPLOT


ALLSEL
TYPE,1
MAT,1
vext,26,30,4,,,0.3       !扶壁（0.3m、MAT1）
vatt,1,,1
VPLOT

ALLSEL
vext,34,38,4,,,6        !扶壁至扶壁段回填土（6.0m、MAT2）
vatt,2,,2
vplot

ALLSEL
vext,42,46,4,,,0.3        !扶壁（0.3m、MAT1）
vatt,1,,1
vplot


ALLSEL
vext,50,54,4,,,3         !扶壁至变形缝段回填土（3.0m、MAT2）
vatt,2,,2
VPLOT

ALLSEL
vext,1,5,1,,,12.6           !墙面板、底板全长（12.6m、MAT1）
vatt,1,,1
VPLOT

ALLSEL
ASEL,s,,,10,13,1
asel,a,,,14,17,1
CM,AREASOIL1,AREA
VEXT,AREASOIL1,AREASOIL1,1,,,12.6           !墙前、后填土（12.6m、MAT2）
vatt,2,,2
VPLOT


ALLSEL
vext,6,9,1,,,12.6         !墙基础夯实填土（12.6m、MAT3）
vatt,3,,2
VPLOT

ALLSEL
vext,18,23,1,,,12.6         !墙开挖周围天然土层（12.6m、MAT4）
vatt,4,,2
VPLOT
ALLSEL


allsel                         !网格划分
ESIZE,0.5
MSHAPE,0,3D
MSHKEY,1
VMESH,ALL
ALLSEL
   
这里就出现个问题，各个体的单元类型和材料属性居然是一样的，是不是之前赋予时，就出现错误。。。
赞一下 回复此楼
10楼2013-03-28 10:24:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1988yyb 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭