平面剛架
matlab
程式語言
:matlab
2003.9.16
2007.2.28
2008.4.1 2009.10.10
變數說明
npoin
nelem
nvfix
nfpoin
nfpres
總結點數
單元數約束個數
受力結點數
非結點力數
coord
lnods
young
結構節點座標陣列
單元定義陣列
彈性模量
fpoin
fpres
force
fixed
結點力陣列，非結點力陣列，總體荷載向量
約束資訊陣列
hkdisp
總體剛度矩陣
結點位移向量
format short e
設定輸出型別
clear
清除記憶體變數
fp1=fopen('6-6.txt','rt')
開啟初始資料檔案
讀入控制資料
nelem=fscanf(fp1,'%d',1);
單元數npoin=fscanf(fp1,'%d',1);
結點數nvfix=fscanf(fp1,'%d',1);
約束數nfpoin=fscanf(fp1,'%d',1);
結點荷載數
nfpres=fscanf(fp1,'%d',1);
非結點荷載數
young=fscanf(fp1,'%f',1);
彈性模量
讀取結構資訊
lnods=fscanf(fp1,'%f',[4,nelem])'
單元定義：
左、右結點號，面積，慣性矩(共計
nelem
組)coord=fscanf(fp1,'%f',[2,npoin])'
座標：x,y
座標(共計
npoin
組)fpoin=fscanf(fp1,'%f',[4,nfpoin])'
節點力(共計
nfpoin
組)：結點號、
x方向力(向右正)
% y方向力(向上正)
m力偶(逆時針正)
fpres=fscanf(fp1,'%f',[4,nfpres])'
均布力(共計
% nfpres
組)：單元號、荷載型別、荷載大小、距離左端長度
fixed=fscanf(fp1,'%f',nvfix)'
約束資訊：約束對應的位移編碼(共計
nvfix
組)hk=zeros(3*npoin,3*npoin);
張成總剛矩陣並清零
force=zeros(3*npoin,1);
張成總荷載向量並清零
形成總剛
for i=1:nelem
對單元個數迴圈
生成區域性單剛
區域性座標
右手座標系
ek=ele_ek(i,lnods,coord,young);
t=zbzh(i,lnods,coord);
座標轉換矩陣
空間剛架matlab 基本平面剛架MATLAB程式
平面剛架 matlab 程式 2003.9.16 2007.2.28 2008.4.1 2009.10 2011.10 2013.9 2014.09 變數說明 npoin nelem nvfix nfpoin nfpres 總結點數 單元數約束個數 受力結點數 非結點力數 coord lnods y...
Matlab基本函式
conndef 建立連通矩陣 imbothat 執行bottom hat濾波 imclearborder 抑制邊緣亮結構 imclose 灰度影象或二值影象的閉運算 imdilate 灰度影象或二值影象的膨脹運算 imerode 灰度影象或二值影象的腐蝕運算 imextendedmax 擴充套件極大...
Matlab基本運算
一.基本矩陣指令 reshape a,m,n 將a變成m n型矩陣 rot90 a 將a逆時針旋轉90度 fliplr a 左右翻轉 flipud a 上下翻轉 size a 求矩陣a的尺寸大小，返回值為a的行數和列數，故前面的返回值應為兩個。另外，a a m,選出a的第m行 b a n 選出a的第...