




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
1、Null-field approach for Laplace problems with circular boundaries using degenerate kernels沈文成 陳正宗 時間:13:20 14:00地點:河工二館 307 室BEM course May 13, 2008 (typicalBVP-L.ppt)1Null-field approach for LaplacOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution an
2、d boundary density Adaptive observer system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions2OutlinesMotivation and literatOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Ada
3、ptive observer system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions3OutlinesMotivation and literatMotivation and literature reviewFictitious BEMBEM/BIEMNull-field approachBump contourLimit process Singular and hypersingularRegularImproper integ
4、ralCPV and HPVIll-posedFictitious boundaryCollocation point4Motivation and literature reviPresent approach1. No principal value2. Well-posedAdvantages of degenerate kernelDegenerate kernelFundamental solutionCPV and HPVNo principal value5Present approach1. No principEngineering problem with arbitrar
5、y geometriesDegenerate boundaryCircular boundaryStraight boundaryElliptic boundary(Fourier series)(Legendre polynomial)(Chebyshev polynomial)(Mathieu function)6Engineering problem with arbitMotivation and literature reviewAnalytical methods for solving Laplace problems with circular holesConformal m
6、appingBipolar coordinateSpecial solutionLimited to doubly connected domainLebedev, Skalskaya and Uyand, 1979, “Work problem in applied mathematics”, Dover PublicationsChen and Weng, 2001, “Torsion of a circular compound bar with imperfect interface”, ASME Journal of Applied MechanicsHonein, Honein a
7、nd Hermann, 1992, “On two circular inclusions in harmonic problem”, Quarterly of Applied Mathematics7Motivation and literature reviFourier series approximationLing (1943) - torsion of a circular tubeCaulk et al. (1983) - steady heat conduction with circular holesBird and Steele (1992) - harmonic and
8、 biharmonic problems with circular holesMogilevskaya et al. (2002) - elasticity problems with circular boundaries8Fourier series approximationLiContribution and goalHowever, they didnt employ the null-field integral equation and degenerate kernels to fully capture the circular boundary, although the
9、y all employed Fourier series expansion.To develop a systematic approach for solving Laplace problems with multiple holes is our goal.9Contribution and goalHowever, OutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Adaptive obser
10、ver system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions10OutlinesMotivation and literatBoundary integral equation and null-field integral equationInterior caseExterior caseNull-field integral equation11Boundary integral equation andOutlinesMot
11、ivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Adaptive observer system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions12OutlinesMotivation and literatExpansions of fundamental soluti
12、on and boundary densityDegenerate kernel - fundamental solutionFourier series expansions - boundary density13Expansions of fundamental soluSeparable form of fundamental solution (1D)Separable propertycontinuousdiscontinuous14Separable form of fundamental Separable form of fundamental solution (2D)15
13、Separable form of fundamental Boundary density discretizationFourier seriesEx . constant elementPresent methodConventional BEM16Boundary density discretizatioOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Adaptive observer sys
14、tem Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions17OutlinesMotivation and literatAdaptive observer systemcollocation point18Adaptive observer systemcollocOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental
15、 solution and boundary density Adaptive observer system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions19OutlinesMotivation and literatVector decomposition technique for potential gradientSpecial case (concentric case) :Non-concentric case:True n
16、ormal direction20Vector decomposition techniqueOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Adaptive observer system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions21Outl
17、inesMotivation and literatLinear algebraic equationwhereColumn vector of Fourier coefficients(Nth routing circle)Index of collocation circleIndex of routing circle 22Linear algebraic equationwhereExplicit form of each submatrix Upk and vector tkFourier coefficientsTruncated terms of Fourier seriesNu
18、mber of collocation points23Explicit form of each submatriFlowchart of present methodPotential of domain pointAnalyticalNumericalAdaptive observer systemDegenerate kernelFourier seriesLinear algebraic equation Collocation point and matching B.C.Fourier coefficientsVector decompositionPotential gradi
19、ent24Flowchart of present methodPotComparisons of conventional BEM and present methodBoundarydensitydiscretizationAuxiliarysystemFormulationObserversystemSingularityConventionalBEMConstant,Linear,QurdratureFundamentalsolutionBoundaryintegralequationFixedobserversystemCPV, RPVand HPVPresentmethodFour
20、ierseriesexpansionDegeneratekernelNull-fieldintegralequationAdaptiveobserversystemNoprincipalvalue25Comparisons of conventional BEOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Adaptive observer system Vector decomposition tec
21、hnique Linear algebraic equationNumerical examplesDegenerate scaleConclusions26OutlinesMotivation and literatNumerical examplesSteady state heat conduction problemsElectrostatic potential of wiresFlow of an ideal fluid pass cylindersA circular bar under torqueAn infinite medium under antiplane shear
22、Half-plane problems27Numerical examplesSteady stateNumerical examplesSteady state heat conduction problemsElectrostatic potential of wiresFlow of an ideal fluid pass cylindersA circular bar under torqueAn infinite medium under antiplane shearHalf-plane problems28Numerical examplesSteady stateSteady
23、state heat conduction problemsCase 1Case 229Steady state heat conduction pSteady state heat conduction problemsCase 3Case 430Steady state heat conduction pCase 1: Isothermal lineExact solution(Carrier and Pearson)BEM-BEPO2D(N=21)FEM-ABAQUS(1854 elements)Present method(M=10)31Case 1: Isothermal lineE
24、xact sRelative error of flux on the small circle32Relative error of flux on the Convergence test - Parsevals sum for Fourier coefficientsParsevals sum33Convergence test - Parsevals Case 2: Isothermal lineCaulks data (1983)IMA Journal of Applied MathematicsPresent method (M=10)FEM-ABAQUS(6502 element
25、s)34Case 2: Isothermal lineCaulksCase 3: Isothermal lineFEM-ABAQUS(8050 elements)Present method (M=10)Caulks data (1983)IMA Journal of Applied Mathematics35Case 3: Isothermal lineFEM-ABACase 4: Isothermal lineFEM-ABAQUS(8050 elements)Present method (M=10)Caulks data (1983)IMA Journal of Applied Math
26、ematics36Case 4: Isothermal lineFEM-ABANumerical examplesSteady state heat conduction problemsElectrostatic potential of wiresFlow of an ideal fluid pass cylindersA circular bar under torqueAn infinite medium under antiplane shearHalf-plane problems37Numerical examplesSteady stateElectrostatic poten
27、tial of wiresHexagonal electrostatic potentialTwo parallel cylinders held positive and negative potentials38Electrostatic potential of wirContour plot of potentialExact solution (Lebedev et al.)Present method (M=10)39Contour plot of potentialExactContour plot of potentialOnishis data (1991)Present m
28、ethod (M=10)40Contour plot of potentialOnishNumerical examplesSteady state heat conduction problemsElectrostatic potential of wiresFlow of an ideal fluid pass cylindersA circular bar under torqueAn infinite medium under antiplane shearHalf-plane problems41Numerical examplesSteady stateFlow of an ide
29、al fluid pass two parallel cylinders is the velocity of flow far from the cylinders is the incident angle42Flow of an ideal fluid pass twVelocity field in different incident anglePresent method (M=10)Present method (M=10)43Velocity field in different iNumerical examplesSteady state heat conduction p
30、roblemsElectrostatic potential of wiresFlow of an ideal fluid pass cylindersA circular bar under torqueAn infinite medium under antiplane shearHalf-plane problems44Numerical examplesSteady stateTorsion bar with circular holes removedThe warping functionBoundary condition whereonTorque45Torsion bar w
31、ith circular holeAxial displacement with two circular holesPresent method (M=10)Caulks data (1983)ASME Journal of Applied MechanicsDashed line: exact solutionSolid line: first-order solution46Axial displacement with two ciAxial displacement with three circular holesPresent method (M=10)Caulks data (
32、1983)ASME Journal of Applied MechanicsDashed line: exact solutionSolid line: first-order solution47Axial displacement with three Axial displacement with four circular holesPresent method (M=10)Caulks data (1983)ASME Journal of Applied MechanicsDashed line: exact solutionSolid line: first-order solut
33、ion48Axial displacement with four cNumerical examplesSteady state heat conduction problemsElectrostatic potential of wiresFlow of an ideal fluid pass cylindersA circular bar under torqueAn infinite medium under antiplane shearHalf-plane problems49Numerical examplesSteady stateInfinite medium under a
34、ntiplane shearThe displacementBoundary conditionTotal displacementon50Infinite medium under antiplanShear stress zq around the hole of radius a1 (x axis)Present method (M=20)Honeins data (1992)Quarterly of Applied Mathematics51Shear stress zq around the hoShear stress zq around the hole of radius a1
35、 (y axis)Present method (M=20)Honeins data (1992)Quarterly of Applied Mathematics52Shear stress zq around the hoShear stress zq around the hole of radius a1 (45 degrees)Present method (M=20)Honeins data (1992)Quarterly of Applied Mathematics53Shear stress zq around the hoShear stress zq around the h
36、ole of radius a1 (Touching)Present methoddiscontinuousdiscontinuousHoneins data (1992)Quarterly of Applied MathematicsGibbs phenomenon54Shear stress zq around the hoTwo equivalent approachesDisplacement approachStress approachPresent methodBird and Steele (1992)ASME Journal of Applied Mechanics55Two
37、 equivalent approachesDisplShear stress zq around the hole of radius a1Present method (M=20)Present method (M=20)Steeles data (1992)Stress approachDisplacement approachHoneins data (1992)5.3485.3494.6475.34513.13%0.02%Analytical0.06%56Shear stress zq around the hoConvergence of stress zq at q=45 deg
38、rees versus R057Convergence of stress zq atThree circular holes with centers on the x axis58Three circular holes with centThree circular holes with centers on the y axis59Three circular holes with centThree circular holes with centers on the line making 45 degrees60Three circular holes with centNume
39、rical examplesSteady state heat conduction problemsElectrostatic potential of wiresFlow of an ideal fluid pass cylindersA circular bar under torqueAn infinite medium under antiplane shearHalf-plane problems61Numerical examplesSteady stateHalf-plane problemsDirichlet boundary condition(Lebedev et al.
40、)Mixed-type boundary condition(Lebedev et al.)62Half-plane problemsDirichlet bDirichlet problemExact solution (Lebedev et al.)Present method (M=10)Isothermal line63Dirichlet problemExact solutioMixed-type problemExact solution (Lebedev et al.)Present method (M=10)Isothermal line64Mixed-type problemE
41、xact solutiOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Adaptive observer system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions65OutlinesMotivation and literatNumerical
42、instability in BEMAnnular caseInterior caseMax errorDegenerate scaleInternational Journal forNumerical Methods in EngineeringEngineering Analysiswith Boundary Elements Matrix singularErrorSingularvalue66Numerical instability in BEMAnDegenerate scale in the multiply connected problema1 =1.0, influenc
43、e matrix U is singular67Degenerate scale in the multipTreatments of degenerate scale problemMethod of adding a rigid body termCHEEF conceptSingularSingularAuxiliary constraint NonsingularCHEEF pointPromote rank68Treatments of degenerate scaleThe minimum singular value versus radius a1Degenerate scal
44、eNumerical failure69The minimum singular value verOutlinesMotivation and literature reviewMathematical formulation Expansions of fundamental solution and boundary density Adaptive observer system Vector decomposition technique Linear algebraic equationNumerical examplesDegenerate scaleConclusions70O
45、utlinesMotivation and literatConclusionsA systematic approach using degenerate kernels, Fourier series and null-field integral equation has been successfully proposed to solve Laplace problems with circular boundaries.Numerical results agree well with available exact solutions, Caulks data, Onishis data and FEM (ABAQUS) for only few terms of Fourier series.71ConclusionsA systematic approaConclusionsMethod of adding a rigid body term and CHEEF approach have been successfully adopted to overcome the degenerate scale for multiply
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 隨筆關(guān)于高一學(xué)生2025年終總結(jié)(9篇)
- 2024年二手車評估師考試復(fù)習(xí)計劃試題及答案
- 2024年小自考行政管理在職人員的必修課試題及答案
- 2025年高分子合成材料項目建議書
- 2024年寵物營養(yǎng)師考試難點分析試題及答案
- 2024年計算機基礎(chǔ)考試基本概念試題及答案
- 2024年計算機基礎(chǔ)重要知識回顧及試題和答案
- 寵物腸道微生物與營養(yǎng)的關(guān)系試題及答案
- 統(tǒng)計圖形與數(shù)據(jù)呈現(xiàn)試題及答案
- 2024年汽車維修工資格考試難度題目及答案
- 合伙開辦教育培訓(xùn)機構(gòu)合同范本
- 嵌入式機器視覺流水線分揀系統(tǒng)設(shè)計
- 江蘇省鹽城市東臺市第一教育聯(lián)盟2024-2025學(xué)年七年級下學(xué)期3月月考英語試題(原卷版+解析版)
- 湖南省2025屆高三九校聯(lián)盟第二次聯(lián)考歷史試卷(含答案解析)
- 2024年全國職業(yè)院校技能大賽(高職組)安徽省集訓(xùn)選拔賽“電子商務(wù)”賽項規(guī)程
- 2025年中考數(shù)學(xué)復(fù)習(xí):翻折問題(含解析)
- (統(tǒng)編版2025新教材)語文七下全冊知識點
- 家具全屋定制的成本核算示例-成本實操
- 第二單元第1課《精彩瞬間》第2課時 課件-七年級美術(shù)下冊(人教版2024)
- 放射科危急值報告流程質(zhì)量控制方案
- 2023-2024學(xué)年江蘇南京江寧區(qū)七年級下冊語文期中試卷及答案A卷
評論
0/150
提交評論