有限元程序設(shè)計改

1、桿系結(jié)構(gòu)有限元程序設(shè)計14265731#include stdafx.h void main()FILE *fp,*wp;int i;float a;float b2;fp=fopen(E:huhut1.txt,r);wp=fopen(E:huhuout.txt,w);for(i=0;i2;i+)fscanf(fp,%fn,&a);bi=a;fprintf(wp,%fn,bi);2for(i=0;iElementsNum;i+)for(j=0;j3;j+)fscanf(fp,%dt,&StartEnd_Eleij); 3有限元計算工作流程一、前期準備工作節(jié)點：編號、坐標單元：編號、節(jié)點組成、材

2、料屬性載荷：作用節(jié)點、方向、大小約束：作用節(jié)點、方向二、中間實戰(zhàn)環(huán)節(jié)單元剛度矩陣：直接選用整體剛度矩陣：總自由度編號、對號入座節(jié)點力列向量：對號入座約束處理：去處奇異性4有限元計算工作流程三、求解及后處理解線性方程組得位移代入約束方程得反力代入單元本構(gòu)方程得應力5有限元計算編程要點一、輸入信息及其計算機語言描述節(jié)點：編號、坐標單元：編號、節(jié)點組成、材料屬性載荷：作用節(jié)點、方向、大小約束：作用節(jié)點、方向1、節(jié)點定義一維數(shù)組：例：含有兩層意思6有限元計算編程要點一、輸入信息及其計算機語言描述2、單元定義兩維數(shù)組：i 表示單元編號，j 包含單元節(jié)點數(shù)、節(jié)點號、材料屬性號、截面特性號等例：7有限元計算

3、編程要點一、輸入信息及其計算機語言描述兩維數(shù)組：i 表示載荷編號，j 包含載荷作用節(jié)點號、自由度號、大小例：3、載荷定義8有限元計算編程要點一、輸入信息及其計算機語言描述兩維數(shù)組：i 表示約束編號，j 包含約束作用節(jié)點號、自由度號、值例：4、約束定義9有限元計算編程要點1、總自由度自動編號程序二、核心程序總自由度編號整體剛度矩陣組裝節(jié)點力列向量形成約束處理i 表示節(jié)點號，j 位置表示節(jié)點自由度順序，數(shù)值表示總自由度號。10目標函數(shù)(輸出信息): 確定每個節(jié)點的每個自由度對應整體自由度的編號已知條件(輸入信息): 節(jié)點總數(shù)nTotalNode 每個節(jié)點自由度數(shù)總數(shù)nNodeFreedomi 表示節(jié)點號，j 位置表示節(jié)點自由度順序，數(shù)值表示總自由度號。整體編碼數(shù)組充“1”11整體編碼數(shù)組生成122、單元剛度矩陣生成程序13單元信息矩陣元素143、總剛度矩陣生成程序總剛置0單剛裝配單元信息單元節(jié)點總自由度編號154、節(jié)點力向量生成程序總節(jié)點力向量