资源描述：
《ansys邓肯-张材料模型》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、ANSYS邓肯-张材料模型楼主给的在ANSYS上实现邓肯－张模型的方法很有用，但其中还有几点需要修正的，这也是楼上的兄弟们有疑问的原因。我把楼主的代码运行了一下，然后对照作了修改，现在上传一下，有问题的兄弟可以仔细对照一下，在这里我对其中几个比较明显的问题说明一下：1.MP命令不能直接给单元加材料，这是对的。在这里，楼主遗漏了一下命令：MPCHG，具体见下面的修改过的代码。2.关于密度的问题。这些要在宏中定义，每修改一种材料(即调用一次邓肯－张子程序)就要修改一次材料的密度，其他有关材料的问题可以类推。3.关于施加重力的问题。要在调用宏后，在同一个循环
2、中重新定义一下重力。以下是我修改过的楼主的代码，希望对兄弟们有所帮助。!用APDL得到初步成果，贴于此供感兴趣的朋友参考，不当之处敬请指正，!欢迎加以完善。!基本思路：!邓肯-张模型的关键点是材料的弹性模量随大小主应力差!及小主应力(围压)的变化而变化，用APDL实现之的基本思路是：!给每个单元定义一个材料号，分级施加荷载，在每个荷载步结束时提取出各!单元的大小主应力，据此计算出下个荷载步的弹性模量Et,修改各单元之MP，!用于下一步计算。!以下是一个简单算例，copy出去可直接运行。!!!常规三轴试验模拟!**********************
3、************************************FINISH/CLEAR/TITLE,NumericalSimulationofthreeaxestestingofsoils/PREP7*dim,SUy,array,50!Settlementrecords*dim,MaxPs,array,120!Maxhistoryp1-p3*dim,MaxDs,array,120!MaxhistoryDs!*dim,EEt,array,50!Etofelememt!!!Duncan-ChangModel!!!Symbols：c-粘滞力，Fai
4、-内摩擦角，Sf-破坏强度(p1-p3)f，!Ds-应力水平，Pa-大气压,P3-围压!***********************************************************************CREATE,Duncan-Chang!CreatMacrofile*afun,deg!Unitofangle*set,Pa,1e5*set,P1,-ArrS3(i)!注意：岩土工程中应力为拉负压正*set,P3,-ArrS1(i)*if,P3,LT,0.1*Pa,thenP3=0.1*Pa !围压最小取值*endifSf
5、0=2*(c0*cos(Fai)+P3*sin(Fai))/(1-sin(Fai))!Mohr-Coulomb破坏强度(p1-p3)fDs=(P1-P3)/Sf0!应力水平，*if,Ds,GT,0.95,thenDs=0.95 !应力水平最大取值*endif!判断加卸荷，如果(P1-P3)小于历史最大值视为卸荷－再加荷过程*if,MaxPs(i),LT,P1-P3,thenEi=k0*Pa*(P3/Pa)**n0Et0=Ei*(1-Rf*Ds)**2 !加荷情况的切线模量MaxPs(i)=P1-P3!保存历史最大应力*elseif,MaxPs(i
6、),GE,P1-P3Et0=Kur*Pa*(P3/Pa)**n0 !卸荷模量*endifmp,ex,i,Et0!修改材料i的Etmp,nuxy,i,Mu0mp,dens,I,1800！重新定义第I种材料的密度*END!***********************单元类型***********************************et,1,42 !平面四节点单元KEYOPT,1,3,2 !平面应变!*********************以下定义材料初始模量*************************mp,ex,1,3.728e7
7、!砂土的弹性模量mp,nuxy,1,0.33mp,dens,1,1800!******************************************************************blc4,0,0,0.08,0.15 !建立几何模型/PNUM,AREA,1/REPLOTaesize,all,0.01!网格划分mat,1amesh,allnsel,s,loc,y,0 !边界条件d,all,Uy!底边界竖向约束nsel,s,loc,x,0d,all,Ux!左侧边界水平向约束nsel,all/replotfini/SOLUtime,
8、0.01!施加围压ACEL,0,9.8,0sfl,all,pres,2e5!200kPasol
