Commit e34a7e96 authored by Mark van Turnhout's avatar Mark van Turnhout
Browse files

supressing history output

parent ab090cea
......@@ -45,8 +45,8 @@ function bas_buildFEM(basil,param)
% For plot parameters and defaults see bas_extractAbaPlot
% param.odbT = 5; (integer) period for output to the abaqus odb-file [s].
% For values <=0, the odb is updated each increment.
% param.indOff = 0; move indenter indOff further from skin
% param.indOff = 0; move indenter indOff [mm] further from skin
%
% See also, bas_prepFEM, bas_MRItoMat, bas_maskBG, bas_segmentMRI,
% bas_mapBone, bas_estCylPos, bas_estIndDisp, bas_extractAbaPlot
......@@ -81,7 +81,7 @@ if size(param.ogden,2) > 3
param.ogden(3) = bas_nu2D(param.ogden(1),param.ogden(4));
end
if ~isfield(param,'odbT'), param.odbT = 5; end % period for out to the abaqus odb-file
if ~isfield(param,'indOff'), param.indOff = 0; end % have abaqus execute the python script (1, default) or not (0)
if ~isfield(param,'indOff'), param.indOff = 0; end % extra offset of indenter from skin [mm]
for b = 1:numel(basil)
......@@ -578,27 +578,7 @@ for b = 1:numel(basil)
fprintf(fid,'myModel.steps[''step-2''].setValues(matrixSolver=DIRECT, matrixStorage=UNSYMMETRIC, solutionTechnique=FULL_NEWTON)\n');
fprintf(fid,'myModel.steps[''step-3''].setValues(matrixSolver=DIRECT, matrixStorage=UNSYMMETRIC, solutionTechnique=FULL_NEWTON)\n');
fprintf(fid,'myModel.steps[''step-4''].setValues(matrixSolver=DIRECT, matrixStorage=UNSYMMETRIC, solutionTechnique=FULL_NEWTON)\n');
% basil(b), diff(bcom(1:2,:),1), bcom(3,:));
%
% fprintf(fid,'myModel.StaticStep(name=''load basil%i'', previous=''Initial'', description=''indent basils leg'', timePeriod=%f, maxNumInc=%i, initialInc=%f, minInc=%f, maxInc=%f, nlgeom=ON)\n',...
% basil(b),param.timePeriod,param.maxNumInc, param.initialInc, param.minInc, param.maxInc);
% if param.fixbone == 0
% fprintf(fid,'myModel.DisplacementBC(name=''moveBone'', createStepName=''load basil%i'', region=region5, u1=%f, u2=%f, u3=%f, ur1=%f, ur2=%f, ur3=%f, amplitude=UNSET, fixed=OFF, distributionType=UNIFORM, fieldName='''', localCsys=None)\n', ...
% basil(b), diff(bcom(1:2,:),1), bcom(3,:));
% else
% fprintf(fid,'myModel.EncastreBC(name=''fixBone'', createStepName=''Initial'', region=region5, localCsys=None)\n');
% end
% fprintf(fid,'region=myModel.rootAssembly.instances[''indenter-1''].sets[''indenterRP'']\n');
% fprintf(fid,'myModel.DisplacementBC(name=''moveIndenter'', createStepName=''load basil%i'', region=region, u1=%f, u2=%f, u3=%f, ur1=0, ur2=0, ur3=0, amplitude=UNSET, fixed=OFF, distributionType=UNIFORM, fieldName='''', localCsys=None)\n', ...
% basil(b), diff(indSE(1:2,:),1));
% if param.cast == 1
% fprintf(fid,'region=myModel.rootAssembly.instances[''cast-1''].sets[''castRP'']\n');
% fprintf(fid,'myModel.EncastreBC(name=''fixCast'', createStepName=''Initial'', region=region, localCsys=None)\n');
% end
% % fprintf(fid,'myModel.steps[''load basil%i''].setValues(reformKernel=4, solutionTechnique=QUASI_NEWTON)\n',basil(b));
if param.discon == 1
fprintf(fid,'myModel.steps[''load basil%i''].control.setValues(allowPropagation=OFF, resetDefaultValues=OFF, discontinuous=ON, lineSearch=(5.0, 1.0, 0.0001, 0.25, 0.1))\n',basil(b));
end
......@@ -611,6 +591,7 @@ for b = 1:numel(basil)
else
fprintf(fid,'myModel.fieldOutputRequests[''F-Output-1''].setValues( variables=(''S'', ''LE'', ''U'', ''RF'', ''COORD'', ''CSTATUS'', ''CFORCE'', ''CDISP'', ''CSTRESS'', ''ELEN'',''ELEDEN'', ''EVOL''))\n');
end
fprintf(fid,'mdb.models[''Basil%i''].historyOutputRequests[''H-Output-1''].suppress()\n',basil(b));
fprintf(fid,'mdb.Job(name=''basil%i'', model=''Basil%i'', description=''Experiment on Basil the Rat'', type=ANALYSIS, numCpus=%i, numDomains=%i)\n', ...
basil(b),basil(b),param.numCpus,param.numDomains);
......
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