G11-MATH6004CM-prj计算方法大作业

于 2023-02-11 21:00:58 发布
阅读量200 收藏
点赞数
CC 4.0 BY-SA版权
文章标签： python 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/joshua_shi_t/article/details/128989041
这段代码展示了使用MATLAB进行2D数据可视化，包括设置坐标轴范围、颜色和标记样式。同时，也包含了OpenFOAM的后处理脚本，用于读取数据并生成2D图形，如2Dcircle.mcr显示了不同变量的轮廓线。另外，showBkg.mcr、showCase.mcr和showObj.mcr分别处理背景、案例和对象的显示。最后，showInt.py是ParaView脚本，用于读取OpenFOAM数据并进行渲染。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

plta1int.m
 

clearvars; clc; % plta1int.m
ya0=importdata("data\xyz_a0int1");
ya0=sortrows(ya0,1);
fig1=figure(1);
set(gca,'xlim',[0.3,0.7]); hold on;
set(gca,'ylim',[-0.1,1.1]); hold on;
x_0=ya0(:,1); y_0=ya0(:,4);
plt0=plot(x_0,y_0); hold on;
plt0.Color='black'; plt0.LineWidth=1.5; plt0.Marker='s'; plt0.MarkerFaceColor='black';
x_1=ya0(:,1); y_1=ya0(:,5);
plt1=plot(x_1,y_1); hold on;
plt1.Color='blue'; plt1.LineWidth=1.5; plt1.Marker='s'; plt1.MarkerFaceColor='blue';
x_2=ya0(:,1); y_2=ya0(:,6);
plt2=plot(x_2,y_2); hold on;
plt2.Color='red'; plt2.LineWidth=1.5; plt2.Marker='s'; plt2.MarkerFaceColor='red';
plt2.LineStyle='-';
leg0=legend([plt0,plt1,plt2]); leg0.String=["ref","invd","prest"];
leg0.Color='none'; leg0.FontSize=15; leg0.Position=[0.8,0.7,0.17,0.19];
saveas(gcf,'plt_1.png'); saveas(gcf,'plt_1.eps','epsc2');

2Dcircle.mcr

#!MC 1410
$!PAGECONTROL CLEAR
$!ReadDataSet  '"StandardSyntax" "1.0" "FEALoaderVersion" "436" "FILENAME_File" "C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\case\system\controlDict" "AutoAssignStrandIDs" "Yes"'
  DataSetReader = 'OpenFOAM (FEA)'
$!PlotType = Cartesian2D
$!TwoDAxis XDetail{RangeMin = 0.0}
$!TwoDAxis XDetail{RangeMax = 1.0}
$!TwoDAxis AxisMode = Independent
$!TwoDAxis YDetail{RangeMin = 0.0}
$!TwoDAxis YDetail{RangeMax = 1.0}
$!TwoDAxis YDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{AxisLine{Show = No}}
$!TwoDAxis XDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{AxisLine{Show = No}}
$!View MakeCurrentViewNice
$!FrameLayout ShowBorder = No
$!TwoDAxis GridArea{DrawBorder = No}
$!FieldLayers ShowMesh = Yes
$!FieldLayers ShowContour = Yes
$!GlobalTime SolutionTime = 0.05
$!FieldMap [4]  Shade{Show = No}
$!FieldMap [4]  Contour{ContourType = Lines}
$!SetContourVar 
  Var = 'alpha0'
  ContourGroup = 1
$!ContourLevels New
  ContourGroup = 1
  RawData
1
0.5
$!FieldMap [4]  Contour{LinePattern = Solid}
$!FieldMap [4]  Contour{PatternLength = 0.5}
$!FieldMap [4]  Contour{LineThickness = 0.5}
$!GlobalContour 1  Legend{Show = No}
$!FieldMap [4]  Mesh{Color = Custom2}
$!FieldMap [4]  Contour{Color = Black}
$!FrameName = 'Frame Background'
$!CreateNewFrame 
  XYPos
    {
    X = 1.0
    Y = 0.25
    }
  Width = 9.0
  Height = 8.0
$!FrameLayout IsTransparent = Yes
$!FrameLayout ShowBorder = No
$!FrameName = 'alpha0_int'
$!PlotType = Cartesian2D
$!TwoDAxis XDetail{RangeMin = 0.0}
$!TwoDAxis XDetail{RangeMax = 1.0}
$!TwoDAxis AxisMode = Independent
$!TwoDAxis YDetail{RangeMin = 0.0}
$!TwoDAxis YDetail{RangeMax = 1.0}
$!View MakeCurrentViewNice
$!TwoDAxis YDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{AxisLine{Show = No}}
$!TwoDAxis XDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{AxisLine{Show = No}}
$!FieldMap [4]  Shade{Show = No}
$!FieldMap [4]  Contour{ContourType = Lines}
$!GlobalTime SolutionTime = 0.05
$!FieldLayers ShowContour = Yes
$!FieldMap [4]  Shade{Show = No}
$!FieldMap [4]  Contour{ContourType = Lines}
$!SetContourVar 
  Var = 'alpha0_int'
  ContourGroup = 1
