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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ansys重启动分析如何改变单元材料属性
131/2返回列表12下一页
查看: 5693  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

匿名

用户注销 (初入文坛)
本帖仅楼主可见

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼2016-01-06 11:06:45
已阅   同方向广播   申请仿真EPI   回复此楼   编辑   查看我的主页
回帖支持 ( 显示支持度最高的前 50 名 )

shgao20

专家顾问 (职业作家)

【答案】应助回帖

如果在载荷步之间，你只是改变了少量单元的材料属性，那么它将只影响这部分单元及其附近单元的结果，对整个模型的结果影响应该不大。
赞一下(1人) 回复此楼
10楼2017-08-12 07:54:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hailx

新虫 (初入文坛)
请问楼主当时解决这个问题没有，使用mpchg命令感觉没有起到预想的效果，和一开始直接设定为改变后的材料号结果是一样的
赞一下 回复此楼
2楼2016-05-17 21:12:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

不要使用重启动，算完一个载荷步，不要离开 /solu 模块，直接在 /solu 中使用 MPCHG 命令修改材料号，然后继续计算就没有你说的问题了。
赞一下 回复此楼
3楼2016-05-18 07:45:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

炮台123

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by shgao20 at 2016-05-18 07:45:34
不要使用重启动，算完一个载荷步，不要离开 /solu 模块，直接在 /solu 中使用 MPCHG 命令修改材料号，然后继续计算就没有你说的问题了。

请问，只要不输入finish，就没有离开/solu模块吧？改变单元材料属性，能否由非线性材料，变为线性材料呢？
赞一下 回复此楼
4楼2017-08-09 17:30:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)
引用回帖:
4楼: Originally posted by 炮台123 at 2017-08-09 17:30:24
请问，只要不输入finish，就没有离开/solu模块吧？改变单元材料属性，能否由非线性材料，变为线性材料呢？...

不输入 finish，也不转换到其它模块 (如 /post1)，就不会离开求解模块。
如果在载荷步之间想将非线性材料改变为线性材料，有两个办法。都不要离开求解模块，然后：
  1  预先定义好需要的线性材料，在载荷步之间，使用 mpchg 命令将单元的非线性材料好改为相应的线性材料号。这种办法可以只修改部分定义的材料属性。
  2  使用 tbdele 命令删除非线性材料定义中的非线性部分，使其变为线性材料。这种办法会将所有该类非线性材料的定义都修改为线性材料。
赞一下 回复此楼
5楼2017-08-10 07:54:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

不输入finish，也没有执行如 /post1 之类的命令转换到其它模块，就不会离开求解模块。
  要改变单元材料属性，由非线性材料，变为线性材料，有两个办法：
    1  预先定义好所需的线性材料，然后在载荷步之间使用 mpchg 命令修改单元的材料编号；
    2  如果要将同一非线性材料的单元全部修改为线性材料，可以使用 tbdele 命令，删除原来材料的非线性部分，使其变为线性材料。
赞一下 回复此楼
6楼2017-08-10 08:00:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

炮台123

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by shgao20 at 2017-08-10 08:00:02
不输入finish，也没有执行如 /post1 之类的命令转换到其它模块，就不会离开求解模块。
  要改变单元材料属性，由非线性材料，变为线性材料，有两个办法：
    1  预先定义好所需的线性材料，然后在载荷步之间使用 ...

大神你好，我按照你提供的方法，在/solu模块中，利用mpchg命令，改变了单元的材料属性，但是在后处理，提取结果时，提示警告：Cumulative iteration 7 may have been solved using different model or boundary condition data than currently stored.POST1 results may be erroneous unless you RESUME from a Jobname.DB file for this substep. 导致无法提取结果，请问如何解决呢？

我做的是地铁隧道开挖模拟，考虑地层损失，在求解时，需要转换单元材料属性，命令流如下：
FINISH
/CLEAR
/filename,sdkw
/prep7  
ET,1,SOLID45            !用于模拟土体的三维实体单元
ET,2,shell181            !用于模拟衬砌的壳单元
r,1,0.35




!输入土层信息
MP,EX,1,3e7                 !输入弹性模量（围岩）
MP,PRXY,1,0.3                  !输入泊松比（围岩）
MP,DENS,1,1830           !输入密度（围岩）
TB,DP,1                  !采用DP准则进行弹塑性分析
TBdata,,22e4,12           !输入凝聚力（围岩）



