function out = model
%
% isoIndentSpherical.m
%
% Model exported on Nov 25 2023, 11:49 by COMSOL 6.1.0.252.

import com.comsol.model.*
import com.comsol.model.util.*

model = ModelUtil.create('Model');

%model.modelPath(['C:\Users\Worker\Desktop\' native2unicode(hex2dec({'65' 'b0'}), 'unicode')  native2unicode(hex2dec({'5e' 'fa'}), 'unicode')  native2unicode(hex2dec({'65' '87'}), 'unicode')  native2unicode(hex2dec({'4e' 'f6'}), 'unicode')  native2unicode(hex2dec({'59' '39'}), 'unicode') ' (3)']);

model.component.create('comp1', true);

model.component('comp1').geom.create('geom1', 2);
model.component('comp1').geom('geom1').axisymmetric(true);

model.component('comp1').mesh.create('mesh1');

model.component('comp1').physics.create('solid', 'SolidMechanics', 'geom1');

model.study.create('std1');
model.study('std1').create('stat', 'Stationary');
model.study('std1').feature('stat').setSolveFor('/physics/solid', true);

model.param.set('Eb', '100');
model.param.set('nu', '0.3');
model.param.set('a', '50');
model.param.set('hb', '20');
model.param.set('ss', 'hb/100');

model.component('comp1').geom('geom1').create('r1', 'Rectangle');
model.component('comp1').geom('geom1').feature('r1').set('size', {'a' 'hb'});
model.component('comp1').geom('geom1').feature('r1').set('pos', {'0' '-hb/2'});
model.component('comp1').geom('geom1').runPre('fin');
model.component('comp1').geom('geom1').run;

model.component('comp1').pair.create('p1', 'Contact');
model.component('comp1').pair.remove('p1');

model.param.set('wi', '-0.001');
model.param.rename('nu', 'nub');
model.param.set('Gt', '7');
model.param.set('Et', '40');
model.param.set('wi', '-0.01');

model.component('comp1').physics('solid').feature('lemm1').set('SolidModel', 'Orthotropic');
model.component('comp1').physics('solid').feature('lemm1').set('TransverseIsotropic', true);
model.component('comp1').physics('solid').feature('lemm1').set('Evect_mat', 'userdef');
model.component('comp1').physics('solid').feature('lemm1').set('Evect', {'Eb' 'Et'});
model.component('comp1').physics('solid').feature('lemm1').set('Evect_mat', 'from_mat');

model.component('comp1').material.create('mat1', 'Common');
model.component('comp1').material.remove('mat1');

model.component('comp1').physics('solid').feature('lemm1').set('TransverseIsotropic', false);
model.component('comp1').physics('solid').feature('lemm1').set('Evector_mat', 'userdef');
model.component('comp1').physics('solid').feature('lemm1').set('Evector', {'Eb' 'Eb' 'Et'});
model.component('comp1').physics('solid').feature('lemm1').set('nuvector_mat', 'userdef');
model.component('comp1').physics('solid').feature('lemm1').set('nuvector', {'nub' '0' '0'});
model.component('comp1').physics('solid').feature('lemm1').set('Gvector_mat', 'userdef');
model.component('comp1').physics('solid').feature('lemm1').set('Gvector', {'Eb/2/(1-nub)' 'Gt' 'Gt'});
model.component('comp1').physics('solid').create('disp1', 'Displacement1', 1);
model.component('comp1').physics('solid').feature('disp1').selection.set([2]);
model.component('comp1').physics('solid').feature('disp1').setIndex('Direction', true, 2);
model.component('comp1').physics('solid').create('wic1', 'WeakInequalityConstraint', 1);

model.param.set('ai', '5');
model.param.set('par', '0.1');

model.component('comp1').physics('solid').feature('wic1').set('constraintExpression', 'ai-sqrt((R+u)^2+(ai+wi*par-w)^2)');
model.component('comp1').physics('solid').feature('wic1').set('method', 'augmentedLagrangian');
model.component('comp1').physics('solid').feature('wic1').set('automaticQuadrature', false);
model.component('comp1').physics('solid').feature('wic1').selection.set([3]);

model.component('comp1').mesh('mesh1').create('map1', 'Map');
model.component('comp1').mesh('mesh1').feature('map1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').selection.set([1]);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').set('numelem', 100);
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').set('type', 'explicit');
model.component('comp1').mesh('mesh1').feature('map1').feature.remove('dis1');
model.component('comp1').mesh('mesh1').feature('map1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('map1').feature.remove('dis1');
model.component('comp1').mesh('mesh1').feature('size').set('custom', true);
model.component('comp1').mesh('mesh1').feature('size').set('hmax', 'ss');
model.component('comp1').mesh('mesh1').run;

model.param.set('ss', 'hb/50');

model.component('comp1').mesh('mesh1').run;

model.study('std1').feature('stat').set('useparam', true);
model.study('std1').feature('stat').setIndex('pname', 'Eb', 0);
model.study('std1').feature('stat').setIndex('plistarr', '', 0);
model.study('std1').feature('stat').setIndex('punit', '', 0);
model.study('std1').feature('stat').setIndex('pname', 'Eb', 0);
model.study('std1').feature('stat').setIndex('plistarr', '', 0);
model.study('std1').feature('stat').setIndex('punit', '', 0);
model.study('std1').feature('stat').setIndex('pname', 'par', 0);
model.study('std1').feature('stat').setIndex('plistarr', 'range(0,1/nsteps,1)', 0);
model.study('std1').feature('stat').setIndex('plistarr', '(range(0,1/nsteps,1))^3', 0);

model.param.set('nsteps', '100');

model.sol.create('sol1');
model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

 %model.result.create('pg1', 'PlotGroup2D');
 %model.result('pg1').set('data', 'dset1');
 %model.result('pg1').set('defaultPlotID', 'stress');
 %model.result('pg1').label('Stress (solid)');
 %model.result('pg1').set('frametype', 'spatial');
 %model.result('pg1').create('surf1', 'Surface');
 %model.result('pg1').feature('surf1').set('expr', {'solid.mises'});
 %model.result('pg1').feature('surf1').set('threshold', 'manual');
 %model.result('pg1').feature('surf1').set('thresholdvalue', 0.2);
 %model.result('pg1').feature('surf1').set('resolution', 'normal');
 %model.result('pg1').feature('surf1').set('colortable', 'Prism');
 %model.result('pg1').feature('surf1').create('def', 'Deform');
 %model.result('pg1').feature('surf1').feature('def').set('scaleactive', true);
 %model.result('pg1').feature('surf1').feature('def').set('scale', '1');
 %model.result('pg1').feature('surf1').feature('def').set('expr', {'u' 'w'});
 %model.result('pg1').feature('surf1').feature('def').set('descr', 'Displacement field');
 %model.result.dataset.create('dset1solidrev', 'Revolve2D');
 %model.result.dataset('dset1solidrev').set('data', 'dset1');
 %model.result.dataset('dset1solidrev').set('revangle', 225);
 %model.result.dataset('dset1solidrev').set('startangle', -90);
 %model.result.dataset('dset1solidrev').set('hasspacevars', true);
 %model.result.create('pg2', 'PlotGroup3D');
 %model.result('pg2').set('data', 'dset1solidrev');
 %model.result('pg2').set('defaultPlotID', 'stress3D');
 %model.result('pg2').label('Stress, 3D (solid)');
 %model.result('pg2').set('frametype', 'spatial');
 %model.result('pg2').create('surf1', 'Surface');
 %model.result('pg2').feature('surf1').set('expr', {'solid.mises'});
 %model.result('pg2').feature('surf1').set('threshold', 'manual');
 %model.result('pg2').feature('surf1').set('thresholdvalue', 0.2);
 %model.result('pg2').feature('surf1').set('colortable', 'Prism');
 %model.result('pg2').feature('surf1').create('def', 'Deform');
 %model.result.dataset('dset1solidrev').set('hasspacevars', true);
 %model.result('pg2').feature('surf1').feature('def').set('revcoordsys', 'cylindrical');
 %model.result('pg2').feature('surf1').feature('def').set('expr', {'u' 'v' 'w'});
 %model.result('pg2').feature('surf1').feature('def').set('descr', 'Displacement field');
 %model.result('pg2').feature('surf1').feature('def').set('descractive', true);
 %model.result('pg2').feature('surf1').feature('def').set('scaleactive', true);
 %model.result('pg2').feature('surf1').feature('def').set('scale', '1');

 % %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;

model.baseSystem('none');

 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 10);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 100);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').create('surf2', 'Surface');
 %model.result('pg1').feature('surf2').set('coloring', 'uniform');
 %model.result('pg1').feature('surf2').set('color', 'black');
 %model.result('pg1').feature('surf2').set('wireframe', true);
 %model.result('pg1').feature('surf2').create('def1', 'Deform');
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').feature('def1').set('scaleactive', true);
 %model.result('pg1').feature('surf2').feature('def1').set('scale', 100);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 1);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').feature('def1').set('scale', 1);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').set('expr', 'w');
 %model.result('pg1').run;

model.param.set('wi', '-0.1');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.param.set('nsteps', '100');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');
 % %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;

model.component('comp1').variable.create('var1');

% To import content from file, use:
% model.component('comp1').variable('var1').loadFile('FILENAME');
model.component('comp1').variable('var1').set('k11', 'wRR');
model.component('comp1').variable('var1').set('k22', 'wR/R');
model.component('comp1').variable('var1').set('m11', 'Dmb*k11+(Dmb-Dgb)*k22');
model.component('comp1').variable('var1').set('m22', 'Dmb*k22+(Dmb-Dgb)*k11');

 %model.result.create('pg3', 'PlotGroup1D');
 %model.result('pg3').run;
 %model.result('pg3').create('ptgr1', 'PointGraph');
 %model.result('pg3').feature('ptgr1').set('markerpos', 'datapoints');
 %model.result('pg3').feature('ptgr1').set('linewidth', 'preference');
 %model.result.dataset.create('arr1', 'Array1D');
 %model.result.dataset('arr1').set('inputmethod', 'manual');
 %model.result.dataset.remove('arr1');
 %model.result.dataset.create('edg1', 'Edge2D');
 %model.result.dataset('edg1').selection.set([3]);
 %model.result.dataset.remove('edg1');
 %model.result.dataset.create('int1', 'Integral');
 %model.result.dataset('int1').set('intsurface', true);
 %model.result.dataset('int1').set('intvolume', true);
 %model.result.dataset('int1').set('level', 'line');
 %model.result('pg3').run;
 %model.result.dataset.remove('int1');

