平面三角形單元常應變單元matlab程序的編制(編輯修改稿)
2025-08-03 14:03 本頁面
【文章內(nèi)容簡介】 單元的節(jié)點編號 for j=1:3 for k=1:3ASTIF((a(j)*21):a(j)*2,(a(k)*21):a(k)*2)=ASTIF((a(j)*21):a(j)*2,(a(k)*21):a(k)*2)+ESTIF(j*21:j*2,k*21:k*2)。 %跟據(jù)節(jié)點編號對應關系將單元剛度分塊疊加到總剛矩陣中 end end endl 說明：FORMSMATX(i)調(diào)用應力子程序，提取i單元的應力矩陣SMATX；a=LNODS(i,:)記錄i單元的三個結點編號；for…end循環(huán)語句表示行從1到3循環(huán)，列從1到3循環(huán)，將單剛中的元素疊加至總剛中：ASTIF((a(j)*21):a(j)*2,(a(k)*21):a(k)*2)表示總剛中第a(j)*21到:a(j)*2行，第a(k)*21到a(k)*2列的元素由單剛中第j*21到j*2行，第k*21到k*2列的元素疊加而得，a(j)*2即將單元中的位移編碼對應到整體的位移編碼。%%將約束信息加入總剛（置0置1法） NUM=1。 %計數(shù)器(當前已分析的節(jié)點數(shù))NVFIX：受約束邊界點數(shù) i=0。 %計數(shù)器(當前已處理的約束數(shù)) tmp(NVFIX)=0。 %臨時存被約束的序號 while iNVFIX for j=1:0 if FIXED(NUM,j+3)==1 %若發(fā)現(xiàn)約束temporary i=i+1。 %計數(shù)器自增 tmp(i)=FIXED(NUM)*2+j。 %求約束序號 end end NUM=NUM+1。 endl 說明： tmp(NVFIX)=0，形成一個元素值均為0的一行，NVFIX列的行向量，執(zhí)行while語句，首先判斷i是否小于控制數(shù)據(jù)NVFIX，若小于則往下進行，若不小于則退出。執(zhí)行for語句，F(xiàn)IXED(NUM,j+3)表示約束信息數(shù)組中第NUM行，第j+3列的元素，j從1到0，即j+3表示2到3列，即約束信息數(shù)組中描述結點x和y方向受約束的情況，判斷FIXED(NUM,j+3)若等于1，則約束數(shù)自增，若不等于1則跳出。FIXED(NUM)表示FIXED(NUM，1)，tmp(i)=FIXED(NUM)*2+j計算整體約束序號，將序號放入tmp行向量中的i列。% for i=1:NVFIX ASTIF(1:2*NPION,tmp(i))=0。 %將一約束序號處的總剛列向量清0 ASTIF(tmp(i),1:2*NPION)=0。 %.將一約束序號處的總剛行向量清0 ASTIF(tmp(i),tmp(i))=1。 %將行列交叉處的元素置為1 endl 說明：后處理法中置0置1法，設（包括），則將總剛中的主元素K jj換為1，j行和j列的其他元素均改為零。%*******************************************************************% 讀取FORMSMATX子程%******************************************************************* function FORMSMATX(ELEMENT) %計算應力矩陣 %引用所需的全局變量 global NPION NELEM COORD LNODS YOUNG POISS global BMATX DMATX SMATX AREA% %生成彈性矩陣D a=YOUNG/(1POISS^2)。 DMATX(1,1)=1*a。DMATX(1,2)=POISS*a。DMATX(2,1)=POISS*a。DMATX(2,2)=1*a。DMATX(3,3)=(1POISS)*a/2。l 說明： 平面應力問題的彈性矩陣,其中，E為彈性模量，為泊松比。% i=ELEMENT。 %i為當前所計算的單元號 %計算當前單元的面積 AREA=det([1 COORD(LNODS(i,1),1) COORD(LNODS(i,1),2)。... 1 COORD(LNODS(i,2),1) COORD(LNODS(i,2),2)。... 1 COORD(LNODS(i,3),1)
環(huán)評公示相關推薦
鄂ICP備17016276號-1