!输入间隙材料信息
MP,EX,2,2e5                 !输入弹性模量（围岩）
MP,PRXY,2,0.3                  !输入泊松比（围岩）
MP,DENS,2,1830           !输入密度（围岩）
TB,DP,2                  !采用DP准则进行弹塑性分析
TBdata,,22e4,12           !输入凝聚力（围岩）



!输入衬砌信息
MP,EX,3,3.45e10         !输入弹性模量
MP,PRXY,3,0.2                  !输入泊松比
MP,DENS,3,2500           !输入密度



!输入隧道开挖土体
MP,EX,4,3e7                 !输入弹性模量（围岩）
MP,PRXY,4,0.3                  !输入泊松比（围岩）
MP,DENS,4,1830           !输入密度（围岩）
TB,DP,4                 !采用DP准则进行弹塑性分析
TBdata,,22e4,12           !输入凝聚力（围岩）


!输入土层损失信息
MP,EX,5,3e7                 !输入弹性模量（围岩）
MP,PRXY,5,0.3                  !输入泊松比（围岩）
MP,DENS,5,1830           !输入密度（围岩）
TB,DP,5                 !采用DP准则进行弹塑性分析
TBdata,,22e4,12           !输入凝聚力（围岩）







!!!!几何建模

csys,0
wpcsys,,0
BLOCK,0,100,0,50,0,50       !建立围岩

allsel
cyl4,50,35,3,,,,50
cyl4,50,35,3.05,,,,50
vptn,all

wpoff,,,25
vsbw,all
wpcsys,-1
wpoff,50,35
wprota,,,-90
vsbw,all
wpoff,,,6
vsbw,all
wpoff,,,-12
vsbw,all
wprota,,90
vsbw,all
wpoff,,,6
vsbw,all
wpoff,,,-12
vsbw,all


csys,0
wpcsys,,0
allsel
vsel,s,loc,y,33.5,36.5
vsel,r,loc,x,48.5,51.5
cm,dwtt,volu
vatt,4,,1


allsel
csys,0
wpcsys,,0
vsel,s,loc,y,32,38
vsel,r,loc,x,37,53
cmsel,u,dwtt
cm,dcss,volu
vatt,5,,1



csys,0
wpcsys,,0
local,12,1,50,35


csys,12
wpcsys,,12
asel,s,loc,x,3
cm,cq,area
aatt,3,1,2  



allsel
csys,0
wpcsys,,0
vsel,all
cmsel,u,dwtt
cmsel,u,dcss
vatt,1,,1


allsel
LSEL,S,LENGTH,,3
LESIZE,ALL,,,4,,,,,1
LSEL,S,LENGTH,,0.05
LESIZE,ALL,,,1,,,,,1
*SET,PI,ACOS(-1)
LSEL,S,LENGTH,,PI*3/2
LESIZE,ALL,,,8,,,,,1


allsel
vsel,all
esize,2
mshape,0
vsweep,all

allsel
mshkey,1
amesh,cq

!!!!边界条件
allsel
csys,0
wpcsys,,0
ACEL,,9.8
NSEL,S,LOC,Y,0
D,ALL,uy
nsel,s,loc,x,0
d,all,ux
nsel,s,loc,x,100
d,all,ux
nsel,s,loc,z,0
d,all,uz
nsel,s,loc,z,50
d,all,uz
allsel,all



   !nummrg,all,,,,low           !合并重复各元素并保留低编号号码
   !numcmp,all                 !压缩各元素编号号码
allsel,all
finish        




!!!!进入求解器
/solu  
antype,static            
deltim,0.1,0.1,0.2         !时间步设置，总长0.1，最小0.1，最大0.2
autots,on                         !使用自动时间步
pred,on                         !打开时间步长预测器
lnsrch,on                         !打开线性搜索
nlgeom,on                         !打开大位移效果
nropt,full                         !设定牛顿-拉普森选项
cnvtol,f,,0.05,2,0.5                   !设定力收敛条件
Save                       !保存数据库
esel,s,mat,,3                         !选择衬砌单元
ekill,all                  !杀死衬砌单元，即在隧道未修建前的自重应力场中不存在衬砌
esel,all                       !选择所有元素
esel,s,live                            !选择活单元，即所有围岩实体单元
nsle,s                                    !选择生单元上的节点
nsel,invert                         !反向选择，即选择了死单元上的节点
d,all,all                 !将死单元上的节点约束所有位移，使其不参与矩阵运算
nsel,all                   !选择所有节点
esel,all                   !选择所有单元
Save                         !保存数据库

