架算例.txt

matlab与有限元算法书的程序

matLab书算例
E=200e6 ;
XYZ1 = [1 ,  0 , 0 , -4] ;
XYZ2 = [2 , -3 , 0 ,  0] ;
XYZ3 = [3 ,  0 , 0 ,  4] ;
XYZ4 = [4 ,  0 , 5 ,  0] ;
K=zeros(12,12) ;
K=SpaceTrussAssemble(K,E,0.001,XYZ1,XYZ4) ;
K=SpaceTrussAssemble(K,E,0.002,XYZ2,XYZ4) ;
K=SpaceTrussAssemble(K,E,0.001,XYZ3,XYZ4) ;
k=K(10:12,10:12) ;
f=[12 0 0] ;
f=f'  ;
u=k\f ;
U=[0 0 0 0 0 0 0 0 0 u'] ;
U= U' ;
F=K*U ;
sigmal1 = SpaceTrussElementStress(E,XYZ1,XYZ4,U)
sigmal2 = SpaceTrussElementStress(E,XYZ2,XYZ4,U)
sigmal3 = SpaceTrussElementStress(E,XYZ3,XYZ4,U)


李老师算例
E=1 ;
XYZ1 = [1 ,  1 , 0.5 , 2] ;
XYZ2 = [2 ,  2 , 0.5 , 2] ;
XYZ3 = [3 ,  0 ,   1 , 0] ;
XYZ4 = [4 ,  0 ,   0 , 0] ;
XYZ5 = [5 ,  3 ,   0 , 0] ;
XYZ6 = [6 ,  3 ,   1 , 0] ;
K=zeros(18,18) ;
K=SpaceTrussAssemble(K,E,1,XYZ1,XYZ2) ; %1
K=SpaceTrussAssemble(K,E,1,XYZ1,XYZ3) ; %2
K=SpaceTrussAssemble(K,E,1,XYZ1,XYZ4) ; %3
K=SpaceTrussAssemble(K,E,1,XYZ2,XYZ5) ; %4
K=SpaceTrussAssemble(K,E,1,XYZ2,XYZ6) ; %5 
K=SpaceTrussAssemble(K,E,1,XYZ1,XYZ6) ; %6  
K=SpaceTrussAssemble(K,E,1,XYZ2,XYZ3) ; %7
K=SpaceTrussAssemble(K,E,1,XYZ2,XYZ4) ; %8 
K=SpaceTrussAssemble(K,E,1,XYZ1,XYZ5) ; %9 
K=SpaceTrussAssemble(K,E,1,XYZ5,XYZ6) ; %10 
K=SpaceTrussAssemble(K,E,1,XYZ3,XYZ6) ; %11
K=SpaceTrussAssemble(K,E,1,XYZ3,XYZ4) ; %12
K=SpaceTrussAssemble(K,E,1,XYZ4,XYZ5) ; %13
k=K(1:6,1:6) ;
f=[10  0 -20  0 -10  0] ;
f=f' ;
u=k\f ;
U=[ u' 0 0 0 0 0 0 0 0 0 0 0 0 ] ;
U= U' ;
F=K*U ;
sigmal1  = SpaceTrussElementStress(E,XYZ1,XYZ2,U)
sigmal2  = SpaceTrussElementStress(E,XYZ1,XYZ3,U)
sigmal3  = SpaceTrussElementStress(E,XYZ1,XYZ4,U)
sigmal4  = SpaceTrussElementStress(E,XYZ2,XYZ5,U)
sigmal5  = SpaceTrussElementStress(E,XYZ2,XYZ6,U)
sigmal6  = SpaceTrussElementStress(E,XYZ1,XYZ6,U)
sigmal7  = SpaceTrussElementStress(E,XYZ2,XYZ3,U)
sigmal8  = SpaceTrussElementStress(E,XYZ2,XYZ4,U)
sigmal9  = SpaceTrussElementStress(E,XYZ1,XYZ5,U)
sigmal10 = SpaceTrussElementStress(E,XYZ5,XYZ6,U)
sigmal11 = SpaceTrussElementStress(E,XYZ3,XYZ6,U)
sigmal12 = SpaceTrussElementStress(E,XYZ3,XYZ4,U)
sigmal13 = SpaceTrussElementStress(E,XYZ4,XYZ5,U)

小塔算例
E=0.206 ;
XYZ1 = [  1   1   1   0 ] ;
XYZ2 = [  2  -1   1   0 ] ;
XYZ3 = [  3   1  -1   0 ] ; 
XYZ4 = [  4  -1  -1   0 ] ;
XYZ5 = [  5   1   1   2 ] ;
XYZ6 = [  6  -1   1   2 ] ;
XYZ7 = [  7   1  -1   2 ] ;
XYZ8 = [  8  -1  -1   2 ] ;
XYZ9 = [  9   1   0   2 ] ;
XYZ10= [ 10  -1   0   2 ] ;
XYZ11= [ 11   0   1   2 ] ;
XYZ12= [ 12   0  -1   2 ] ;
XYZ13= [ 13   2   2   6 ] ;
XYZ14= [ 14  -2   2   6 ] ;
XYZ15= [ 15   2  -2   6 ] ;
XYZ16= [ 16  -2  -2   6 ] ;
K=zeros(48,48) ;
K = SpaceTrussAssemble( K, E , 235.9, XYZ1 ,  XYZ5)   ; %1
K = SpaceTrussAssemble( K, E , 235.9, XYZ2 ,  XYZ6)   ; %2
K = SpaceTrussAssemble( K, E , 235.9, XYZ3 ,  XYZ7)   ; %3
K = SpaceTrussAssemble( K, E , 235.9, XYZ4 ,  XYZ8)   ; %4
K = SpaceTrussAssemble( K, E , 235.9, XYZ1 ,  XYZ6)   ; %5
K = SpaceTrussAssemble( K, E , 235.9, XYZ2 ,  XYZ5)   ; %6
K = SpaceTrussAssemble( K, E , 235.9, XYZ3 ,  XYZ8)   ; %7
K = SpaceTrussAssemble( K, E , 235.9, XYZ4 ,  XYZ7)   ; %8
K = SpaceTrussAssemble( K, E , 235.9, XYZ1 ,  XYZ7)   ; %9
K = SpaceTrussAssemble( K, E , 235.9, XYZ2 ,  XYZ8)   ; %10
K = SpaceTrussAssemble( K, E , 235.9, XYZ3 ,  XYZ5)   ; %11
K = SpaceTrussAssemble( K, E , 235.9, XYZ4 ,  XYZ6)   ; %12
K = SpaceTrussAssemble( K, E , 235.9, XYZ9 ,  XYZ11)  ; %13
K = SpaceTrussAssemble( K, E , 235.9, XYZ11,  XYZ13)  ; %14
K = SpaceTrussAssemble( K, E , 235.9, XYZ9 ,  XYZ13)  ; %15
K = SpaceTrussAssemble( K, E , 235.9, XYZ5 ,  XYZ11)  ; %16
K = SpaceTrussAssemble( K, E , 235.9, XYZ5 ,  XYZ13)  ; %17
K = SpaceTrussAssemble( K, E , 235.9, XYZ5 ,  XYZ9)   ; %18
K = SpaceTrussAssemble( K, E , 235.9, XYZ10,  XYZ11)  ; %19
K = SpaceTrussAssemble( K, E , 235.9, XYZ11,  XYZ14)  ; %20
K = SpaceTrussAssemble( K, E , 235.9, XYZ10,  XYZ14)  ; %21
K = SpaceTrussAssemble( K, E , 235.9, XYZ6 ,  XYZ10)  ; %22
K = SpaceTrussAssemble( K, E , 235.9, XYZ6 ,  XYZ14)  ; %23
K = SpaceTrussAssemble( K, E , 235.9, XYZ6 ,  XYZ11)  ; %24
K = SpaceTrussAssemble( K, E , 235.9, XYZ9 ,  XYZ12)  ; %25
K = SpaceTrussAssemble( K, E , 235.9, XYZ12,  XYZ15)  ; %26
K = SpaceTrussAssemble( K, E , 235.9, XYZ9 ,  XYZ15)  ; %27
K = SpaceTrussAssemble( K, E , 235.9, XYZ7 ,  XYZ9)   ; %28
K = SpaceTrussAssemble( K, E , 235.9, XYZ7 ,  XYZ15)  ; %29
K = SpaceTrussAssemble( K, E , 235.9, XYZ7 ,  XYZ12)  ; %30
K = SpaceTrussAssemble( K, E , 235.9, XYZ10,  XYZ12)  ; %31
K = SpaceTrussAssemble( K, E , 235.9, XYZ12,  XYZ16)  ; %32
K = SpaceTrussAssemble( K, E , 235.9, XYZ10,  XYZ16)  ; %33
K = SpaceTrussAssemble( K, E , 235.9, XYZ8 ,  XYZ12)  ; %34
K = SpaceTrussAssemble( K, E , 235.9, XYZ8 ,  XYZ16)  ; %35
K = SpaceTrussAssemble( K, E , 235.9, XYZ8 ,  XYZ10)  ; %36
K = SpaceTrussAssemble( K, E , 235.9, XYZ9 ,  XYZ10)  ; %37
k = K(1 : 36 , 1: 36 ) ;
%  0                 1                   2                   3            
%  1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6
f = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6.92 ] ;
f = f'  ;
u = k\f ;
%  0                 1                   2                   3            
%  1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6
U = [u' 0 0 0 0 0 0 0 0 0 0 0 0 ] ; 
U = U' ;
F = K*U ;
sigmal_10_41_01_05  = SpaceTrussElementStress( E ,  XYZ1  ,  XYZ5  , U)    %1
sigmal_11_41_02_06  = SpaceTrussElementStress( E ,  XYZ2  ,  XYZ6  , U)    %2
sigmal_12_42_03_07  = SpaceTrussElementStress( E ,  XYZ3  ,  XYZ7  , U)    %3
sigmal_13_43_04_08  = SpaceTrussElementStress( E ,  XYZ4  ,  XYZ8  , U)    %4
sigmal_10_41_01_06  = SpaceTrussElementStress( E ,  XYZ1  ,  XYZ6  , U)    %5
sigmal_11_40_02_05  = SpaceTrussElementStress( E ,  XYZ2  ,  XYZ5  , U)    %6
sigmal_12_43_03_08  = SpaceTrussElementStress( E ,  XYZ3  ,  XYZ8  , U)    %7
sigmal_13_42_04_07  = SpaceTrussElementStress( E ,  XYZ4  ,  XYZ7  , U)    %8
sigmal_10_42_01_07  = SpaceTrussElementStress( E ,  XYZ1  ,  XYZ7  , U)    %9
sigmal_11_43_02_08  = SpaceTrussElementStress( E ,  XYZ2  ,  XYZ8  , U)    %10
sigmal_12_40_03_05  = SpaceTrussElementStress( E ,  XYZ3  ,  XYZ5  , U)    %11
sigmal_13_41_04_06  = SpaceTrussElementStress( E ,  XYZ4  ,  XYZ6  , U)    %12
sigmal_40_60_05_11  = SpaceTrussElementStress( E ,  XYZ5  ,  XYZ11 , U)    %16
sigmal_41_60_06_11  = SpaceTrussElementStress( E ,  XYZ6  ,  XYZ11 , U)    %24
sigmal_42_62_07_12  = SpaceTrussElementStress( E ,  XYZ7  ,  XYZ12 , U)    %30
sigmal_43_62_08_12  = SpaceTrussElementStress( E ,  XYZ8  ,  XYZ12 , U)    %34
sigmal_40_50_05_09  = SpaceTrussElementStress( E ,  XYZ5  ,  XYZ9  , U)    %18
sigmal_41_51_06_10  = SpaceTrussElementStress( E ,  XYZ6  ,  XYZ10 , U)    %22
sigmal_42_50_07_09  = SpaceTrussElementStress( E ,  XYZ7  ,  XYZ9  , U)    %28
sigmal_43_51_08_10  = SpaceTrussElementStress( E ,  XYZ8  ,  XYZ10 , U)    %36
sigmal_50_60_09_11  = SpaceTrussElementStress( E ,  XYZ9  ,  XYZ11 , U)    %13
sigmal_51_60_10_11  = SpaceTrussElementStress( E ,  XYZ10 ,  XYZ11 , U)    %19
sigmal_50_62_09_12  = SpaceTrussElementStress( E ,  XYZ9  ,  XYZ12 , U)    %25
sigmal_51_62_10_12  = SpaceTrussElementStress( E ,  XYZ10 ,  XYZ12 , U)    %31
sigmal_50_51_09_10  = SpaceTrussElementStress( E ,  XYZ9  ,  XYZ10 , U)    %37
sigmal_60_70_11_13  = SpaceTrussElementStress( E ,  XYZ11 ,  XYZ13 , U)    %14
sigmal_50_70_09_13  = SpaceTrussElementStress( E ,  XYZ9  ,  XYZ13 , U)    %15
sigmal_40_70_05_13  = SpaceTrussElementStress( E ,  XYZ5  ,  XYZ13 , U)    %17
sigmal_60_71_11_14  = SpaceTrussElementStress( E ,  XYZ11 ,  XYZ14 , U)    %20
sigmal_51_71_10_14  = SpaceTrussElementStress( E ,  XYZ10 ,  XYZ14 , U)    %21
sigmal_41_71_06_14  = SpaceTrussElementStress( E ,  XYZ6  ,  XYZ14 , U)    %23
sigmal_62_72_12_15  = SpaceTrussElementStress( E ,  XYZ12 ,  XYZ15 , U)    %26
sigmal_50_72_09_15  = SpaceTrussElementStress( E ,  XYZ9  ,  XYZ15 , U)    %27
sigmal_42_72_07_15  = SpaceTrussElementStress( E ,  XYZ7  ,  XYZ15 , U)    %29
sigmal_62_73_12_16  = SpaceTrussElementStress( E ,  XYZ12 ,  XYZ16 , U)    %32
sigmal_51_73_10_16  = SpaceTrussElementStress( E ,  XYZ10 ,  XYZ16 , U)    %33
sigmal_43_73_08_16  = SpaceTrussElementStress( E ,  XYZ8  ,  XYZ16 , U)    %35