基于matlab的三维杆系有限元软件

本有限元软件是用 matlab 语言不依靠任何现有程序开发的。采用空间三次欧拉梁单元， 每个结点 6 个自由度，可根据需要实现杆与杆之间的铰接和固接。可在全局坐标系内对每个 杆的结点做 x,y,z 方向的平动约束和 x,y,z 方向的转动约束。载荷的施加方式支持在结点上施 加集中力和力矩，支持在单元上施加均匀分布力和力矩，还支持在单元上施加线性分布的力 和力矩。 模型的输入方式采用数组输入。

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        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\