iswrite,on                      !打开初应力文件生成开关
solve                                 !进行自重地应力场模拟计算，第一荷载步
Save

/post1
/solu

isfile,read,sdkw,ist,,2        !读入初应力文件，作为荷载施加
     !isfile,list                     !查看施加的初应力荷载
solve                           !进行地应力平衡计算（此时的荷载为重力荷载和初应力荷载）






!!!开挖土体
!!!以下进行土体的开挖操作，先选择每次开挖的围岩单元，然后将其赋予“死属性”

!第一个荷载开挖步
esel,s,mat,,4
    nsle,s
    nsel,r,loc,z,0,25
    esln,r,1   
    ekill,all


!以下进行支护结构的施加操作，先选择支护结构壳单元，然后将其赋予“生属性”
esel,s,mat,,3               
    nsle,s
    nsel,r,loc,z,0,25
    esln,r,1
    ealive,all  
    nsle,s
    ddele,all,all                 !将定义的约束条件删除     


esel,s,mat,,5
    nsle,s
    nsel,r,loc,z,0,25
    esln,r,1
    mpchg,2,all      !将材料号为5的单元变为材料号2



!选择生单元，即包括支护结构壳单元和未开挖部分围岩实体单元
    esel,all
    esel,s,live             
    nsle,s
!反向选择，并将死单元上的节点约束所有自由度
    nsel,invert          
    d,all,all
    nsel,all                              !选择所有节点
    esel,all                              !选择所有单元
    solve



!第二个荷载开挖步
esel,s,mat,,4
    nsle,s
    nsel,r,loc,z,25,50
    esln,r,1   
    ekill,all


!以下进行支护结构的施加操作，先选择支护结构壳单元，然后将其赋予“生属性”
esel,s,mat,,3               
    nsle,s
    nsel,r,loc,z,25,50
    esln,r,1
    ealive,all  
    nsle,s
    ddele,all,all                 !将定义的约束条件删除     


esel,s,mat,,5
    nsle,s
    nsel,r,loc,z,25,50
    esln,r,1
    mpchg,2,all



!选择生单元，即包括支护结构壳单元和未开挖部分围岩实体单元
    esel,all
    esel,s,live             
    nsle,s
!反向选择，并将死单元上的节点约束所有自由度
    nsel,invert          
    d,all,all
    nsel,all                              !选择所有节点
    esel,all                              !选择所有单元
    solve





/post1
/DEVICE,VECTOR,0
SET,1,LAST,1
path,dbcj,2
ppath,1,,50,50,0
ppath,2,,50,50,50
pdef,dbcj,u,y
prpath,dbcj               !提取隧道顶部地表竖向位移
赞一下 回复此楼
7楼2017-08-10 15:57:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

你遇到的警告信息，是因为在载荷步之间修改了单元的材料属性，因此提示你不同载荷步使用的计算模型是不一样的。这个警告信息应该不影响对结果的后处理。
  如果你忽略该警告信息后确实无法提取结果，可以在 post1 中先点击 result summary 看一下是否确实得到了结果。如果确实有结果，那么该警告信息是不会影响结果的提取的。
赞一下 回复此楼
8楼2017-08-11 07:56:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

炮台123

新虫 (初入文坛)
引用回帖:
8楼: Originally posted by shgao20 at 2017-08-11 07:56:09
你遇到的警告信息，是因为在载荷步之间修改了单元的材料属性，因此提示你不同载荷步使用的计算模型是不一样的。这个警告信息应该不影响对结果的后处理。
  如果你忽略该警告信息后确实无法提取结果，可以在 post1  ...

结果好像可以提取了，但是改变材料属性计算出来的结果，跟不改变材料属性计算的结果差别很小。。。。但是我每算一步，都会看一下单元的材料号，材料号确实发生了改变。。。。
赞一下 回复此楼
9楼2017-08-11 10:07:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xujiangshan 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