[26502] | 1 | %Test Name: SquareSheetConstrainedStressMLHO3d
|
---|
| 2 | md=triangle(model(),'../Exp/Square.exp',200000.);
|
---|
| 3 | md=setmask(md,'','');
|
---|
| 4 | md=parameterize(md,'../Par/SquareSheetConstrained.par');
|
---|
| 5 | md=setflowequation(md,'MLHO','all');
|
---|
| 6 | md = extrude(md, 5, 1);
|
---|
| 7 | md.cluster=generic('name',oshostname(),'np',3);
|
---|
| 8 | md.stressbalance.requested_outputs={'default','VxSurface','VySurface','VxShear','VyShear','VxBase','VyBase'};
|
---|
| 9 | md=solve(md,'Stressbalance');
|
---|
| 10 |
|
---|
| 11 | %Fields and tolerances to track changes
|
---|
| 12 | %Fields and tolerances to track changes
|
---|
| 13 | field_names ={'Vx','Vy','Vz','Vel','Pressure','VxSurface','VySurface','VxShear','VyShear','VxBase','VyBase'};
|
---|
| 14 | field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
|
---|
| 15 | field_values={...
|
---|
| 16 | (md.results.StressbalanceSolution.Vx),...
|
---|
| 17 | (md.results.StressbalanceSolution.Vy),...
|
---|
| 18 | (md.results.StressbalanceSolution.Vz),...
|
---|
| 19 | (md.results.StressbalanceSolution.Vel),...
|
---|
| 20 | (md.results.StressbalanceSolution.Pressure),...
|
---|
| 21 | (md.results.StressbalanceSolution.VxSurface),...
|
---|
| 22 | (md.results.StressbalanceSolution.VySurface),...
|
---|
| 23 | (md.results.StressbalanceSolution.VxShear),...
|
---|
| 24 | (md.results.StressbalanceSolution.VyShear),...
|
---|
| 25 | (md.results.StressbalanceSolution.VxBase),...
|
---|
| 26 | (md.results.StressbalanceSolution.VyBase),...
|
---|
| 27 | };
|
---|
| 28 |
|
---|