找回密码
 注册
Simdroid-非首页
查看: 256|回复: 7

[新手疑问区] 跪求:基于邓肯-张模型的尾矿坝的 APDL二次开发源程序范例!

[复制链接]
发表于 2007-7-3 13:23:01 | 显示全部楼层 |阅读模式 来自 广西桂林
正在写论文,可是没有相关的参考资料,自己试了几次,可是全都不能收敛,请有相关成功经验的朋友们能割爱,千千万万!
 楼主| 发表于 2007-7-11 17:23:53 | 显示全部楼层 来自 广西桂林

我的这个命令流能否解决坝的问题

Simdroid开发平台
我在做ansys二次开发中邓肯张模型的应用,现碰到问题,模型的建立考虑有限单元法,逐个单元赋予弹性模量和泊松比,然后求解。问题就在于怎么给每个单元赋予弹性模量和泊松比,土体单元就一种42号单元,而材料却随单元数目的增加而增加,我是在求解了地应力后修改每个单元的弹性模量和泊松比的,修改后在求解,大家可以看看或帮忙改进一下,也给我提点建议,因为模型不够通用,不能象ansys中的的D-P模型那样通用,也希望大家一起来讨论,怎么建立一个通用的D-C模型,顺便说下,本模型在求解过程中还存在问题,我没找出来,还有很多问题要解决,希望有高人继续补充,开发出可以通用的模型。谢谢
FINISHokok.org
/CLEARokok.org
/FILNAME,PIT
/TITLE,PIT EXCAVATIONokok.org
/PREP7okok.org
*dim,Str_max,array,180
*dim,S_max,Array,180
*dim,Settlement,array,5
!*dim,M_et,array,5
!*dim,M_u,array,5

*create,Duncan-Chang
   *afun,degokok.org
   *set,Pa,1e5
   *set,P1,-ArrS3(i)
   *set,P3,-ArrS1(i)
        *IF,P3,LT,0.1*Pa,then
            P3=0.1*Pa
        *ENDIF
   Str=2*(c*cos(Fai)+P3*sin(Fai))/(1-sin(Fai))okok.org
   S=(P1-P3)/Str
        *if,S,GT,0.95,then
            S=0.95
        *endif

        *if,Str_max(i),gt,P1-P3,and,S_max(i),gt,S,then
            Et=Kur*Pa*(P3/Pa)**Nur
        *elseif,Str_max(i),gt,P1-P3,and,S_max(i),le,S,then
            Ei=k*Pa*(P3/Pa)**n
            Et=Ei*(1-Rf*S)**2
            S_max(i)=Sokok.org
        *elseif,Str_max(i),le,P1-P3,and,S_max(i),gt,S,then
            Ei=k*Pa*(P3/Pa)**n
            Et=Ei*(1-Rf*S)**2
            Str_max(i)=P1-P3okok.org
        *elseif,Str_max(i),le,P1-P3,and,S_max(i),le,S,then
            Ei=k*Pa*(P3/Pa)**n
            Et=Ei*(1-Rf*S)**2
            Str_max(i)=P1-P3okok.org
            S_max(i)=Sokok.org
        *endifokok.org
   !----------------------------------------------------------------------
   Ei=k*Pa*(P3/Pa)**n
   A=(p1-p3)*d/(Ei*(1-Rf*S))
   Mu=(G-F*log10(P3/Pa))/(1-A)**2
        *if,Mu,GE,0.49,then
            Mu=0.49
        *endif
   mp,ex,i,Etokok.org
   mp,nuxy,i,Mu
   !mpchg,i,
*end

*SET,T1,0.5
*SET,R1,25E-3
*SET,E1,4E7okok.org
*SET,NU1,0.32okok.org
*SET,DENS1,1800
*SET,E2,1.85E10
*SET,NU2,0.2
*SET,E3,2.1E11okok.org
*SET,NU3,0.3
okok.org
ET,1,42
KEYOPT,1,3,2
ET,2,3okok.org
KEYOPT,2,6,1
ET,3,1okok.org

*DO,I,1,3,1
MP,EX,I,E1okok.org
MP,PRXY,I,NU1okok.org
MP,DENS,I,DENS1
!TB,DP,I
!TBDATA,1,C_R,FI_Rokok.org
*ENDDO
      
MP,EX,4,E2okok.org
MP,PRXY,4,NU2okok.org
okok.org
*DO,I,5,6,1okok.org
MP,EX,I,E3okok.org
MP,PRXY,I,NU3okok.org
*ENDDO
   
R,1,T1,T1*T1*T1/12,T1okok.org
R,2,3.14*(R1/2)**2

