讲稿3-6(3学时).pdfVIP

3.5 外力功 节点等效载荷 作用在物体上的外载荷，按其作用形式可分为体积力、表面 力和集中力。对于表面力和集中力，约定： 集中力作用处均布置节点； 表面力作用在一个或若干个单元的完整边界上； 尽可能避免使用集中力。 不建议使用的有限元划分方式。 P 作用在有限元离散体系上的各种形式的外载荷，可通过适当 方式转化成节点等效载荷。 一. 体积力的节点等效载荷 体积力是指外载荷作用在体系的各个微元上。重力(惯性力) 是最常见的体积力 设单位体积中的体积力为f =[f x f y ]T -2 -2 f 的量纲为kg ·m ·s ，即密度与加速度的乘积。 体积力f 在一个单元上做功 T T T T B u f d V N f d V P  e  e e   e e PeB N T f d  e [N f N f  N f N f ]T d  1 x 1 y 8 x 8 y  e 上式表明，作用在单元内每个质点上的体积力，在(虚)功等 效的基础上，可转化成作用在单元各个节点上等效载荷。 体积力f 在所有单元上的总功为 VT P B VT A T P B VT A T P B VT P B  e e  e e  e e e e e 定义体积力的节点等效载荷PB P B A T P B  e e e 求和运算的实质，是将按不同单元计算出的同一自由度上的 体积力等效载荷，进行叠加的过程。 轴对称情况下，离心力的等效载荷 1 1  2  B T T  x Pe N f d 2π  N  xdetJdd  1 1  f y  e   其中，是物体密度，是物体绕y 轴作定轴转动的角速度。 计算体积力节点等效载荷的程序实现 void CacuBF(double Density, double Thick, double ax,double ay, double *xo, double *yo, double *Load, int *Ele,int Iopt, int nGauss) /*参数说明： Density—单元密度； Thick—单元厚度，仅用于平面应力情况下； ax—在x方向的惯性加速度，若Iopt==2 ，则为转动角速度；