forums.silverfrost.com

[yuvadon]

[code:1:702e135d02]c#######################################################################
c
c SPH: Smoothed Particle Hydrodynamics
c
c#######################################################################
c
program SPH
c
implicit real*8(a-h,o-z)
c
include 'parameter.f'
c
parameter(
& TOTAL_TIME=10.0D-6, ! ƒVƒ~ƒ…ƒŒ�[ƒVƒ‡ƒ“‚̍�‘�ŽžŠ�
& MAXLOOP=100000) ! ƒVƒ~ƒ…ƒŒ�[ƒVƒ‡ƒ“‚̍�‘�ŒJ‚�•�‚��”
c
common/model/ DISA,XMAX,XMIN,XSIZE,YMAX,YMIN,YSIZE,
& NUMX1,NUMY1,NUMX2,NUMY2,NUMP,NP1,NP2
c
common/particle/ X(NDIM,MAXP),X0(NDIM,MAXP)
c
common/property/ RHO1,RHO2,CS1,CS2,G1,G2,PK1,PK2
& RHO(MAXP),PMASS(MAXP),CS(MAXP),PK(MAXP)
c
common/composite/ CMAT(3,3,2),MATER(MAXP)
c
common/boundary/
& V_BC_GIVEN(NDIM,MAXBC),V_LOAD_GIVEN(NDIM,MAXBC),
& I_BC_GIVEN(NDIM,MAXBC),I_LOAD_GIVEN(NDIM,MAXBC),
& N_BC_GIVEN(NDIM),N_LOAD_GIVEN(NDIM),SPEED
c
common/calc/ V(NDIM,MAXP),VISC(NEIGHBOR,MAXP),P(MAXP),E(MAXP)
c
common/param/ H0,H,DT,PDT
c
common/neighbor_table/
& NEIGH(NEIGHBOR,MAXP),NNEIGH(MAXP),
& NEIGH0(NEIGHBOR,MAXP),NNEIGH0(MAXP),JCONNECT0(NEIGHBOR,MAXP),
& RIJ(NDIM,NEIGHBOR,MAXP),RIJ1(NEIGHBOR,MAXP),RIJ2(NEIGHBOR,MAXP),
& RIJ0(NEIGHBOR,MAXP)
c
common/kernel/ GRADW(NDIM,NEIGHBOR,MAXP)
c
common/stress/
& SIGM(3,MAXP),HENSAS(3,MAXP),HENSASOLD(3,MAXP),STRNDOT(3,MAXP),STRN(3,MAXP),STRNOLD(3,MAXP)
c
common/vardot/ DRHODT(MAXP),DVDT(NDIM,MAXP),DEDT(MAXP)
c
common/rkgill/ DX(NDIM,MAXP),DV(NDIM,MAXP),DR(MAXP),DE(MAXP)
c
common/yield/ SGMY(MAXP),EPSP(MAXP),T(MAXP),JYIELD(MAXP)
c
common/damage/ DMG(MAXP)

c
common/phase/ MELT(MAXP)
c
c==== output files
open(2,file='rst.txt')
open(3,file='out_damage.plt')
open(4,file='rst_damage.csv')
open(7,file='out_stress.plt')
open(100,file='check_coordinates.csv')
c
c==== input
call inp(angle_inc)
c
call set_param
c
c==== initialize
do ip=1,NUMP
X(1:NDIM,ip) = X0(1:NDIM,ip) ! ��•W
E(ip) = 0.0D0 ! ƒGƒlƒ‹ƒM�[
do k=1,3
SIGM(k,ip) = 0.0D0 ! ‰ž—�
HENSAS(k,ip) = 0.0D0 ! •΍�‰ž—�
HENSASOLD(k,ip) = 0.0D0
STRNDOT(k,ip) = 0.0D0 ! ‚�‚�‚�‘�“x
STRN(k,ip) = 0.0D0
STRNOLD(k,ip) = 0.0D0

enddo
c
DX(1:NDIM,ip) = 0.0D0 ! ƒ‹ƒ“ƒQ�EƒNƒbƒ^–@
DV(1:NDIM,ip) = 0.0D0 ! ‚�‚�ŽžŠԐ�•�‚�Žg‚�
DR(ip) = 0.0D0
DE(ip) = 0.0D0
c
SGMY(ip) = YIELD0 ! �~•š‰ž—�
EPSP(ip) = 0.0D0 ! ‘Š“–‘Y��‚�‚�‚�
T(ip) = 300.0D0 ! ‰�“x
c
DMG(ip) = 0.0D0 ! ‘���ƒpƒ‰ƒ��[ƒ^
MELT(ip) = 0.0D0
enddo

[DanRRight]

LOL Ian already helped you with that...in August

http://forums.silverfrost.com/viewtopic.php?t=2627&highlight=

Say thank you