shape1d.h

这是等参单元的有限元程序

//#include <stdlib.h>
//#include <iomanip.h>
//#include <fstream.h>
////////////////////////////////////////////////////////////////////////////////
//
//	shape 函数：一维lagrange插值函数
//  shape_a函数：对一维lagrange插值基函数局部坐标求导
//  
//////////////////////////////////////////////////////////////////////////////

#ifndef __SHAPE1D_
#define __SHAPE1D_
void shape(double A,double *N,int *Ele);
void shape_a(double A,double *NA,int *Ele);

void shape(double A,double *N,int *Ele)
/* 参数说明如下：
    A—局部坐标；
    N—长度为3的数组的指针，存放三个插值基函数的值；
    Ele—长度为3的数组的指针，存放插值区间(单元)的节点
              构成；如果某元素为零，则意味该节点不存在。
*/
{   if(Ele[2]) N[2]=1.0-A*A; else N[2]=0.0;

      N[0]=0.5*(1.0-A)-0.5*N[2];
      N[1]=0.5*(1.0+A)-0.5*N[2];
   return;
}

void shape_a(double A,double *NA,int *Ele)
/* 参数说明如下：
    A—局部坐标；
    NA—长度为3的数组的指针，存放三个插值基函数的导数值；
    Ele—长度为3的数组的指针，存放插值区间(单元)的节点
              构成；如果某元素为零，则意味该节点不存在。
*/
{ if(Ele[2]) NA[2]=-2*A; else NA[2]=0.0;

   NA[0]=-0.5-0.5*NA[2];
   NA[1]=0.5-0.5*NA[2];
   return;
}
#endif