Commit 36e73b00 authored by matecellab's avatar matecellab
Browse files

updated help of extractAbaPlot

parent ebf9e11a
function bas_extractAbaPlot(basil,param)
% function bas_extractAbaPlot(basil,param)
%
% Writes a python file to be fed at Abaqus (tested on v6.14.1) with
% instructions to build (and subsequently run) a complete FE model:
% geometry, mesh, boundary conditions, contact, and the likes.
% Writes a python file to be fed at Abaqus (tested on v6.14-1) with
% instructions to plot (and subsequently export) results from the Abaqus
% odb-file.
%
% output: basil<basil>_abaPlot.py;
% output: basil<basilid>_abaPlot.py;
%
% and then e.g.:
% system('abaqus cae noGUI=basil<basil>_job.py')
% will run your job/model.
% system('abaqus cae noGUI=basil<basilid>_abaPlot.py')
% will call Abaqus to run the script (and plot and export your images).
%
% As of yet, the following options (= default) are supported:
% param.basilhome = '~/tue/WillekeBasil/'; root folder for
......@@ -22,7 +22,8 @@ function bas_extractAbaPlot(basil,param)
% param.pstate = 1; plot contours on deformed mesh (1, default) or on the
% undeformed mesh (0).
% param.pvar = 'ESEDEN'; Abaqus variable to plot, a case-sensitive string.
% default: 'ESEDEN'
% default: 'ESEDEN' (Element Strain Energy
% DENsity)
% param.pz = []; plot xy-view cuts at z = pz-bcom(3), i.e. define pz
% compared to the mesh contours before the model is offset
% to the bone's centre of mass: -9.5 <= pz <= 9.5 (20 slices).
......@@ -43,7 +44,7 @@ function bas_extractAbaPlot(basil,param)
if nargin == 1, param.default = 1; end
% set defaults
if ~isfield(param,'basilhome'); param.basilhome = bas_getDataRoot; end % root folder for basil<xxxxxx>/data
if ~isfield(param,'basilhome'), param.basilhome = bas_getDataRoot; end % root folder for basil<xxxxxx>/data
if ~isfield(param,'prange'), param.prange = [NaN NaN]; end % set plotting value range [min max], NaN = auto compute
if ~isfield(param,'pcast'), param.pcast = 0; end % plot the cast. 0 = no (default), 1 = yes.
if ~isfield(param,'pind'), param.pind = 1; end % plot the indenter. 0 = no, 1 = yes (default).
......@@ -183,9 +184,10 @@ for b = 1:numel(basil)
for pz = 1:numel(param.pz)
writez(fid,basil(b),param,param.pz(pz)-bcom(1,3),pz);
end
end
end
fclose(fid);
% do we extract the plot NOW?
if param.pexec > 0
system(sprintf('abaqus cae noGUI=basil%i_abaPlot.py',basil(b)));
end
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment