MohamedGaber
Civil/Environmental
- Jun 29, 2015
- 1
Hi all,
I am using Ansys Mechanical APDL16 in modelling a concrete Dam in which it is subjected to varying pressure load.
The code shall be as following :
Code:
/prep7
!Element type
ET,1,user300
!Keyopt,1,1,1
!USRELEM, NNODES, NDIM, KeyShape, NREAL, NSAVEVARS, NRSLTVAR, KEYANSMAT, NINTPNTS, KESTRESS, KEYSYM
USRELEM,8,3,Brick,8,168,112,0,8,3,1
USRDOF,define,ux,uy,uz,pres
R, 1, 25000, .2,300, .96, .0002, .5 !, 2.4e-9
RMORE, 10, 2
ET,2,Plane182
!geometry
TYPE,2
!Keypoints for triangular shapes
K,1,46,0
K,2,72,0
K,3,46,33
K,4,20,33
K,5,20,66
K,6,14,66
K,7,14,102
l,1,2
l,1,3
l,2,3
l,3,4
l,3,5
l,4,5
l,5,6
l,6,7
l,5,7
!Create the triangles
lsel,s,line,,1,3,1
Al,all
lsel,s,line,,4,6,1
AL,all
lsel,s,line,,7,9,1
Al,all
!creating rectangles
BLC4,0,0,14,33
BLC4,0,33,14,33
BLC4,0,66,14,36
BLC4,14,0,6,33
BLC4,14,33,6,33
BLC4,20,0,26,33
ASEL,ALL
AGLUE, ALL
!THEN I DISCRETED THE MODEL!
TYPE,1
EEXTRUDE, PLANE, 1, 1, 100
!Creating a line for the pressure load in z=100
KNODE, 275,275
KNODE, 308,308
L,275,308
Esel,All
ASEL,none
EMODIF, All, TYPE, 1 ! change to user300
ASEL,none
!solution
/SOLU
nSEL, None
nSEL, S, LOC, y,0
d,all,all ! displacement constraints
nsel,all
!surface load on y location from 66 to 102 (z=0)
sfl,40,pres,225.830,225830+(1001.000)*(36/102)
!surface load on y location from 66 to 102 (z=100)
sfl,17,pres,-(225.830),-(225.830+(1001.000)*(36/102))
nsel,all
SFTRAN
SBCTRAN
AUTOTS,on
NROPT,UNSYM
OUTRES,ALL,1 ! WRITE ALL OUTPUT
nsel,all
SOLVE
FINISH
The model acts like there is no load. I think there is a problem in surface load applying method.
Thank you in advance,
_________________
Mohamed Gaber
Student at the Technical University of Dresden
I am using Ansys Mechanical APDL16 in modelling a concrete Dam in which it is subjected to varying pressure load.
The code shall be as following :
Code:
/prep7
!Element type
ET,1,user300
!Keyopt,1,1,1
!USRELEM, NNODES, NDIM, KeyShape, NREAL, NSAVEVARS, NRSLTVAR, KEYANSMAT, NINTPNTS, KESTRESS, KEYSYM
USRELEM,8,3,Brick,8,168,112,0,8,3,1
USRDOF,define,ux,uy,uz,pres
R, 1, 25000, .2,300, .96, .0002, .5 !, 2.4e-9
RMORE, 10, 2
ET,2,Plane182
!geometry
TYPE,2
!Keypoints for triangular shapes
K,1,46,0
K,2,72,0
K,3,46,33
K,4,20,33
K,5,20,66
K,6,14,66
K,7,14,102
l,1,2
l,1,3
l,2,3
l,3,4
l,3,5
l,4,5
l,5,6
l,6,7
l,5,7
!Create the triangles
lsel,s,line,,1,3,1
Al,all
lsel,s,line,,4,6,1
AL,all
lsel,s,line,,7,9,1
Al,all
!creating rectangles
BLC4,0,0,14,33
BLC4,0,33,14,33
BLC4,0,66,14,36
BLC4,14,0,6,33
BLC4,14,33,6,33
BLC4,20,0,26,33
ASEL,ALL
AGLUE, ALL
!THEN I DISCRETED THE MODEL!
TYPE,1
EEXTRUDE, PLANE, 1, 1, 100
!Creating a line for the pressure load in z=100
KNODE, 275,275
KNODE, 308,308
L,275,308
Esel,All
ASEL,none
EMODIF, All, TYPE, 1 ! change to user300
ASEL,none
!solution
/SOLU
nSEL, None
nSEL, S, LOC, y,0
d,all,all ! displacement constraints
nsel,all
!surface load on y location from 66 to 102 (z=0)
sfl,40,pres,225.830,225830+(1001.000)*(36/102)
!surface load on y location from 66 to 102 (z=100)
sfl,17,pres,-(225.830),-(225.830+(1001.000)*(36/102))
nsel,all
SFTRAN
SBCTRAN
AUTOTS,on
NROPT,UNSYM
OUTRES,ALL,1 ! WRITE ALL OUTPUT
nsel,all
SOLVE
FINISH
The model acts like there is no load. I think there is a problem in surface load applying method.
Thank you in advance,
_________________
Mohamed Gaber
Student at the Technical University of Dresden