[Mohamed.gaber]

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