FEAP User Forum

FEAP => Input File Issues => Topic started by: JoGe on August 13, 2020, 06:06:15 AM

Title: proportionally applied rotational body force
Post by: JoGe on August 13, 2020, 06:06:15 AM
Dear Prof. Taylor, dear FEAP users,

I am trying to solve a problem in a quasi-static manner in parallel using FEAP 8.5 where a body which has one side clamped experiences pressure on its "free" surfaces as well as a rotational body force.
As you see below the part is rotating around its z-axis:

GLOB
OMEG,,100.0
OMEG,node,0.0,0.0,0.0
OMEG,vector,0.0,0.0,-100.0

The pressure is applied using the pressure surface element along with the corresponding material set:
MATE,2
  PRES
    LOAD,-1.0,1
    FOLL

Unfortunately the <prop,,1> command in my solve file does only capture the pressure load - not the rotational body force!

Is there any way of applying a rotational body force proportionally? I have tried many different ways yet without success.

I would appreciate your support!

Kind regards,
JoGe
Title: Re: proportionally applied rotational body force
Post by: Prof. R.L. Taylor on August 13, 2020, 08:22:59 AM
The OMEGa is an angular velocity, do you want to change the value during the solution?  If so, the code needs to be modified.  the loading and its tangent for finite deformation problems are computed in pbodyw.f in the directory ./elements/solid3d.  Addionally, you could set a proportional load number in ./program/pglobal.f to be used.  This routine computes inertial loading due to the spin velocity.

If you are just doing a finite rotation there is the SPIN command -- see Sect 5.8 of current user manual.