K,,,,okok.org
K,,4
K,,10okok.org
L,1,2okok.org
L,2,3okok.org
LGEN,2,1,2,1,,-1.5okok.org
LGEN,2,3,4,1,,-1.5okok.org
LGEN,2,5,6,1,,-6

L,1,4okok.org
L,4,7okok.org
L,7,10okok.org
LGEN,2,9,11,1,4
LGEN,2,12,14,1,6

NUMMRG,KP, , , ,LOW

LSEL,S,LINE,,1,8,1okok.org
LESIZE,ALL,1
LSEL,S,LINE,,9,15,3
LSEL,A,LINE,,10,16,3
LESIZE,ALL,0.5okok.org
LSEL,S,LINE,,11,17,3
LESIZE,ALL,0.5
ALLSEL

A,1,2,5,4okok.org
A,2,3,6,5okok.org
A,4,5,8,7okok.org
A,5,6,9,8okok.org
A,7,8,11,10
A,8,9,12,11

AATT,1,,1okok.org
AMESH,1
AATT,2,,1okok.org
AMESH,3
AATT,3,,1okok.org
AMESH,2
AMESH,4,6,1
ALLSEL

TYPE,2
MAT,4
REAL,1
E,2,7
E,7,8
E,8,6
E,6,22
E,22,23
E,23,21
E,21,79
E,79,80
E,80,81
E,81,82

TYPE,3
MAT,5
REAL,2
E,8,46
E,46,71
E,71,62

TYPE,3okok.org
MAT,6okok.org
REAL,2
E,23,64
E,64,158
E,158,140

NSEL,S,LOC,X,0okok.org
NSEL,A,LOC,X,10
D,ALL,UX
NSEL,S,LOC,Y,-9
D,ALL,UY
ALLSEL

FINISHokok.org
/SOLUokok.org
ANTYPE,STATICokok.org
AUTOTS,ONokok.org
PRED,ONokok.org
LNSRCH,ONokok.org
NLGEOM,ONokok.org
NROPT,FULLokok.org
OUTRES,ALL,ALLokok.org
CNVTOL,F,,0.02,2,0.5okok.org
NSUBST,1okok.org

TIME,1
ACEL,,9.8
ESEL,S,MAT,,4,6,1
EKILL,ALL
ESEL,ALL
ESEL,S,LIVE
NSLE,S
NSEL,INVE
D,ALL,ALL

ALLSEL
SOLVE
SAVE,S1,DB

TIME,2okok.org
ESEL,S,MAT,,1okok.org
EKILL,ALL
ESEL,S,MAT,,4,5,1
EALIVE,ALLokok.org
NSLE,Sokok.org
DDELE,ALL,ALLokok.org
ESEL,ALL
ESEL,S,LIVE
NSLE,Sokok.org
NSEL,INVEokok.org
D,ALL,ALL

/post1
etable,etabs1,s,1okok.org
etable,etabs3,s,3okok.org
*dim,ArrS1,array,180
*dim,ArrS3,array,180
  *do,i,1,180okok.org
    *get,ArrS1(i),elem,i,etab,etabs1okok.org
    *get,ArrS3(i),elem,i,etab,etabs3
  *enddo

/prep7
  *do,i,13,180,1okok.org
    c=16$Fai=37$Rf=0.71$k=26.8$n=2.01$Kur=26.8$Nur=2.01$G=0.3
    $F=0.04$d=0.8
    *use,Duncan-Chang,c,Fai,Rf,k,n,Kur,Nur,G,F,d
  *enddo
        
M_et=Et
M_u=Mu
okok.org
/solu
alls
antype,,new,
solveokok.org
SAVE,S2,DB

TIME,3okok.org
ESEL,S,MAT,,2okok.org
EKILL,ALL
ESEL,S,MAT,,6
EALIVE,ALLokok.org
NSLE,Sokok.org
DDELE,ALL,ALLokok.org
ESEL,ALL
ESEL,S,LIVE
NSLE,Sokok.org
NSEL,INVEokok.org
D,ALL,ALL

/post1
etable,etabs1,s,1okok.org
etable,etabs3,s,3okok.org
*dim,ArrS1,array,180
*dim,ArrS3,array,180
  *do,i,1,180okok.org
    *get,ArrS1(i),elem,i,etab,etabs1okok.org
    *get,ArrS3(i),elem,i,etab,etabs3
  *enddo

/prep7
  *do,i,25,180,1okok.org
    c=16$Fai=37$Rf=0.71$k=26.8$n=2.01$Kur=26.8$Nur=2.01$G=0.3
    $F=0.04$d=0.8
    *use,Duncan-Chang,c,Fai,Rf,k,n,Kur,Nur,G,F,d
  *enddo
        
