C I have just sent the plotv user subroutine tailored to calculate the cord
C angle in the deform
C ed configuration. The analysis of cord-rubber composites is carried out in the
C total lagrangia
C n framework and hence only the initial cord angle info. is available.
C  This can be used for the following products : tires, air springs, hoses,
C bladders,
C robotic pickers and grippers (in general any cord-rubber composite structure).
C Please let me know if there are any questions ....... GOOD LUCK !!
C                                                                Sanjay
C Choudhry,
                                                                (415) 329-6836.


c********************************************************************************
c
c      Calculation of the Updated Cord Angles in the Deformed Geometry
c      due to the Pantographing Mechanism in Inflated Cord-Rubber Composites
c
c********************************************************************************
c
      subroutine plotv(v,s,sp,etot,eplas,ecreep,t,m,idum0,layer,idum1,
     *                 idum2,jpltcd)
c* * * * * *
c
c     select a variable contour plotting (user subroutine).
c
c     v          variable
c     s (idss)   stress array
c     sp         stresses in preferred direction
c     etot       total strain (generalized)
c     eplas      total plastic strain
c     ecreep     total creep strain
c     t          current temperature
c     m          element number
c     nn         integration point number
c     layer      layer number
c     ndi (3)    number of direct stress components
c     nshear (3) number of shear stress components
c     lsub       control variable for the part of element assembly function
c
c* * * * * *
c
      implicit real*8 (a-h,o-z)
                                              dp
      dimension s(1),etot(1),eplas(1),ecreep(1),sp(1)
      dimension m(2),oa(2,2),dispt(12)
      include '/disk1/sanjay/common/blnk'
      include '/disk1/sanjay/common/deveub'
      include '/disk1/sanjay/common/dimen'
      include '/disk1/sanjay/common/elmcom'
      include '/disk1/sanjay/common/pival'
      include '/disk1/sanjay/common/space'
      include '/disk1/sanjay/common/parpts'
      include '/disk1/sanjay/common/concom'
      include '/disk1/sanjay/common/heat'
      include '/disk1/sanjay/common/arrays'
      include '/disk1/sanjay/common/array5'
      include '/disk1/sanjay/common/lass'
      include '/disk1/sanjay/common/array2'
c
c
c Input file must include LARGE DISP option (but no Update flag !!!!)
c
c The basis for calculation is : n = F.N
c
c Calculate the initial angle (undeformed geometry)
c
      lofr=(n-1)*nelstr
      theta=vars(iangle+lofr)*pi/180.d0
      ct=dcos(theta)
      st=dsin(theta)
      itel=2
      v=0.d0
      if(jtype.ne.18) return
      if(nn.eq.1) then
      call scla(dircos,0.d0,3,3,0)
      do 761 ijk=1,3
761   dircos(ijk,ijk)=1.d0
      call zdirc(vars(idrcos),4)
      lovl=4
c
c Transformation matrix
c
      lsub=7
      call ellib
c
c Wt. functions
c
      lsub=5
c
      call ellib
      if(jparel.eq.1) then
        call scla(vars(ibp),0.d0,ngenel,n1*4,0)
        if(ismall.eq.0) call scla(vars(idumpp),0.d0,nldsrc,n1*4,0)
        if(ismall.eq.0) call scla(vars(idump1p),0.d0,nldsrc,n1*4,0)
      endif
      endif
c
c Linear beta matrix
c
      lsub=2
      if(jparel.eq.0) then
        call scla(vars(idump),0.d0,n1,1,0)
        lsub=4
      endif
      call ellib
      call elvec(dispt,nnode,vars(idsxt),ndegmx,ndegel,numnp,5,lm)
      if(jparel.eq.1) then
        call dispdp(dispt,oa,vars(idump1p),itel,1.d0,ndi,ityp,n1,nn)
      else
        call dispdr(dispt,oa,vars(idump1),itel,1.d0,ndi,ityp,n1,nn)
      endif
c
c oa has deformation gradient
c
      a1=oa(1,1)*ct+oa(1,2)*st
      a2=oa(2,1)*ct+oa(2,2)*st
      cost=a1/dsqrt(a1**2+a2**2)
      tnew=dacos(cost)*180.d0/pi
c
c Plot variables : initial and final angles
c
      if(jpltcd.eq.1) v=theta
      if(jpltcd.eq.2) v=tnew
c
      if((m(1).eq.430).or.(m(1).eq.470).or.(m(1).eq.510)) then
       write(49,*) 'increment=',inc,' iteration=',ncycle
       write(49,*) 'element  =',m(1),' int pt=',nn,' n=',n
       write(49,*) 'init angle=',theta*180.d0/pi,' final angle=',tnew
      endif
c
      return
      end