Thank you for your reply, David.As youadvised, I edited that part in "RB_construction.C" and ranmy code again.As a result, I could obtain the correct result.I was able to solve the problem much simpler than I expected thanks to your help.Best regards,Kang------------------------------------------------------------ ShinseongKang GraduateStudent PusanNationalUniversity, South KoreaTel.: +82-051-510-3052 H.P.:+82-010-9770-6595 E-mail:ss.k...@pusan.ac.kr
----------------------------------------------------------------- Original Message -----From : David Knezevic <david.kneze...@akselos.com>To : "Shinseong Kang" <ss.k...@pusan.ac.kr>Cc : "Libmesh user group" <libmesh-users@lists.sourceforge.net>Sent : 2019-10-08 21:53:36Subject : Re: [RB] Assembly for a Interior (neighbor) side between elementsHello,I think your approach is basically correct, butissue is that by default the code skips assembly on internal sides for the sake of efficiency (since usually we do not need to do assembly on internal sides). In order to not skip internal sides you can set "impose_internal_fluxes = true" on your RBConstruction object. That is false by default and if you look atRBConstruction::add_scaled_matrix_and_vector() you will see that you have to set it to true in order to loop over non-boundary sides.Let me know if that works for you.Best,DavidOn Tue, Oct 8, 2019 at 6:57 AM Shinseong Kang <ss.k...@pusan.ac.kr> wrote:Hello, all.I triedto definean output in a simply supported 3-D beam problem.My output is an average displacement in the middle of the beam.For simplicity, I only used two brick elements, i.e, ■|■, where"■" is an element, and"|" is a location I want to compute the average output.To do this, I utilized the codes in RB example 5.First, Ideclared "ThetaOutput" (parameter-dependent term) and alsodeclared "AssemblyOutput" (parameter-independent term) as follows:=======================================voidAssemblyOutput::boundary_assembly(FEMContext & c)}=======================================As you can see, The above code is very similar to"Assemblyp1" in "assembly.C" ofRB example 5.However, I could not obtain the correct output because the above code cannot deal with interior side between elements.In other words, theboundary information of FEMContext only include outer side.I tried to find another way, but I could not find the right solution.Therefore,I'd like your help.To summarize my questions areas follows:1. Is there any way to assemble an interior (neighbor) sidebetween elements?2. How to compute anaverage displacement for the middleof a beam in terms of the RB codes?Thank you.Best regards,Kang------------------------------------------------------------ ShinseongKang GraduateStudent PusanNationalUniversity, South KoreaTel.: +82-051-510-3052 H.P.:+82-010-9770-6595 E-mail:ss.k...@pusan.ac.kr ------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