M_et=Et
M_u=Mu
okok.org
/solu
alls
antype,,new,
solveokok.org
SAVE,S3,DB
模型简介:这个是简单的基坑开挖问题,分两次开挖
回复 不支持

使用道具 举报

发表于 2011-10-15 18:17:48 | 显示全部楼层 来自 贵州黔西南州兴义
楼主你好!想问下你的命令流中,出现了很多okok.org,这是什么意思啊。还有,我也看过邓肯-张模型的理论,与楼主写的命令流中的模型有点不相同,想问问楼主你用的参考资料是什么书啊!求教了,谢谢
回复 不支持

使用道具 举报

发表于 2014-7-20 01:53:01 | 显示全部楼层 来自 云南昆明
FINISH
/CLEAR
/FILNAME,PIT
/TITLE,PIT EXCAVATION
/PREP7
*dim,Str_max,array,180
*dim,S_max,Array,180
*dim,Settlement,array,5
!*dim,M_et,array,5
!*dim,M_u,array,5

*create,Duncan-Chang
*afun,deg
*set,Pa,1e5
*set,P1,-ArrS3(i)
*set,P3,-ArrS1(i)
*IF,P3,LT,0.1*Pa,then
P3=0.1*Pa
*ENDIF
Str=2*(c*cos(Fai)+P3*sin(Fai))/(1-sin(Fai))
S=(P1-P3)/Str
*if,S,GT,0.95,then
S=0.95
*endif

*if,Str_max(i),gt,P1-P3,and,S_max(i),gt,S,then
Et=Kur*Pa*(P3/Pa)**Nur
*elseif,Str_max(i),gt,P1-P3,and,S_max(i),le,S,then
Ei=k*Pa*(P3/Pa)**n
Et=Ei*(1-Rf*S)**2
S_max(i)=S
*elseif,Str_max(i),le,P1-P3,and,S_max(i),gt,S,then
Ei=k*Pa*(P3/Pa)**n
Et=Ei*(1-Rf*S)**2
Str_max(i)=P1-P3
*elseif,Str_max(i),le,P1-P3,and,S_max(i),le,S,then
Ei=k*Pa*(P3/Pa)**n
Et=Ei*(1-Rf*S)**2
Str_max(i)=P1-P3
S_max(i)=S
*endif
!----------------------------------------------------------------------
Ei=k*Pa*(P3/Pa)**n
A=(p1-p3)*d/(Ei*(1-Rf*S))
Mu=(G-F*log10(P3/Pa))/(1-A)**2
*if,Mu,GE,0.49,then
Mu=0.49
*endif
mp,ex,i,Et
mp,nuxy,i,Mu
!mpchg,i,
*end

*SET,T1,0.5
*SET,R1,25E-3
*SET,E1,4E7
*SET,NU1,0.32
*SET,DENS1,1800
*SET,E2,1.85E10
*SET,NU2,0.2
*SET,E3,2.1E11
*SET,NU3,0.3

ET,1,42
KEYOPT,1,3,2
ET,2,3
KEYOPT,2,6,1
ET,3,1

*DO,I,1,3,1
MP,EX,I,E1
MP,PRXY,I,NU1
MP,DENS,I,DENS1
!TB,DP,I
!TBDATA,1,C_R,FI_R
*ENDDO

MP,EX,4,E2
MP,PRXY,4,NU2

*DO,I,5,6,1
MP,EX,I,E3
MP,PRXY,I,NU3
*ENDDO

R,1,T1,T1*T1*T1/12,T1
R,2,3.14*(R1/2)**2

K,,,,
K,,4
K,,10
L,1,2
L,2,3
LGEN,2,1,2,1,,-1.5
LGEN,2,3,4,1,,-1.5
LGEN,2,5,6,1,,-6

L,1,4
L,4,7
L,7,10
LGEN,2,9,11,1,4
LGEN,2,12,14,1,6

NUMMRG,KP, , , ,LOW

LSEL,S,LINE,,1,8,1
LESIZE,ALL,1
LSEL,S,LINE,,9,15,3
LSEL,A,LINE,,10,16,3
LESIZE,ALL,0.5
LSEL,S,LINE,,11,17,3
LESIZE,ALL,0.5
ALLSEL

A,1,2,5,4
A,2,3,6,5
A,4,5,8,7
A,5,6,9,8
A,7,8,11,10
A,8,9,12,11

AATT,1,,1
AMESH,1
AATT,2,,1
AMESH,3
AATT,3,,1
AMESH,2
AMESH,4,6,1
ALLSEL

TYPE,2
MAT,4
REAL,1
E,2,7
E,7,8
E,8,6
E,6,22
E,22,23
E,23,21
E,21,79
E,79,80
E,80,81
E,81,82

