|
回复 4楼 的帖子
subroutine forcem(press,th1,th2,nn,n)
include 'D:/MSC.Software/MSC.Marc/2005/marc2005/common/implicit'
common/lpres3/prnorm(3)
!dimension n(7)
c* * * * * *
c
c defined non-uniformed distributed force on an element.
c
c press distributed load increment magnitude
c if follower force then give total magnitude
c th1 coordinate
c th2 coordinate
c nn integration point number
c n(1) user element number
c n(2) parameter identifying the type of load
c n(3) is the integration point number
c n(4) not used
c n(5) is the distributed load index
c for 2005+ input format this is the boundary condition id
c n(6) =0 if conventional pressure
c =1 if user returns 2 or 3 components for pressure in global direction
c n(7) is the internal element number
c
c for distributed load in a given direction
c prnorm is the direction cosine of the direction of the load
c with respect to the global system
c
c* * * * * *
dimension th1(3),th2(3),n(7)
include 'D:/MSC.Software/MSC.Marc/2005/marc2005/common/creeps'
b=1.0d0
v=sqrt(2.0d0)
distan=v*(cptim+timinc)
xc=0.5d0*sqrt(2.0d0)*distan
yc=0.5d0*sqrt(2.0d0)*distan
xmin=xc-b/2.0d0
xmax=xc+b/2.0d0
ymin=yc-b/2.0d0
ymax=yc+b/2.0d0
press=0.0d0
if (th1(1).le.xmax.and.th1(1).ge.xmin.and.
$ th1(2).le.ymax.and.th1(2).ge.ymin) then
press=2.0d0
end if
th2(1)=0.0d0
th2(2)=0.0d0
th2(3)=-1.0d0
open(1,file="exam_run_FORfile.txt",form="unformatted")
write(1)"this subroutine has been called"
close(1)
return
end subroutine forcem |
评分
-
1
查看全部评分
-
|