$!ContourLevels New
  ContourGroup = 1
  RawData
1
0.5
$!FieldMap [4]  Contour{LinePattern = Solid}
$!FieldMap [4]  Contour{PatternLength = 0.5}
$!FieldMap [4]  Contour{LineThickness = 0.5}
$!GlobalContour 1  Legend{Show = No}
$!FieldMap [4]  Mesh{Color = Custom2}
$!FieldMap [4]  Contour{Color = Blue}
$!CreateNewFrame 
  XYPos
    {
    X = 1.0
    Y = 0.25
    }
  Width = 9.0
  Height = 8.0
$!FrameLayout IsTransparent = Yes
$!FrameLayout ShowBorder = No
$!FrameName = 'alpha1'
$!PlotType = Cartesian2D
$!TwoDAxis XDetail{RangeMin = 0.0}
$!TwoDAxis XDetail{RangeMax = 1.0}
$!TwoDAxis AxisMode = Independent
$!TwoDAxis YDetail{RangeMin = 0.0}
$!TwoDAxis YDetail{RangeMax = 1.0}
$!View MakeCurrentViewNice
$!TwoDAxis YDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{AxisLine{Show = No}}
$!TwoDAxis XDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{AxisLine{Show = No}}
$!FieldMap [4]  Shade{Show = No}
$!FieldMap [4]  Contour{ContourType = Lines}
$!GlobalTime SolutionTime = 0.05
$!FieldLayers ShowContour = Yes
$!FieldMap [4]  Shade{Show = No}
$!FieldMap [4]  Contour{ContourType = Lines}
$!SetContourVar 
  Var = 'alpha1'
  ContourGroup = 1
$!ContourLevels New
  ContourGroup = 1
  RawData
1
0.5
$!FieldMap [4]  Contour{LinePattern = Solid}
$!FieldMap [4]  Contour{PatternLength = 0.5}
$!FieldMap [4]  Contour{LineThickness = 0.5}
$!GlobalContour 1  Legend{Show = No}
$!FieldMap [4]  Mesh{Color = Custom2}
$!FieldMap [4]  Contour{Color = Red}
$!AttachGeom 
  GeomType = Rectangle
  AnchorPos
    {
    X = 0.3122263836870571
    Y = 0.2832336203122713
    }
  Color = Custom4
  LineThickness = 0.4
  RawData
0.376862257719 -0.0656646713614 
$!CreateNewFrame 
  XYPos
    {
    X = 2.4
    Y = 1.55
    }
  Width = 6.3
  Height = 5.5
$!FrameLayout ShowBorder = No
$!FrameLayout IsTransparent = Yes
$!FrameName = 'Frame Front'
$!ExportSetup ExportFormat = JPEG
$!ExportSetup ExportRegion = CurrentFrame
$!ExportSetup ImageWidth = 1000
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\2Dcircle.jpg'
$!Export 
  ExportRegion = CurrentFrame
$!ExportSetup ExportFormat = EPS
$!ExportSetup ExportRegion = CurrentFrame
$!ExportSetup PrintRenderType = Image
$!ExportSetup ImageWidth = 1000
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\2Dcircle.eps'
$!Export 
  ExportRegion = CurrentFrame

showBkg.mcr

#!MC 1410
$!PAGECONTROL CLEAR
$!ReadDataSet  '"StandardSyntax" "1.0" "FEALoaderVersion" "436" "FILENAME_File" "C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\case\bkg-bak\system\controlDict" "AutoAssignStrandIDs" "Yes"'
  DataSetReader = 'OpenFOAM (FEA)'
$!PlotType = Cartesian2D
$!TwoDAxis XDetail{RangeMin = 0.0}
$!TwoDAxis XDetail{RangeMax = 1.0}
$!TwoDAxis AxisMode = Independent
$!TwoDAxis YDetail{RangeMin = 0.0}
$!TwoDAxis YDetail{RangeMax = 1.0}
$!TwoDAxis YDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{AxisLine{Show = No}}
$!TwoDAxis XDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{AxisLine{Show = No}}
$!View MakeCurrentViewNice
$!FrameLayout ShowBorder = No
$!TwoDAxis GridArea{DrawBorder = No}
$!FieldLayers ShowMesh = Yes
$!FieldMap [2]  Shade{Show = No}
$!RedrawAll
$!ExportSetup ExportFormat = EPS
$!ExportSetup ExportRegion = AllFrames
$!ExportSetup ImageWidth = 1000
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\showBkg.eps'
$!Export 
  ExportRegion = AllFrames
