Write ‘inp’ file#

Consider an FEM problem of a plate with a circular hole on the edge (see figure below).

../../_images/example_mesh.svg

An example of ‘inp’ file with some output commands is given below :

****calcul                    % beginning of a FEM problem (note no indent)
 ***mesh small_deformation    % defines the element formulation
 ***resolution newton         % resolution method, and its options
  **sequence                  % loading sequence
    ...
 ***bc                        % definition of boundary conditions
  **impose_nodal_dof          % creation of BCs
    top      U3 3. time       % impose U3=3*time on 'top' nset
    x=0      U1 0.0           % ...
    y=0      U2 0.0           %
    z=0      U3 0.0           %

 ***output                    % specifies output options
  **save_material_rotations
  **verbose                   % verbose output to see the front size etc
  **contour                   % contour output for extrapolated integ->nodal
  **contour_by_element        % contour values of integ values element by element
  **value_at_integration      % these are the REAL stresses and strains.
  **curve                     % curve output into default file plate.test
   *node_var node2 U2         % print U2 values for node2 nset
   *nset_var top  U2          % print the applied load (F2) over the surface
                              % defined by the node set 'top'

 ***material                  % material definition
  *this_file                  % the behavior is defined within the same file
****return