TYPE,3
MAT,5
REAL,2
E,8,46
E,46,71
E,71,62

TYPE,3
MAT,6
REAL,2
E,23,64
E,64,158
E,158,140

NSEL,S,LOC,X,0
NSEL,A,LOC,X,10
D,ALL,UX
NSEL,S,LOC,Y,-9
D,ALL,UY
ALLSEL

FINISH
/SOLU
ANTYPE,STATIC
AUTOTS,ON
PRED,ON
LNSRCH,ON
NLGEOM,ON
NROPT,FULL
OUTRES,ALL,ALL
CNVTOL,F,,0.02,2,0.5
NSUBST,1

TIME,1
ACEL,,9.8
ESEL,S,MAT,,4,6,1
EKILL,ALL
ESEL,ALL
ESEL,S,LIVE
NSLE,S
NSEL,INVE
D,ALL,ALL

ALLSEL
SOLVE
SAVE,S1,DB

TIME,2
ESEL,S,MAT,,1
EKILL,ALL
ESEL,S,MAT,,4,5,1
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ESEL,ALL
ESEL,S,LIVE
NSLE,S
NSEL,INVE
D,ALL,ALL

/post1
etable,etabs1,s,1
etable,etabs3,s,3
*dim,ArrS1,array,180
*dim,ArrS3,array,180
*do,i,1,180
*get,ArrS1(i),elem,i,etab,etabs1
*get,ArrS3(i),elem,i,etab,etabs3
*enddo

/prep7
*do,i,13,180,1
c=16$Fai=37$Rf=0.71$k=26.8$n=2.01$Kur=26.8$Nur=2.01$G=0.3
$F=0.04$d=0.8
*use,Duncan-Chang,c,Fai,Rf,k,n,Kur,Nur,G,F,d
*enddo

M_et=Et
M_u=Mu

/solu
alls
antype,,new,
solve
SAVE,S2,DB

TIME,3
ESEL,S,MAT,,2
EKILL,ALL
ESEL,S,MAT,,6
EALIVE,ALL
NSLE,S
DDELE,ALL,ALL
ESEL,ALL
ESEL,S,LIVE
NSLE,S
NSEL,INVE
D,ALL,ALL

/post1
etable,etabs1,s,1
etable,etabs3,s,3
*dim,ArrS1,array,180
*dim,ArrS3,array,180
*do,i,1,180
*get,ArrS1(i),elem,i,etab,etabs1
*get,ArrS3(i),elem,i,etab,etabs3
*enddo


*do,i,25,180,1
c=16$Fai=37$Rf=0.71$k=26.8$n=2.01$Kur=26.8$Nur=2.01$G=0.3
$F=0.04$d=0.8
*use,Duncan-Chang,c,Fai,Rf,k,n,Kur,Nur,G,F,d
*enddo

M_et=Et
M_u=Mu

/solu
alls
antype,,new,
回复 不支持

使用道具 举报

发表于 2016-3-23 10:45:02 | 显示全部楼层 来自 大连理工大学
inkey 发表于 2014-7-20 01:53
FINISH
/CLEAR
/FILNAME,PIT

这是你编好的邓肯-张模型吗?只用生死单元没用重启动分析?生死单元允许进入后处理吗?我最近也在编这个,被卡住了。。。
回复 不支持

使用道具 举报

发表于 2016-3-23 10:47:58 | 显示全部楼层 来自 大连理工大学
虽然时间相隔甚远,但是还是想问一下楼主现在还有邓肯-张模型的程序吗?我已经被这个问题卡了很久了,如果还有希望楼主可以慷慨相赠借我参考下,之所以请求相赠绝非自己偷懒,而是自己愚笨已经研究很久了依然无果,不敢再继续耗费时间了,我的邮箱:2623144801@qq.com非常感谢!
回复 不支持

使用道具 举报

发表于 2017-2-16 02:14:40 | 显示全部楼层 来自 云南
仿真计算,重在真这字上呀,ANSYS结构计算应是一流的了,ABQUS,FLAC,ADINA,各软都有强项,我属笨点的行例,多数用ANSYS,别的玩不转啊!有空我发在实例里,找找破电脑里的东东,不过这些都是世界一流工程,有小湾,锦屏的部分仿真计算,也有混凝土面板坝工程,(计算的重点是要参数,参数要可XXX准);P;P算出来的就是八九不离十了。:P
回复 不支持

使用道具 举报

发表于 2019-12-25 20:53:43 | 显示全部楼层 来自 湖北武汉
跪求,可以看看我的命令流哪里出现问题了吗?
难道inistate命令不可用嘛

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 09:31 , Processed in 0.044847 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表