程序示例说明：
利用matlab处理有限元（或光滑有限元）的边界条件施加示例，因为项目测试所需，
需要从剖分好的网格数据文件中，找到所需关系，按照规则，施加边界条件并保存为新的文件以供后续使用。
主要涉及到的有：索引规律迭代，matlab文件读写，etc.
---------------------------------------
Note：好久没写博客了，趁着一个大程序在跑，
等待结果的过程闲来无事，写个小内容如下。
---------------------------------------
直接上代码吧，代码里相应注释比较清晰：
Code 1 ------------- 主程序脚本 TmpScriptforboundary.m -------------------------
%% Tmp_Apply_BoundaryConstraints %% import clear all clc struct = tmp_importfile('coord of point.out'); clear textdata %% --------simple description------------ % third edge % * * * * * * * * * * * * * * % 4th * * 2nd % edg * * edge % * * * * * * * * * * * * * * % 1st edge % #ID x y boundmarker % 1 0 0 1 4 % 2 1 0 1 % 3 2 0 1 2 % 4 2 .5 2 % 5 2 1 2 3 % 6 1 1 3 % 7 0 1 3 4 % 8 0 .5 4 % ---------------------------- %% export % ess_bound.out , the fourth boundary edge ess_val = 0; % ess_boundValue ux = ess_val*cos(rand); uy = ess_val*sin(rand); nodeid = find(data(:,2) == 0); fid = fopen('ess_bound.out', 'wt'); fprintf(fid, '#ess_bound\tID\tnodid\tux\tuy\n'); for i = 1:length(nodeid) %ess_bound(i,:) = [i nodeid(i) ux uy]; fprintf(fid, '%d\t%d\t%d\t%d\n', i, nodeid(i), ux, uy); end fclose(fid); clear ess_val ux uy nodeid fid % tract_bound.out, the second boundary edge tract_val = 1000; % tract_boundValue tx = tract_val*cos(0); ty = tract_val*sin(0); nodeid = find(data(:,2) == max(data(:,2))); fid = fopen('tract_bound.out', 'wt'); fprintf(fid, '#tract_bound\tID\tnodid\ttx\tty\n'); for i = 1:length(nodeid) %tract_bound(i,:) = [i nodeid(i) tx ty]; fprintf(fid, '%d\t%d\t%d\t%d\n', i, nodeid(i), tx, ty); end fclose(fid); clear tract_val tx ty nodeid fid %clear all
Code 2 ------------- 子函数，通过matlab导入数据自动生成的函数 tmp_importfile.m--------------
function [newData1,vars]=tmp_importfile(fileToRead1) %IMPORTFILE(FILETOREAD1) % Imports data from the specified file % FILETOREAD1: file to read % Auto-generated by MATLAB on 27-Apr-2015 15:16:10 % Import the file newData1 = importdata(fileToRead1); % Create new variables in the base workspace from those fields. vars = fieldnames(newData1); for i = 1:length(vars) assignin('base', vars{i}, newData1.(vars{i})); end
东西比较简单，完结；
若为所需，仅供参考。