model.component('comp1').cpl.create('intop1', 'Integration');
model.component('comp1').cpl('intop1').set('axisym', true);
model.component('comp1').cpl('intop1').selection.geom('geom1', 1);
model.component('comp1').cpl('intop1').selection.set([3]);

model.sol('sol1').updateSolution;

 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result('pg3').set('xlabelactive', true);
 %model.result('pg3').set('xlabel', 'wi');
 %model.result('pg3').set('ylabelactive', true);
 %model.result('pg3').set('ylabel', 'force');
 %model.result('pg3').run;
 %model.result('pg3').feature('ptgr1').selection.all;
 %model.result('pg3').run;
 %model.result.remove('pg3');
 %model.result('pg2').run;
 %model.result.param.clear;
 %model.result.create('pg3', 'PlotGroup1D');
 %model.result('pg3').run;
 %model.result.table.create('tbl1', 'Table');
 %model.result.table.remove('tbl1');
 %model.result.numerical.create('gpev1', 'EvalPointMatrix');
 %model.result.numerical.remove('gpev1');
 %model.result.numerical.create('int1', 'IntLine');
 %model.result.numerical('int1').set('intsurface', true);
 %model.result.numerical('int1').selection.set([3]);
 %model.result.numerical('int1').set('expr', {'solid.wic1.fc'});
 %model.result.numerical('int1').set('descr', {'Constraint force'});
 %model.result.numerical('int1').set('unit', {''});
 %model.result.table.create('tbl1', 'Table');
 %model.result.table('tbl1').comments('Line Integration 1');
 %model.result.numerical('int1').set('table', 'tbl1');
 %model.result.numerical('int1').setResult;
 %model.result('pg3').run;
 %model.result('pg3').create('tblp1', 'Table');
 %model.result('pg3').feature('tblp1').set('markerpos', 'datapoints');
 %model.result('pg3').feature('tblp1').set('linewidth', 'preference');
 %model.result('pg3').feature('tblp1').set('xaxisdata', 1);
 %model.result('pg3').feature('tblp1').set('plotcolumninput', 'manual');
 %model.result('pg3').feature('tblp1').set('plotcolumns', [2]);
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp1').set('linemarker', 'point');
 %model.result('pg1').run;
 %model.result.table('tbl1').setIndex('headers', 'par*', 0, 1);
 %model.result.table('tbl1').setIndex('headers', '-par*wi', 0, 1);
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result.table('tbl1').setIndex('headers', 'par', 0, 1);
 %model.result.table('tbl1').setIndex('headers', '-par*wi', 0, 1);
 %model.result.table('tbl1').setIndex('headers', '-par*0.1', 0, 1);
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp1').set('rowbased', false);
 %model.result('pg3').run;
 %model.result.numerical('int1').setIndex('expr', 'solid.wic1.fc*2*pi*R', 0);
 %model.result('pg3').run;
 %model.result.numerical('int1').set('table', 'tbl1');
 %model.result.numerical('int1').appendResult;
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp1').set('plotcolumns', [3]);
 %model.result('pg3').run;