$!ExportSetup ExportFormat = JPEG
$!ExportSetup ExportRegion = AllFrames
$!ExportSetup ImageWidth = 1000
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\showBkg.jpg'
$!Export 
  ExportRegion = AllFrames

showCase.mcr

#!MC 1410
$!PAGECONTROL CLEAR
$!ReadDataSet  '"StandardSyntax" "1.0" "FEALoaderVersion" "436" "FILENAME_File" "C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\case\system\controlDict" "AutoAssignStrandIDs" "Yes"'
  DataSetReader = 'OpenFOAM (FEA)'
$!PlotType = Cartesian2D
$!TwoDAxis XDetail{RangeMin = 0.0}
$!TwoDAxis XDetail{RangeMax = 1.0}
$!TwoDAxis AxisMode = Independent
$!TwoDAxis YDetail{RangeMin = 0.0}
$!TwoDAxis YDetail{RangeMax = 1.0}
$!TwoDAxis YDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{AxisLine{Show = No}}
$!TwoDAxis XDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{AxisLine{Show = No}}
$!View MakeCurrentViewNice
$!FrameLayout ShowBorder = No
$!TwoDAxis GridArea{DrawBorder = No}
$!FieldLayers ShowContour = Yes
$!GlobalTime SolutionTime = 0.05
$!FieldMap [4]  Shade{Show = No}
$!SetContourVar 
  Var = 'alpha0'
  ContourGroup = 1
$!ContourLevels New
  ContourGroup = 1
  RawData
1
0.5
$!GlobalContour 1  Legend{Show = No}
$!CreateNewFrame 
  XYPos
    {
    X = 2.4
    Y = 1.55
    }
  Width = 6.3
  Height = 5.5
$!FrameLayout ShowBorder = No
$!FrameLayout IsTransparent = Yes
$!FrameName = 'Frame Front'
$!ExportSetup ExportFormat = JPEG
$!ExportSetup ExportRegion = CurrentFrame
$!ExportSetup ImageWidth = 1000
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\showCase.jpg'
$!Export 
  ExportRegion = CurrentFrame
$!ExportSetup ExportFormat = EPS
$!ExportSetup ExportRegion = CurrentFrame
$!ExportSetup ImageWidth = 1000
$!ExportSetup PrintRenderType = Image
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\showCase.eps'
$!Export 
  ExportRegion = CurrentFrame

showObj.mcr

#!MC 1410
$!PAGECONTROL CLEAR
$!ReadDataSet  '"StandardSyntax" "1.0" "FEALoaderVersion" "436" "FILENAME_File" "C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\case\obj0\system\controlDict" "AutoAssignStrandIDs" "Yes"'
  DataSetReader = 'OpenFOAM (FEA)'
