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

explicitely fix indenter until it oughta move

parent e34a7e96
......@@ -546,10 +546,14 @@ for b = 1:numel(basil)
fprintf(fid,'\n# create loading steps, BNC''s\n');
fprintf(fid,'import step\n');
% first step, rotate bone over z-axis
fprintf(fid,'# first step, rotate bone over z-axis\n');
fprintf(fid,'Iregion=myModel.rootAssembly.instances[''indenter-1''].sets[''indenterRP'']\n');
fprintf(fid,'Bregion=myModel.rootAssembly.instances[''leg-1''].sets[''boneRP'']\n');
fprintf(fid,'# first step, rotate bone over z-axis\n');
fprintf(fid,'myModel.StaticStep(name=''step-1'', previous=''Initial'', description=''Rotate bone over z-axis'', timePeriod=1, maxNumInc=100, initialInc=1, minInc=1e-6, maxInc=1, nlgeom=ON)\n');
fprintf(fid,'myModel.DisplacementBC(name=''moveBone'', createStepName=''step-1'', region=region5, u1=0, u2=0, u3=0, amplitude=UNSET, fixed=OFF, distributionType=UNIFORM, fieldName='''', localCsys=None)\n');
fprintf(fid,'myModel.VelocityBC(name=''rotateBone'', createStepName=''step-1'', region=region5, vr1=0, vr2=0, vr3=%f, amplitude=UNSET, localCsys=None, distributionType=UNIFORM, fieldName='''')\n',bcom(3,3));
fprintf(fid,'myModel.DisplacementBC(name=''moveIndenter'', createStepName=''step-1'', region=Iregion, u1=0, u2=0, u3=0, amplitude=UNSET, fixed=OFF, distributionType=UNIFORM, fieldName='''', localCsys=None)\n');
fprintf(fid,'myModel.DisplacementBC(name=''moveBone'', createStepName=''step-1'', region=Bregion, u1=0, u2=0, u3=0, amplitude=UNSET, fixed=OFF, distributionType=UNIFORM, fieldName='''', localCsys=None)\n');
fprintf(fid,'myModel.VelocityBC(name=''rotateBone'', createStepName=''step-1'', region=Bregion, vr1=0, vr2=0, vr3=%f, amplitude=UNSET, localCsys=None, distributionType=UNIFORM, fieldName='''')\n',bcom(3,3));
% second step, rotate bone over y-axis
fprintf(fid,'# second step, rotate bone over y-axis\n');
fprintf(fid,'myModel.StaticStep(name=''step-2'', previous=''step-1'', description=''Rotate bone over y-axis'', timePeriod=1, maxNumInc=100, initialInc=1, minInc=1e-6, maxInc=1, nlgeom=ON)\n');
......@@ -563,15 +567,13 @@ for b = 1:numel(basil)
fprintf(fid,'myModel.StaticStep(name=''step-4'', previous=''step-3'', description=''indent basil%is 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);
fprintf(fid,'myModel.boundaryConditions[''rotateBone''].setValuesInStep(stepName=''step-4'', vr1=0)\n');
fprintf(fid,'myModel.boundaryConditions[''moveBone''].setValuesInStep(stepName=''step-4'', u1=%f, u2=%f, u3=%f)\n', diff(bcom(1:2,:),1));
fprintf(fid,'region=myModel.rootAssembly.instances[''indenter-1''].sets[''indenterRP'']\n');
fprintf(fid,'myModel.DisplacementBC(name=''moveIndenter'', createStepName=''step-4'', region=region, u1=%f, u2=%f, u3=%f, ur1=0, ur2=0, ur3=0, amplitude=UNSET, fixed=OFF, distributionType=UNIFORM, fieldName='''', localCsys=None)\n', ...
diff(indSE(1:2,:),1));
fprintf(fid,'myModel.boundaryConditions[''moveBone''].setValuesInStep(stepName=''step-4'', u1=%f, u2=%f, u3=%f)\n', diff(bcom(1:2,:),1));
fprintf(fid,'myModel.boundaryConditions[''moveIndenter''].setValuesInStep(stepName=''step-4'', u1=%f, u2=%f, u3=%f)\n', diff(indSE(1:2,:),1));
% fix cast (in initial step)
if param.cast == 1
fprintf(fid,'# fix cast (in initial step)\n');
fprintf(fid,'region=myModel.rootAssembly.instances[''cast-1''].sets[''castRP'']\n');
fprintf(fid,'myModel.EncastreBC(name=''fixCast'', createStepName=''Initial'', region=region, localCsys=None)\n');
fprintf(fid,'Cregion=myModel.rootAssembly.instances[''cast-1''].sets[''castRP'']\n');
fprintf(fid,'myModel.EncastreBC(name=''fixCast'', createStepName=''Initial'', region=Cregion, localCsys=None)\n');
end
fprintf(fid,'# set solution technique\n');
fprintf(fid,'myModel.steps[''step-1''].setValues(matrixSolver=DIRECT, matrixStorage=UNSYMMETRIC, solutionTechnique=FULL_NEWTON)\n');
......
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