%?“三角形單元“求解“平面應力問題“——處理器部分
%?說明：在運行本程序之前，請做好以下準備工作-預處理
%?1、預處理
%?1.0?實常數(shù)定義
%?????彈性模量?EX
%?????泊松比???PRXY
%?????板厚?????t
%?注意：以上各量單位封閉，自成體系????
%?1.1?三角形單元位移模式?u=a1+a2*x+a3*y;
%???????????????????????v=a4+a5*x+a6*y;
%?1.2?彈性體剖分，劃分三角形單元???1、2...T?（T?表示三角形的個數(shù)）
%????????????????對節(jié)點進行標號???1、2...n?（n?表示節(jié)點的個數(shù)）
%????????????????給出節(jié)點坐標信息?Nodes[i1]（i=1..n）?表示第i個節(jié)點的“x“坐標
%????????????????????????????????Nodes[i2]（i=1..n）?表示第i個節(jié)點的“y“坐標
%????????????????寫出單元連接信息?Element（ki）?（k=1..T）?表示第“k“個三角形的“i“節(jié)點
%????????????????????????????????Element（kj）?（k=1..T）?表示第“k“個三角形的“j“節(jié)點
%????????????????????????????????Element（km）?（k=1..T）?表示第“k“個三角形的“m“節(jié)點
%?1.3?外載荷等效節(jié)點載荷?Load
%???????????????????????Load（ixy）（i=1..n）?表示i節(jié)點的等效載荷沿?“x“、“y“方向的分力

%?2、處理器
%?2.1?計算單元剛度矩陣???
%?????（1）、[K]=（66i）（i=1..T）表示第i個三角形的剛度矩陣?
%?????（2）、[K]=t*AREA（i）*[B]‘[D][B]
%??????????其中?[D]為平面應力問題中的物理矩陣
%???????????????[B]為幾何矩陣，可通過下式求出
%?2.2?組裝形成整體剛度矩陣??[k]=（2*n2*n）即大小為2倍的節(jié)點數(shù)的方陣
%?2.3?引入位移邊界條件??
%?????Hold（2*節(jié)點個數(shù)1）=1?表示該位移被限制，?Hold（2*節(jié)點個數(shù)11）=0?表示為自由位移
%?2.4?模型的求解
%?????（1）消元法????即“劃行劃列法“
%?????（2）置大數(shù)法
%?????說明，本程序中采用的?前一種方法

%?3、后處理
%?3.1?求應力?????????Thigma?????最大應力?Max_Thigma
%?3.2?繪制變形后網(wǎng)絡
%?3.3?繪制出應力場圖
%?3.4?導出圖表、網(wǎng)頁等

%界面初始化
clc;
clear;

%加載初始信息
load?Indata.mat;
%?~~EX_彈性模量
%?~~Element_單元連接信息
%?~~Hold_節(jié)點位移約束
%?~~Nodes_節(jié)點坐標
%?~~PRXY_泊松比
%?~~t_板厚

%?計算物理矩陣?D（33）;
D=EX/（1-PRXY^2）*[1?PRXY?0;PRXY?1?0;0?0?（1-PRXY）/2];

%?計算每個三角形的面積?AREA
??AREA=zeros（size（Element1）1）;??????%?size（Element2）?測量劃分三角形的個數(shù)
??for?i=1:?size（Element1）
??????AREA（i）=1/2*det（...
????????????[?1Nodes（Element（i1）1）Nodes（Element（i1）2）;?...
??????????????1Nodes（Element（i2）1）Nodes（Element（i2）2）;?...
??????????????1Nodes（Element（i3）1）Nodes（Element（i3）2）]）;?
??end????
%?定義單元剛度矩陣?大小為6*6
??K=zeros（66size（Element1））;
%?定義單元幾何矩陣
??B=zeros（36size（Element1））;
??for?i=1:size（Element1）
??????delta=...
??????[?1Nodes（Element（i1）1）Nodes（Element（i1）2）;?...
????????1Nodes（Element（i2）1）Nodes（Element（i2）2）;?...
????????1Nodes（Element（i3）1）Nodes（Element（i3）2）];?
??????for?j=1:3
??????????%其中?Math（deltaji）表示求delta?矩陣j行i列的代數(shù)余子式
??????????B（12*j-1i）=Math（deltaj2）*1/2/AREA（i）;
??????????B（32*ji）=Math（deltaj2）*1/2/AREA（i）;
??????????B（22*ji）=Math（deltaj3）*1/2/AREA（i）;
??????????B（32*j-1i）=Math（deltaj3）*1/2/AREA（i）;
??????end
??end
%?計算單元剛度矩陣
??for?i=1:size（Element1）
??????K（::i）=t*AREA（i）*B（::i）‘*D*B（::i）;
??end
%?疊加得整體剛度矩陣
Node=0;????????????????????%?標記節(jié)點序號
%?定義整體剛度矩陣
k=zeros（2*size（Nodes1）2*size（Nodes1））;
%單元剛度矩陣的投放
for?i=1:size（Element1）
????for?j=1:3
????Node=Element（ij）;
????k（2*Node-1:2*Node2*Node-1:2*Node）=k（2*Node-1:2*Node2*Node-1:2*Node）+K（2*j-1:2*j2*j-1:2*ji）;
????end
????k（2*Element（i1）-1:2*Element（i1）2*Element（i2）-1:2*Element（i2））=...
????????k（2*Element（i1）-1:2*Element（i1）2*Element（i2）-1:2*Element（i2））+K（1:23:4i）;
????k（2*Element（i1）-1:2*Element（i1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????8275??2010-01-09?19:43??In.txt

?????文件???????1290??2010-01-10?15:43??Indata.mat

?????文件???????5487??2010-01-10?16:15??Main.m

?????文件???????5488??2010-01-10?15:54??Main.asv

?????文件????????289??2009-12-26?12:35??Math.m

?????文件????????230??2009-12-26?11:31??Math.asv

?????文件???????4217??2010-01-09?19:41??Out.txt

?????文件?????142177??2010-01-09?19:28??實驗報告.pdf

-----------?---------??----------?-----??----

???????????????167453????????????????????8