$!PlotType = Cartesian2D
$!TwoDAxis XDetail{RangeMin = 0.0}
$!TwoDAxis XDetail{RangeMax = 1.0}
$!TwoDAxis AxisMode = Independent
$!TwoDAxis YDetail{RangeMin = 0.0}
$!TwoDAxis YDetail{RangeMax = 1.0}
$!TwoDAxis YDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis YDetail{AxisLine{Show = No}}
$!TwoDAxis XDetail{Ticks{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{TickLabel{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{Title{ShowOnAxisLine = No}}
$!TwoDAxis XDetail{AxisLine{Show = No}}
$!View MakeCurrentViewNice
$!FrameLayout ShowBorder = No
$!TwoDAxis GridArea{DrawBorder = No}
$!FieldLayers ShowMesh = Yes
$!FieldMap [2]  Shade{Show = No}
$!RedrawAll
$!ExportSetup ExportFormat = EPS
$!ExportSetup ExportRegion = AllFrames
$!ExportSetup ImageWidth = 1000
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\showObj.eps'
$!Export 
  ExportRegion = AllFrames
$!ExportSetup ExportFormat = JPEG
$!ExportSetup ExportRegion = AllFrames
$!ExportSetup ImageWidth = 1000
$!ExportSetup ExportFName = 'C:\myPC\myFile\SJTU_G11\Sundry\research\postPro\showObj.jpg'
$!Export 
  ExportRegion = AllFrames

showInt.py

# trace generated using paraview version 5.10.0
#import paraview
#paraview.compatibility.major = 5
#paraview.compatibility.minor = 10

#### import the simple module from the paraview
from paraview.simple import *
#### disable automatic camera reset on 'Show'
paraview.simple._DisableFirstRenderCameraReset()

# create a new 'OpenFOAMReader'
casefoam = OpenFOAMReader(registrationName='case.foam', FileName='C:\\myPC\\myFile\\SJTU_G11\\Sundry\\research\\postPro\\case\\case.foam')
casefoam.MeshRegions = ['internalMesh']
casefoam.CellArrays = ['U', 'alpha.water', 'alpha0', 'alpha0_int', 'alpha1', 'cellTypes', 'div', 'fl_vof_out', 'myProcNo_out', 'p_rgh', 'stencilOut', 'stencilOutWeights', 'thincPolyCoeffs0', 'thincPolyCoeffs1', 'thincPolyCoeffs2', 'tmpData', 'zoneID']

# get animation scene
animationScene1 = GetAnimationScene()

# update animation scene based on data timesteps
animationScene1.UpdateAnimationUsingDataTimeSteps()

# get active view
renderView1 = GetActiveViewOrCreate('RenderView')

# show data in view
casefoamDisplay = Show(casefoam, renderView1, 'UnstructuredGridRepresentation')

# trace defaults for the display properties.
casefoamDisplay.Representation = 'Surface'
casefoamDisplay.ColorArrayName = [None, '']
casefoamDisplay.SelectTCoordArray = 'None'
casefoamDisplay.SelectNormalArray = 'None'
casefoamDisplay.SelectTangentArray = 'None'
casefoamDisplay.OSPRayScaleArray = 'U'
casefoamDisplay.OSPRayScaleFunction = 'PiecewiseFunction'
casefoamDisplay.SelectOrientationVectors = 'U'
casefoamDisplay.ScaleFactor = 0.1
casefoamDisplay.SelectScaleArray = 'None'
casefoamDisplay.GlyphType = 'Arrow'
casefoamDisplay.GlyphTableIndexArray = 'None'
casefoamDisplay.GaussianRadius = 0.005
casefoamDisplay.SetScaleArray = ['POINTS', 'U']
casefoamDisplay.ScaleTransferFunction = 'PiecewiseFunction'
casefoamDisplay.OpacityArray = ['POINTS', 'U']
casefoamDisplay.OpacityTransferFunction = 'PiecewiseFunction'
casefoamDisplay.DataAxesGrid = 'GridAxesRepresentation'
casefoamDisplay.PolarAxes = 'PolarAxesRepresentation'
casefoamDisplay.ScalarOpacityUnitDistance = 0.23829607324923416
casefoamDisplay.OpacityArrayName = ['POINTS', 'U']

# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
casefoamDisplay.ScaleTransferFunction.Points = [0.0, 0.0, 0.5, 0.0, 1.1757813367477812e-38, 1.0, 0.5, 0.0]

# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
casefoamDisplay.OpacityTransferFunction.Points = [0.0, 0.0, 0.5, 0.0, 1.1757813367477812e-38, 1.0, 0.5, 0.0]

# reset view to fit data
renderView1.ResetCamera(False)

# get the material library
materialLibrary1 = GetMaterialLibrary()

# update the view to ensure updated data information
renderView1.Update()

# change representation type
casefoamDisplay.SetRepresentationType('Surface With Edges')

# set scalar coloring
ColorBy(casefoamDisplay, ('CELLS', 'cellTypes'))

# rescale color and/or opacity maps used to include current data range
casefoamDisplay.RescaleTransferFunctionToDataRange(False, False)

# show color bar/color legend
casefoamDisplay.SetScalarBarVisibility(renderView1, False)

# get color transfer function/color map for 'cellTypes'
cellTypesLUT = GetColorTransferFunction('cellTypes')

# get opacity transfer function/opacity map for 'cellTypes'
cellTypesPWF = GetOpacityTransferFunction('cellTypes')

#================================================================
# addendum: following script captures some of the application
# state to faithfully reproduce the visualization during playback
#================================================================
layout1 = GetLayout()
# get layout
layout1.SetSize(500, 500)

#-----------------------------------
# saving camera placements for views

# current camera placement for renderView1
renderView1.CameraPosition = [0.5, 0.5, 3.8460652149512318]
renderView1.CameraFocalPoint = [0.5, 0.5, 0.5]
renderView1.CameraParallelScale = 0.5644739300537772
renderView1.CameraParallelProjection = 1
renderView1.OrientationAxesVisibility = 0
renderView1.CameraParallelProjection = 1
renderView1.UseColorPaletteForBackground = 0
renderView1.Background = [1.0, 1.0, 1.0]
#--------------------------------------------
# uncomment the following to render all views
# RenderAllViews()
# alternatively, if you want to write images, you can use SaveScreenshot(...).
ExportView('C:/myPC/myFile/SJTU_G11/Sundry/research/postPro/showInt.eps', view=renderView1)

jmsyh
关注
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

jmsyh

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值