-
大小: 268KB文件類型: .zip金幣: 2下載: 1 次發(fā)布日期: 2021-05-13
- 語(yǔ)言: Matlab
- 標(biāo)簽:
資源簡(jiǎn)介
本有限元軟件是用 matlab 語(yǔ)言不依靠任何現(xiàn)有程序開(kāi)發(fā)的。采用空間三次歐拉梁?jiǎn)卧?每個(gè)結(jié)點(diǎn) 6 個(gè)自由度,可根據(jù)需要實(shí)現(xiàn)桿與桿之間的鉸接和固接。可在全局坐標(biāo)系內(nèi)對(duì)每個(gè)
桿的結(jié)點(diǎn)做 x,y,z 方向的平動(dòng)約束和 x,y,z 方向的轉(zhuǎn)動(dòng)約束。載荷的施加方式支持在結(jié)點(diǎn)上施
加集中力和力矩,支持在單元上施加均勻分布力和力矩,還支持在單元上施加線性分布的力
和力矩。 模型的輸入方式采用數(shù)組輸入。

代碼片段和文件信息
function?position?=?inner_node(?karnode1node2Ue_general?)
point_num=8;
T=TT(kar);
Ue_local=T*Ue_general;
position=zeros(3point_num+2);
L=norm(node1-node2);
for?ii=1:point_num+2
????ksi=-1+(ii-1)*2./(point_num+1);
????N=[
??1/2?-?ksi/2???????????0?????????????????????????0?????????????????????????????????0??????????????????????????0??????????????????????????????????0?ksi/2?+?1/2???????????0???????????????????????????0?????????????????????????????????0???????????????????????????0??????????????????????????????????0;
????????????0?1/2?-?ksi/2?????????????????????????0?????????????????????????????????0??????????????????????????0??????????????????????????????????0???????????0?ksi/2?+?1/2???????????????????????????0?????????????????????????????????0???????????????????????????0??????????????????????????????????0;
????????????0???????????0?ksi^3/4?-?(3*ksi)/4?+?1/2?????????????????????????????????0??????????????????????????0?-(L*(-?ksi^3?+?ksi^2?+?ksi?-?1))/8???????????0???????????0?-?ksi^3/4?+?(3*ksi)/4?+?1/2?????????????????????????????????0???????????????????????????0?-(L*(-?ksi^3?-?ksi^2?+?ksi?+?1))/8;
????????????0???????????0?????????????????????????0?????????(3*ksi^2)/4?-?ksi/2?-?1/4?-(2*((3*ksi^2)/4?-?3/4))/L??????????????????????????????????0???????????0???????????0???????????????????????????0?????????(3*ksi^2)/4?+?ksi/2?-?1/4???(2*((3*ksi^2)/4?-?3/4))/L??????????????????????????????????0;
????????????0???????????0?????????????????????????0?(L*(-?ksi^3?+?ksi^2?+?ksi?-?1))/8??ksi^3/4?-?(3*ksi)/4?+?1/2??????????????????????????????????0???????????0???????????0???????????????????????????0?(L*(-?ksi^3?-?ksi^2?+?ksi?+?1))/8?-?ksi^3/4?+?(3*ksi)/4?+?1/2??????????????????????????????????0;
????????????0???????????0?(2*((3*ksi^2)/4?-?3/4))/L?????????????????????????????????0??????????????????????????0??????????(3*ksi^2)/4?-?ksi/2?-?1/4???????????0???????????0??-(2*((3*ksi^2)/4?-?3/4))/L?????????????????????????????????0???????????????????????????0??????????(3*ksi^2)/4?+?ksi/2?-?1/4];
?
????Ue_inner_local=N*Ue_local;
????Ue_inner_general=T(1:61:6)‘*Ue_inner_local;
????
????position(:ii)=node1+(node2-node1)*(ii-1)./(point_num+1)+...
????[Ue_inner_general(1);Ue_inner_general(3);Ue_inner_general(5)];
end???
end
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????2189??2019-01-26?22:29??finite_element_of_trussing\K_beam.m
?????文件????????2186??2019-01-28?09:45??finite_element_of_trussing\TT.m
?????文件????????2345??2019-01-28?01:07??finite_element_of_trussing\inner_node.m
?????文件?????????219??2019-01-28?09:56??finite_element_of_trussing\ka_r.m
?????文件???????11717??2019-01-28?10:13??finite_element_of_trussing\prepare.m
?????文件????????3620??2019-01-28?01:07??finite_element_of_trussing\solve.m
?????文件??????322484??2019-01-28?10:19??finite_element_of_trussing\三維桿系有限元軟件第一版.pdf
?????目錄???????????0??2019-01-28?10:26??finite_element_of_trussing\
評(píng)論
共有 條評(píng)論