21.01.2015 Views

COMSOL Multiphysics™

COMSOL Multiphysics™

COMSOL Multiphysics™

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

femlin<br />

Purpose<br />

femlin<br />

Solve linear or linearized stationary PDE problem.<br />

Syntax<br />

Description<br />

fem.sol = femlin(fem,...)<br />

fem = femlin(fem,'Out', {'fem'},...)<br />

[Ke,Le,Null,ud] = femlin(fem,...)<br />

[Kl,Ll,Nnp] = femlin(fem,...)<br />

[Ks,Ls] = femlin(fem,...)<br />

fem.sol = femlin('In',{'K' K 'L' L 'M' M 'N' N},...)<br />

fem.sol = femlin(fem) solves a linear or linearized stationary PDE problem<br />

described by the (possibly extended) FEM structure fem. See femstruct for details<br />

on the FEM structure.<br />

fem.sol = femlin(fem, 'pname','P', 'plist',list,...) solves a linear or<br />

linearized stationary PDE problem for several values of the parameter P. The values<br />

of the parameter P are given in the vector list.<br />

fem = femlin(fem,'out',{'fem'}) modifies the FEM structure to include the<br />

solution structure, fem.sol.<br />

[Ke,Le,Null,ud] = femlin(fem) partially solves the PDE problem by<br />

eliminating the constraints. The solution of PDE problem can be obtained by the<br />

scripting command u = u0+Null*(Ke\Le)+ud, where u0 is the linearization point.<br />

[Kl,Ll,Nnp] = femlin(fem) partially solves the PDE problem by using the<br />

Lagrange method. The solution can then be obtained by u = Kl\Ll, and then u =<br />

u0+u(1:Nnp).<br />

[Ks,Ls] = femlin(fem) partially solves the PDE problem by approximating the<br />

constraints with stiff springs. The solution to the PDE problem is u = u0+Ks\Ls.<br />

fem.sol = femlin('in',{'K' K 'N' N 'L' L 'M' M}) solves a pre-assembled<br />

PDE problem.<br />

u = femlin('in’,{'K' K 'L' L},’out’,’u’) is equivalent to solving the linear<br />

system using u = K\L, with the important difference that you have access to all<br />

linear system solvers (except Geometric multigrid) using the Linsolver property.<br />

Consider the FEM discretization of a stationary PDE problem:<br />

0 = L–<br />

N T Λ ,<br />

M<br />

where L, N, and M depend on the solution vector U. femlin solves the linearized<br />

form of this problem:<br />

124 | CHAPTER 1: COMMAND REFERENCE

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!