fixed.h

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

//////////////////////////////////////////////////////////////////////////////
//
//	 Fixed 函数：	
//
//		功能：处理位移约束
//
//		输入变量：	k:  表示在第k个自由度上加位移约束		
//
//		            u： 加上的约束的位移值
//
//                  jd: 指示第k个自由度在一维总刚度矩阵中的位置
//
//                  zk： 一维总刚度矩阵 zk[0]未用
//
//                  p：  载荷向量,长度为NF+1，P[0]未用
//
//                  NF:  结构的自由度总数,值为NodeNumber*ND
//
//////////////////////////////////////////////////////////////////////////////

#include "skdd.h"
#ifndef __FIXED_
#define __FIXED_
void Fixed(int k,double u,double *zk, double*P,int*jd,int NF);

void Fixed(int k,double u,double *zk, double*P,int*jd,int NF)

{    int i;
     for(i=1;i<=NF;i++)
	   if(i!=k)
		 P[i]-=GetStiff(i,k,zk,jd)*u;
	        else 
				P[i]=u;                      //修改等效节点载荷;
    for(i=1;i<=NF;i++)
	{ SetStiff(k,i,0,zk,jd);
      SetStiff(i,k,0,zk,jd);
	}
    SetStiff(k,k,1.0,zk,jd);
}                                           //修改刚度矩阵k行k列;
#endif