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))
Входные параметры
Выходные аргументы
Смотрите также
StructuralModel | StaticStructuralResults | interpolateDisplacement | interpolateStress | interpolateStrain | interpolateVonMisesStress | evaluateReaction | evaluatePrincipalStress