空間剛架matlab 基本平面剛架MATLAB程式

2021-10-12 15:06:35 字數 1642 閱讀 9048

平面剛架

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

young

結構節點座標陣列

單元定義陣列

彈性模量

afpoin

fpres

force

fixed

結點力陣列,非結點力陣列,總體荷載向量

約束資訊陣列

hkdisp

總體剛度矩陣,結點位移向量

format short e

設定輸出型別

clear

清除記憶體變數

fp1=fopen('xx.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

組):單元號、荷載型別、荷載大小、距離左端長度

荷載型別

1-均布荷載

2-橫向集中力

3-縱向集中力

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);

空間剛架matlab 基本平面剛架程式matlab

平面剛架 matlab 程式語言 matlab 2003.9.16 2007.2.28 2008.4.1 2009.10.10 變數說明 npoin nelem nvfix nfpoin nfpres 總結點數 單元數約束個數 受力結點數 非結點力數 coord lnods young 結構節點座標...

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的第...