evaluatePrincipalStrain
Оценка основной деформации в узловых точках
Описание
pStrain = evaluatePrincipalStrain(structuralresults)structuralresults. Для структурных моделей с переходной и частотными характеристиками, evaluatePrincipalStrain оценивает основную деформацию для всех временных или частотных шагов, соответственно.
Примеры
Октаэдрическая деформация сдвига для биметаллического кабеля под натяжением
Решить статическую структурную модель, представляющую биметаллический кабель под натяжением, и вычислить октаэдрическую деформацию сдвига.
Создайте несущую модель.
structuralmodel = createpde('structural','static-solid');
Создайте геометрию и включите ее в модель. Постройте график геометрии.
gm = multicylinder([0.01,0.015],0.05); structuralmodel.Geometry = gm; pdegplot(structuralmodel,'FaceLabels','on', ... 'CellLabels','on', ... 'FaceAlpha',0.5)
Задайте модуль Янга и коэффициент Пуассона для каждого металла.
structuralProperties(structuralmodel,'Cell',1,'YoungsModulus',110E9, ... 'PoissonsRatio',0.28); structuralProperties(structuralmodel,'Cell',2,'YoungsModulus',210E9, ... 'PoissonsRatio',0.3);
Задайте, что грани 1 и 4 являются фиксированными контурами.
structuralBC(structuralmodel,'Face',[1,4],'Constraint','fixed');
Задайте поверхностную тягу для граней 2 и 5.
structuralBoundaryLoad(structuralmodel,'Face',[2,5], ... 'SurfaceTraction',[0;0;100]);
Сгенерируйте mesh и решите проблему.
generateMesh(structuralmodel); structuralresults = solve(structuralmodel)
structuralresults =
StaticStructuralResults with properties:
Displacement: [1x1 FEStruct]
Strain: [1x1 FEStruct]
Stress: [1x1 FEStruct]
VonMisesStress: [22281x1 double]
Mesh: [1x1 FEMesh]
Оцените основную деформацию в узловых местах.
pStrain = evaluatePrincipalStrain(structuralresults);
Используйте основной штамм для оценки первого и второго инвариантов штамма.
I1 = pStrain.e1 + pStrain.e2 + pStrain.e3;
I2 = pStrain.e1.*pStrain.e2 + pStrain.e2.*pStrain.e3 + pStrain.e3.*pStrain.e1;
tauOct = sqrt(2*(I1.^2 -3*I2))/3;
pdeplot3D(structuralmodel,'ColorMapData',tauOct)
Основная деформация для 3-D структурной динамической задачи
Оцените основную деформацию и октаэдрическую деформацию сдвига в луче при гармоническом возбуждении.
Создайте переходную динамическую модель для 3-D задачи.
structuralmodel = createpde('structural','transient-solid');
Создайте геометрию и включите ее в модель. Постройте график геометрии.
gm = multicuboid(0.06,0.005,0.01); structuralmodel.Geometry = gm; pdegplot(structuralmodel,'FaceLabels','on','FaceAlpha',0.5) view(50,20)
Задайте модуль Юнга, отношение Пуассона и массовую плотность материала.
structuralProperties(structuralmodel,'YoungsModulus',210E9, ... 'PoissonsRatio',0.3, ... 'MassDensity',7800);
Закрепите один конец балки.
structuralBC(structuralmodel,'Face',5,'Constraint','fixed');
Применить синусоидальное перемещение по оси Y на конце, противоположном фиксированному концу балки.
structuralBC(structuralmodel,'Face',3,'YDisplacement',1E-4,'Frequency',50);
Сгенерируйте mesh.
generateMesh(structuralmodel,'Hmax',0.01);Задайте нулевое начальное перемещение и скорость.
structuralIC(structuralmodel,'Displacement',[0;0;0],'Velocity',[0;0;0]);
Решить модель.
tlist = 0:0.002:0.2; structuralresults = solve(structuralmodel,tlist);
Оцените основную деформацию в балке.
pStrain = evaluatePrincipalStrain(structuralresults);
Используйте основной штамм для оценки первого и второго инвариантов.
I1 = pStrain.e1 + pStrain.e2 + pStrain.e3; I2 = pStrain.e1.*pStrain.e2 + pStrain.e2.*pStrain.e3 + pStrain.e3.*pStrain.e1;
Используйте инварианты напряжения, чтобы вычислить октаэдрическую деформацию сдвига.
tauOct = sqrt(2*(I1.^2 -3*I2))/3;
Постройте график результатов.
figure
pdeplot3D(structuralmodel,'ColorMapData',tauOct(:,end))
Входные параметры
Выходные аргументы
См. также
evaluatePrincipalStress | evaluateReaction | interpolateDisplacement | interpolateStrain | interpolateStress | interpolateVonMisesStress | StaticStructuralResults | StructuralModel