Tang Laoya
Reputation
Top tag
Next privilege 250 Rep.
 Jul 22 awarded Popular Question Aug 24 awarded Popular Question Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? By the way, could you please send me the code to output the whole 6*6 integrated matrix K? Thanks Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? Thanks for your kindly reply. He studied the high order problem, which is much more difficult than one order problem. Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? Thanks for your kindly reply. There are many papers about tetrahedron, for example: dspace.uta.edu/bitstream/handle/10106/933/… Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? Could you please also give me the code to output the whole matrix K by integrating? Thanks Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? Thanks for your so beautiful code:) Feb 10 revised Is it possible to give the closed-form of the stiffness matrix of triangular prism element? added 59 characters in body Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? In fact, the formula of $$N_i$$is just a convenient expression. We can write as follows for clearness:$$N_1=\frac{1}{2}(1-\zeta)(1-\xi-\eta)\\N_4=\frac{1}{2}(1+\zeta)(1-\xi-‌​\eta)$$etc. Thanks Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? Hi DumpsterDoofus, thanks for your kindly reply. The $$x_i, y_i, z_i$$ are coordinates of vertices of original triangular prism. The x and y-coordinates of the vertices are used in the integration. Yes, $$J^{-T}$$ is the shorthand for $$(J^{-1})^{T}$$Thanks Feb 10 comment Is it possible to give the closed-form of the stiffness matrix of triangular prism element? oh, sorry, $$\zeta_i$$ are the z-coordinate of standard triangular prism, i.e., $$\zeta_i=-1, i=1,2,3, \\ \zeta_i=1, i=4,5,6.$$ Thanks Feb 9 comment How to generate closed-form stiff matrix of triangular prism element by mathematica OK, thanks Rahul Narain for your kindly remind. Feb 9 comment How to generate closed-form stiff matrix of triangular prism element by mathematica Hi Bill, thanks for your kindly reply. I have posted a new question here: mathematica.stackexchange.com/questions/41979/… could you please help me to take a look at it? Feb 9 asked Is it possible to give the closed-form of the stiffness matrix of triangular prism element? Feb 9 awarded Teacher Feb 9 asked Output most efficient Fortran code by mathematica Feb 9 comment Latest Mathematica 9.0 can't use Format.m to generate optimized Fortran code Hi rm, thanks for your process. The old problem is solved by used an updated Format.m and placed it to the correct location. But I have some new problems as I described. I think that they are the same kind of problems so I placed it here. Thanks. Feb 9 revised Latest Mathematica 9.0 can't use Format.m to generate optimized Fortran code added 2653 characters in body Feb 9 comment How to generate closed-form stiff matrix of triangular prism element by mathematica Hi Bill, thanks for your kindly reply. I need to: 1) create an interpolation function in triangular prism just like V=a+bx+cy+d*z for tetrahedron (4 vertices and 4 unknowns); 2) create the shape functions (for tetrahedron it is the volume coordinate but I don't know how to express them for triangular prism); 3) calculate the stiffness matrix, for example, the term Int(d^2 N/dx^2+d^2 N/dy^2+d^2 N/dz^2) dxdydz, where N is the shape functions 4) give the closed-form express by mathematica. Thanks Feb 9 asked How to generate closed-form stiff matrix of triangular prism element by mathematica