1、桿系結(jié)構(gòu)有限元程序設(shè)計(jì)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有限元計(jì)算工作流程一、前期準(zhǔn)備工作節(jié)點(diǎn)：編號(hào)、坐標(biāo)單元：編號(hào)、節(jié)點(diǎn)組成、材

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

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