有没有人做个陈火红二次开发上移动载荷那个例子？

hanfeifei

有没有人做个陈火红二次开发上移动载荷那个例子？用户子程序该怎么写呢？有人做通过么？讨论下啊？或者牛人出来指点下啊！

liebebilder

我做通了，有什么问题？

hanfeifei

终于等到高人回复了，请问你能否把做通的程序发我看看，我想确定我是因为程序的问题，还是因为几个软件连接的问题而做不通。我的邮箱：aldan@163.com,如盟赐教，感激不尽。

hanfeifei

我给您发短消息了，请查收，急切需要您的帮助。

liebebilder

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

hanfeifei

真是十分感谢！你真是大好人啊！

liebebilder

我发到你邮箱了，不知道你收到了吗，祝顺利

hanfeifei

原帖由 liebebilder 于 2008-8-5 17:01 发表
我发到你邮箱了，不知道你收到了吗，祝顺利

感谢！我回复邮件给你了！谢谢！

hanfeifei

原帖由 liebebilder 于 2008-8-5 17:01 发表
我发到你邮箱了，不知道你收到了吗，祝顺利

我的问题基本上已经解决了，非常感谢你的帮助！程序的安装版本和顺序确实很重要，我按照marc操作手册上的说明进行了环境变量的设置，这样把问题解决了。谢谢你!