我在做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
模型简介:这个是简单的基坑开挖问题,分两次开挖 |