杆刚度校核c语言程序,c语言求解单元刚度矩阵

运用c语言求解有限元单元刚度矩阵

/*函数功能：计算单元刚度矩阵

函数入口参数：实型数组nodeCoo[8][3]，存储单元的8个结点的总体自由度编号实型数组elementK[24][24]，存储单元刚度矩阵

函数返回值：无

void Getnumber(double nodeCoo[8][3],double elementK[24][24])

{

double Gauss[2]={0.57735026919,-0.57735026919}; //高斯点

double WeiG[2]={1.0,1.0}; //权重

int local[3][8]={{1,1,-1,-1,1,1,-1,-1},{-1,1,1,-1,-1,1,1,-1},{-1,-1,-1,-1,1,1,1,1}}; //结点局部坐标

double diffN[3][8]; //插值函数关于局部坐标偏导数

double diffNInverJ[3][8]={0}; //插值函数关于整体坐标偏导数

double DJ;

double J[3][3];

double *MatFirstadd;

double InverJ[3][3]; //存储J的逆矩阵 double B[6][24]={0}; //存储B矩阵 double BInver[24][6]; //存储B的转置 double *MatrixInverFirstadd;

double BInverD[24][6]={0}; //存储B的转置与D的乘积

int i,j,ksi,eta,zeta,m,n;

double Gaussksi,Gausseta,Gausszeta;

double G=E/(2*(1+v)),lam=E*v/((1+v)*(1-2*v)); //已知常数E，v

double

D[6][6]={{lam+2*G,lam,lam,0,0,0},{lam,lam+2*G,lam,0,0,0},{lam,lam,lam+2*G,0,0,0},{0,0,0,G,0,0},{0,0,0,0,G,0},{0,0,0,0,0,G}}; //D矩阵

double BTDBJ[24][24]={0};

for(ksi=0;ksi<2;ksi++) //以ksi方向高斯点计算积分

{

Gaussksi=Gauss[ksi];

for(eta=0;eta<2;eta++) //以eta方向高斯点计算积分

{

Gausseta=Gauss[eta];

for(zeta=0;zeta<2;zeta++) //以zeta方向高斯点计算积分

{

Gausszeta=Gauss[zeta];

for(i=0;i<3;i++) //以坐标进行循环

{

if (i==0)

{

for(j=0;j<8;j++) //以节点进行循环

{

杆刚度校核c语言程序,c语言求解单元刚度矩阵相关推荐

棱形旋转c语言程序_C 语言时隔 5 年重回巅峰，这 20 个热门项目拿去练手！
在上个月的 TIOBE 编程语言排名中,C 语言和 Java 的差距只有 0.01%.在近日 TIOBE 公布的 2020 年 5 月编程语言排行榜中,C 语言成功超越了 Java,重返第一的王者宝座 ...
蜂鸣器发出7种音阶c语言程序_C语言编程新思路
第一章单元测试 1.单选题: 关于一个C语言程序执行的起点和终点,以下选项正确的是 ( ). 选项: A: main 任意 B: main 最后一个 C: main main D: 第一个最后一个 ...
文件的记录c语言程序,c语言程序学生籍贯信息记录簿设计.docx
c 语言程序学生籍贯信息记录簿设计学生籍贯信息记录簿课程设计报告书班级: 方 0909-1 学号:姓名: 苑小叶指导教师 : 康亚男石家庄铁道大学四方学院 2010年 07月 ...
学生实验平台搭建c语言程序,c语言程序设计实验学生用.doc
c语言程序设计实验学生用 C语言程序设计实验指导 (学生用) 计算机基础教研室 <C语言程序设计>课程组 2012年9月前言 <C语言程序设计>是计算机科学技术系面向全校 ...
插入法排序c语言程序,C语言之插入排序算法
一.什么是直接插入算法? 直接插入排序是一种简单的插入排序法,其基本思想是:把待排序的纪录按其关键码值的大小逐个插入到一个已经排好序的有序序列中,直到所有的纪录插入完为止,得到一个新的有序序列. 选择 ...
青花瓷音乐的单片机c语言程序,c语言曲谱_单片机c语言音乐简谱代码
51单片机曲谱编写音乐程序,我看不懂曲谱,谁能教下我! //<世上只有妈妈好>51单片机C语言程序和音乐采灯程序 //此程序在硬件上调试通过 //本程序的单片机晶振采用11．0592MHZ ...
框图c语言程序,C语言程序设计框图
<C语言程序设计框图>由会员分享,可在线阅读,更多相关<C语言程序设计框图(86页珍藏版)>请在人人文库网上搜索. 1.第三章控制结构,返回总目录,目录,3.1节目结构框,3. ...
简单谱子C语言程序,c语言曲谱_单片机c语言音乐简谱代码
51单片机曲谱编写音乐程序,我看不懂曲谱,谁能教下我! //<世上只有妈妈好>51单片机C语言程序和音乐采灯程序 //此程序在硬件上调试通过 //本程序的单片机晶振采用11．0592MHZ ...
用直接分解法求方程组的C语言程序,c语言编程求解线性方程组论文
计算机编程求解线性方程组第一章绪论在自然科学.工程技术.经济和医学各领域中产生的许多实际问题都可以通过数学语言描述为数学问题,也就是说,由实际问题建立数学模型,然后应用各种数学方法和技巧来求解 ...
记得每天锻炼身体c语言程序,c语言程序
一实验名称计算出1000以内10个最大素数之和二.实验目的1.熟练掌握if.if-else.if-else if语句和witch语句格式及使用方法,掌握if语句中的嵌套关系和匹配原则,利用if语句和s ...

杆刚度校核c语言程序,c语言求解单元刚度矩阵

杆刚度校核c语言程序,c语言求解单元刚度矩阵相关推荐

最新文章

热门文章

杆刚度校核c语言程序,c语言 求解单元刚度矩阵

杆刚度校核c语言程序,c语言 求解单元刚度矩阵相关推荐

最新文章

热门文章

杆刚度校核c语言程序,c语言求解单元刚度矩阵

杆刚度校核c语言程序,c语言求解单元刚度矩阵相关推荐