Issue with high number of DOFs

[PK]

Hello,
 
I am a new user of FEAP. My problem has 19 degrees of freedom and hence FEAP is not taking input- EBOUndary, EDISP, EPRO- of all the dofs. Is there any way to make FEAP take more than 16  items of data in these records?  Any suggestion will be greatly appreciated.

Thanks,
Prasenjit

[FEAP_Admin]

wrapped the additional data onto the next line; i.e. start the next line with item 17 18 etc.  Then go to the following line for
the next node.

[Colin McAuliffe]

Hi Prasenjit, in addition to wrapping the data I also had to modify program/plinka.f around line 100 from

  40  if(.not.pcomp(yyy(i:i+2),fend,3)) then
        iclink = iclink + 1
        go to 10
      else
        close(iop)
      endif
      return

to

  40  if(.not.pcomp(yyy(i:i+2),fend,3)) then
        iclink = iclink + 1
        go to 10
      else
          !add extra blank line before end of file
        write(iop,*)
        close(iop)
      endif
      return

This adds an extra blank line to the temporary files generated by commands like ebou/edisp etc. and prevents a premature end of file when data is contained on multiple lines.

[PK]

Thank you so much. The issue is resolved now.

[mubeen]

wrapped the additional data onto the next line; i.e. start the next line with item 17 18 etc.  Then go to the following line for
the next node.

Hi FEAP Admin,
I am using FEAP 8.2 and I have wrapped the additional values in INITial DISPlacement record on next line, i.e.

Code: [Select]

BATCh
INITial DISPlacement
END
1   0 0.0 0.0 0.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0
 3.0
27  0 0.0 0.0 0.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0
 3.0

and when I submitted the input file to FEAP, the result was FEAP crash, with an error in Output file:

Code: [Select]

     Material  Number   1 for User Element Type 11
                         
     Degree of Freedom Assignments    Local    Global
                                     Number    Number
                                        1         0
                                        2         0
                                        3         0
                                        4         0
                                        5         0
                                        6         0
                                        7         0
                                        8         0
                                        9         0
                                       10         0
                                       11         0
                                       12         0
                                       13         0
                                       14    210840
                                       15    130420                                                                         

 *ERROR* ENDCLR: End of file encountered in TINPUT
         Unit Number =  15


I have then modified the program/plinka.f (as suggested by Colin, above), by adding

Code: [Select]


  40  if(.not.pcomp(yyy(i:i+2),fend,3)) then
        iclink = iclink + 1
        go to 10
      else
        write(iop,*) !extra line as suggested by Colin McAuliffe
        close(iop)
      endif

but the error persists. I guess the solution suggested by Colin is for EBOUN, EDISP records, and perhaps it does not work with record inside BATCh - END.

Is there any possibility to specify initial values for more than 15 DOF/displacements in FEAP input file?

I have attached the input file for this problem.

[Prof. R.L. Taylor]

I am traveling so cannot check much.

You need to debug.  Input should call routine pinitl.f which then calls routine genvec.f.

Ver 8.2 is old and has many issues related to dimensioning. You should update to new version if possible.

Will check your file when I finish travels

[Prof. R.L. Taylor]

The attached file excutes correctly in the current release of feap.  I can not test your version still.

Note the places I modified to get to read -- I do not know if 8.2 has all of them and since I do not have your user element do not know if that will work properly for all the analysis in either version of feap