matlab有限元工具箱国外资源

从国外一个权威fem分析网站上下载的工具箱，希望大家能好好利用

function [SIGsys] = BRICK_stress （in_data resp）



fj = size（in_data.EL）;
SIGsys = zeros（size（in_data.ND1）1）;

for i=1:fj（1） 
   if in_data.EL（i2）==6
      node1 = find（in_data.ND（:1）==in_data.EL（i3））;
      node2 = find（in_data.ND（:1）==in_data.EL（i4））;
      node3 = find（in_data.ND（:1）==in_data.EL（i5））;
      node4 = find（in_data.ND（:1）==in_data.EL（i6））;
      node5 = find（in_data.ND（:1）==in_data.EL（i7））;
      node6 = find（in_data.ND（:1）==in_data.EL（i8））;
      node7 = find（in_data.ND（:1）==in_data.EL（i9））;
      node8 = find（in_data.ND（:1）==in_data.EL（i10））;
      Em    = in_data.EL（i11）;
      miu_1 = in_data.EL（i12）;
      rho_  = in_data.EL（i13）;
      
      [KlcB_lcE_lcV] = D3_LB （in_data.ND（node12）in_data.ND（node13）in_data.ND（node14）...
         in_data.ND（node22）in_data.ND（node23）in_data.ND（node24）in_data.ND（node32）in_data.ND（node33）...
         in_data.ND（node34）in_data.ND（node42）in_data.ND（node43）in_data.ND（node44）in_data.ND（node52）...
         in_data.ND（node53）in_data.ND（node54）in_data.ND（node62）in_data.ND（node63）in_data.ND（node64）...
         in_data.ND（node72）in_data.ND（node73）in_data.ND（node74）in_data.ND（node82）in_data.ND（node83）...
         in_data.ND（node84） Em miu_1）;
      Dlocal = [resp.static.D（node1*3-2） resp.static.D（node1*3-1） resp.static.D（node1*3） ...
              resp.static.D（node2*3-2） resp.static.D（node2*3-1） resp.static.D（node2*3）  ...
              resp.static.D（node3*3-2） resp.static.D（node3*3-1） resp.static.D（node3*3） ...
              resp.static.D（node4*3-2） resp.static.D（node4*3-1） resp.static.D（node4*3） ...
              resp.static.D（node5*3-2） resp.static.D（node5*3-1） resp.static.D（node5*3） ...
              resp.static.D（node6*3-2） resp.static.D（node6*3-1） resp.static.D（node6*3） ...
              resp.static.D（node7*3-2） resp.static.D（node7*3-1） resp.static.D（node7*3） ...
              resp.static.D（node8*3-2） resp.static.D（node8*3-1） resp.static.D（node8*3）];
      
      SIGlocal_c = E_lc*B_lc（::1）*Dlocal‘;
      SIGlocal_1 = E_lc*B_lc（::2）*Dlocal‘;
      SIGlocal_2 = E_lc*B_lc（::3）*Dlocal‘;
      SIGlocal_3 = E_lc*B_lc（::4）*Dlocal‘;
      SIGlocal_4 = E_lc*B_lc（::5）*Dlocal‘;
      SIGlocal_5 = E_lc*B_lc（::6）*Dlocal‘;
      SIGlocal_6 = E_lc*B_lc（::7）*Dlocal‘;
      SIGlocal_7 = E_lc*B_lc（::8）*Dlocal‘;
      SIGlocal_8 = E_lc*B_lc（::9）*Dlocal‘;
      
      SIGsys（node11:6） = SIGlocal_1‘;
      SIGsys（node21:6） = SIGlocal_2‘;
      SIGsys（node31:6） = SIGlocal_3‘;
      SIGsys（node41:6） = SIGlocal_4‘;
      SIGsys（node51:6） = SIGlocal_5‘;
      SIGsys（node61:6） = SIGlocal_6‘;
      SIGsys（node71:6） = SIGlocal_7‘;
      SIGsys（node81:6） = SIGlocal_8‘;
  end;   
end;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2814  2008-09-14 15:04  fem_toolbox\fem\assemble\BRICK_stress.m

     文件        610  2008-09-14 15:05  fem_toolbox\fem\assemble\BRICK_stress_main.m

     文件       1498  2008-09-14 15:07  fem_toolbox\fem\assemble\build_object_KM.m

     文件       2791  2008-09-14 15:09  fem_toolbox\fem\assemble\CST_CSQ_stress.m

     文件       1646  2008-09-14 15:11  fem_toolbox\fem\assemble\DBC_solve.m

     文件        906  2008-09-14 15:12  fem_toolbox\fem\assemble\frame2d_uniform_loads.m

     文件      16586  2008-09-14 15:20  fem_toolbox\fem\assemble\K_assembly.m

     文件        353  2008-09-14 15:21  fem_toolbox\fem\assemble\K_CON.m

     文件       1675  2008-09-14 15:22  fem_toolbox\fem\assemble\K_CST.m

     文件       1403  2008-09-14 15:25  fem_toolbox\fem\assemble\K_pre_allocate.m

     文件        247  2008-09-14 15:25  fem_toolbox\fem\assemble\KM_info.m

     文件        267  2008-09-14 15:26  fem_toolbox\fem\assemble\LOAD_V.m

     文件        476  2008-09-14 15:26  fem_toolbox\fem\assemble\M_2D_beam_assembly.m

     文件        705  2008-09-14 15:27  fem_toolbox\fem\assemble\M_3D_beam_assembly.m

     文件        249  2008-09-14 15:27  fem_toolbox\fem\assemble\matrx_info.m

     文件        556  2008-09-14 15:28  fem_toolbox\fem\assemble\MFC_master_slave.m

     文件        659  2008-09-14 15:29  fem_toolbox\fem\assemble\type_of_elem.m

     文件      12528  2004-10-27 20:03  fem_toolbox\fem\dynamics\bedr.txt

     文件        308  2008-09-14 15:30  fem_toolbox\fem\dynamics\dxdtU.m

     文件        380  2008-09-14 15:31  fem_toolbox\fem\dynamics\dxxdt.m

     文件        376  2008-09-14 15:31  fem_toolbox\fem\dynamics\dxxdt_p.m

     文件       5006  2008-09-14 19:55  fem_toolbox\fem\dynamics\dyn_FEM.m

     文件        600  2008-09-14 15:39  fem_toolbox\fem\dynamics\eigFEM.m

     文件        777  2008-09-14 15:40  fem_toolbox\fem\dynamics\eigFEM_s.m

     文件       1076  2008-09-14 15:41  fem_toolbox\fem\dynamics\FFTsig.m

     文件        803  2008-09-14 15:41  fem_toolbox\fem\dynamics\freq_damp_FEM.m

     文件        312  2008-09-14 15:42  fem_toolbox\fem\dynamics\FT_freq.m

     文件        683  2008-09-14 15:42  fem_toolbox\fem\dynamics\modal_analysis.m

     文件       3098  2008-09-14 15:46  fem_toolbox\fem\dynamics\odeRK6.m

     文件       2002  2006-06-04 08:21  fem_toolbox\fem\dynamics\pulse_a2.txt

............此处省略83个文件信息