Commit afe3fc09 authored by matecellab's avatar matecellab
Browse files

bug fixing with Willeke

parent d9bbfdc0
......@@ -286,7 +286,7 @@ for b = 1:numel(basil)
fprintf(fid,'\n# assist disc to help partitioning of skin for contact\n');
fprintf(fid,'myModel.Part(name=''assistDisc'', dimensionality=THREE_D, type=ANALYTIC_RIGID_SURFACE)\n');
fprintf(fid,'s1 = myModel.ConstrainedSketch(name=''discSketch'', sheetSize=10.0)\n');
fprintf(fid,'s1.Line(point1=(%1.1f, %f), point2=(%1.1f, %f))\n',param.ri*[4 -1 4.5 1]);
fprintf(fid,'s1.Line(point1=(%1.1f, %f), point2=(%1.1f, %f))\n',param.ri*[4 1 4.5 1]);
fprintf(fid,'s1.ConstructionLine(point1=(0.0, -%1.1f), point2=(0.0, 0.0))\n',param.ri);
fprintf(fid,'p = myModel.parts[''assistDisc'']\n');
fprintf(fid,'p.AnalyticRigidSurfRevolve(sketch=s1)\n');
......
......@@ -66,6 +66,7 @@ for b = 1:numel(basil)
kidx = Rs <= param.ri+1 & Vr(:,3) <= param.ri ;
% find highest vertex, use as end point
D = max(Vr(kidx,3));
D = D(1);
% set indenter end position
indSE(2,1:3) = tip + D*idir';
......@@ -74,11 +75,11 @@ for b = 1:numel(basil)
if param.figs % some graphical feedback...
% load undeformed skin data
load(fullfile(num2str(basil(b)),'contour_bone_preload-ncc.mat'));
load(fullfile(param.basilhome, num2str(basil(b)),'contour_bone_preload-ncc.mat'));
% create patch data
Vbu = bas_cylSurfPatch(Vcs,150);
% load deformed skin data
load(fullfile(num2str(basil(b)),'contour_bone_loading-ncc.mat'));
load(fullfile(param.basilhome, num2str(basil(b)),'contour_bone_loading-ncc.mat'));
% create patch data
Vbd = bas_cylSurfPatch(Vcs,150);
......
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