151 reputation
7
bio website
location
age
visits member for 1 year, 3 months
seen Mar 31 at 3:47

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+b*x+c*y+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
Feb
9
answered Latest Mathematica 9.0 can't use Format.m to generate optimized Fortran code
Feb
9
comment Latest Mathematica 9.0 can't use Format.m to generate optimized Fortran code
thanks. I have replaced ALL 'FilterOptions' to 'Sequence@@FilterRules' and the Format.m can be run successfully. Then I changed the beginning of my program from << format`Format` to << format`Format.m` the following errors still displayed: Get::noopen: Cannot open format`Format.m`. >> How to modify the program?