model.component('comp1').physics('solid').create('weak1', 'WeakContribution', 2);
model.component('comp1').physics('solid').feature('weak1').selection.set([1]);
model.component('comp1').physics('solid').feature('weak1').set('weakExpression', '-2*pi*R*(test(k11)*m11+test(k22)*m22)');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.param.set('t', '0.3201');
model.param.descr('t', 'thickness per layer');
model.param.set('Dm', '0.14');
model.param.descr('Dm', 'bending modulus');
model.param.set('Dmb', 'Dm/t');
model.param.descr('Dmb', 'inner homo bending modulus');
model.param.set('Dg', '0.097286');
model.param.descr('Dg', '');
model.param.set('Dgb', 'Dg/t');
model.param.descr('Dgb', '');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').physics('solid').prop('ShapeProperty').set('order_displacement', '3s');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');
 %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result('pg2').run;
 %model.result('pg3').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepFirst(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result('pg1').run;
 %model.result('pg1').setIndex('looplevel', 13, 0);
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result.table('tbl1').setIndex('headers', '-(par)^3*0.1', 0, 1);
 %model.result('pg3').run;
 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result.table('tbl1').setIndex('headers', 'par', 0, 1);
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg1').run;

model.label('SHomo_inden_axisy_SolidMechanics_iec.mph');

 %model.result('pg3').run;
 %model.result('pg2').run;
 %model.result('pg3').run;
 %model.result('pg2').run;
 %model.result('pg1').run;
 %model.result('pg1').stepFirst(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result('pg2').run;
 %model.result('pg2').run;
 %model.result('pg2').run;
 %model.result('pg2').run;
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg1').run;
 %model.result('pg1').stepFirst(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepLast(0);
 %model.result('pg1').run;
 %model.result('pg2').run;
 %model.result('pg2').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 10);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 1);
 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').stepFirst(0);
 %model.result('pg1').run;
 %model.result('pg1').setIndex('looplevel', 37, 0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 10);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 100);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 10);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').setIndex('looplevel', 37, 0);
 %model.result('pg1').run;
 %model.result.create('pg4', 'PlotGroup1D');
 %model.result('pg4').run;
 %model.result('pg4').label('1D Plot Group 4--contact pressure');
 %model.result('pg4').create('lngr1', 'LineGraph');
 %model.result('pg4').feature('lngr1').set('markerpos', 'datapoints');
 %model.result('pg4').feature('lngr1').set('linewidth', 'preference');
 %model.result('pg4').feature('lngr1').selection.set([3]);
 %model.result('pg4').feature('lngr1').set('expr', 'solid.wic1.fc');
 %model.result('pg4').feature('lngr1').set('descr', 'Constraint force');
 %model.result('pg4').feature('lngr1').set('xdata', 'expr');
 %model.result('pg4').feature('lngr1').set('xdataexpr', 'R');
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').setIndex('looplevelinput', 'manual', 0);
 %model.result('pg4').setIndex('looplevel', [1], 0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result.table.remove('tbl1');
 %model.result.table.create('tbl1', 'Table');
 %model.result.table('tbl1').comments('Line Integration 1');
 %model.result.numerical('int1').set('table', 'tbl1');
 %model.result.numerical('int1').setResult;
 %model.result('pg3').run;
 %model.result('pg3').create('tblp1', 'Table');
 %model.result('pg3').feature('tblp1').set('markerpos', 'datapoints');
 %model.result('pg3').feature('tblp1').set('linewidth', 'preference');
 %model.result('pg3').feature('tblp1').set('xaxisdata', 1);
 %model.result('pg3').feature('tblp1').set('plotcolumninput', 'manual');
 %model.result('pg3').feature('tblp1').set('plotcolumns', [2]);
 %model.result('pg3').feature('tblp1').set('preprocx', 'none');
 %model.result('pg3').feature('tblp1').set('linemarker', 'point');
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').set('xlabelactive', true);
 %model.result('pg3').set('ylabelactive', true);
 %model.result('pg3').set('ylabel', 'Constraint Force (nN)');
 %model.result('pg3').set('xlabel', 'h/h_{max}');
 %model.result('pg3').run;
 %model.result('pg3').set('xlabel', 'h/hmax');
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg4').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg4').run;
 %model.result('pg4').setIndex('looplevel', [12], 0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg1').run;
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').set('xdataexpr', 'R/0.4');
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').stepFirst(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepLast(0);
 %model.result('pg4').run;
 %model.result('pg4').stepFirst(0);
 %model.result('pg4').run;
 %model.result('pg4').setIndex('looplevel', [33], 0);
 %model.result('pg4').run;
 %model.result('pg4').set('axislimits', true);
 %model.result('pg4').set('xmin', -1);
 %model.result('pg4').set('xmax', 4);
 %model.result('pg4').set('ymin', -0.1);
 %model.result('pg4').set('ymax', '1e-2');
 %model.result('pg4').run;
 %model.result('pg4').set('ymin', '-1e-3');
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepFirst(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').set('xmin', -1);
 %model.result('pg4').set('xmax', 4);
 %model.result('pg4').set('ymin', -0.1);
 %model.result('pg4').set('ymax', 0.5);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').set('ymax', 1);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').set('ymax', 2);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').set('ymax', 5);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;

model.component('comp1').mesh('mesh1').feature.remove('map1');
model.component('comp1').mesh('mesh1').create('fq1', 'FreeQuad');
model.component('comp1').mesh('mesh1').feature('fq1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis1').selection.set([1]);
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis1').set('numelem', 'ceil(hb/(abs(wi)/10))');
model.component('comp1').mesh('mesh1').feature('fq1').create('dis2', 'Distribution');
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis2').selection.set([4]);
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis2').set('numelem', 50);
model.component('comp1').mesh('mesh1').feature('fq1').create('cr1', 'CornerRefinement');
model.component('comp1').mesh('mesh1').feature('fq1').feature('cr1').selection.geom('geom1', 2);
model.component('comp1').mesh('mesh1').feature('fq1').feature('cr1').selection.geom('geom1');
model.component('comp1').mesh('mesh1').feature('fq1').feature.remove('cr1');
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis2').selection.set([3]);
model.component('comp1').mesh('mesh1').feature('size').set('hmin', 0.0015);
model.component('comp1').mesh('mesh1').feature.remove('fq1');
model.component('comp1').mesh.remove('mesh1');
model.component('comp1').mesh.create('mesh1');
model.component('comp1').mesh('mesh1').create('fq1', 'FreeQuad');
model.component('comp1').mesh('mesh1').feature('fq1').create('se1', 'SizeExpression');
model.component('comp1').mesh('mesh1').feature('fq1').feature('se1').set('evaltype', 'initialexpression');
model.component('comp1').mesh('mesh1').feature('fq1').feature.remove('se1');
model.component('comp1').mesh('mesh1').feature('fq1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis1').selection.set([3]);
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis1').set('numelem', 'ceil(hb/(abs(wi)/10))');
model.component('comp1').mesh('mesh1').feature('fq1').feature.duplicate('dis2', 'dis1');
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis2').selection.set([]);
model.component('comp1').mesh('mesh1').feature('fq1').feature.remove('dis2');
model.component('comp1').mesh('mesh1').feature('fq1').feature('dis1').selection.set([1 3]);
model.component('comp1').mesh('mesh1').feature.remove('fq1');
model.component('comp1').mesh('mesh1').create('ftri1', 'FreeTri');
model.component('comp1').mesh('mesh1').feature('ftri1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('ftri1').feature('dis1').selection.set([1 3]);
model.component('comp1').mesh('mesh1').feature('ftri1').feature('dis1').set('numelem', 'ceil(hb/(abs(wi)/10))');
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('ftri1').feature('dis1').selection.set([3]);
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('size').set('hmax', 'ss');
model.component('comp1').mesh('mesh1').run;

model.component('comp1').geom('geom1').run('r1');
model.component('comp1').geom('geom1').create('r2', 'Rectangle');
model.component('comp1').geom('geom1').feature('r2').set('size', {'ai*2' 'ai*2'});
model.component('comp1').geom('geom1').feature('r2').set('pos', {'0' 'hb/2-ai*2'});
model.component('comp1').geom('geom1').run('r2');
model.component('comp1').geom('geom1').run('r2');
model.component('comp1').geom('geom1').create('pard1', 'PartitionDomains');
model.component('comp1').geom('geom1').feature('pard1').selection('domain').set('r1', 1);
model.component('comp1').geom('geom1').feature('pard1').set('partitionwith', 'objects');
model.component('comp1').geom('geom1').feature('pard1').selection('object').set({'r2'});
model.component('comp1').geom('geom1').run('fin');
model.component('comp1').geom('geom1').feature.remove('r2');
model.component('comp1').geom('geom1').run('r1');
model.component('comp1').geom('geom1').feature.remove('pard1');
model.component('comp1').geom('geom1').run('r1');
model.component('comp1').geom('geom1').create('ls1', 'LineSegment');
model.component('comp1').geom('geom1').feature('ls1').set('specify1', 'coord');
model.component('comp1').geom('geom1').feature('ls1').set('coord1', {'0' '5*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls1').set('specify2', 'coord');
model.component('comp1').geom('geom1').feature('ls1').set('coord2', {'a' '5*abs(wi)'});
model.component('comp1').geom('geom1').feature.duplicate('ls2', 'ls1');
model.component('comp1').geom('geom1').feature('ls2').set('coord1', {'5*abs(wi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls2').set('coord2', {'5*abs(wi)' 'hb/2'});
model.component('comp1').geom('geom1').runPre('fin');
model.component('comp1').geom('geom1').feature('ls1').set('coord1', {'0' 'hb/2-5*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls1').set('coord2', {'a' 'hb/2-5*abs(wi)'});
model.component('comp1').geom('geom1').runPre('fin');
model.component('comp1').geom('geom1').feature.duplicate('ls3', 'ls1');
model.component('comp1').geom('geom1').feature.duplicate('ls4', 'ls2');
model.component('comp1').geom('geom1').feature('ls3').set('coord1', {'0' 'hb/2-10*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls3').set('coord2', {'a' 'hb/2-10*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls4').set('coord1', {'10*abs(wi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls4').set('coord2', {'10*abs(wi)' 'hb/2'});
model.component('comp1').geom('geom1').runPre('fin');

model.param.set('wi', '-1');

model.component('comp1').geom('geom1').run('fin');
model.component('comp1').geom('geom1').run('fin');
model.component('comp1').geom('geom1').create('pard1', 'PartitionDomains');
model.component('comp1').geom('geom1').feature('pard1').selection('domain').set('fin', 7);
model.component('comp1').geom('geom1').feature.remove('pard1');

model.component('comp1').mesh('mesh1').feature.remove('ftri1');

model.param.set('ss', 'hb/20');

model.component('comp1').mesh('mesh1').create('map1', 'Map');
model.component('comp1').mesh('mesh1').feature('map1').selection.geom('geom1', 2);
model.component('comp1').mesh('mesh1').feature('map1').selection.set([3 7]);
model.component('comp1').mesh('mesh1').feature('map1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').selection.set([5 6 7 12]);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').set('numelem', 'ceil(abs(wi)/10)');

model.param.set('s1', '5');
model.param.set('s2', '10');

model.component('comp1').geom('geom1').feature('ls1').set('coord1', {'0' 'hb/2-s1*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls1').set('coord2', {'a' 'hb/2-s1*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls2').set('coord1', {'s1*abs(wi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls2').set('coord2', {'s1*abs(wi)' 'hb/2'});
model.component('comp1').geom('geom1').feature('ls3').set('coord1', {'0' 'hb/2-s2*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls3').set('coord2', {'a' 'hb/2-s2*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls4').set('coord1', {'s2*abs(wi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls4').set('coord2', {'s2*abs(wi)' 'hb/2'});
model.component('comp1').geom('geom1').runPre('fin');
model.component('comp1').geom('geom1').run('ls4');
model.component('comp1').geom('geom1').create('pard1', 'PartitionDomains');
model.component('comp1').geom('geom1').feature('pard1').selection('domain').set('r1', 1);
model.component('comp1').geom('geom1').feature('pard1').set('partitionwith', 'objects');
model.component('comp1').geom('geom1').feature('pard1').selection('object').set({'ls1' 'ls2' 'ls3' 'ls4'});
model.component('comp1').geom('geom1').feature('pard1').set('keepobject', false);
model.component('comp1').geom('geom1').run('pard1');
model.component('comp1').geom('geom1').run;
model.component('comp1').geom('geom1').runPre('fin');

model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').set('numelem', 's1*10');
model.component('comp1').mesh('mesh1').feature('map1').feature.duplicate('dis2', 'dis1');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis2').set('numelem', 's2*10');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis2').selection.set([15 16 18 22]);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis2').set('type', 'explicit');
model.component('comp1').mesh('mesh1').feature('map1').feature.remove('dis2');
model.component('comp1').mesh('mesh1').feature('map1').create('size1', 'Size');
model.component('comp1').mesh('mesh1').feature('map1').feature('size1').selection.set([7]);
model.component('comp1').mesh('mesh1').feature('map1').feature('size1').set('custom', true);
model.component('comp1').mesh('mesh1').feature('map1').feature('size1').set('hmaxactive', true);
model.component('comp1').mesh('mesh1').feature('map1').feature('size1').set('hmax', 'ss');
model.component('comp1').mesh('mesh1').create('fq1', 'FreeQuad');
model.component('comp1').mesh('mesh1').feature('fq1').selection.geom('geom1', 2);
model.component('comp1').mesh('mesh1').feature('fq1').selection.set([1 2 4 5 6 9]);
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('fq1').selection.set([1 2 4 5 6 8 9]);
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('fq1').selection.set([2 5 6]);
model.component('comp1').mesh('mesh1').feature('map1').selection.set([1 3 4 7 8 9]);
model.component('comp1').mesh('mesh1').run;

model.study('std1').feature('stat').setIndex('plistarr', '(range(0,1/nsteps,1))', 0);

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.param.rename('s1', 'ss1');
model.param.rename('s2', 'ss2');

model.component('comp1').geom('geom1').feature('ls1').set('coord1', {'0' 'hb/2-ss1*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls1').set('coord2', {'a' 'hb/2-ss1*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls2').set('coord1', {'ss1*abs(wi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls2').set('coord2', {'ss1*abs(wi)' 'hb/2'});
model.component('comp1').geom('geom1').feature('ls3').set('coord1', {'0' 'hb/2-ss2*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls3').set('coord2', {'a' 'hb/2-ss2*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls4').set('coord1', {'ss2*abs(wi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls4').set('coord2', {'ss2*abs(wi)' 'hb/2'});
model.component('comp1').geom('geom1').run('fin');

model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').set('numelem', 'ss1*10');
model.component('comp1').mesh('mesh1').run;

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').material.create('mat1', 'Common');
model.component('comp1').material('mat1').label('Material 1 --iso');
model.component('comp1').material('mat1').propertyGroup.create('Enu', 'Young''s_modulus_and_Poisson''s_ratio');
model.component('comp1').material('mat1').propertyGroup('Enu').set('E', {'200'});
model.component('comp1').material('mat1').propertyGroup('Enu').set('nu', {'0.3'});
model.component('comp1').material('mat1').propertyGroup('def').set('density', {'1e-3'});
model.component('comp1').material.create('mat2', 'Common');
model.component('comp1').material('mat2').propertyGroup.create('OrthotropicVoGrp', 'Orthotropic_Voigt_notation');
model.component('comp1').material('mat2').propertyGroup('OrthotropicVoGrp').set('Evector', {'Eb' 'Eb' 'Et'});
model.component('comp1').material('mat2').propertyGroup('OrthotropicVoGrp').set('nuvector', {'nub' '0' '0'});
model.component('comp1').material('mat2').propertyGroup('OrthotropicVoGrp').set('GvectorVo', {'Eb/2/(1-nub)' 'Gt' 'Gt'});
model.component('comp1').material('mat2').selection.all;
model.component('comp1').material('mat2').propertyGroup('def').set('density', {'1e-3'});
model.component('comp1').material('mat2').label('Material 2 --trans');

model.component('comp1').physics('solid').feature('lemm1').set('Evector_mat', 'from_mat');
model.component('comp1').physics('solid').feature('lemm1').set('nuvector_mat', 'from_mat');
model.component('comp1').physics('solid').feature('lemm1').set('Gvector_mat', 'from_mat');
model.component('comp1').physics('solid').create('lemm2', 'LinearElasticModel', 2);
model.component('comp1').physics('solid').feature('lemm2').selection.all;
model.component('comp1').physics('solid').feature('lemm2').set('SolidModel', 'Orthotropic');
model.component('comp1').physics('solid').feature.remove('lemm2');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').physics('solid').feature('lemm1').set('SolidModel', 'Isotropic');

model.component('comp1').material('mat2').active(false);

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

 %model.result('pg1').run;
 %model.result('pg1').setIndex('looplevel', 1, 0);

model.study('std1').feature('stat').setIndex('plistarr', '(range(0,1/nsteps,1))^3', 0);

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

 %model.result('pg1').run;
 %model.result('pg1').setIndex('looplevel', 8, 0);
 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result.numerical('int1').set('table', 'tbl1');

model.component('comp1').physics('solid').feature('lemm1').set('ForceLinearStrainRes', true);

 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').active(false);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').stepPrevious(0);
 %model.result('pg1').feature('surf2').active(true);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').active(false);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;

model.component('comp1').physics('solid').feature('lemm1').set('E_mat', 'userdef');
model.component('comp1').physics('solid').feature('lemm1').set('E', 200);
model.component('comp1').physics('solid').feature('lemm1').set('nu_mat', 'userdef');
model.component('comp1').physics('solid').feature('lemm1').set('nu', 0.3);
model.component('comp1').physics('solid').feature('lemm1').set('ForceLinearStrainRes', false);

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').mesh('mesh1').feature('size').set('hgrad', 1.2);
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('size').set('hgrad', 1.15);
model.component('comp1').mesh('mesh1').run;

 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg2').run;
 %model.result('pg2').run;
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepPrevious(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepFirst(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepNext(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepNext(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepNext(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepNext(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepNext(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepNext(0);
 %model.result('pg2').run;
 %model.result('pg2').feature('surf1').stepNext(0);
 %model.result('pg2').run;

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').mesh('mesh1').feature('size').set('hgrad', 1.1);
model.component('comp1').mesh('mesh1').run;

model.component('comp1').physics('solid').feature('weak1').active(false);

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').material('mat1').active(false);
model.component('comp1').material('mat2').active(true);

model.component('comp1').physics('solid').feature('lemm1').set('SolidModel', 'Orthotropic');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_wic1_lm'});
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').setIndex('looplevel', [1], 0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').stepPrevious(0);
 %model.result('pg4').run;
 %model.result('pg4').stepNext(0);
 %model.result('pg4').run;

model.param.set('nsteps', '200');

model.component('comp1').material('mat2').active(false);

model.component('comp1').physics('solid').feature('lemm1').set('SolidModel', 'Isotropic');

model.title('SHomo inden axisy SolidMechanics iec');

model.label('SHomo inden axisy SolidMechanics iec');

model.description(['spherical indenter,\ngradient mesh']);

model.component('comp1').physics('solid').feature('wic1').selection.set([7]);

 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;

% To import content from file, use:
% model.param.loadFile('FILENAME');
model.param.descr('a', 'half length');
model.param.set('Gt', '7.677');
model.param.descr('Gt', 'Shear modulus');
model.param.set('Et', '37.3499');
model.param.descr('Et', 'interlayer Young''s modulus');
model.param.set('ai', '0.5');
model.param.descr('ai', 'radius of the indenter');
model.param.set('nsteps', '500');
model.param.set('ss', '0.1');
model.param.descr('ss', 'global seed size');
model.param.set('n', '50', 'num of layers');
model.param.set('ht', '(n-1)*t', 'total thickness');
model.param.set('nu', '0.3051', 'possison''s ratio');
model.param.set('E2D', '247.3221', '2D Young''s modulus');
model.param.set('E2Db', 'E2D/t', 'inner homo Young''s modulus');
model.param.set('hin', '(n-2)*t', 'inner thickness');
model.param.set('ui', '0.01');
model.param.set('sp', '0.2');
model.param.set('ks', '10');
model.param.remove('Eb');
model.param.remove('nub');
model.param.remove('a');
model.param.remove('hb');
model.param.remove('wi');
model.param.remove('Gt');
model.param.remove('Et');
model.param.remove('ai');
model.param.remove('par');
model.param.remove('nsteps');
model.param.remove('t');
model.param.remove('Dm');
model.param.remove('Dmb');
model.param.remove('Dg');
model.param.remove('Dgb');
model.param.remove('ss');
model.param.remove('ss1');
model.param.remove('ss2');
model.param.remove('n');
model.param.remove('ht');
model.param.remove('nu');
model.param.remove('E2D');
model.param.remove('E2Db');
model.param.remove('hin');
model.param.remove('ui');
model.param.remove('sp');
model.param.remove('ks');

% To import content from file, use:
% model.param.loadFile('FILENAME');
model.param.set('a', '50', 'half length');
model.param.set('n', '50', 'num of layers');
model.param.set('t', '0.3201', 'thickness per layer');
model.param.set('ht', '(n-1)*t', 'total thickness');
model.param.set('nu', '0.3051', 'possison''s ratio');
model.param.set('E2D', '247.3221', '2D Young''s modulus');
model.param.set('Gt', '7.677', 'Shear modulus');
model.param.set('Dm', '0.14', 'bending modulus');
model.param.set('Et', '37.3499', 'interlayer Young''s modulus');
model.param.set('E2Db', 'E2D/t', 'inner homo Young''s modulus');
model.param.set('Dmb', 'Dm/t', 'inner homo bending modulus');
model.param.set('Dg', '0.097286');
model.param.set('Dgb', 'Dg/t');
model.param.set('hin', '(n-2)*t', 'inner thickness');
model.param.set('ss', '0.1', 'global seed size');
model.param.set('ai', '0.5', 'radius of the indenter');
model.param.set('ui', '0.01');
model.param.set('nsteps', '500');
model.param.set('sp', '0.2');
model.param.set('ks', '10');
model.param.set('ss1', '5');
model.param.set('ss2', '10');

model.component('comp1').geom('geom1').feature('r1').set('size', {'a' 'ht'});
model.component('comp1').geom('geom1').feature('r1').set('pos', {'0' '-ht/2'});
model.component('comp1').geom('geom1').feature('ls1').set('coord1', {'0' 'ht/2-ss1*abs(wi)'});
model.component('comp1').geom('geom1').feature('ls1').set('coord2', {'a' 'ht/2-ss1*abs(wi)'});

model.param.rename('ui', 'hi');
model.param.set('hi', '1');

model.component('comp1').geom('geom1').run('r1');

model.param.set('nsteps', '100');
model.param.rename('sp', 'par');

model.component('comp1').geom('geom1').feature('ls1').set('coord1', {'0' 'ht/2-ss1*abs(hi)'});
model.component('comp1').geom('geom1').feature('ls1').set('coord2', {'a' 'ht/2-ss1*abs(hi)'});
model.component('comp1').geom('geom1').feature('ls2').set('coord1', {'ss1*abs(hi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls2').setIndex('coord1', '-ht/2', 1);
model.component('comp1').geom('geom1').feature('ls2').set('coord2', {'ss1*abs(hi)' 'hb/2'});
model.component('comp1').geom('geom1').feature('ls2').setIndex('coord2', 'ht/2', 1);
model.component('comp1').geom('geom1').feature('ls3').set('coord1', {'0' 'ht/2-ss2*abs(hi)'});
model.component('comp1').geom('geom1').feature('ls3').set('coord2', {'a' 'ht/2-ss2*abs(hi)'});
model.component('comp1').geom('geom1').feature('ls4').set('coord1', {'ss2*abs(hi)' '-hb/2'});
model.component('comp1').geom('geom1').feature('ls4').setIndex('coord1', '-ht/2', 1);
model.component('comp1').geom('geom1').feature('ls4').set('coord2', {'ss2*abs(hi)' 'hb/2'});
model.component('comp1').geom('geom1').feature('ls4').setIndex('coord2', 'ht/2', 1);
model.component('comp1').geom('geom1').run('fin');

model.component('comp1').material('mat1').active(true);

model.component('comp1').physics('solid').feature('wic1').set('constraintExpression', 'ai-sqrt((R+u)^2+(ai-hi*par-w)^2)');
model.component('comp1').physics('solid').feature('wic1').set('ks', 'ks');

model.component('comp1').mesh('mesh1').run;

model.param.set('ss', 'ht/20');

model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('size').set('hgrad', 1.05);
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('size').set('hgrad', 1.1);
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('size').set('hmin', 'hi/10*0.8');
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('size').set('hmin', 'hi/10');
model.component('comp1').mesh('mesh1').run;

model.component('comp1').physics('solid').feature('wic1').selection.set([7 14]);

model.study('std1').feature('stat').setIndex('pname', 'a', 0);
model.study('std1').feature('stat').setIndex('plistarr', '', 0);
model.study('std1').feature('stat').setIndex('punit', '', 0);
model.study('std1').feature('stat').setIndex('pname', 'a', 0);
model.study('std1').feature('stat').setIndex('plistarr', '', 0);
model.study('std1').feature('stat').setIndex('punit', '', 0);
model.study('std1').feature('stat').setIndex('pname', 'par', 0);
model.study('std1').feature('stat').setIndex('plistarr', 'range(0,1/nsteps,1)', 0);
model.study('std1').feature('stat').setIndex('plistarr', '(range(0,1/nsteps,1))^3', 0);

model.component('comp1').geom('geom1').run('pard1');
model.component('comp1').geom('geom1').create('c1', 'Circle');
model.component('comp1').geom('geom1').feature.move('c1', 0);
model.component('comp1').geom('geom1').feature('c1').set('r', 'ai');
model.component('comp1').geom('geom1').feature('c1').set('angle', 90);
model.component('comp1').geom('geom1').feature('c1').set('pos', {'0' 'ai+ht/2'});
model.component('comp1').geom('geom1').run('c1');
model.component('comp1').geom('geom1').feature('c1').set('rot', -90);
model.component('comp1').geom('geom1').run('c1');
model.component('comp1').geom('geom1').runPre('fin');
model.component('comp1').geom('geom1').run('fin');
model.component('comp1').geom('geom1').feature.move('c1', 6);
model.component('comp1').geom('geom1').feature('fin').set('action', 'assembly');
model.component('comp1').geom('geom1').feature('fin').set('pairtype', 'contact');
model.component('comp1').geom('geom1').feature('fin').set('splitpairs', false);
model.component('comp1').geom('geom1').run('fin');

model.param.set('ai', '5');

model.component('comp1').geom('geom1').run('fin');

model.component('comp1').physics('solid').selection.set([1 2 3 4 5 6 7 8 9 10]);
model.component('comp1').physics('solid').create('lemm2', 'LinearElasticModel', 2);
model.component('comp1').physics('solid').feature('lemm2').selection.set([]);
model.component('comp1').physics('solid').feature('lemm2').selection.all;
model.component('comp1').physics('solid').feature('lemm2').selection.set([1 2 3 4 5 6 7 8 9]);
model.component('comp1').physics('solid').create('rd1', 'RigidDomain', 2);
model.component('comp1').physics('solid').feature('rd1').selection.set([10]);
model.component('comp1').physics('solid').feature('rd1').feature('init1').set('w', 'hi*par');
model.component('comp1').physics('solid').feature.remove('wic1');

model.component('comp1').pair.create('p1', 'Contact');
model.component('comp1').pair('p1').source.set([27]);
model.component('comp1').pair('p1').destination.set([7 14]);

model.component('comp1').physics('solid').create('cnt1', 'SolidContact', 1);
model.component('comp1').physics('solid').feature('cnt1').set('pairs', {'p1'});
model.component('comp1').physics('solid').feature('cnt1').set('ContactMethodCtrl', 'AugmentedLagrange');

model.component('comp1').mesh('mesh1').run;

model.param.set('a', '25');

model.component('comp1').geom('geom1').run('fin');

model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').run;

model.component('comp1').physics('solid').feature('rd1').feature('init1').set('w', '-hi*par');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').physics('solid').feature('rd1').feature('init1').set('w', 0);
model.component('comp1').physics('solid').create('disp2', 'Displacement1', 1);
model.component('comp1').physics('solid').feature.remove('disp2');
model.component('comp1').physics('solid').create('disp2', 'Displacement0', 0);
model.component('comp1').physics('solid').feature.remove('disp2');
model.component('comp1').physics('solid').create('disp2', 'Displacement2', 2);
model.component('comp1').physics('solid').feature('rd1').create('pdr1', 'PrescribedDispRot', -1);
model.component('comp1').physics('solid').feature('rd1').feature('pdr1').set('WeakConstraints', true);
model.component('comp1').physics('solid').feature('rd1').feature('pdr1').set('W0', '-hi*par');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');
 %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').active(true);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf2').stepNext(0);
 %model.result('pg1').run;
 %model.result.numerical('int1').selection.set([7 14]);
 %model.result.table('tbl1').clearTableData;
 %model.result.numerical('int1').set('table', 'tbl1');
 %model.result.numerical('int1').setIndex('expr', '', 0);
 %model.result.numerical.remove('int1');
 %model.result.numerical.create('pev1', 'EvalPoint');
 %model.result.numerical('pev1').set('expr', {'solid.RFz'});
 %model.result.numerical('pev1').set('descr', {'Reaction force, z-component'});
 %model.result.numerical('pev1').set('unit', {''});
 %model.result.table.create('tbl2', 'Table');
 %model.result.table('tbl2').comments('Point Evaluation 1');
 %model.result.numerical('pev1').set('table', 'tbl2');
 %model.result.numerical('pev1').setResult;

model.nodeGroup.create('grp1', 'Results');
model.nodeGroup('grp1').set('type', 'table');
model.nodeGroup('grp1').add('table', 'tbl1');
model.nodeGroup.remove('grp1');

 %model.result.numerical.remove('pev1');
 %model.result.table.remove('tbl2');
 %model.result.numerical.create('gev1', 'EvalGlobal');
 %model.result.numerical('gev1').set('expr', {'solid.RFtotalz'});
 %model.result.numerical('gev1').set('descr', {'Total reaction force, z-component'});
 %model.result.numerical('gev1').set('unit', {''});
 %model.result.table.create('tbl1', 'Table');
 %model.result.table('tbl1').comments('Global Evaluation 1');
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').setResult;
 %model.result.numerical('gev1').remove('unit', 0);
 %model.result.numerical('gev1').remove('descr', 0);
 %model.result.numerical('gev1').remove('expr', [0]);
 %model.result.numerical('gev1').set('descr', {});
 %model.result.numerical('gev1').set('unit', {});
 %model.result.numerical('gev1').set('expr', {'solid.RFtotalr' 'solid.RFtotalphi' 'solid.RFtotalz'});
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').appendResult;
 %model.result.numerical.remove('gev1');
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepLast(0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').selection.set([7 14]);
 %model.result('pg4').feature('lngr1').set('expr', 'solid.Tnz');
 %model.result('pg4').feature('lngr1').set('descr', 'Contact pressure, z-component');
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepLast(0);
 %model.result('pg4').run;

model.component('comp1').physics('solid').feature('cnt1').create('fric1', 'Friction', 1);
model.component('comp1').physics('solid').feature('cnt1').feature('fric1').set('penaltyCtrlAuglag', 'FromParent');
model.component('comp1').physics('solid').feature('cnt1').feature.remove('fric1');
model.component('comp1').physics('solid').feature.remove('disp2');
model.component('comp1').physics('solid').feature('cnt1').create('adh1', 'Adhesion', 1);
model.component('comp1').physics('solid').feature('cnt1').feature.remove('adh1');
model.component('comp1').physics('solid').feature('rd1').feature('pdr1').set('UnidirectionalConstraints', true);

model.component('comp1').mesh('mesh1').create('ftri1', 'FreeTri');
model.component('comp1').mesh('mesh1').feature('ftri1').selection.geom('geom1', 2);
model.component('comp1').mesh('mesh1').feature('ftri1').selection.set([10]);
model.component('comp1').mesh('mesh1').feature('ftri1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('ftri1').feature('dis1').selection.set([27]);
model.component('comp1').mesh('mesh1').feature('ftri1').feature('dis1').set('numelem', 'ss1*10');
model.component('comp1').mesh('mesh1').run('ftri1');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').mesh('mesh1').feature('ftri1').feature('dis1').set('numelem', 'ss1*20');
model.component('comp1').mesh('mesh1').run;

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').physics('solid').feature('lemm2').active(false);

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.component('comp1').physics('solid').feature('lemm2').active(true);

 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepPrevious(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg4').run;
 %model.result('pg1').run;

model.component('comp1').physics('solid').feature('cnt1').set('tunedFor', 'Speed');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');
 %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').feature('surf1').feature('def').set('scale', 1);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').set('expr', 'solid.Tn');
 %model.result('pg4').feature('lngr1').set('descr', 'Contact pressure');
 %model.result('pg4').feature('lngr1').set('xdataexpr', 'R');
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepNext(0);
 %model.result('pg4').run;
 %model.result('pg4').feature('lngr1').stepLast(0);
 %model.result('pg4').run;
 %model.result.table('tbl1').clearTableData;
 %model.result.numerical.create('gev1', 'EvalGlobal');
 %model.result.numerical('gev1').set('descr', {});
 %model.result.numerical('gev1').set('unit', {});
 %model.result.numerical('gev1').set('expr', {'solid.RFtotalr' 'solid.RFtotalphi' 'solid.RFtotalz'});
 %model.result.table.create('tbl2', 'Table');
 %model.result.table('tbl2').comments('Global Evaluation 1');
 %model.result.numerical('gev1').set('table', 'tbl2');
 %model.result.numerical('gev1').setResult;
 %model.result.table.remove('tbl1');
 %model.result.table.remove('tbl2');
 %model.result.table.create('tbl1', 'Table');
 %model.result.table('tbl1').comments('Global Evaluation 1');
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').setResult;
 %model.result('pg3').run;
 %model.result('pg3').create('tblp1', 'Table');
 %model.result('pg3').feature('tblp1').set('markerpos', 'datapoints');
 %model.result('pg3').feature('tblp1').set('linewidth', 'preference');
 %model.result('pg3').feature('tblp1').set('xaxisdata', 1);
 %model.result('pg3').feature('tblp1').set('plotcolumninput', 'manual');
 %model.result('pg3').feature('tblp1').set('plotcolumns', [4]);
 %model.result('pg3').feature('tblp1').set('preprocx', 'linear');
 %model.result('pg3').feature('tblp1').set('scalingx', 'hi');
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').set('xlabel', 'hi (nm)');
 %model.result('pg3').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg3').run;

model.study('std1').feature('stat').setIndex('plistarr', '(range(0,1/nsteps,1))^3 range(1,1/nsteps*5,-1/nsteps*5)', 0);
model.study('std1').feature('stat').setIndex('plistarr', 'range(1,1/nsteps*5,-1/nsteps*5)', 0);

model.func.create('an1', 'Analytic');
model.func.remove('an1');
model.func.create('pw1', 'Piecewise');

model.study('std1').feature('stat').setIndex('plistarr', '(range(0,1/nsteps,1))^3', 0);

model.func('pw1').setIndex('pieces', 0, 0, 0);
model.func('pw1').setIndex('pieces', 1, 0, 1);
model.func('pw1').setIndex('pieces', 'x^3', 0, 2);
model.func('pw1').setIndex('pieces', 1, 1, 0);
model.func('pw1').setIndex('pieces', 2, 1, 1);
model.func('pw1').setIndex('pieces', '2-x', 1, 2);

model.study('std1').feature('stat').setIndex('plistarr', 'pw1({range(0,2/(nsteps-1),2)})', 0);

model.func('pw1').set('smooth', 'cont');

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*32.447882642939895');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.func('pw1').set('smooth', 'none');

model.component('comp1').physics('solid').feature('rd1').feature('pdr1').set('W0', '-hi*pw1(par)');

model.study('std1').feature('stat').setIndex('plistarr', 'range(0,2/(nsteps-1),2)', 0);

 %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').stepFirst(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result('pg1').stepNext(0);
 %model.result('pg1').run;
 %model.result.export.create('anim1', 'Animation');
 %model.result.export('anim1').set('target', 'player');
 %model.result.export('anim1').set('plotgroup', 'pg1');
 %model.result.export('anim1').run;
 %model.result.export('anim1').showFrame;
 %model.result.export('anim1').showFrame;
 %model.result.export('anim1').run;
 %model.result.export('anim1').showFrame;
 %model.result.export('anim1').set('showframe', 1);
 %model.result.export('anim1').set('framesel', 'all');
 %model.result.export('anim1').showFrame;
 %model.result.export('anim1').set('maxframes', 100);
 %model.result.export('anim1').run;
 %model.result.export('anim1').run;
 %model.result.export('anim1').run;
 %model.result.table.remove('tbl1');
 %model.result.table.create('tbl1', 'Table');
 %model.result.table('tbl1').comments('Global Evaluation 1');
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').setResult;
 %model.result('pg4').run;
 %model.result('pg3').run;
 %model.result('pg3').create('tblp1', 'Table');
 %model.result('pg3').feature('tblp1').set('markerpos', 'datapoints');
 %model.result('pg3').feature('tblp1').set('linewidth', 'preference');
 %model.result('pg3').feature('tblp1').set('xaxisdata', 1);
 %model.result('pg3').feature('tblp1').set('plotcolumninput', 'manual');
 %model.result('pg3').feature('tblp1').set('plotcolumns', [4]);
 %model.result('pg3').feature('tblp1').set('preprocx', 'none');
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result.numerical('gev1').setIndex('looplevelinput', 'manual', 0);
 %model.result.numerical('gev1').setIndex('looplevel', [1], 0);
 %model.result.numerical('gev1').setIndex('looplevel', [1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50], 0);

model.study('std1').feature('stat').setIndex('plistarr', 'range(0,2/(nsteps),2)', 0);

 %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result('pg3').run;
 %model.result.table.remove('tbl1');
 %model.result.numerical('gev1').setIndex('looplevel', [1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51], 0);
 %model.result.numerical('gev1').remove('unit', 0);
 %model.result.numerical('gev1').remove('descr', 0);
 %model.result.numerical('gev1').remove('expr', [0]);
 %model.result.numerical('gev1').remove('unit', 0);
 %model.result.numerical('gev1').remove('descr', 0);
 %model.result.numerical('gev1').remove('expr', [0]);
 %model.result.table.create('tbl1', 'Table');
 %model.result.table('tbl1').comments('Global Evaluation 1');
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').setResult;
 %model.result.numerical.duplicate('gev2', 'gev1');
 %model.result.numerical('gev2').setIndex('looplevel', [51], 0);
 %model.result.numerical('gev2').setIndex('looplevel', [51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101], 0);
 %model.result.numerical('gev2').set('expr', {'solid.RFtotalz' 'par'});
 %model.result.numerical('gev2').set('descr', {'Total reaction force, z-component' ''});
 %model.result.numerical('gev2').setIndex('expr', '2-par', 1);
 %model.result.numerical('gev2').set('allowdescrupdate', false);
 %model.result.numerical('gev2').move('expr', [1], -1);
 %model.result.numerical('gev2').move('unit', [1], -1);
 %model.result.numerical('gev2').move('descr', [1], -1);
 %model.result.numerical('gev2').set('allowdescrupdate', true);
 %model.result.numerical('gev1').set('expr', {'solid.RFtotalz' 'par'});
 %model.result.numerical('gev1').set('descr', {'Total reaction force, z-component' ''});
 %model.result.numerical('gev1').setIndex('expr', '(par)^3', 1);
 %model.result.numerical('gev1').set('allowdescrupdate', false);
 %model.result.numerical('gev1').move('expr', [1], -1);
 %model.result.numerical('gev1').move('unit', [1], -1);
 %model.result.numerical('gev1').move('descr', [1], -1);
 %model.result.numerical('gev1').set('allowdescrupdate', true);
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').appendResult;
 %model.result.table('tbl1').clearTableData;
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').setResult;
 %model.result.numerical('gev2').set('table', 'tbl1');
 %model.result.numerical('gev2').appendResult;
 %model.result.table('tbl1').clearTableData;
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').setResult;
 %model.result.table.create('tbl2', 'Table');
 %model.result.numerical('gev2').set('table', 'tbl2');
 %model.result.numerical('gev2').setResult;
 %model.result('pg3').run;
 %model.result('pg3').create('tblp1', 'Table');
 %model.result('pg3').feature('tblp1').set('markerpos', 'datapoints');
 %model.result('pg3').feature('tblp1').set('linewidth', 'preference');
 %model.result('pg3').feature('tblp1').set('xaxisdata', 2);
 %model.result('pg3').feature('tblp1').set('plotcolumninput', 'manual');
 %model.result('pg3').feature('tblp1').set('plotcolumns', [3]);
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').create('tblp2', 'Table');
 %model.result('pg3').feature('tblp2').set('markerpos', 'datapoints');
 %model.result('pg3').feature('tblp2').set('linewidth', 'preference');
 %model.result('pg3').feature('tblp2').set('table', 'tbl2');
 %model.result('pg3').feature('tblp2').set('xaxisdata', 2);
 %model.result('pg3').feature('tblp2').set('plotcolumninput', 'manual');
 %model.result('pg3').feature('tblp2').set('plotcolumns', [3]);
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp2').set('linestyle', 'dashed');
 %model.result('pg3').feature('tblp2').set('linewidth', 3);
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp1').set('linewidth', 3);
 %model.result('pg3').feature('tblp1').set('legendmethod', 'manual');
 %model.result('pg3').feature('tblp1').setIndex('legends', 'Loading', 0);
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp2').set('legendmethod', 'manual');
 %model.result('pg3').feature('tblp2').setIndex('legends', 'Unloading', 0);
 %model.result('pg3').feature('tblp2').set('legend', true);
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp1').set('legend', true);
 %model.result('pg3').run;
 %model.result('pg3').set('legendpos', 'lowerright');
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result.numerical.duplicate('gev3', 'gev1');
 %model.result.numerical('gev1').setIndex('expr', '(par)^3*hi', 0);
 %model.result.numerical('gev2').setIndex('expr', '(2-par)*hi', 0);
 %model.result.numerical('gev3').setIndex('expr', '(par)^3*hi', 0);
 %model.result.numerical('gev3').setIndex('expr', '3*solid.RFtotalz*(1-nu^2)/4/sqrt(ai*((par)^3*hi)^3)', 1);
 %model.result.table.create('tbl3', 'Table');
 %model.result.table('tbl3').comments('Global Evaluation 3');
 %model.result.numerical('gev3').set('table', 'tbl3');
 %model.result.numerical('gev3').setResult;
 %model.result.numerical('gev1').label('Global Evaluation 1--loading');
 %model.result.numerical('gev2').label('Global Evaluation 2--unloading');
 %model.result.numerical('gev3').label('Global Evaluation 3--Young''s Modulus');
 %model.result.table('tbl1').label('Table 1--loading');
 %model.result.table('tbl2').label('Table 2--unloading');
 %model.result.table('tbl3').label('Table 3--Young''s modulus');
 %model.result.table('tbl1').setIndex('headers', 'hi', 1, 1);
 %model.result.table('tbl2').setIndex('headers', 'hi', 1, 1);
 %model.result.table('tbl3').setIndex('headers', 'hi', 1, 1);
 %model.result.table('tbl3').setIndex('headers', 'E', 2, 1);
 %model.result('pg3').run;
 %model.result('pg3').label('1D Plot Group 3--loading/unloading curve');
 %model.result.create('pg5', 'PlotGroup1D');
 %model.result('pg5').run;
 %model.result('pg5').label('1D Plot Group 5--Young''s modulus');
 %model.result('pg5').create('tblp1', 'Table');
 %model.result('pg5').feature('tblp1').set('markerpos', 'datapoints');
 %model.result('pg5').feature('tblp1').set('linewidth', 'preference');
 %model.result('pg5').feature('tblp1').set('table', 'tbl3');
 %model.result('pg5').feature('tblp1').set('xaxisdata', 2);
 %model.result('pg5').feature('tblp1').set('plotcolumninput', 'manual');
 %model.result('pg5').feature('tblp1').set('plotcolumns', [3]);
 %model.result('pg5').run;
 %model.result('pg5').feature('tblp1').set('legendmethod', 'manual');
 %model.result('pg5').feature('tblp1').set('legend', true);
 %model.result('pg5').feature('tblp1').setIndex('legends', 'Infered by Indentation', 0);
 %model.result('pg5').run;
 %model.result('pg5').create('lngr1', 'LineGraph');
 %model.result('pg5').feature('lngr1').set('markerpos', 'datapoints');
 %model.result('pg5').feature('lngr1').set('linewidth', 'preference');
 %model.result('pg5').feature('lngr1').selection.set([7 14]);
 %model.result('pg5').feature('lngr1').set('expr', '200');
 %model.result.numerical('gev3').setIndex('expr', '3*solid.RFtotalz*(1-0.3^2)/4/sqrt(ai*((par)^3*hi)^3)', 1);
 %model.result.table('tbl3').clearTableData;
 %model.result.numerical('gev3').set('table', 'tbl3');
 %model.result.numerical('gev3').setResult;
 %model.result('pg5').run;
 %model.result('pg5').run;
 %model.result('pg5').feature('lngr1').set('xdata', 'expr');
 %model.result('pg5').feature('lngr1').set('xdataexpr', 'R');
 %model.result('pg5').run;
 %model.result('pg5').feature('tblp1').set('xaxisdata', 2);
 %model.result('pg5').feature('tblp1').set('plotcolumns', [3]);
 %model.result('pg5').run;
 %model.result('pg5').run;
 %model.result('pg5').feature('lngr1').set('data', 'dset1');
 %model.result('pg5').feature('lngr1').set('xdataexpr', 'par');
 %model.result('pg5').run;
 %model.result('pg5').feature('lngr1').set('legend', true);
 %model.result('pg5').feature('lngr1').set('legendmethod', 'manual');
 %model.result('pg5').feature.remove('lngr1');
 %model.result('pg5').run;
 %model.result('pg5').run;
 %model.result('pg5').set('xlabelactive', true);
 %model.result('pg5').set('xlabel', 'hi');
 %model.result('pg5').set('ylabelactive', true);
 %model.result('pg5').set('ylabel', 'E (GPa)');
 %model.result('pg5').set('xlabel', 'hi (nm)');
 %model.result('pg4').run;
 %model.result('pg5').run;
 %model.result('pg5').run;
 %model.result('pg4').run;
 %model.result('pg4').run;
 %model.result('pg5').run;
 %model.result('pg5').run;
 %model.result('pg5').run;

model.component('comp1').pair('p1').mapping('initial');

 %model.result('pg1').run;

model.component('comp1').physics('solid').feature('lemm1').set('ForceLinearStrainRes', true);
model.component('comp1').physics('solid').prop('ShapeProperty').set('order_displacement', '2s');
model.component('comp1').physics('solid').feature('lemm2').set('ForceLinearStrainRes', true);

model.component('comp1').geom('geom1').run('c1');
model.component('comp1').geom('geom1').create('ls5', 'LineSegment');
model.component('comp1').geom('geom1').feature.move('ls5', 5);
model.component('comp1').geom('geom1').feature('ls5').set('specify1', 'coord');
model.component('comp1').geom('geom1').feature('ls5').set('coord1', {'0' 'ht/2-t/2'});
model.component('comp1').geom('geom1').feature('ls5').set('specify2', 'coord');
model.component('comp1').geom('geom1').feature('ls5').set('coord2', {'a' 'ht/2-t/2'});
model.component('comp1').geom('geom1').feature.duplicate('ls6', 'ls5');
model.component('comp1').geom('geom1').feature('ls6').set('coord1', {'0' 'ht/2-t'});
model.component('comp1').geom('geom1').feature('ls6').set('coord2', {'a' 'ht/2-t'});
model.component('comp1').geom('geom1').run('ls6');
model.component('comp1').geom('geom1').feature('pard1').selection('object').set({'ls1' 'ls2' 'ls3' 'ls4' 'ls5' 'ls6'});
model.component('comp1').geom('geom1').runPre('fin');

model.param.set('aci', 'sqrt(ai^2-(ai-hi)^2)');
model.param.set('ss1r', '1.1');
model.param.rename('ss1', 'ss1z');
model.param.set('ss1z', '1');
model.param.set('ss2', '2');

model.component('comp1').geom('geom1').feature('ls1').set('coord1', {'0' 'ht/2-ss1z*aci'});
model.component('comp1').geom('geom1').feature('ls1').set('coord2', {'a' 'ht/2-ss1z*aci'});
model.component('comp1').geom('geom1').feature('ls2').set('coord1', {'ss1r*aci' '-ht/2'});
model.component('comp1').geom('geom1').feature('ls2').set('coord2', {'ss1r*aci' 'ht/2'});
model.component('comp1').geom('geom1').feature('ls3').set('coord1', {'0' 'ht/2-ss2*aci'});
model.component('comp1').geom('geom1').feature('ls3').set('coord2', {'a' 'ht/2-ss2*aci'});
model.component('comp1').geom('geom1').feature('ls4').set('coord1', {'ss2*aci' '-ht/2'});
model.component('comp1').geom('geom1').feature('ls4').set('coord2', {'ss2*aci' 'ht/2'});
model.component('comp1').geom('geom1').runPre('fin');
model.component('comp1').geom('geom1').run('c1');
model.component('comp1').geom('geom1').create('sq1', 'Square');
model.component('comp1').geom('geom1').feature('sq1').set('size', 'ai+2');
model.component('comp1').geom('geom1').feature('sq1').set('pos', {'-1' 'ht/2+hi*1.1'});
model.component('comp1').geom('geom1').run('sq1');
model.component('comp1').geom('geom1').create('dif1', 'Difference');
model.component('comp1').geom('geom1').feature('dif1').selection('input').set({'c1'});
model.component('comp1').geom('geom1').feature('dif1').selection('input2').set({'sq1'});
model.component('comp1').geom('geom1').runPre('fin');
model.component('comp1').geom('geom1').run;

model.param.set('sss', 'hi/100');
model.param.set('a', '20');
model.param.set('n', '60');

model.component('comp1').mesh('mesh1').feature('size').set('hmin', 'sss');

model.component('comp1').geom('geom1').run;
model.component('comp1').geom('geom1').run('fin');

model.component('comp1').mesh('mesh1').feature('map1').selection.set([1 3 5 6 11 12 13 14 15]);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').selection.set([10 11]);
model.component('comp1').mesh('mesh1').feature('map1').create('size2', 'Size');
model.component('comp1').mesh('mesh1').feature('map1').feature.move('size2', 0);
model.component('comp1').mesh('mesh1').feature('map1').feature.move('size1', 0);
model.component('comp1').mesh('mesh1').feature('map1').feature('size1').selection.set([5]);
model.component('comp1').mesh('mesh1').feature('map1').feature('size1').set('hmax', 'sss');
model.component('comp1').mesh('mesh1').feature('map1').feature.duplicate('size3', 'size2');
model.component('comp1').mesh('mesh1').feature('map1').feature.move('size3', 2);
model.component('comp1').mesh('mesh1').feature('map1').feature('size2').selection.set([3]);
model.component('comp1').mesh('mesh1').feature('map1').feature('size2').set('custom', true);
model.component('comp1').mesh('mesh1').feature('map1').feature('size2').set('hmaxactive', true);
model.component('comp1').mesh('mesh1').feature('map1').feature('size2').set('hmax', '10*sss');
model.component('comp1').mesh('mesh1').feature('map1').feature('size2').set('hminactive', false);
model.component('comp1').mesh('mesh1').feature('map1').feature.remove('size2');
model.component('comp1').mesh('mesh1').feature('map1').feature.remove('size3');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').selection.set([5 16]);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').set('numelem', 'ceil((ss1z*aci-t)/sss/10)');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').selection.set([5 16]);
model.component('comp1').mesh('mesh1').feature('map1').create('dis2', 'Distribution');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis2').selection.set([6 8]);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis2').set('numelem', 'ceil((ss1r*aci/sss/10))');
model.component('comp1').mesh('mesh1').feature.remove('fq1');
model.component('comp1').mesh('mesh1').feature('ftri1').selection.set([2 4 7 8 9 10 16]);
model.component('comp1').mesh('mesh1').feature.duplicate('ftri2', 'ftri1');
model.component('comp1').mesh('mesh1').feature('ftri1').selection.set([2 4 7 8 9 10]);
model.component('comp1').mesh('mesh1').feature('ftri1').feature.remove('dis1');
model.component('comp1').mesh('mesh1').feature('ftri2').selection.set([16]);
model.component('comp1').mesh('mesh1').create('edg1', 'Edge');
model.component('comp1').mesh('mesh1').feature.move('edg1', 3);
model.component('comp1').mesh('mesh1').feature('edg1').selection.set([41]);
model.component('comp1').mesh('mesh1').feature('edg1').create('size1', 'Size');
model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('custom', true);
model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('hmaxactive', true);
model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('hmax', 'sss*2');
model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('hgradactive', true);
model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('hgrad', 2.5);
model.component('comp1').mesh('mesh1').feature('ftri2').feature.remove('dis1');
model.component('comp1').mesh('mesh1').run;
model.component('comp1').mesh('mesh1').feature('ftri2').create('size1', 'Size');
model.component('comp1').mesh('mesh1').feature('ftri2').feature('size1').set('custom', true);
model.component('comp1').mesh('mesh1').feature('ftri2').feature('size1').set('hgradactive', true);
model.component('comp1').mesh('mesh1').feature('ftri2').feature('size1').set('hgrad', 2.5);
model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('hgradactive', false);
model.component('comp1').mesh('mesh1').run;

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

 %model.result('pg1').run;
 %model.result('pg1').setIndex('looplevel', 51, 0);
 %model.result('pg1').run;
 %model.result('pg1').run;

model.param.remove({'nu' 'E2D' 'Gt' 'Dm' 'Et' 'E2Db' 'Dmb' 'Dg' 'Dgb'});

model.component('comp1').physics('solid').feature.remove('weak1');

 %model.result.table('tbl1').clearTableData;
 %model.result.numerical('gev1').set('table', 'tbl1');
 %model.result.numerical('gev1').setResult;
 %model.result.table('tbl3').clearTableData;
 %model.result.numerical('gev3').set('table', 'tbl3');
 %model.result.numerical('gev3').setResult;
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp1').set('xaxisdata', 2);
 %model.result('pg3').feature('tblp1').set('plotcolumns', [3]);
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp2').active(false);
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp2').set('linemarker', 'point');
 %model.result('pg3').feature('tblp2').active(true);
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp2').active(false);
 %model.result('pg3').run;
 %model.result('pg3').run;
 %model.result('pg3').feature('tblp1').set('linemarker', 'point');
 %model.result('pg4').run;
 %model.result('pg5').run;
 %model.result('pg5').feature('tblp1').set('xaxisdata', 2);
 %model.result('pg5').feature('tblp1').set('plotcolumns', [3]);
 %model.result('pg5').run;
 %model.result('pg5').feature('tblp1').set('linemarker', 'point');
 %model.result('pg5').feature('tblp1').set('linewidth', 3);
 %model.result('pg5').run;
 %model.result('pg5').set('legendpos', 'lowerright');
 %model.result('pg5').run;
 %model.result('pg5').run;

model.label('isoIndentSpherical.mph');

 %model.result('pg5').run;
 %model.result('pg1').run;

model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('hmax', 'sss/4');
model.component('comp1').mesh('mesh1').run;

model.study('std1').feature('stat').setIndex('plistarr', 'range(0,2/(nsteps),1)', 0);

model.sol('sol1').study('std1');

model.study('std1').feature('stat').set('notlistsolnum', 1);
model.study('std1').feature('stat').set('notsolnum', 'auto');
model.study('std1').feature('stat').set('listsolnum', 1);
model.study('std1').feature('stat').set('solnum', 'auto');

model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').create('v1', 'Variables');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol1').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol1').feature('v1').feature('comp1_u').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').create('s1', 'Stationary');
model.sol('sol1').feature('s1').create('p1', 'Parametric');
model.sol('sol1').feature('s1').feature.remove('pDef');
model.sol('sol1').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol1').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').feature('s1').create('se1', 'Segregated');
model.sol('sol1').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol1').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol1').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol1').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol1').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol1').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol1').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');
 %model.sol('sol1').runAll;

 %model.result('pg1').run;
 %model.result.table.create('evl2', 'Table');
 %model.result.table('evl2').comments('Interactive 2D values');
 %model.result.table('evl2').label('Evaluation 2D');
 %model.result.table('evl2').addRow([0.6517249345779419 8.535236358642578 -1], [0 0 0]);

model.component('comp1').mesh('mesh1').feature('edg1').feature('size1').set('hmax', 'sss/2');
model.component('comp1').mesh('mesh1').run;

model.component('comp1').pair('p1').source.set([11 41]);
model.component('comp1').pair('p1').destination.set([11 22 41]);
model.component('comp1').pair('p1').source.set([11 41]);
model.component('comp1').pair('p1').destination.set([11 41]);

 %model.result('pg1').run;

model.study.create('std2');
model.study('std2').create('stat', 'Stationary');
model.study('std2').feature('stat').setSolveFor('/physics/solid', true);
model.study('std2').label('Study 2--');
model.study('std2').label('Study 2--as both source and destination');
model.study('std2').feature('stat').set('useparam', true);
model.study('std2').feature('stat').setIndex('pname', 'a', 0);
model.study('std2').feature('stat').setIndex('plistarr', '', 0);
model.study('std2').feature('stat').setIndex('punit', '', 0);
model.study('std2').feature('stat').setIndex('pname', 'a', 0);
model.study('std2').feature('stat').setIndex('plistarr', '', 0);
model.study('std2').feature('stat').setIndex('punit', '', 0);
model.study('std2').feature('stat').setIndex('pname', 'par', 0);
model.study('std2').feature('stat').setIndex('plistarr', 'range(0,2/(nsteps),1)', 0);

model.sol.create('sol2');
model.sol('sol2').study('std2');
model.sol('sol2').attach('std2');

model.component('comp1').pair.duplicate('p2', 'p1');
model.component('comp1').pair('p1').source.set([41]);
model.component('comp1').pair('p1').destination.set([11]);
model.component('comp1').pair('p2').source.set([11]);
model.component('comp1').pair('p2').destination.set([41]);

model.component('comp1').physics('solid').feature('cnt1').set('pairs', {'p1' 'p2'});

model.sol('sol2').study('std2');
model.sol('sol2').attach('std2');

model.component('comp1').pair.remove('p2');

model.study('std2').label('Study 2--increase scale factor of contact pressure');

model.sol('sol2').study('std2');

model.study('std2').feature('stat').set('notlistsolnum', 1);
model.study('std2').feature('stat').set('notsolnum', 'auto');
model.study('std2').feature('stat').set('listsolnum', 1);
model.study('std2').feature('stat').set('solnum', 'auto');

model.sol('sol2').create('st1', 'StudyStep');
model.sol('sol2').feature('st1').set('study', 'std2');
model.sol('sol2').feature('st1').set('studystep', 'stat');
model.sol('sol2').create('v1', 'Variables');
model.sol('sol2').feature('v1').feature('comp1_solid_Tn_p1').set('scalemethod', 'manual');
model.sol('sol2').feature('v1').feature('comp1_solid_rd1_w').set('scalemethod', 'manual');
model.sol('sol2').feature('v1').feature('comp1_u').set('scalemethod', 'manual');
model.sol('sol2').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '100000000');
model.sol('sol2').feature('v1').feature('comp1_solid_rd1_w').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol2').feature('v1').feature('comp1_u').set('scaleval', '1e-2*28.274302799715503');
model.sol('sol2').feature('v1').set('control', 'stat');
model.sol('sol2').create('s1', 'Stationary');
model.sol('sol2').feature('s1').create('p1', 'Parametric');
model.sol('sol2').feature('s1').feature.remove('pDef');
model.sol('sol2').feature('s1').feature('p1').set('porder', 'constant');
model.sol('sol2').feature('s1').feature('p1').set('control', 'stat');
model.sol('sol2').feature('s1').set('control', 'stat');
model.sol('sol2').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol2').feature('s1').create('se1', 'Segregated');
model.sol('sol2').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol2').feature('s1').feature('se1').create('ss1', 'SegregatedStep');
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('segvar', {'comp1_u' 'comp1_solid_rd1_w' 'comp1_solid_rd1_RFz'});
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('subdtech', 'ddog');
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('subtermauto', 'itertol');
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('subiter', 7);
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('subntolfact', 1);
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('linsolver', 'dDef');
model.sol('sol2').feature('s1').feature('se1').feature('ss1').label('Solid Mechanics');
model.sol('sol2').feature('s1').feature('se1').create('ls1', 'LumpedStep');
model.sol('sol2').feature('s1').feature('se1').feature('ls1').set('segvar', {'comp1_solid_Tn_p1'});
model.sol('sol2').feature('s1').feature('se1').set('maxsegiter', 15);
model.sol('sol2').feature('s1').feature.remove('fcDef');
model.sol('sol2').attach('std2');
model.sol('sol2').feature('s1').feature('dDef').set('linsolver', 'pardiso');
model.sol('sol2').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '1e9');

model.component('comp1').mesh('mesh1').run;

model.study('std2').feature('stat').set('plot', true);

 %model.result.create('pg6', 'PlotGroup2D');
 %model.result('pg6').set('data', 'dset2');
 %model.result('pg6').set('defaultPlotID', 'stress');
 %model.result('pg6').label('Stress (solid) 1');
 %model.result('pg6').set('frametype', 'spatial');
 %model.result('pg6').create('surf1', 'Surface');
 %model.result('pg6').feature('surf1').set('expr', {'solid.mises'});
 %model.result('pg6').feature('surf1').set('threshold', 'manual');
 %model.result('pg6').feature('surf1').set('thresholdvalue', 0.2);
 %model.result('pg6').feature('surf1').set('resolution', 'normal');
 %model.result('pg6').feature('surf1').set('colortable', 'Prism');
 %model.result('pg6').feature('surf1').create('def', 'Deform');
 %model.result('pg6').feature('surf1').feature('def').set('scaleactive', true);
 %model.result('pg6').feature('surf1').feature('def').set('scale', '1');
 %model.result('pg6').feature('surf1').feature('def').set('expr', {'u' 'w'});
 %model.result('pg6').feature('surf1').feature('def').set('descr', 'Displacement field');
 %model.result.dataset.create('dset2solidrev', 'Revolve2D');
 %model.result.dataset('dset2solidrev').set('data', 'dset2');
 %model.result.dataset('dset2solidrev').set('revangle', 225);
 %model.result.dataset('dset2solidrev').set('startangle', -90);
 %model.result.dataset('dset2solidrev').set('hasspacevars', true);
 %model.result.create('pg7', 'PlotGroup3D');
 %model.result('pg7').set('data', 'dset2solidrev');
 %model.result('pg7').set('defaultPlotID', 'stress3D');
 %model.result('pg7').label('Stress, 3D (solid) 1');
 %model.result('pg7').set('frametype', 'spatial');
 %model.result('pg7').create('surf1', 'Surface');
 %model.result('pg7').feature('surf1').set('expr', {'solid.mises'});
 %model.result('pg7').feature('surf1').set('threshold', 'manual');
 %model.result('pg7').feature('surf1').set('thresholdvalue', 0.2);
 %model.result('pg7').feature('surf1').set('colortable', 'Prism');
 %model.result('pg7').feature('surf1').create('def', 'Deform');
 %model.result.dataset('dset2solidrev').set('hasspacevars', true);
 %model.result('pg7').feature('surf1').feature('def').set('revcoordsys', 'cylindrical');
 %model.result('pg7').feature('surf1').feature('def').set('expr', {'u' 'v' 'w'});
 %model.result('pg7').feature('surf1').feature('def').set('descr', 'Displacement field');
 %model.result('pg7').feature('surf1').feature('def').set('descractive', true);
 %model.result('pg7').feature('surf1').feature('def').set('scaleactive', true);
 %model.result('pg7').feature('surf1').feature('def').set('scale', '1');
 %model.result.remove('pg7');
 %model.result.remove('pg6');
 %model.result.dataset.remove('dset2solidrev');
 %model.result('pg1').run;
 %model.result('pg1').set('data', 'dset2');

 %model.sol('sol2').runAll;

 %model.result('pg1').run;

model.sol('sol2').feature('v1').feature('comp1_solid_Tn_p1').set('scaleval', '1e14');
 %model.sol('sol2').runAll;

 %model.result('pg1').run;

out = model;