Оцените главное напряжение в узловых местоположениях
pStress = evaluatePrincipalStress(structuralresults)structuralresults. Для переходного процесса и частотной характеристики структурные модели, evaluatePrincipalStress оценивает главное напряжение навсегда и шаги частоты, соответственно.
Решите статическую структурную модель, представляющую биметаллический кабель под силой, и вычислите восьмигранное напряжение сдвига.
Создайте структурную модель.
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]
Оцените главное напряжение в узловых местоположениях.
pStress = evaluatePrincipalStress(structuralresults);
Используйте главное напряжение, чтобы оценить первый и второй инвариант напряжения.
I1 = pStress.s1 + pStress.s2 + pStress.s3; I2 = pStress.s1.*pStress.s2 + ... pStress.s2.*pStress.s3 + ... pStress.s3.*pStress.s1; tauOct = sqrt(2*(I1.^2 -3*I2))/3; pdeplot3D(structuralmodel,'ColorMapData',tauOct)
Оцените главное напряжение и восьмигранное напряжение сдвига в луче при гармоническом возбуждении.
Создайте переходную динамическую модель для 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);
Оцените главное напряжение в луче.
pStress = evaluatePrincipalStress(structuralresults);
Используйте главное напряжение, чтобы оценить первые и вторые инварианты.
I1 = pStress.s1 + pStress.s2 + pStress.s3; I2 = pStress.s1.*pStress.s2 + ... pStress.s2.*pStress.s3 + ... pStress.s3.*pStress.s1;
Используйте инварианты напряжения, чтобы вычислить восьмигранное напряжение сдвига.
tauOct = sqrt(2*(I1.^2 -3*I2))/3;
Постройте график результатов.
figure
pdeplot3D(structuralmodel,'ColorMapData',tauOct(:,end))
StructuralModel | StaticStructuralResults | interpolateDisplacement | interpolateStress | interpolateStrain | interpolateVonMisesStress | evaluateReaction | evaluatePrincipalStrain