drN
Reputation
1,954
Top tag
Next privilege 2,000 Rep.
 Aug 8 revised FEM Solution desired for “Plate with orifice” deflection: Application of Boundary Conditions and use of Regions added 339 characters in body Aug 8 comment FEM Solution desired for “Plate with orifice” deflection: Application of Boundary Conditions and use of Regions No it doesn't: Encountered non-numerical value for a derivative at r == 0.005., but I suppose I'll amend my question: I would like to use a domain as specified by a region difference. the idea is to plot the deflections on the domain to show the utility of FEM to solve such problems. I hope that makes sense. The link you provide is relevant. Aug 8 comment FEM Solution desired for “Plate with orifice” deflection: Application of Boundary Conditions and use of Regions @chris Likely you are correct. I probably haven't used this correctly. However, I am also concerned about the other double and triple derivative conditions. Is there some way I can appoint these boundary conditions on the domain $\Omega$ without needing to resort to the (seemingly cryptic ;)') DirichletCondition and NeumannValue types. Aug 8 asked FEM Solution desired for “Plate with orifice” deflection: Application of Boundary Conditions and use of Regions Jul 30 awarded Popular Question Jul 29 accepted Potential flow over cylinder with FEM: Pressure field is flipped Jul 29 comment Potential flow over cylinder with FEM: Pressure field is flipped Is there some way I could choose the range of data represented by the PlotLegend? Thank you for your answer; I'll try this very soon. Jul 29 comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) @user21 You make a good point. I'll try this with MeshOrder Jul 28 comment Potential flow over cylinder with FEM: Pressure field is flipped @marcob yikes my mistake. The code I commented out may be treated as "old version". Sorry about that. Jul 28 asked Potential flow over cylinder with FEM: Pressure field is flipped Jul 18 awarded Popular Question Jul 10 comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) Yes, I just realized that. Now I am on the hunt for a generalized Hookes law that describes a cantilever with transverse load (original mma example). Would you mind putting the final statement in your answer "The results are unchanged from those produced by the code in the question but differ from that of the test case, because the differential equations are different." in bold? This is vital. I was under the impression that the plane stress operator in full would still solve the linear deflection situation while it solves only the transverse load condition. Jul 10 accepted Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) Jul 10 comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) I suppose I am unable to understand what the x and y directions are. Is x along the thickness of this bar element and y along its span? If so, should the load be an +xload? (of course the difference between - and + loads would be the sign of the deflection alone) Jul 9 comment How to solve heat equation with Robin type conditions with NDSolve? Are you sure of your boundary conditions? This example has Robin boundary conditions. (it is a link to a pdf file) Jul 9 comment Displaying a matrix with integral entries What I have noticed is m.n would "work" to produce a meaningful result for matrices m and n but m//MatrixForm . n//MatrixForm` does not compute. This would seem to be not explicit in the documentation? Jul 9 comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) @bbgodfrey yes. But in a purely axial load setting, there will be deflection only in axial x direction (u here) Jul 9 revised Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) cleaning up of code Jul 9 revised Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) deleted 529 characters in body Jul 9 comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation) @bbgodfrey But with inactive, the operators do execute. I still do get an interpolating function result. The real question for me is to understand why the example problem in the mma documentation uses inactive if it is of no use for this situation. I'll clean up the code like you say.